car.pref <- read.table("C:\\Kurssit\\Glim\\Glim02\\Datat\\Car.preferences.txt", header = TRUE, sep = ",") car.pref.d <- read.table("C:\\Kurssit\\Glim\\Glim02\\Datat\\Car.pref.d.txt", header = TRUE, sep = ",") names(car.pref) library(nnet) pref.d<- multinom(response ~ x1 + x2 + x3,weight=frequency, data = car.pref.d) pref<- multinom(response ~ sex + age ,weight=frequency, data = car.pref) px <- expand.grid(x3 = c(0,1), x2= c(0,1), x1= c(0,1)) predict(pref.d,px,type="prob") pref.d<- multinom(response ~ x1 + age,weight=frequency, data = car.pref.d)