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) # 写入文件
}
发表评论