0
我寫了一個循環,在該循環中遍歷給定.csv的列並運行anova和posthoc測試。然後我將每個結果合併到一個數據框中並將其導出到一個.csv文件。但是,我無法獲得rbind()來構建我的data.frame。對此有何幫助?下面是該腳本:rbind()只迭代最後3個anova結果?
setwd("~/School/Lab/mice/sugar_study_2015/MG-RAST and Metagenassist/Trimmed/R. CSV")
#Save your Datasheet into variable X
x <- read.csv("T0_B_Class_Anova.csv")
x = x[1:9,]
x[is.na(x)] <- 0
DF.Anova <- data.frame()
DF.Tukey <- data.frame()
#Counts through the columns
for(i in 2:(ncol(x)-1)){
columns <- names(x[i])
##Runs an ANOVA - 'Group' being a grouping factor
anovaresult <- anova(aov(x[,i]~Group,data=x))
DF.Anova <- rbind(DF.Anova, anovaresult)
##fix anova into data frame
Famall = colnames(x)
Famall = as.data.frame(Famall)
Famall = Famall[2:83,]
Famall = as.data.frame(Famall)
DFanovanames = rep(Famall, each = 2)
DFanovanames = as.data.frame(DFanovanames)
#install.packages("tidyr")
library(tidyr)
anovanames = data.frame(Names=unlist(DFanovanames, use.names = FALSE))
o.anovanames = dplyr::arrange(anovanames, Names)
###dont forget to change this**************************
finalanova_BFT0 = cbind(rn = rownames(DF.Anova), DF.Anova, o.anovanames)
##Runs Tukeys Post-hoc test on Anova
posthocresult <- TukeyHSD(aov(x[,i]~Group,data=x))
DF.Tukey <- rbind(DF.Tukey, posthocresult$Group)
##fix tukey into data frame
Famname = colnames(x)
Famname = as.data.frame(Famname)
Famname = Famname[2:83,]
Famname = as.data.frame(Famname)
DFposthocnames = rep(Famname, each = 3)
DFposthocnames = data.frame(DFposthocnames)
#install.packages("tidyr")
library(tidyr)
library(dplyr)
posthocnames = data.frame(Names=unlist(DFposthocnames, use.names = FALSE))
o.posthocnames = dplyr::arrange(posthocnames, Names)
###dont forget to change this****************************
finalposthoc_BFT0 = cbind(rn = rownames(DF.Tukey), DF.Tukey, o.posthocnames)
##Prints posthoc results into txt file
print(columns)
print(anovaresult)
print(posthocresult)
}
write.csv(finalanova_BFT0, file="testfinalanova_BCT0")
write.csv(finalposthoc_BFT0, file="finalposthoc_BCT0")
你可以找到示例.csv here
你確定您的數據是否正確讀取?你的'setwd'調用看起來不正確 – tbradley
如果你提供了一個[最小可重現的例子],回答你的問題會更容易(https://stackoverflow.com/questions/5963269/how-to-make-a-great-r-可重現的例子),而不是鏈接到數據。 –
我已編輯數據集以包含3列。我不想從腳本中刪除任何內容,因爲我不確定它在哪裏出錯。我盯着這件事一兩個小時,無法弄清楚。 – Haley