{
data(scotvote)
attach(scotvote)
Y <- PerYesParl
X <- cbind(PrivateHousingStarts,CouncilTax,Percentage5to15,
PrimaryPTRatio,PerBirthsOut,PerClaimantFemale)
detach(scotvote)
glmdm.linear(Y,X,num.reps=5000,log=FALSE)
}Run the code above in your browser using DataLab