if age < 40 else '>45')
if params.excluded_trials:
db = db[~db['Trial_Name'].str.contains('|'.join(params.
excluded_trials))]
db['Trial'].apply(pre_fetch)
# relevant parameters
output_path = Path(params.output_dir)
# logging
fileConfig(config_dir / 'logging.ini', disable_existing_loggers=False)
log = logging.getLogger(params.logger_name)
db_elev = db.loc[db['Trial_Name'].str.contains('_CA_|_SA_|_FE_')].copy()
prepare_db(db_elev, params.torso_def, params.scap_lateral,
params.dtheta_fine, params.dtheta_coarse,
[params.min_elev, params.max_elev])
#%%
if bool(distutils.util.strtobool(params.parametric)):
spm_test = partial(spm1d.stats.ttest2, equal_var=False)
infer_params = {}
else:
spm_test = spm1d.stats.nonparam.ttest2
infer_params = {'force_iterations': True}
alpha = 0.05
color_map = plt.get_cmap('Dark2')
markers = ['^', 'o', 's', '*']
act_row = {'ca': 0, 'sa': 1, 'fe': 2}
BiplaneViconSubject,
include_anthro=True)
db['age_group'] = db['Age'].map(lambda age: '<35' if age < 40 else '>45')
if params.excluded_trials:
db = db[~db['Trial_Name'].str.contains('|'.join(params.excluded_trials)
)]
db = db.loc[[params.trial_name]]
# logging
fileConfig(config_dir / 'logging.ini', disable_existing_loggers=False)
log = logging.getLogger(params.logger_name)
prepare_db(db,
params.torso_def,
params.scap_lateral,
params.dtheta_fine,
params.dtheta_coarse, [params.min_elev, params.max_elev],
should_fill=False,
should_clean=False)
trial = db.loc[params.trial_name]
#%%
# calculate st true axial rotation "manually"
st_angvel_proj = extended_dot(trial.st.ang_vel, trial.ht.rot_matrix[:, :,
1])
st_axial_rot_manual = cumtrapz(st_angvel_proj, dx=trial.st.dt, initial=0)
st_axial_rot = trial.ht.true_axial_rot - trial.gh.true_axial_rot
# plot
start_idx = trial.up_down_analysis.max_run_up_start_idx
end_idx = trial.up_down_analysis.max_run_up_end_idx
