I have a problem with "Error in comp[[i]] : subscript out of bounds"
I want to indicate compact letter display in my analysis but it did not work. I have tried to find some solution in the internet but it failed.
Could anyone help me?
this is my data
and this is the script
library(datasets) library(ggplot2) library(multcompView) library(dplyr) library(datasets) library(tidyverse) library(multcomp) Data = read.csv("data.csv", h= TRUE) qplot(x = species, y = Ratio, geom = "point", data = Data) + facet_grid(.~Strain) # creating a variable as factor for the ANOVA Data$Strain <- as.factor(Data$Strain) Data$species <- as.factor(Data$species) str(Data) # analysis of variance anova <- aov(Ratio ~ Strain*factor(species), data = Data) summary(anova) # table with factors, means and standard deviation data_summary <- group_by(Data, Strain, species) %>% summarise(mean=mean(Ratio), sd=sd(Ratio)) %>% arrange(desc(mean)) print(data_summary) # Tukey's test tukey <- TukeyHSD(anova) print(tukey) # compact letter display coba = multcompLetters4(anova, tukey) print(coba) # creating the compact letter display tukey.cld <- multcompLetters4(anova, tukey) print(tukey.cld)
It will be easier for me if I can distinguish the significant of data based on the letters.
I really need help here as soon as possible.