Hola Anthony,
Mi recomendación sería que empleases plsr() en lugar de mrv(), por ejemplo
siguiendo los pasos descritos aquí:
https://rubenfcasal.github.io/aprendizaje_estadistico/pca-pls.html#regresi%C3%B3n-por-m%C3%ADnimos-cuadrados-parciales-plsr
La idea es no fijar el número de componentes a priori y seleccionarlas a
posteriori.
Aprovecho para comentar que en el cuatrimestre pasado impartimos por
primera vez una asignatura de aprendizaje estadístico (en el master
interuniversitario "gallego" de estadística) y preparamos unos apuntes
en
formato bookdown: https://rubenfcasal.github.io/aprendizaje_estadistico.
Cualquier sugerencia será bien recibida...
Un saludo, Rubén.
El mié, 10 mar 2021 a las 11:28, Emilio L. Cano (<emilopezcano en
gmail.com>)
escribió:
> Anthony,
>
> ¿Si no usas validación cruzada también te da el error? ¿Cuáles son las
> dimensiones de moTest?
>
> Yo probaría sin validación cruzada, o con LOO en vez de CV. Sin datos con
> qué reproducirlo no puedo decir mucho más.
>
> Un saludo,
>
> Emilio L. Cano
> http://emilio.lcano.com
>
>
>
>
> > El 10 mar 2021, a las 11:17, Anthony Gabourel <anthogab14 en
gmail.com>
> escribió:
> >
> > Estimados/as:
> >
> > Estoy trabajando con el paquete "pls" para generar modelos
predictivos.
> Sin
> > embargo, me encuentro con un problema a la hora de establecer el
número
> de
> > componentes, ya que me aparece el siguiente error y no me deja
trabajar
> con
> > el número de componentes que deseo (solo funciona con 1):
> >
> >
> > *Error in pls::mvr(MO ~ s350:s2500, ncomp = 10, data = moTest,
> validation > > "CV", : Invalid number of components,
ncomp*
> >
> > El script en cuestión es el siguiente:
> >
> > library(pls)
> > library(readxl)
> > library(dplyr)
> >
> > pls.options(plsralg="oscorespls")
> >
> > db_spectra <-
read_excel("D:\\Espectros\\db_spectra.xlsx")
> > head(db_spectra)
> >
> > moTest <- select(db_spectra, MO, s350:s2500)
> >
> > mo1 <- plsr(MO ~ s350:s2500, ncomp=10, data=moTest,
validation="CV")
> >
> > summary(mo1)
> >
> > plot(mo1, ncomp=1, asp=1, line=TRUE)
> >
> > Si alguien me pudiera echar una mano para solucionar el problema, lo
> > agradecería mucho.
> >
> > Un cordial saludo,
> >
> > Anthony Gabourel
> >
> > [[alternative HTML version deleted]]
> >
> > _______________________________________________
> > R-help-es mailing list
> > R-help-es en r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-help-es
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
--
Rubén Fernández Casal
https://rubenfcasal.github.io
Dep. Matemáticas, Universidade da Coruña
Centro de Investigación en Tecnologías de la Información y las
Comunicaciones (CITIC)
[[alternative HTML version deleted]]