glmdm (version 2.01)

glmdm.probit: glmdm.probit

Description

R code for simulation of GLMDM

Usage

glmdm.probit(Y, X, num.reps=1000, a1=3, b1=2, d=0.25, MM=15,VV=30,...)

Arguments

Y

Y

X

X

num.reps

num.reps

a1

a1

b1

b1

d

d

MM

MM

VV

VV

...

..

Examples

{
data(ssas)
attach(ssas)
Y <- ssas[,18]
X <- cbind(ssas$househld,ssas$rsex,ssas$rage,ssas$relgsums,ssas$ptyallgs,ssas$idlosem,
                 ssas$marrmus,ssas$ukintnat,ssas$natinnat,ssas$voiceuk3,ssas$nhssat,ssas$hincdif2,
                 ssas$unionsa,ssas$whrbrn,ssas$hedqual2)
detach(ssas)
glmdm.out <- glmdm.probit(Y,X,num.reps=10000)


data(asia)
Y <- asia[,1]
X <- as.matrix(asia[,-1]) 
glmdm.out <- glmdm.probit(Y,X,num.reps=10000,log=TRUE) 
}