center_of_mass_i += (pt_ij[0] + pt_ij[1] + pt_ij[2])
center_of_mass_i = center_of_mass_i / float(nPt)
pt_mat_i = np.zeros([2, nPt])
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
center_of_mass_i += (pt_ij[0] + pt_ij[1] + pt_ij[2])
center_of_mass_i = center_of_mass_i / float(nPt)
pt_mat_i = np.zeros([2, nPt])
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)
shape_i.SetPoint(pt_mat_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
center_of_mass_i += (pt_ij[0] + pt_ij[1] + pt_ij[2])
center_of_mass_i = center_of_mass_i / float(nPt)
pt_mat_i = np.zeros([2, nPt])
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)
shape_i.SetPoint(pt_mat_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
print(center_of_mass_i)
pt_mat_i = np.zeros([2, nPt])
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_list.append(pt_mat_i)
shape_0 = manifolds.kendall2D(nPt)
print(len(shape_list[0]))
print(nPt)
shape_0.SetPoint(shape_list[0])
print(shape_0.pt[0, 1])
# Log Map Check
shape0 = manifolds.kendall2D(nPt)
shape0.SetPoint(shape_list[0])
shape1 = manifolds.kendall2D(nPt)
shape1.SetPoint(shape_list[1])
print("Construction Done")
tVec_check = shape0.LogMap(shape1)