# NOT RUN {
# }
# NOT RUN {
your_data=data.frame(mvrnorm(n=20,mu=c(0,0,0,0,0),
Sigma=matrix(c(4,2,2,2,2,
2,4,2,2,2,
2,2,4,2,2,
2,2,2,4,2,
2,2,2,2,4),
nrow=5, ncol=5)))
colnames(your_data)=c("x1","x2","x3","x4","x5")
mod='tau=~x1+x2+x3+x4+x5'
fit=bsem(mod,data=your_data)
bomega(K=5,mod=fit,alpha=3.51,beta=1.75,CI=0.95)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab