#%% 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)):