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 <[email protected]> > 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 (<[email protected] > <mailto:[email protected]>>) 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 <[email protected] > > <mailto:[email protected]>> 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 > > [email protected] <mailto:[email protected]> > > https://stat.ethz.ch/mailman/listinfo/r-help-es > > <https://stat.ethz.ch/mailman/listinfo/r-help-es> > > _______________________________________________ > R-help-es mailing list > [email protected] <mailto:[email protected]> > 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]] _______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
