ANOVA в R с использованием сводных данных

можно ли запустить ANOVA в r только со средними значениями, стандартным отклонением и n-значением? Вот мой фрейм данных:

q2data.mean <- c(90,85,92,100,102,106)
q2data.sd <- c(9.035613,11.479667,9.760268,7.662572,9.830258,9.111457)
q2data.n <- c(9,9,9,9,9,9)
q2data.frame <- data.frame(q2data.mean,q2data.sq,q2data.n)

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

Любая помощь могла бы быть полезна! :)


person y3trgfhsfgr    schedule 02.10.2014    source источник
comment
Это был первый результат Google для дисперсии R по среднему и стандартному отклонению: r.789695.n4.nabble.com/   -  person thelatemail    schedule 03.10.2014
comment
В интернете тоже ничего не нашел, поэтому и задал вопрос. :)   -  person y3trgfhsfgr    schedule 03.10.2014
comment
Кроме того, поиск в Google вашего точного названия этого вопроса дает ссылку на stats.stackexchange.com/questions/95949/, который указывает на: inside-r.org/packages/cran/rpsychi/docs/ind.oneway.second   -  person thelatemail    schedule 03.10.2014
comment
О, вторая ссылка выглядит многообещающе; Я отправлю ответ, если он работает. :)   -  person y3trgfhsfgr    schedule 03.10.2014


Ответы (1)


Вот и все, используя ind.oneway.second из пакета rspychi:

library(rpsychi)
with(q2data.frame, ind.oneway.second(q2data.mean,q2data.sd,q2data.n) )

#$anova.table
#                SS df     MS     F
#Between (A) 2923.5  5 584.70 6.413
#Within      4376.4 48  91.18      
#Total       7299.9 53   
# etc etc

В качестве несвязанного примечания, ваши данные могут быть переименованы. q2data.frame это data.frame, не надо указывать это в заголовке. Кроме того, нет необходимости указывать q2data.mean внутри q2data.frame — наверняка mean будет достаточно. Это просто означает, что вы получите сложный код, например:

q2data.frame$q2data.mean

когда:

q2$mean

даст вам всю необходимую информацию.

person thelatemail    schedule 02.10.2014
comment
Вы - находка, спасибо! Также: для тех, у кого нет библиотеки rpsychi, перед приведенным выше сценарием был включен следующий код: install.packages(rpsychi) install.packages(gtools) - person y3trgfhsfgr; 03.10.2014
comment
@BrettCochrane - install.packages должен обрабатывать все зависимости - нет необходимости устанавливать gtools отдельно. - person thelatemail; 03.10.2014