## Not run:
#
# # create list of parameters
# Item1 <- c(1,-2,c(-0.5,0.3,0.2),c(-0.5,-0.3,0.8))
# names(Item1) <- c("a","b",paste("zeta1",1:3,sep=""),paste("lamb",1:3,sep=""))
#
# Item2 <- c(1,-1,c(-0.5,-0.3,0.8),c(-0.5,0.3,0.2))
# names(Item2) <- c("a","b",paste("zeta1",1:3,sep=""),paste("lamb",1:3,sep=""))
#
# Item3 <- c(1,0,c(-0.5,-0.3,0.8),c(-0.5,0.3,0.2))
# names(Item3) <- c("a","b",paste("zeta1",1:3,sep=""),paste("lamb",1:3,sep=""))
#
# Item4 <- c(1,1,c(-0.5,-0.3,0.8),c(-0.5,0.3,0.2))
# names(Item4) <- c("a","b",paste("zeta1",1:3,sep=""),paste("lamb",1:3,sep=""))
#
# Item5 <- c(1,2,c(-0.5,-0.3,0.8),c(-0.5,0.3,0.2))
# names(Item5) <- c("a","b",paste("zeta1",1:3,sep=""),paste("lamb",1:3,sep=""))
#
# ParList <- list(Item1=Item1,Item2=Item2,Item3=Item3,Item4=Item4,Item5=Item5)
#
# # simulate
# perp1 <- rnorm(1000,0,1)
# simdat1 <- NLM.sim(ParList,perp1)
#
# # reshape
# reshOBJ <- reshMG(simdat1,items=1:5,groups=NA,correct=rep(0,5),design="nodif",echo=TRUE,TYPE="NLM")
#
# # estimate a nested logit model
# res.nlm <- nelm(reshOBJ=reshOBJ)
#
# # plot the estimated data
# plot(res.nlm)
#
#
# ## End(Not run)
Run the code above in your browser using DataLab