Как я могу сделать HSD-тест Тьюки на сплит-дизайн? - PullRequest
Новая
       96

Как я могу сделать HSD-тест Тьюки на сплит-дизайн?

0 голосов
/ 21 октября 2019

Я пытаюсь проанализировать данные из дизайна разделенного графика. Я думаю, что у меня недостаточно опыта, чтобы самому написать HSD-тест Тьюки, хотя я смог сделать это для других экспериментальных проектов. Я искал переполнение стека, а также в Google, но мне не удалось получить рабочий набор кода.

Вот глава моего набора данных. Я могу поделиться больше, если это необходимо: 

 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")

Мне интересно увидеть разные категории букв для разных курсов. Я хотел бы поместить их на график, который имеет гистограмму для каждого сорта (всего шесть) с totalCWTacCulls на оси y и 5 курсами на каждой оси x

Вот анова, которую я используюдля этих данных: 

require(agricole)
model <- with(rawData, sp.plot(Rep, Rate, Variety, totalCWTacCulls))

Я пытаюсь использовать функцию HSD.test, но у меня недостаточно четкого понимания того, что я делаю, чтобы это работало. Вот моя текущая попытка: 

HSD <- HSD.test(model, "Rate", alpha = 0.05, group=TRUE)
print(HSD$groups)

Однако я получаю ответ о том, что мой MSerror отсутствует. У меня не было этой проблемы, когда я использовал этот метод для анализа своих экспериментов ранее, но это первый анализ, который я сделал для дизайна с разделением графиков.

Любая помощь или направление будут высоко оценены. Спасибо!

...