for j in range(nPt):
pt_ij = output_i.GetPoint(j)
pt_ij_norm = [(pt_ij[0] - center_of_mass_i) / size_i,
(pt_ij[1] - center_of_mass_i) / size_i]
pt_mat_i[0, j] = pt_ij_norm[0]
pt_mat_i[1, j] = pt_ij_norm[1]
shape_i = manifolds.kendall2D(nPt)
scale_i = manifolds.pos_real(1)
scale_i.SetPoint(scale_obs_list[i])
shape_i.SetPoint(pt_mat_i)
scale_shape_i = manifolds.scale_kendall2D(nPt)
scale_shape_i.SetPoint([scale_i, shape_i])
# Covariates
if i in male_list:
s_i = 1
else:
s_i = 0
# if i in autism_list:
# autism_i = 1
# else:
# autism_i = 0
# if i in autismspec_list:
# autismspec_i = 1
for j in range(nPt):
pt_ij = output_i.GetPoint(j)
pt_ij_norm = [(pt_ij[0] - center_of_mass_i) / size_i,
(pt_ij[1] - center_of_mass_i) / size_i]
pt_mat_i[0, j] = pt_ij_norm[0]
pt_mat_i[1, j] = pt_ij_norm[1]
shape_i = manifolds.kendall2D(nPt)
scale_i = manifolds.pos_real(1)
scale_i.SetPoint(scale_obs_list[i])
shape_i.SetPoint(pt_mat_i)
scale_shape_i = manifolds.scale_kendall2D(nPt)
scale_shape_i.SetPoint([scale_i, shape_i])
# Covariates
if i in male_list:
s_i = 1
else:
s_i = 0
# if i in autism_list:
# autism_i = 1
# else:
# autism_i = 0
# if i in autismspec_list:
# autismspec_i = 1