[R] Help please with error from nnet::multinom
Lars Bishop
lars52r at gmail.com
Sun Jun 26 20:32:41 CEST 2016
Thanks Bert.
But I it doesn't complain when predict is used on X instead of X_new
(using nnet_7.3-12), which is even more puzzling to me:
pred <- predict(fit, X, type = "probs")
head(pred)
ysim1 ysim2 ysim3
1 0.3059421 0.3063284 0.3877295
2 0.3200219 0.3202551 0.3597230
3 0.3452414 0.3451460 0.3096125
4 0.3827077 0.3819603 0.2353320
5 0.2973288 0.2977994 0.4048718
6 0.3817027 0.3809759 0.2373214
Thanks again,
Lars.
On Sun, Jun 26, 2016 at 1:05 PM, Bert Gunter <bgunter.4567 at gmail.com> wrote:
> Well, for one thing, there is no "probs" method for predict.nnet, at
> least in my version: nnet_7.3-12
>
> Cheers,
> Bert
>
>
> Bert Gunter
>
> "The trouble with having an open mind is that people keep coming along
> and sticking things into it."
> -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
>
>
> On Sun, Jun 26, 2016 at 9:27 AM, Lars Bishop <lars52r at gmail.com> wrote:
> > Hello,
> >
> > I'd appreciate your help in spotting the reason for the error and warning
> > messages below.
> >
> > library(nnet)
> > set.seed(1)
> > ysim <- gl(3, 100)
> > y <- model.matrix(~ysim -1)
> > X <- matrix( 3 * runif(length(ysim)), nrow = 300, ncol = 3)
> > X_new <- matrix( 3 * runif(length(ysim)), nrow = 200, ncol = 3)
> >
> > fit <- multinom(y ~ X, trace = FALSE)
> > pred <- predict(fit, X_new, type = "probs")
> >
> > Error in predict.multinom(fit, X_new, type = "probs") :
> > NAs are not allowed in subscripted assignments
> > In addition: Warning message:
> > 'newdata' had 200 rows but variables found have 300 rows
> >
> > Thanks,
> > Lars.
> >
> > [[alternative HTML version deleted]]
> >
> > ______________________________________________
> > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> > and provide commented, minimal, self-contained, reproducible code.
>
[[alternative HTML version deleted]]
More information about the R-help
mailing list