data(we8there)
y <- normalize(we8thereRatings$Overall)
F <- normalize(freq(we8thereCounts))
fit <- pls(F, y)
plot(fit, y=as.factor(we8thereRatings$Overall), col=c(2,2,2,3,3))
## Uncomment the following to get multinomial inverse regression.
## The MAP estimation can take up to 5 minutes, depending on your machine.
# fit <- mnlm(we8thereCounts, y)
# par(mfrow=c(1,2))
# plot(fit, y= as.factor(we8thereRatings$Overall), col=c(2,2,2,3,3))
# plot(fit, cat=2842, delta=0.05)Run the code above in your browser using DataLab