Displaying 1 result from an estimated 1 matches for "npergroup".
Did you mean:
apegroup
2010 Mar 22
0
using lmer weights argument to represent heteroskedasticity
...any advice.
J.R.
###############################################################################
## EXAMPLE: generate balanced nested data with heteroskedastic errors of
## variance 0.5, 1, 2, or 4
###############################################################################
ngroup <- 50
npergroup <- 20
n <- ngroup*npergroup
set.seed(5046)
d <- data.frame(unitid = 1:n, groupid = rep(1:ngroup,
each=npergroup), verror = sample(c(0.5, 1, 2, 4), size=n, replace=T),
x = 0.1*rnorm(n))
groupeffx <- data.frame(groupid = 1:ngroup, theta = rnorm(ngroup, sd = 0.25))
d <...