Displaying 1 result from an estimated 1 matches for "w_gik".
Did you mean:
v_gik
2005 Dec 01
1
Simulate Correlated data from complex sample
..., sqrt(.99))
# Organize into a dataframe
data <- data.frame(schid = rep(1:J, rps), stuid = 1:N, cbind(v_gk,
v_gik), psi = v_gk + v_gik + a_g[1])
# Now create theta
B_g <- .95 # This is correlation between within-grade trait and
vertical trait
w_gk <- 0 # fixed at zero for now
data$w_gik <-rnorm(N, 0, sqrt(.0975))
data$theta <- (B_g * data$psi) + w_gk + data$w_gik
[[alternative HTML version deleted]]