# # gamma1 * s t
# est_gamma1 = tangent_slope_arr[ 0 ]
# est_gamma1.Write( '/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/HMG_AutismSpecOnly_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_AutismSpecOnly_ADOS_New/est_gamma2.tVec')
# Read
est_beta0 = manifolds.kendall2D(nPt)
est_beta0.Read(
'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/HMG_AutismSpecOnly_ADOS_New/est_beta0.pt'
)
est_beta1 = manifolds.kendall2D_tVec(nPt)
est_beta1.Read(
'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/HMG_AutismSpecOnly_ADOS_New/est_beta1.tVec'
)
est_beta2 = manifolds.kendall2D_tVec(nPt)
est_beta2.Read(
'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/HMG_AutismSpecOnly_ADOS_New/est_beta2.tVec'
)
est_gamma0 = manifolds.kendall2D_tVec(nPt)
est_gamma0.Read(
'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/HMG_AutismSpecOnly_ADOS_New/est_gamma0.tVec'
)
est_gamma1 = manifolds.kendall2D_tVec(nPt)
mean_shape.Read(
'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/MLMG_LRHR_Autism/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]
beta_i = manifolds.kendall2D_tVec(nPt)
for k in range(2):
for j in range(nPt):
beta_i.tVector[k, j] = (est_beta1.tVector[k, j] * s_i +
est_beta2.tVector[k, j] * hr_i +
est_beta3.tVector[k, j] * 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)
gamma3_p1_i = est_beta0.ParallelTranslateAtoB(est_beta0, p1_i, est_gamma3)
gamma_i = manifolds.kendall2D_tVec(nPt)
# est_gamma5.Write( '/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/MLMG/est_gamma5.tVec')
# # gamma6 * healthy t
# est_gamma6 = tangent_slope_arr[ 5 ]
# est_gamma6.Write( '/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/MLMG/est_gamma6.tVec')
print("# of Subjects")
print(len(shape_list_a))
# Read
est_beta0 = manifolds.kendall2D(nPt)
est_beta0.Read(
'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/MLMG/est_beta0.pt'
)
est_beta1 = manifolds.kendall2D_tVec(nPt)
est_beta1.Read(
'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/MLMG/est_beta1.tVec'
)
est_beta2 = manifolds.kendall2D_tVec(nPt)
est_beta2.Read(
'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/MLMG/est_beta2.tVec'
)
est_beta3 = manifolds.kendall2D_tVec(nPt)
est_beta3.Read(
'/media/shong/IntHard1/Projects/4DShapeAnalysis/Data/ACE_IBIS_CC/Data/MLMG/est_beta3.tVec'
)
est_beta4 = manifolds.kendall2D_tVec(nPt)