## Not run:
# data(simICCdata)
# require(lme4)
#
# #computing icca
# vy = var(simICCdata$Y)
# lmm0 <- lmer(Y ~ (1|l2id),data=simICCdata,REML=F)
# VarCorr(lmm0)$l2id[1,1]/vy
#
# #Estimating random slopes model
# lmm1 <- lmer(Y~I(X1-m_X1)+I(X2-m_X2) +(I(X1-m_X1)+I(X2-m_X2)|l2id),data=simICCdata2,REML=F)
# X = model.matrix(lmm1)
# p=ncol(X)
# T1 = VarCorr(lmm1) $l2id[1:p,1:p]
# #computing iccb
# #Notice '+1' because icc_beta assumes l2ids are from 1 to 30.
# icc_beta(X,simICCdata2$l2id+1,T1,vy)$rho_beta
# ## End(Not run)
Run the code above in your browser using DataLab