Каков наилучший способ использовать агролы для создания ANOVA на сплит-проекте?

Question

Я пытаюсь запустить ANOVA для данных из эксперимента с разделением графиков, в идеале используя пакет agricolae. Прошло много времени с тех пор, как я взял урок статистики и хотел быть уверенным, что правильно анализирую эти данные, поэтому я провел поиск в Интернете и не смог найти последовательности в том, как люди анализируют свои эксперименты с раздельными графиками. Какой лучший способ для меня это сделать?

Вот заголовок моих данных:

dput(head(rawData))
structure(list(ï..Plot = 2111:2116, Variety = structure(c(5L, 
4L, 3L, 6L, 1L, 2L), .Label = c("Burbank", "Hodag", "Lamoka", 
"Norkotah", "Silverton", "Snowden"), class = "factor"), Rate = c(4L, 
4L, 4L, 4L, 4L, 4L), Rep = c(1L, 1L, 1L, 1L, 1L, 1L), totalTubers = c(594L, 
605L, 656L, 729L, 694L, 548L), totalOzNoCulls = c(2544.18, 2382.07, 
2140.69, 2401.56, 2440.56, 2503.5), totalCWTacNoCulls = c(461.76867, 
432.345705, 388.535235, 435.88314, 442.96164, 454.38525), avgLWratio = c(1.260615419, 
1.287949374, 1.111981583, 1.08647584, 1.350686661, 1.107173509
), Hollow = c(14L, 15L, 22L, 25L, 14L, 13L), Double = c(10L, 
13L, 15L, 22L, 11L, 9L), Knob = c(86L, 80L, 139L, 156L, 77L, 
126L), Researcher = structure(c(1L, 1L, 1L, 1L, 1L, 1L), .Label = "Wang", class = "factor"), 
    CullsPounds = c(1.75, 1.15, 4.7, 1.85, 0.8, 5.55), CullsOz = c(28, 
    18.4, 75.2, 29.6, 12.8, 88.8), totalOz = c(2572.18, 2400.47, 
    2215.89, 2431.16, 2453.36, 2592.3), totalCWTacCulls = c(466.85067, 
    435.685305, 402.184035, 441.25554, 445.28484, 470.50245)), row.names = c(NA, 
6L), class = "data.frame")

Для этих данных весь график равен Rate, разделенный график равен Variety, блок Rep, и для обсуждения здесь мы можем посмотреть totalCWTacNoCulls как ответ.

Любая помощь будет очень признательна! У меня все еще есть опыт переполнения стека, поэтому, если я допустил какие-либо ошибки или поделился своими данными неправильно, пожалуйста, дайте мне знать, и я изменю это. Спасибо!

Answer 1

Вы можете сделать это, используя agricolae пакет следующим образом

library(agricolae)
attach(rawData)
Rate = factor(Rate)
Variety = factor(Variety)
Rep = factor(Rep)
sp.plot(Rep, Rate, Variety, totalCWTacNoCulls)

Использование в соответствии с agricolae Пакет

sp.plot(block, pplot, splot, Y)

где, блок - репликации, pplot - main-Коэффициент сюжета, сплот - это коэффициент подзаговора и переменная ответа Y