3.MEMA_2grps.txt 2.3 KB

123456789101112131415161718192021222324252627282930313233343536
  1. # 3 genetypic groups: Val/Val (12), Val/Met (10), Met/Met (9)
  2. # Unequal number of subjects across groups
  3. # 1 effect estimate per subject
  4. # By default xxi3dMEMA++ can only compare two groups
  5. # Group comparison: B-A
  6. 3dMEMA -prefix output.MEMA_2grps -jobs 2 \
  7. -groups VV VM -missing_data 0 \
  8. -covariates cov_2grps.1D -covariates_center FD = 0 0 \
  9. -covariates_model center=same slope=different \
  10. -HKtest -model_outliers -residual_Z \
  11. -set VV \
  12. as data/as+tlrc'[HO_resvnr_mvn#0_Coef]' data/as+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  13. cg data/cg+tlrc'[HO_resvnr_mvn#0_Coef]' data/cg+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  14. cm data/cm+tlrc'[HO_resvnr_mvn#0_Coef]' data/cm+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  15. db data/db+tlrc'[HO_resvnr_mvn#0_Coef]' data/db+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  16. dh data/dh+tlrc'[HO_resvnr_mvn#0_Coef]' data/dh+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  17. dt data/dt+tlrc'[HO_resvnr_mvn#0_Coef]' data/dt+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  18. gm data/gm+tlrc'[HO_resvnr_mvn#0_Coef]' data/gm+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  19. hb data/hb+tlrc'[HO_resvnr_mvn#0_Coef]' data/hb+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  20. js data/js+tlrc'[HO_resvnr_mvn#0_Coef]' data/js+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  21. kb data/kb+tlrc'[HO_resvnr_mvn#0_Coef]' data/kb+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  22. mg data/mg+tlrc'[HO_resvnr_mvn#0_Coef]' data/mg+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  23. ph data/ph+tlrc'[HO_resvnr_mvn#0_Coef]' data/ph+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  24. -set VM \
  25. pl data/pl+tlrc'[HO_resvnr_mvn#0_Coef]' data/pl+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  26. ro data/ro+tlrc'[HO_resvnr_mvn#0_Coef]' data/ro+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  27. te data/te+tlrc'[HO_resvnr_mvn#0_Coef]' data/te+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  28. tf data/tf+tlrc'[HO_resvnr_mvn#0_Coef]' data/tf+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  29. zy data/zy+tlrc'[HO_resvnr_mvn#0_Coef]' data/zy+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  30. ee data/ee+tlrc'[HO_resvnr_mvn#0_Coef]' data/ee+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  31. jl data/jl+tlrc'[HO_resvnr_mvn#0_Coef]' data/jl+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  32. jp data/jp+tlrc'[HO_resvnr_mvn#0_Coef]' data/jp+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  33. kc data/kc+tlrc'[HO_resvnr_mvn#0_Coef]' data/kc+tlrc'[HO_resvnr_mvn#0_Tstat]' \
  34. rb data/rb+tlrc'[HO_resvnr_mvn#0_Coef]' data/rb+tlrc'[HO_resvnr_mvn#0_Tstat]' \