Displaying 1 result from an estimated 1 matches for "numsimulations".
2012 Apr 26
2
Lambert (1992) simulation
...am trying to recreate Table 2. But my estimates for Gamma
are not working out. Any ideas why? Please cc me on a reply!
Thanks,
Chris
## ZIP simulations based on Lambert (1992)'s conditions
# Empty workspace
rm(list=ls())
# Load zero-inflation package
library(pscl)
# Simulation conditions #3
NumSimulations=2000
X=seq(from=0,to=1,length.out=N)
Model.Matrix=cbind(rep(1,length(X)),X)
Gamma=c(-1.5,2)
Beta=c(1.5,-2)
P=exp(Model.Matrix%*%Gamma)/exp(1+Model.Matrix%*%Gamma)
Lambda=exp(Model.Matrix%*%Beta)
CoefSimulations=matrix(nrow=NumSimulations,ncol=2*dim(Model.Matrix)[2])
for(i in 1 : NumSimulations){
L...