# "Simulate" data following random and fixed effect specifications # Chris Adolph # June 2011 rm(list=ls()) library(simcf) library(tile) library(RColorBrewer) set.seed(123456) # "Random effects" model # concept: grade vs hours of study on a set of assignments # Create covariates n <- 100 t <- 5 beta = 0.75 sigmaAlpha <- 0.7 sigma <- 0.2 rXm <- c(-0.7,.7) sdX <- 0.3 collist <- brewer.pal(9,"Set1") while (length(collist)