[R-es] paquete pls

[Emilio L. Cano]

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 using 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 using gmail.com <mailto:emilopezcano using 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 using gmail.com <mailto:anthogab14 using 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 using r-project.org <mailto:R-help-es using 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 using r-project.org <mailto:R-help-es using 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]]