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)