'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/HMG_ScaleKendall_HR_ADOS_New/mean_shape.pt'
)
for i in range(len(shape_list_all)):
shape_i = shape_list_all[i]
t_i = t_list_all[i]
s_i = s_list_all[i]
hr_i = hr_list_all[i]
lr_i = lr_list_all[i]
autism_i = autism_list_all[i]
autismSpec_i = autismspec_list_all[i]
healthy_i = none_list_all[i]
ADOS_i = ADOS_list_all[i]
beta_i = manifolds.scale_kendall2D_tVec(nPt)
for k in range(2):
for j in range(nPt):
beta_i.tVector[1].tVector[k, j] = (
est_beta1.tVector[1].tVector[k, j] * s_i +
est_beta2.tVector[1].tVector[k, j] * autism_i)
beta_i.tVector[0].tVector[0] = (est_beta1.tVector[0].tVector[0] * s_i +
est_beta2.tVector[0].tVector[0] * autism_i)
p1_i = est_beta0.ExponentialMap(beta_i)
gamma0_p1_i = est_beta0.ParallelTranslateAtoB(est_beta0, p1_i, est_gamma0)
gamma1_p1_i = est_beta0.ParallelTranslateAtoB(est_beta0, p1_i, est_gamma1)
gamma2_p1_i = est_beta0.ParallelTranslateAtoB(est_beta0, p1_i, est_gamma2)
'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/HMG_ScaleKendall_HR_ADOS_New/est_gamma1.tVec'
)
# gamma2 * ADOS t
est_gamma2 = tangent_slope_arr[1]
est_gamma2.Write(
'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/HMG_ScaleKendall_HR_ADOS_New/est_gamma2.tVec'
)
# Read
est_beta0 = manifolds.scale_kendall2D(nPt)
est_beta0.Read(
'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/HMG_ScaleKendall_HR_ADOS_New/est_beta0.pt'
)
est_beta1 = manifolds.scale_kendall2D_tVec(nPt)
est_beta1.Read(
'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/HMG_ScaleKendall_HR_ADOS_New/est_beta1.tVec'
)
est_beta2 = manifolds.scale_kendall2D_tVec(nPt)
est_beta2.Read(
'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/HMG_ScaleKendall_HR_ADOS_New/est_beta2.tVec'
)
est_gamma0 = manifolds.scale_kendall2D_tVec(nPt)
est_gamma0.Read(
'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/HMG_ScaleKendall_HR_ADOS_New/est_gamma0.tVec'
)
est_gamma1 = manifolds.scale_kendall2D_tVec(nPt)