#%%
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}
x = db_elev.iloc[0]['traj_interp'].common_ht_range_fine
init_graphing(params.backend)
plt.close('all')
fig = plt.figure(figsize=(190 / 25.4, 190 / 25.4))
axs = fig.subplots(3, 3)
# set axes limits
ax_limits = [(-35, 10), (-40, 20), (-20, 40)]
for i in range(3):
for j in range(3):
axs[i, j].set_ylim(ax_limits[i][0], ax_limits[i][1])
axs[i, j].yaxis.set_major_locator(ticker.MultipleLocator(10))
style_axes(axs[i, j],
'Humerothoracic Elevation (Deg)' if i == 2 else None,
'Axial Rotation (Deg)' if j == 0 else None)
zip(*db_elev.apply(shr_compute, axis=1))
#%%
if bool(distutils.util.strtobool(params.parametric)):
spm_test = spm1d.stats.ttest_paired
infer_params = {}
else:
spm_test = spm1d.stats.nonparam.ttest_paired
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}
plot_utils.init_graphing(params.backend)
plt.close('all')
x = np.arange(0, 100 + 0.25, 0.25)
fig = plt.figure(figsize=(110 / 25.4, 190 / 25.4))
axs = fig.subplots(3, 1)
for i in range(3):
style_axes(axs[i],
'Humerothoracic Elevation (Deg)' if i == 2 else None,
'Elevation (deg)')
axs[i].set_ylim(0, 12)
leg_mean = []
for idx, (activity, activity_df) in enumerate(
db_elev.groupby('Activity', observed=True)):
