glmdm (version 1.02)

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(househld,rsex,rage,relgsums,ptyallgs,idlosem,marrmus,ukintnat,natinnat,voiceuk3,nhssat,hincdif2,unionsa,whrbrn,hedqual2)
detach(ssas)
glmdm.out <- glmdm.probit(Y,X,num.reps=10000)


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