123456789101112131415161718192021222324252627282930313233343536 |
- # 3 genetypic groups: Val/Val (12), Val/Met (10), Met/Met (9)
- # Unequal number of subjects across groups
- # 1 effect estimate per subject
- # By default xxi3dMEMA++ can only compare two groups
- # Group comparison: B-A
- 3dMEMA -prefix output.MEMA_2grps -jobs 2 \
- -groups VV VM -missing_data 0 \
- -covariates cov_2grps.1D -covariates_center FD = 0 0 \
- -covariates_model center=same slope=different \
- -HKtest -model_outliers -residual_Z \
- -set VV \
- as data/as+tlrc'[HO_resvnr_mvn#0_Coef]' data/as+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- cg data/cg+tlrc'[HO_resvnr_mvn#0_Coef]' data/cg+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- cm data/cm+tlrc'[HO_resvnr_mvn#0_Coef]' data/cm+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- db data/db+tlrc'[HO_resvnr_mvn#0_Coef]' data/db+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- dh data/dh+tlrc'[HO_resvnr_mvn#0_Coef]' data/dh+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- dt data/dt+tlrc'[HO_resvnr_mvn#0_Coef]' data/dt+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- gm data/gm+tlrc'[HO_resvnr_mvn#0_Coef]' data/gm+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- hb data/hb+tlrc'[HO_resvnr_mvn#0_Coef]' data/hb+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- js data/js+tlrc'[HO_resvnr_mvn#0_Coef]' data/js+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- kb data/kb+tlrc'[HO_resvnr_mvn#0_Coef]' data/kb+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- mg data/mg+tlrc'[HO_resvnr_mvn#0_Coef]' data/mg+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- ph data/ph+tlrc'[HO_resvnr_mvn#0_Coef]' data/ph+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- -set VM \
- pl data/pl+tlrc'[HO_resvnr_mvn#0_Coef]' data/pl+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- ro data/ro+tlrc'[HO_resvnr_mvn#0_Coef]' data/ro+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- te data/te+tlrc'[HO_resvnr_mvn#0_Coef]' data/te+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- tf data/tf+tlrc'[HO_resvnr_mvn#0_Coef]' data/tf+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- zy data/zy+tlrc'[HO_resvnr_mvn#0_Coef]' data/zy+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- ee data/ee+tlrc'[HO_resvnr_mvn#0_Coef]' data/ee+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- jl data/jl+tlrc'[HO_resvnr_mvn#0_Coef]' data/jl+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- jp data/jp+tlrc'[HO_resvnr_mvn#0_Coef]' data/jp+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- kc data/kc+tlrc'[HO_resvnr_mvn#0_Coef]' data/kc+tlrc'[HO_resvnr_mvn#0_Tstat]' \
- rb data/rb+tlrc'[HO_resvnr_mvn#0_Coef]' data/rb+tlrc'[HO_resvnr_mvn#0_Tstat]' \
|