我有一个简单的方差分析模拟功能。然而,我想知道为什么我的
replicate()
命令不生成数字表?
fun <- function(eta.sq = .25, groups = 4, n = 10, bet.var = 10){
with.var = bet.var*(1/eta.sq - 1)
N = groups*n
sim.means = rnorm(n = groups, mean = 0, sd = sqrt(bet.var))
sim.data = data.frame(group = gl(groups, 1, length = N),
response = rnorm(N, sim.means, sqrt(with.var)))
sim.anova = anova(aov(response ~ group, sim.data))
}
# Problem part:
t(replicate(5, fun())) # Why I don't get a numeric table?
Df Sum Sq Mean Sq F value Pr(>F)
[1,] Integer,2 Numeric,2 Numeric,2 Numeric,2 Numeric,2
[2,] Integer,2 Numeric,2 Numeric,2 Numeric,2 Numeric,2
[3,] Integer,2 Numeric,2 Numeric,2 Numeric,2 Numeric,2
[4,] Integer,2 Numeric,2 Numeric,2 Numeric,2 Numeric,2
[5,] Integer,2 Numeric,2 Numeric,2 Numeric,2 Numeric,2
