Блог

Главная - Вопросы по программированию - Выполнение одностороннего ANOVA с использованием rpsychi (R) и последующих последовательных парных сравнений

Выполнение одностороннего ANOVA с использованием rpsychi (R) и последующих последовательных парных сравнений

#r #anova

#r #anova

Вопрос:

У меня есть 4 выборки, представленные их средними значениями и стандартными отклонениями. Если я считаю размер выборки для всех из них равным 100, как я могу выполнить односторонний ANOVA, используя библиотеку rpsychi в R?

Все выборки нормально распределены и независимы.

Это их значения:

S1: среднее значение 1 = 5,5, sd1 = 0,6

S2: среднее значение 2 = 5,6, sd2 = 0,5

S3: среднее значение 3 = 5,9, sd3 = 0,5

S4: среднее значение 4 = 6,1, sd4 = 0,6

Как я могу впоследствии выполнить последующее парное сравнение (Tukey-HSD), также используя R??

Комментарии:

1. Вы смотрели справку в rpsychi? Специально для ind.oneway.second ? Кроме того, этот вопрос, вероятно, слишком похож на ваш другой, который был задан и на который был дан ответ. stats.stackexchange.com/q/102703/601 . Ваш расчет статистики теста Tukey HSD прост, и его можно посмотреть в Интернете. Все, что вам нужно знать в R, это то, что p-значение из изученного распределения диапазона ( q ) может быть получено с ptukey помощью .

2. Извините за задержку. Джон, да, я проверил ind.oneway.second в PDF-файле справки pfsychi. Тем не менее, я несколько раз безуспешно пытался ввести свои данные. Я не знаю, как написать идеальный код для выполнения этого ANOVA. Итак, я спрашиваю, знает ли кто-нибудь точный код моего примера, и таким образом, изучите его. Спасибо

3. Добавьте то, что вы пробовали тогда. Сайт предназначен не только для того, чтобы дать ответ, но и для того, чтобы помочь людям учиться. Если мы не знаем, как вы терпите неудачу, вы не научитесь. И если вы не покажете им свои усилия, люди часто неохотно будут помогать. Кроме того, если речь идет только о вводе команд R, то это неправильный форум, и он должен быть в StackOverflow с тегом R.

Ответ №1:

Недавно я обнаружил сообщение на другом форуме с некоторыми моими собственными проблемами: http://www.wenda.io/questions/2570717/anova-in-r-using-summary-data.html

Это частично решило то, что мне было нужно. Применение кода, опубликованного в этом интернете, к моим данным:

 mean <- c(5.5,5.6,5.9,6.1)
sd <- c(0.6,0.5,0.5,0.6)
n <- c(100,100,100,100)
q2data.frame <- data.frame(mean,sd,n)
library(rpsychi)
with(q2data.frame, ind.oneway.second(mean,sd,n))

Я получаю следующую таблицу ANOVA:

 $anova.table
                SS  df    MS      F
Between (A)  22.75   3 7.583 24.863
Within      120.78 396 0.305       
Total       143.53 399             

$omnibus.es
      etasq etasq.lower etasq.upper 
      0.159       0.094       0.219 

$raw.contrasts
    mean.diff  lower  upper   std
1-2      -0.1 -0.254  0.054 0.078
1-3      -0.4 -0.554 -0.246 0.078
1-4      -0.6 -0.754 -0.446 0.078
2-3      -0.3 -0.454 -0.146 0.078
2-4      -0.5 -0.654 -0.346 0.078
3-4      -0.2 -0.354 -0.046 0.078

$standardized.contrasts
        es  lower  upper   std
1-2 -0.181 -0.459  0.097 0.141
1-3 -0.724 -1.002 -0.446 0.141
1-4 -1.086 -1.364 -0.808 0.141
2-3 -0.543 -0.821 -0.265 0.141
2-4 -0.905 -1.183 -0.627 0.141
3-4 -0.362 -0.640 -0.084 0.141

$power
 small medium  large 
 0.363  0.993  1.000

Кто-нибудь знает, как я могу получить p-значение анализа ANOVA? Более того, как я могу выполнить TukeyHSD?

Ответ №2:

Используйте библиотеку (asbio). Это даст вам t (s), средние различия, верхнюю и нижнюю границы, HLSD Diff Нижнее верхнее решение Adj. решение о значении p и скорректированное значение p.

  #pairwise comparisons
 library(asbio)
 bonf <- pairw.anova(data$dv, data$group, method="tukey") #also try "bonf" or "lsd"
 print(tukey)
 #plot(tukey) #can plot the CFs
Метки: