#R version 3.3.2 
  
data(iris)
str(iris)
attach(iris)
library(nnet)
names(iris)
levels(Species)

#Creating the Model
DV=relevel(Species, ref="setosa")
Mod=multinom(DV~Sepal.Length+Sepal.Width+Petal.Length+Petal.Width)
summary(Mod)

#Classification by the model
Predict=predict(Mod, iris)
Prob=predict(Mod, iris, type="prob")

#Check for misclassification
Tab=table(predict(Mod), iris$Species)
Tab