R 循环的写输出
我在数千个文件中运行一个循环来执行Pagel LambdaR 循环的写输出,r,loops,R,Loops,我在数千个文件中运行一个循环来执行Pagel Lambda files <- list.files(pattern=".txt") column_names <- data.frame(files = "files", V1 = "V1", Lambda = "lamda", LogL = "logL", LogL10 = "logL10", pvalue = "P") write.table(column_names, file = "output.csv", row.names
files <- list.files(pattern=".txt")
column_names <- data.frame(files = "files", V1 = "V1", Lambda = "lamda", LogL = "logL", LogL10 = "logL10", pvalue = "P")
write.table(column_names, file = "output.csv", row.names = FALSE,
append = FALSE, col.names = FALSE, sep = ", ", quote = TRUE)
for(i in files){
graph <- read.table(i,header=F)
row.names(graph)<-graph[,1]
trait<-graph[,-1]
names(trait)<-row.names(graph)
result <- data.frame(t(c(files[i], phylosig(mytree, trait, method="lambda", test=TRUE, nsim=999))))
write.table(result, file = "output.csv", row.names = FALSE,
append = TRUE, col.names = FALSE, sep = ", ")
}
我需要这些值出现在输出文件和输入文件的名称中(我还没有在循环中设置它,我不知道怎么做)
因此,我希望输出如下所示:
[1] "some data in x given as 'NA', dropping corresponding species from tree"
V1 lambda logL logL0 P
1 EOG7B0F7J.txt 6.610696e-05 -59.62533 -59.62509 1
FileA V1 lambda logL logL0 P
FileB V1 lambda logL logL0 P
FileC V1 lambda logL logL0 P
或者类似的事情——如果有人能帮我解决这个问题,我会非常高兴
非常感谢我想出来了,如果有人感兴趣,下面是答案:
column_names <- data.frame(i = "files", Lambda = "lamda", LogL = "logL", LogL10 = "logL10", pvalue = "P")
write.table(column_names, file = "output.csv", row.names = FALSE,
append = FALSE, col.names = FALSE, sep = ", ", quote = TRUE)
for(i in files){
graph <- read.table(i,header=F)
row.names(graph)<-graph[,1]
trait<-graph[,-1]
names(trait)<-row.names(graph)
result <- data.frame(t(c(i, phylosig(mytree, trait, method="lambda", test=TRUE, nsim=999))))
my.df <- data.frame(lapply(result, as.character), stringsAsFactors=FALSE)
write.table(my.df, file = "output.csv", row.names = FALSE, append = TRUE, col.names = FALSE, sep = ", ")
}
column\u name您的data.frame
行中有一个输入错误。是的,我已经纠正了!谢谢