表型按列分文件forGAPIT

myFile <- choose.files()
setwd(dirname(myFile))
print(getwd())
readTabler <- read.table(myFile,header=T,sep=",")
for (i in 2:ncol(readTabler)) {  
  dir.create(names(readTabler)[i], showWarnings = FALSE)  # 创建目录,避免警告  
  temp <- cbind(readTabler[1], readTabler[i])  # 合并第一列和当前列  
  write.table(temp, paste0(names(readTabler)[i], "/GAPIT.Pheno.txt"), quote = FALSE, sep = "\t", row.names = FALSE)  # 写入文件  
}  

评论

发表评论

了解 数据控|突破是我们的每一步 的更多信息

立即订阅以继续阅读并访问完整档案。

继续阅读