if (FALSE) {
## Model-based hierarchical classification of items from simulated data
# Setup
r = 6 # number of items
n = 1000 # sample size
bev = rep(0,r)
k = r/2
multi = rbind(1:(r/2),(r/2+1):r)
L = chol(matrix(c(1,0.6,0.6,1),2,2))
data = matrix(0,n,r)
model = 1
# Create data
Th = matrix(rnorm(2*n),n,2)
for(i in 1:n) for(j in 1:r){
if(j<=r/2){
pc = exp(Th[i,1]-bev[j]); pc = pc/(1+pc)
}else{
pc = exp(Th[i,2]-bev[j]); pc = pc/(1+pc)
}
data[i,j] = runif(1)
Run the code above in your browser using DataLab