Я пытаюсь рассчитать размеры эффекта, но все, что я могу найти в Интернете, касается только величины эффекта разницы между двумя группами.Что меня интересует, так это размер эффекта различия между групповыми различиями.
Не очень четкое объяснение, я знаю ... Я буду использовать набор данных IRIS в качестве примера:
Iхотел бы знать, если
разница между лепестком. Длина видов setosa и лепесток. Длина видов лишай (Δsetosa-versicolor)
больше или меньше, чем
разница между лепестками. Длина видов setosa и лепестками. Длина видов viriginica (Δsetosa-virginica).
Поэтому я хотел бы рассчитать величину эффекта между Δsetosa-versicolor и Δsetosa-virginica.
Для моего набора данных я уже выполнил GLM и последующий пост-специальный тест с использованием пользовательских кальстатов, чтобы оценить, действительно ли существует различие между группами (что соответствует значительному различию между сетозой и разноцветной длиной лепестка, а также между сетозой идлина лепестка virginica в наборе данных IRIS).
Некоторые попытки, которые я предпринял:
используйте функцию dist () для вычисления матрицы разностей Petal.Length между двумя видами.Но я теряю информацию о видах, поэтому не очень полезна.Более того, я не уверен, что я получаю ... все возможные парные сравнения?Только те, между группами?последующие сравнения?
as.numeric(dist(iris$Petal.Length))
abs(apply(combn(iris$Petal.Length,2), 2, diff))
Использовать функцию mes () из пакета compute.es.Но таким образом я могу получить только величину эффекта, например, между сетозой и версиколором, а не то, что мне интересно.
mes(mean(iris$Petal.Length[iris$Species=="setosa"],na.rm=T), mean(iris$Petal.Length[iris$Species=="versicolor"],na.rm=T),
sd(iris$Petal.Length[iris$Species=="setosa"],na.rm=T), sd(iris$Petal.Length[iris$Species=="versicolor"],na.rm=T),
length(iris$Petal.Length[iris$Species=="setosa"]), length(iris$Petal.Length[iris$Species=="versicolor"]))
Что я могу сделать?Я действительно заблокирован здесь
Большое спасибо за вашу помощь!=) * * Тысяча двадцать-семь