if (FALSE) {
data(italy2006)
election <- set.data(italy2006 , shape="wide", choice="vote", varying=c(5:14), sep="_")
m <- mlogit(vote~prox+partyID | gov_perf+sex+age+education, election, reflevel = "UL")
true.pos <- list(FI=7.59, UL=3.50, RC=1.95, AN=8.08, UDC=5.66)
true.votes <- list(FI=.24, UL=.40, RC=.10, AN=.18, UDC=.08)
# comparison against true votes and party positions
nash.eq <- equilibrium(model=m, data=election, pos=true.pos, votes=true.votes)
nash.eq
par(mfrow=c(3,1))
plot(nash.eq)
# bootstrap confidence intervals
nash.eq <- equilibrium(model=m, data=election, boot=10)
nash.eq
plot(nash.eq)
par(mfrow=c(1,1))
}
Run the code above in your browser using DataLab