Postprocessing Code for 'Falling Heads': investigating reflexive responses to head-neck perturbationsdoi:10.18419/darus-2526DaRUS2022-04-011Wochner, Isabell; Schmitt, Syn, 2022, "Postprocessing Code for 'Falling Heads': investigating reflexive responses to head-neck perturbations", https://doi.org/10.18419/darus-2526, DaRUS, V1Postprocessing Code for 'Falling Heads': investigating reflexive responses to head-neck perturbationsdoi:10.18419/darus-2526Wochner, IsabellSchmitt, Syn390740016info:eu-repo/grantAgreement/EC/H2020/768947DaRUSWochner, IsabellSchmitt, SynWochner, Isabell2022-02-17Computer and Information ScienceEngineeringMedicine, Health and Life SciencesPhysicsReflex BehaviourHead-Neck PerturbationsMotor Control<p>This dataset contains the code necessary to postprocess the experimental trajectory data of the Falling Heads experiment and perform the inverse dynamics analysis as described in the paper "’Falling Heads’: investigating reflexive responses to head-neck perturbations" by Isabell Wochner, Lennart V. Nölle, Oleksandr V. Martynenko and Syn Schmitt. Always cite the paper together with this dataset because this dataset is not self-explanatory.</p>
<p>The motion analyses were performed with custom software written in Matlab (Mathworks, Natick, MA) based on the recorded marker positions. From these recorded trajectories, we calculated displacements, velocities and accelerations and processed the signals with a Butterworth low-pass filter (cut-off frequency 15 Hz; fourth-order).</p>The folder contains 14 matlab codes, which you can and should run in their consecutive order to perform the postprocessing of the experimental trajectories (starting from A0.. up to A13). <br>
Furthermore, the matlab data files named "kinematic_analysis_expdata" contains the experimental data which is needed to run the scripts (they should be in the same folder).Wochner, Isabell, 2022, "Trajectory Data for Falling Heads: investigating reflexive responses to head-neck perturbations", <a href="https://doi.org/10.18419/darus-1132">https://doi.org/10.18419/darus-1132</a>, DaRUS.<br>Wochner, Isabell, 2022, "Neck muscle data for Falling Heads: investigating reflexive responses to head-neck perturbations", <a href="https://doi.org/10.18419/darus-1145">https://doi.org/10.18419/darus-1145</a>, DaRUS.<br>Wochner, Isabell; Schmitt, Syn, 2022, "EMG data for Falling Heads: investigating reflexive responses in active Human Body Models", <a href="https://doi.org/10.18419/darus-1038">https://doi.org/10.18419/darus-1038</a>, DaRUS.Wochner, I., Nölle, L. V., Martynenko, O. V., & Schmitt, S. (2022). ‘Falling heads’: investigating reflexive responses to head-neck perturbations. Biomedical engineering online, 21(1), article no. 25.10.1186/s12938-022-00994-9Wochner, I., Nölle, L. V., Martynenko, O. V., & Schmitt, S. (2022). ‘Falling heads’: investigating reflexive responses to head-neck perturbations. Biomedical engineering online, 21(1), article no. 25.A0_LoadTrajectories.mThis script is needed to load the experimental trajectories. Furthermore it saves the x and y displacement of the centre of gravity (head), the x-y displacement of a second node, and the x-y displacement of the chin and neck marker as "node_traj.mat" file. text/plainA10_plot_fallinghead_verticaldisp_shaded.mThis script plots the differences in vertical displacement for different ages and sexes; run this script to reproduce Figure 3 of the paper. text/plainA11_CalculateStiffness_onlyk.mThis script plots the effective stiffness for all participants and all trials; run this script to reproduce Figure 2 of the paper. text/plainA12_PlotDisplacement.mThis scripts plots the vertical displacmement of the 'falling heads' experiment for all participants and all trials; run this script to reproduce Figure 1 of the paper. text/plainA13_PlotAccelerations.mtext/plainA1_Calculate_LinearAcceleration_x.mThis script calculates and saves the linear acceleration in x direction of the centre of gravity (head). It creates two figures as output (for the supine and prone case) and saves the calculated data as "linearacc_x.mat" file. text/plainA2_Calculate_LinearAcceleration_y.mThis script calculates and saves the linear acceleration in y direction of the centre of gravity (head). It creates two figures as output (for the supine and prone case) and saves the calculated data as "linearacc_y.mat" file.
Furthermore, it produces the first half of the latex tables Table 1 and Table 3 (linear peak accelerations) of the paper. text/plainA3_CalculateAngle.mThis script calculates the angle (rotational displacement) of the head. It creates two figures as output (for the supine and prone case) and saves the calculated data as "Angle.mat" file. text/plainA4_Calculate_RotationalAcceleration.mThis script calculates and saves the rotational acceleration of the centre of gravity (head). It creates two figures as output (for the supine and prone case) and saves the calculated data as "rotationalpeakacc.mat" file.
Furthermore, it produces the second half of the latex tables Table 1 and Table 3 (rotational peak accelerations) of the paper. text/plainA5_CalculateHeadmassAndMomentOfInertia.mThis script calculates and saves the headmass and the moment of inertia of all participants based on their size and height according to description given in the paper (section 5.5.). text/plainA6_Calculate_rx_ry.mThis script calculates and saves the moment arms r_x and r_y of the headneck joint necessary to calculate the inverse dynamics. text/plainA7_CalculateFxy.mThis script calculates the force acting in the sagittal x-y plane on the rigid head-neck segment. The force in x-direction can be calculated according to equation (6) in the paper, the force in y-direction according to equation (7). text/plainA8_CalculateMz.mThis script calculates the joint torque M_net exerted on the head-neck segment by the trunk segment at the connecting joint. This can be calculated according to equation (8) in the paper. text/plainA9_PlotInvDynamikAnalysis.mThis script plots the results of the inverse dynamic analysis. With the current settings, it plots Fig (6). If you want to reproduce Figure 7, you need to switch the "caseflagreproduce_figure7_exptraj" to on, and the other caseflagreproduces off; same applied to reproduce figure S1 of the supplementary material. text/plainkinematic_analysis_expdata.matThis data file is needed to run A0_LoadTrajectories and contains the trajectory data of the three markers. application/matlab-matkinematic_analysis_expdata_chinmarker.matThis data file is needed to run A0_LoadTrajectories and contains the trajectory data of the chin marker. application/matlab-matkinematic_analysis_expdata_neckmarker.matThis data file is needed to run A0_LoadTrajectories and contains the trajectory data of the neck marker. application/matlab-matmatrix2latex.mtext/plain