glmdm (version 2.60)

glmdm: Generalized Linear Mixed Dirichlet Model

Description

R code for simulation of GLMDM

Usage

glmdm(formula, family=gaussian, data, num.reps=1000, a1=3, b1=2, d=0.25, MM=15, VV=30, ...)

Arguments

formula
a symbolic description of the model to be fit.
family
a descreption of the error distribution and link function to be used in the model.
data
an optional data frame, list or environment containing the variables in the model.
num.reps
num.reps
a1
a1
b1
b1
d
d
MM
MM
VV
VV
...
..

Examples

## Not run: 
# 
# data(scotvote)
# glmdm.linear.out <- glmdm(PerYesParl ~ PrivateHousingStarts + CouncilTax 
# + Percentage5to15 + PrimaryPTRatio + PerBirthsOut + PerClaimantFemale, 
# data=scotvote, num.reps=5000)
# 
# data(ssas)
# glmdm.probit.ssas <- glmdm(scotpar2 ~ househld + rsex + rage + relgsums + ptyallgs
# + idlosem + marrmus + ukintnat + natinnat + voiceuk3 + nhssat, data=ssas, 
# family=binomial(link="probit"), num.reps=10000, log=TRUE)
# 
# data(asia)
# glmdm.probit.asia <- glmdm(ATT ~ DEM + FED + SYS + AUT, data=asia, 
# family=binomial(link="probit"), num.reps=10000, log=TRUE)
# 
# ## End(Not run)