所有的材料准备就绪,正式开始分析 这里我是用了一个循环的shell脚本进行多个表型的同时分析 #!/bin/sh for i in `ls /public200T/lintao_cau/lintao/workspace/yangtailai/tri_filter_emmax/phenotype/emmax*.txt` do emmax-intel64 -v -d 10 -t /public200T/
emmax.cov.txt 2.4 表型数据 emmax.trait.**.txt emmax中一次只能输入一个性状,文件格式与.tfam类似,需要在R中进行转换。 >tfam<-read.table("emmax_in.tfam",header=F,stringsAsFactors=F)>tr<-read.table("trait.txt",header=T,check.names=F,stringsAsFactors=F)>head(tr)<Trait>TLSAADNR1146.8610.100.69...