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