Я все еще довольно новый в R и новый в stackoverflow.
Я хочу выполнить групповые post-hoc тесты (такие как Tukey, SNK) после выполнения ANOVA для различных групп в кадре данных без поднабора данных и повторенияодно и то же несколько раз.Лучший способ отобразить результаты - это снова кадр данных с переменными «by_group», «group», «letter».Чтобы отобразить букву, я использую пакет Agricola
. Я пробовал разные вещи, такие как pipe operator или group_by, но все не удалось.Я надеюсь, что стало ясно, чего я хочу, и кто-нибудь может мне помочь.
# create dataframe with 10 diffrent experiments (that should analysed seperately)
Exp.No<-10 #number of experiments
s.size<- 6 # sample size per experiment play with it
k<- 5 # Number of categories
pop.mean<-370 ;pop.sd=38 #define population parameters
#sample data form population
dv<-rnorm(n=Exp.No*s.size*k,mean=pop.mean,sd=pop.sd)
#construct dataframe for the experiments
my.df<- data.frame(expand.grid(k=as.factor(1:k),rep=(1:s.size),Exp=(1:Exp.No)),y=dv)
View(my.df)
#perform ANOVA for each experiment and evaluate result
#define by group create a "grouped dataframe"
by.Exp<-group_by(my.df,Exp);str(by.Exp)
anova.by.exp<-by.Exp%>% do(tidy(Anova(aov(y~k,data=.))))#%>% filter(term=="k")
#for a single experiment I could use the aov model and perform a post hoc test like this
modExp1<- aov(y~k,data=subset(my.df, Exp==1))
library(agricolae)
lines<-HSD.test(modExp1,"k");lines$groups
#I would like some thing like this but with an additional variable for the Experiment