График взаимодействия - ошибка кодирования для двусторонней ановы в R

Я пытаюсь проанализировать и представить результаты соотношения численности брюхоногих моллюсков к высоте на берегу и видов в двухфакторном дисперсионном анализе.

Это набор данных http://dropcanvas.com/ercb6.

Данные были преобразованы sqrt, и тест ANOVA прошел успешно. В приведенном ниже примере кода я смог построить график взаимодействия обычных данных Gastropod. Однако я хочу создать график взаимодействия с использованием преобразованных данных, но я не знаю, как изменить код, чтобы это разрешить. Предложения?

Gastropods = read.csv(file = "MaroubraZones.csv", header = TRUE)
boxplot(Abundance ~ Zone*Species,data = Gastropods, names = c("A.high", "A.mid", "A.low", "C.high", "C.mid", "C.low", "N.high", "N.mid", "N.low"))
Gastropods.ANOVA = aov(Abundance ~ Zone * Species, data = Gastropods)
hist(Gastropods.ANOVA$residuals)
plot(Gastropods.ANOVA)
summary(Gastropods.ANOVA)
Gastropods$sqrtAbundance = sqrt(Gastropods$Abundance +1)
Gastropods.aov = aov(Gastropods$sqrtAbundance ~ Zone + Species + Zone:Species, data = Gastropods)
summary(Gastropods.aov)

interaction.plot(Gastropods$Zone, Gastropods$Species, Gastropods$Abundance, main= "Gastropod Interaction Plot", xlab = "Gastropod Zone", ylab= "Mean of Gastropod Abundance", legend = TRUE)

r interaction transformation anova sqrt
person JoshM8    schedule 01.06.2015    source источник
comment
Это не работает? interaction.plot(Gastropods$Zone, Gastropods$Species, Gastropods$sqrtAbundance, main= "Gastropod Interaction Plot", xlab = "Gastropod Zone", ylab= "Mean of SQRT Gastropod Abundance", legend = TRUE)   -  person Robert    schedule 01.06.2015
comment
фейспалм Вау. Клянусь, я попробовал это и получил сообщение об ошибке... Большое спасибо за помощь с моим крайним нубизмом. Кроме того, мне было интересно, не могли бы вы объяснить, как маркировать полосы ошибок на гистограмме? Я искал везде и ничего не придумал. Я полный новичок в этом. Ссылка на соответствующий источник даже была бы фантастической. Еще раз спасибо.   -  person JoshM8    schedule 01.06.2015
comment
Используйте yaxt или xaxt =n внутри графика или гистограммы. Whit axis() вы можете пометить ось как хотите. Вводные материалы могут быть такими: introductoryr.co.uk/R_Resources_for_Beginners.html.   -  person Robert    schedule 01.06.2015
comment
Я был здесь всего день или два, но сообщество было потрясающим. Спасибо за помощь и ссылку. Мне официально разрешено голосовать, поэтому я сделал это для вас. Спасибо еще раз.   -  person JoshM8    schedule 01.06.2015

Ответы (1)


arrow_upward
1
arrow_downward

Просто чтобы закрыть вопрос:

interaction.plot(Gastropods$Zone, Gastropods$Species, Gastropods$sqrtAbundance, main= "Gastropod Interaction Plot", xlab = "Gastropod Zone", ylab= "Mean of SQRT Gastropod Abundance", legend = TRUE)
person Robert    schedule 01.06.2015