Post date: May 13, 2014 2:15:36 AM
I am running the cross-validation analyses for survival and adult weight for the combined experiments, Ac, Ms, GLA, SLA, and all. The conditions are as described previously for the linear model, here is an example for weight:
cd /home/A01963476/projects/lycaeides_hostplant/melGemma/
sleep 5
gemma -g geno_combSla.txt -p pheno_trainWgtSla.txt -bslmm 1 -n 10 -o wgtTrain_combSlaSet10Ch1 -rpace 40 -w 1000000 -s 5000000
gemma -g geno_combSla.txt -p pheno_trainWgtSla.txt -n 10 -gk 1 -o wgtTrain_combSlaSet10Ch1
gemma -g geno_combSla.txt -p pheno_trainWgtSla.txt -n 10 -epm ./output/wgtTrain_combSlaSet10Ch1.param.txt -emu ./output/wgtTrain_combSlaSet10Ch1.log.txt -ebv ./output/wgtTrain_combSlaSet10Ch1.bv.txt -k ./output/wgtTrain_combSlaSet10Ch1.cXX.txt -predict 1 -o wgtTrain_combSlaSet10Ch1k
I need to do the same thing for family-level cross validation.