Rubén,
Vi los apuntes el otro día: espectacular el trabajo que habéis hecho. Muchas
gracias por compartir (pronto lo haremos en otro foro como sabes)
El caso es que el error que pone @Anthony es de mvr, pero en el código usa plsr,
no sé si se equivocó al pegar, o es que plsr llama a mvr.
Un saludo a todos,
Emilio L. Cano
Tel: 665 676 225
http://emilio.lcano.com
> El 10 mar 2021, a las 16:21, Rubén Fernández Casal <rubenfcasal en
gmail.com> escribió:
>
> 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
<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
<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 <mailto: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 <http://emilio.lcano.com/>
>
>
>
>
> > El 10 mar 2021, a las 11:17, Anthony Gabourel <anthogab14 en
gmail.com <mailto: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 <mailto:R-help-es en r-project.org>
> > https://stat.ethz.ch/mailman/listinfo/r-help-es
<https://stat.ethz.ch/mailman/listinfo/r-help-es>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org <mailto:R-help-es en r-project.org>
> https://stat.ethz.ch/mailman/listinfo/r-help-es
<https://stat.ethz.ch/mailman/listinfo/r-help-es>
>
>
> --
>
> Rubén Fernández Casal
> https://rubenfcasal.github.io <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]]