注意表型首行使用标签<Trait>
脚本mrmlm.r
#!/public/home/caisl/biosoft/miniconda3/envs/r4.1/bin/R
library(‘mrMLM‘)
mrMLM(fileGen=‘184.hmp.txt‘,filePhe=‘pheno‘,Genformat=‘Hmp‘,method=c("mrMLM","FASTmrMLM","FASTmrEMMA"),trait=1:8,CriLOD=3,DrawPlot=TRUE,Plotformat=‘jpeg‘,dir="/public/home/caisl/lee/gwas/cold184/mrmlm/result/")
#首行为R的路径
qsub脚本
#PBS -N mrMLM
#PBS -l nodes=1:ppn=8
#PBS -q middle
cd $PBS_O_WORKDIR
NP=`cat $PBS_NODEFILE|wc -l`
/public/home/caisl/biosoft/miniconda3/envs/r4.1/bin/Rscript mrmlm.r