[R-es] paquete pls

[Anthony Gabourel]

Hola, Emilio y Rubén:

Gracias por sus respuestas.

Parece que ya he solucionado el problema, sustituyendo las variables
"s350:s2500" por ".", de manera que ahora sí se seleccionan los datos
correctamente y me es posible establecer el número de componentes que deseo.

Por otro lado, en el código uso plsr() pero en el mensaje de error sí que
aparecía mvr().

Saludos cordiales,

Anthony Gabourel






El mié, 10 mar 2021 a las 17:01, Emilio L. Cano (<emilopezcano using gmail.com>)
escribió:

> 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
> 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 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
>>
>>
>>
>>
>> > El 10 mar 2021, a las 11:17, Anthony Gabourel <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
>> > https://stat.ethz.ch/mailman/listinfo/r-help-es
>>
>> _______________________________________________
>> R-help-es mailing list
>> R-help-es using 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]]