library(nnet)
library(foreign)
ml <- read.dta("http://www.ats.ucla.edu/stat/data/hsbdemo.dta")
ml$prog2 <- relevel(ml$prog, ref = "academic")
# Multinomial model
m <- multinom(prog2 ~ ses + write, data = ml)
summary(m)
logitgof(ml$prog2, fitted(m))
logitgof(ml$prog2, fitted(m))$exp
logitgof(ml$prog2, fitted(m), g = 5)$exp
logitgof(ml$prog2, fitted(m), g = 5)
# Binary logistic model
# Outcome with factor names
m2 <- glm(female ~ ses + write, data=ml, family="binomial")
logitgof(ml$female, fitted(m2))
# Outcome with 1/0 coding
ml$f2 <- 1
ml[ml$female=="male", ]$f2 <- 0
m3 <- glm(f2 ~ ses + write, data=ml, family="binomial")
logitgof(ml$f2, fitted(m3))
Run the code above in your browser using DataLab