R软件的方差分析，没有原始数据怎么办？

已知均数、标准差、样本量，未知原始数据。

用rpsychi包实现，ind.oneway.second是关键命令，缺陷是不能提供p值。

## 示例1：

install.packages(“rpsychi”)#安装rpsychi包

library(rpsychi)#加载rpsychi包

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)#把均数标准差样本量整合到一个名为q2data.frame的数据框中

with(q2data.frame, ind.oneway.second(mean,sd,n,sig.level = 0.05))

## 原始示例2

##Kline (2004) Table 6.3

dat <- data.frame(y = c(9,12,13,15,16,

8,12,11,10,14,

10,11,13,11,15),

x = rep(factor(c(“a”,”b”,”c”)), each=5)

)

##contrast 1: a – c, contrast 2: 1/2(a + c) – b

my.cont <- matrix(c(1,0,-1,1/2,-1,1/2), ncol=3, nrow=2, byrow=TRUE)

ind.oneway.second(m = tapply(dat$y, dat$x, mean),

sd = tapply(dat$y, dat$x, sd),

n= tapply(dat$y, dat$x, length))

ind.oneway.second(m = tapply(dat$y, dat$x, mean),

sd = tapply(dat$y, dat$x, sd),

n= tapply(dat$y, dat$x, length),

contr = my.cont)