R にはたくさんのパッケージがあるので

同じことをするにも，何通りかある。そして，パッケージによっては出力される結果の種類に違いがある（答だって，微妙に違うこともあるけど(^_^)

さて，「多値のロジスティック回帰を復習しておこう&Rで実装しておこう!」で（実装という言葉の意味を間違えているだろうけど）

> 多項ロジスティック回帰には nnet パッケージにある multinom() 関数を使います。
 :
> この関数、問題はp値が算出されないのが問題ですね。

とあるが，当然の所，各群の所属確率なんだろうね。predict で type="probs" を指定すれば，出るよ!

 http://rpubs.com/r-de-r/1651 を参照