Eva Prieto Castro
2012-Jul-01 17:28 UTC
[R-es] Regresión lineal múltiple: modelo polinómico de grado 3 superpuesto a componentes cosenoidales
Hola: Tengo un modelo de regresión lineal en el cual las componentes son cosenoidales, y lo construyo del siguiente modo: modelo = "y ~ I(t) + I(t^2) + I(t^3) + x1[, 1] + x2[, 1]" x1[, 1] = cos(2 * pi * t / periods[1]) x2[, 1] = sin(2 * pi * t / periods[1]) for (i in 2:nComp) { x1[, i] = cos(2 * pi * t / periods[i]) x2[, i] = sin(2 * pi * t / periods[i]) modelo = paste(modelo, " + x1[,", i, "] + x2[,", i, "]", sep="") } res = lm (as.formula(modelo)) Este modelo deriva en resultados correctos de su análisis correspondiente (summary(res)). Sin embargo, si añado regresión polinómica, el summary no me devuelve los estimadores correctos para las componentes lineal, cuadrática y cúbica. La línea del modelo la modifico del siguiente modo: modelo = "y ~ I(t) + I(t^2) + I(t^3) + x1[, 1] + x2[, 1]" He probado también así: modelo = "y ~ t + t^2 + t^3 + x1[, 1] + x2[, 1]" ¡Pero no hay manera! ¿Me podéis orientar?. Gracias Eva [[alternative HTML version deleted]]
Eva Prieto Castro
2012-Jul-01 17:29 UTC
[R-es] Regresión lineal múltiple: modelo polinómico de grado 3 superpuesto a componentes cosenoidales
Corrección: en el primer caso (sólo componentes cosenoidales) no utilizo el I(t) + I(t^2)... Ahí lo puse pero evidentemente no lo lleva. Saludos. --- El dom, 1/7/12, Eva Prieto Castro <evapcastro@yahoo.es> escribió: De: Eva Prieto Castro <evapcastro@yahoo.es> Asunto: Regresión lineal múltiple: modelo polinómico de grado 3 superpuesto a componentes cosenoidales Para: "R" <r-help-es@r-project.org> Fecha: domingo, 1 de julio, 2012 19:28 Hola: Tengo un modelo de regresión lineal en el cual las componentes son cosenoidales, y lo construyo del siguiente modo: modelo = "y ~ I(t) + I(t^2) + I(t^3) + x1[, 1] + x2[, 1]" x1[, 1] = cos(2 * pi * t / periods[1]) x2[, 1] = sin(2 * pi * t / periods[1]) for (i in 2:nComp) { x1[, i] = cos(2 * pi * t / periods[i]) x2[, i] = sin(2 * pi * t / periods[i]) modelo = paste(modelo, " + x1[,", i, "] + x2[,", i, "]", sep="") } res = lm (as.formula(modelo)) Este modelo deriva en resultados correctos de su análisis correspondiente (summary(res)). Sin embargo, si añado regresión polinómica, el summary no me devuelve los estimadores correctos para las componentes lineal, cuadrática y cúbica. La línea del modelo la modifico del siguiente modo: modelo = "y ~ I(t) + I(t^2) + I(t^3) + x1[, 1] + x2[, 1]" He probado también así: modelo = "y ~ t + t^2 + t^3 + x1[, 1] + x2[, 1]" ¡Pero no hay manera! ¿Me podéis orientar?. Gracias Eva [[alternative HTML version deleted]]
Marcuzzi, Javier Rubén
2012-Jul-01 21:06 UTC
[R-es] Regresión lineal múltiple: modelo polinómico de grado 3 superpuesto a componentes cosenoidales
Estimada Eva Prietro En el libro: Pinheiro, J. C. and Bates, D. M. (2000), \emph{Mixed-Effects Models in S and S-PLUS}, Springer, New York. (Appendix A.18) hay un ejemplo donde utilizan seno, coseno, ... utilizan el tamaño de folículos. De la ayuda de nlme copio el siguiente código: # Pinheiro and Bates use in nlme: # from p. 240 needed on p. 396 fm1Ovar.lme <- lme(follicles ~ sin(2*pi*Time) + cos(2*pi*Time), data = Ovary, random = pdDiag(~sin(2*pi*Time))) fm5Ovar.lme <- update(fm1Ovar.lme, corr = corARMA(p = 1, q = 1)) # p. 396 fm1Ovar.nlme <- nlme(follicles~ A+B*sin(2*pi*w*Time)+C*cos(2*pi*w*Time), data=Ovary, fixed=A+B+C+w~1, random=pdDiag(A+B+w~1), start=c(fixef(fm5Ovar.lme), 1) ) # p. 397 fm2Ovar.nlme <- update(fm1Ovar.nlme, corr=corAR1(0.311) ) No sabría responder el resto de su pregunta. Javier -----Mensaje original----- From: Eva Prieto Castro Sent: Sunday, July 01, 2012 2:29 PM To: R Subject: Re: [R-es] Regresión lineal múltiple: modelo polinómico de grado 3 superpuesto a componentes cosenoidales Corrección: en el primer caso (sólo componentes cosenoidales) no utilizo el I(t) + I(t^2)... Ahí lo puse pero evidentemente no lo lleva. Saludos. --- El dom, 1/7/12, Eva Prieto Castro <evapcastro en yahoo.es> escribió: De: Eva Prieto Castro <evapcastro en yahoo.es> Asunto: Regresión lineal múltiple: modelo polinómico de grado 3 superpuesto a componentes cosenoidales Para: "R" <r-help-es en r-project.org> Fecha: domingo, 1 de julio, 2012 19:28 Hola: Tengo un modelo de regresión lineal en el cual las componentes son cosenoidales, y lo construyo del siguiente modo: modelo = "y ~ I(t) + I(t^2) + I(t^3) + x1[, 1] + x2[, 1]" x1[, 1] = cos(2 * pi * t / periods[1]) x2[, 1] = sin(2 * pi * t / periods[1]) for (i in 2:nComp) { x1[, i] = cos(2 * pi * t / periods[i]) x2[, i] = sin(2 * pi * t / periods[i]) modelo = paste(modelo, " + x1[,", i, "] + x2[,", i, "]", sep="") } res = lm (as.formula(modelo)) Este modelo deriva en resultados correctos de su análisis correspondiente (summary(res)). Sin embargo, si añado regresión polinómica, el summary no me devuelve los estimadores correctos para las componentes lineal, cuadrática y cúbica. La línea del modelo la modifico del siguiente modo: modelo = "y ~ I(t) + I(t^2) + I(t^3) + x1[, 1] + x2[, 1]" He probado también así: modelo = "y ~ t + t^2 + t^3 + x1[, 1] + x2[, 1]" ¡Pero no hay manera! ¿Me podéis orientar?. Gracias Eva [[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
Carlos Ortega
2012-Jul-02 09:39 UTC
[R-es] Regresión lineal múltiple: modelo polinómico de grado 3 superpuesto a componentes cosenoidales
Hola Eva, Cuando dices que no devuelve los estimadores correctos: ¿quieres decir que no los devuelve? ¿que no son correctos comparados con qué otra forma de estimar?... En cuanto a la forma de especificar e modelo, efectivamente es con el operado I(). En el siguiente ejemplo, el modelo sí que estima la parte polinomial y la trigonométrica:> x <- c(32,64,96,118,126,144,152.5,158) > y <- c(99.5,104.8,108.5,100,86,64,35.3,15) > > mod.fit<-lm(y ~ x + I(x^2) + I(x^3) + cos(x) + sin(x)) > > summary(mod.fit)Call: lm(formula = y ~ x + I(x^2) + I(x^3) + cos(x) + sin(x)) Residuals: 1 2 3 4 5 6 7 8 -0.2225 0.4749 1.0493 -2.5407 -1.4982 4.4489 1.0239 -2.7356 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.180e+02 2.630e+01 4.485 0.0463 * x -1.017e+00 1.083e+00 -0.939 0.4467 I(x^2) 2.107e-02 1.262e-02 1.669 0.2370 I(x^3) -1.170e-04 4.381e-05 -2.670 0.1163 cos(x) -1.556e+00 4.513e+00 -0.345 0.7631 sin(x) -3.850e+00 3.219e+00 -1.196 0.3543 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 4.382 on 2 degrees of freedom Multiple R-squared: 0.9956, Adjusted R-squared: 0.9844 F-statistic: 89.52 on 5 and 2 DF, p-value: 0.01108 Vaya, como otras veces, si envías un conjunto de datos reproducibles, o el tipo de error... mejor. Saludos, Carlos Ortega www.qualityexcellence.es El 1 de julio de 2012 19:28, Eva Prieto Castro <evapcastro@yahoo.es>escribió:> Hola: > > Tengo un modelo de regresión lineal en el cual las componentes son > cosenoidales, y lo construyo del siguiente modo: > > > modelo = "y ~ I(t) + I(t^2) + I(t^3) + x1[, 1] + x2[, 1]" > x1[, 1] = cos(2 * pi * t / periods[1]) > x2[, 1] = sin(2 * pi * t / periods[1]) > for (i in 2:nComp) { > x1[, i] = cos(2 * pi * t / periods[i]) > x2[, i] = sin(2 * pi * t / periods[i]) > modelo = paste(modelo, " + x1[,", i, "] + x2[,", i, "]", sep="") > } > res = lm (as.formula(modelo)) > > Este modelo deriva en resultados correctos de su análisis correspondiente > (summary(res)). Sin embargo, si añado regresión polinómica, el summary no > me devuelve los estimadores correctos para las componentes lineal, > cuadrática y cúbica. > > La línea del modelo la modifico del siguiente modo: > > modelo = "y ~ I(t) + I(t^2) + I(t^3) + x1[, 1] + x2[, 1]" > > He probado también así: > > modelo = "y ~ t + t^2 + t^3 + x1[, 1] + x2[, 1]" > > ¡Pero no hay manera! > > ¿Me podéis orientar?. > > Gracias > > Eva > [[alternative HTML version deleted]] > > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > >-- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]]
Hola colegas: Necesito instalar el paquete ggplot2 y me está dando este problemas con el repositorio. W: Error de GPG: http://cran.stat.sfu.ca oneiric/ Release: Las firmas siguientes no se pudieron verificar porque su clave pública no está disponible: NO_PUBKEY 51716619E084DAB9 Ya he cambiado varias veces el repositorio y me continua sucediendo lo mismo. Alguna ayuda??? Un saludo Leonardo ------------ próxima parte ------------ An embedded and charset-unspecified text was scrubbed... Name: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20120702/0dccbb34/attachment.pl>
Hola, ¿Has probado a bajarte el paquete (.zip o .tar.gz) directamente e instalarlo desde local? Yo acabo de bajarme el paquete para Windows sin ningún problema. Vaya, si usas la versión de R para Windows esta podría ser una alternativa inmediata. Saludos, Carlos Ortega www.qualityexcellence.es El 2 de julio de 2012 15:02, Leonardo Hernández Pérez < leonardo.hernandez@etecsa.cu> escribió:> > Hola colegas: Necesito instalar el paquete ggplot2 y me está dando este > problemas con el repositorio. > > W: Error de GPG: http://cran.stat.sfu.ca oneiric/ Release: Las firmas > siguientes no se pudieron verificar porque su clave pública no está > disponible: NO_PUBKEY 51716619E084DAB9 > > Ya he cambiado varias veces el repositorio y me continua sucediendo lo > mismo. > > Alguna ayuda??? > > Un saludo > > Leonardo > --- > This message was processed by Kaspersky Mail Gateway 5.6.28/RELEASE > running at host imx2.etecsa.cu > Visit our web-site: <http://www.kaspersky.com>, <http://www.viruslist.com> > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > >-- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]]
Hola, ¿qué tal? Es una pregunta específica de Ubuntu que han respondido otros aquí: http://ubuntuforums.org/showthread.php?t=1694202 Un saludo, Carlos J. Gil Bellosta http://www.datanalytics.com> El 2 de julio de 2012 15:02, Leonardo Hernández Pérez < > leonardo.hernandez en etecsa.cu> escribió: > >> >> Hola colegas: Necesito instalar el paquete ggplot2 y me está dando este >> problemas con el repositorio. >> >> W: Error de GPG: http://cran.stat.sfu.ca oneiric/ Release: Las firmas >> siguientes no se pudieron verificar porque su clave pública no está >> disponible: NO_PUBKEY 51716619E084DAB9 >> >> Ya he cambiado varias veces el repositorio y me continua sucediendo lo >> mismo. >> >> Alguna ayuda??? >> >> Un saludo >> >> Leonardo >> --- >> This message was processed by Kaspersky Mail Gateway 5.6.28/RELEASE >> running at host imx2.etecsa.cu >> Visit our web-site: <http://www.kaspersky.com>, <http://www.viruslist.com> >> >> _______________________________________________ >> R-help-es mailing list >> R-help-es en r-project.org >> https://stat.ethz.ch/mailman/listinfo/r-help-es >> >> > > > -- > Saludos, > Carlos Ortega > www.qualityexcellence.es > > [[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 >
> Hola colegas: > > Lamentablemente este problema no he podido resolverlo, he intentado > dar todas las soluciones que me han indicado pro no he podido lograrlo. > > Aquí les resumo para si alguien tiene alguna otra me lo diga. > > El primer problema es que no puedo conectarme al repositorio, dándome > este error. > > W: Error de GPG: http://cran.stat.sfu.ca oneiric/ Release: Las firmas > siguientes no se pudieron verificar porque su clave pública no está > disponible: NO_PUBKEY 51716619E084DAB9 > > Intento resolverlo Come me sugiere Carlos J. Gil Bellosta en: > http://ubuntuforums.org/showthread.php?t=1694202 > > gpg --keyserver keyserver.ubuntu.com --recv 51716619E084DAB9 > gpg --export --armor 06F90DE5381BA480 | sudo apt-key add - > > Y al dar el primer comando me devuelve: > > leonardo en TI16:~$ gpg --keyserver keyserver.ubuntu.com --recv > 51716619E084DAB9 > gpg: solicitando clave E084DAB9 de hkp servidor keyserver.ubuntu.com > ?: keyserver.ubuntu.com: Host not found > gpgkeys: HTTP fetch error 7: couldn't connect: No such file or directory > gpg: no se han encontrados datos OpenPGP válidos > gpg: Cantidad total procesada: 0 > > también pruebo otra solución que es similar pero incluye el proxy: > > wget --user=leonardo.hernandez --password=xxxxxxxxx -q > "http://keyserver.ubuntu.com:11371/pks/lookup?op=get&search=0x51716619E084DAB9" > -O- | sudo apt-key add - > > y me devuelve: gpg: no se han encontrados datos OpenPGP válidos > > También intenté instalar el paquete manualmente como sugiere Carlos > Ortega pero tiene alrededor de 5 dependencias que a su vez tiene otras > más, es un proceso muy engorroso que comúnmente da errores (Estoy > instalándolo sobre Ubuntu Linux). > > Además reinstalé la PC como sugiere Luis Gutiérrez ahora tengo la > version 2.14 de R pero los problemas para instalar ggplot2 son los mismos. > > Alguna otra ayuda????????? > > Un Saludos, > > Leonardo------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20120704/c9db46f2/attachment.html> ------------ próxima parte ------------ An embedded and charset-unspecified text was scrubbed... Name: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20120704/c9db46f2/attachment.pl>
Hola Fernando: Increíblemente acabo de instalarlo, tuve que descargar manualmente todas las dependencias y comenzar a instalar a medida que fueron necesitándose. Estas son los paquetes y dependencias. colorspace_1.1-1 dichromat_1.2-4 digest_0.5.2 ggplot2_0.9.1 labeling_0.1 memoise_0.1 munsell_0.3 plyr_1.7.1 RColorBrewer_1.0-5 reshape2_1.2.1 scales_0.2.1 stringr_0.6 Saludos Leonardo On 04/07/12 15:06, Fernando Macedo wrote:> Sé que existe de forma global para la sesión y seguramente será > posible configurarlo solamente para el terminal. Lo que no recuerdo es > como hacerlo desde un terminal y dependiendo de la versión que tengas > (K/L/X/Ubuntu) tenes herramientas gráficas, aunque en el momento no lo > recuerdo (lo hice en Ubuntu con Gnome 2 hace mucho tiempo). > > Quedará en la ayuda de otro o por tu cuenta y google buscar la forma. > > Busca por configurar proxy en ubuntu. > > Saludos y disculpa no poder ayudarte mejor en este tema. > > Fernando > > El 04/07/12 15:08, Leonardo Hernández Pérez escribió: >> La limitación que tengo es que salgo a través de un proxy y me >> imagino que por eso no puedo llegar a esa dirección. >> >> Dime si existe la posibilidad de declarar el proxy en consola. >> >> >> Saludos, >> >> Leonardo >> >> On 04/07/12 13:30, Fernando Macedo wrote: >>> Estas conectándote en dónde? >>> >>> Lo probé acá y no tuvo ningún problema en conectarse. De repente >>> estas en una red con limitaciones. >>> >>> Saludos. >>> >>> Fernando Macedo >>> >>> El 04/07/12 13:34, Leonardo Hernández Pérez escribió: >>>> Hola Fernando: >>>> >>>> También lo hice y no me funcionó. >>>> >>>> Saludos, >>>> >>>> Leonardo >>>> >>>> On 04/07/12 12:17, Fernando Macedo wrote: >>>>> En realidad no va el número completo sino solamente los 8 >>>>> caracteres finales. >>>>> >>>>> prueba así: >>>>> >>>>> *gpg --keyserver keyserver.ubuntu.com --recv E084DAB9 >>>>> >>>>> gpg --export --armor E084DAB9 | sudo apt-key add - * >>>>> >>>>> Eso debería solucionar el problema de los repositorios (si es >>>>> Ubuntu o derivados y me parece no recuerdo bien en Debian) >>>>> >>>>> Luego intenta la instalación de los paquetes. >>>>> >>>>> Saludos >>>>> >>>>> Fernando Macedo >>>>> >>>>> >>>>> El 04/07/12 12:19, Leonardo Hernández Pérez escribió: >>>>>>> Hola colegas: >>>>>>> >>>>>>> Lamentablemente este problema no he podido resolverlo, he >>>>>>> intentado dar todas las soluciones que me han indicado pro no he >>>>>>> podido lograrlo. >>>>>>> >>>>>>> Aquí les resumo para si alguien tiene alguna otra me lo diga. >>>>>>> >>>>>>> El primer problema es que no puedo conectarme al repositorio, >>>>>>> dándome este error. >>>>>>> >>>>>>> W: Error de GPG: http://cran.stat.sfu.ca oneiric/ Release: Las >>>>>>> firmas >>>>>>> siguientes no se pudieron verificar porque su clave pública no está >>>>>>> disponible: NO_PUBKEY 51716619E084DAB9 >>>>>>> >>>>>>> Intento resolverlo Come me sugiere Carlos J. Gil Bellosta en: >>>>>>> http://ubuntuforums.org/showthread.php?t=1694202 >>>>>>> >>>>>>> gpg --keyserver keyserver.ubuntu.com --recv 51716619E084DAB9 >>>>>>> gpg --export --armor 06F90DE5381BA480 | sudo apt-key add - >>>>>>> >>>>>>> Y al dar el primer comando me devuelve: >>>>>>> >>>>>>> leonardo en TI16:~$ gpg --keyserver keyserver.ubuntu.com --recv >>>>>>> 51716619E084DAB9 >>>>>>> gpg: solicitando clave E084DAB9 de hkp servidor keyserver.ubuntu.com >>>>>>> ?: keyserver.ubuntu.com: Host not found >>>>>>> gpgkeys: HTTP fetch error 7: couldn't connect: No such file or >>>>>>> directory >>>>>>> gpg: no se han encontrados datos OpenPGP válidos >>>>>>> gpg: Cantidad total procesada: 0 >>>>>>> >>>>>>> también pruebo otra solución que es similar pero incluye el proxy: >>>>>>> >>>>>>> wget --user=leonardo.hernandez --password=xxxxxxxxx -q >>>>>>> "http://keyserver.ubuntu.com:11371/pks/lookup?op=get&search=0x51716619E084DAB9" >>>>>>> -O- | sudo apt-key add - >>>>>>> >>>>>>> y me devuelve: gpg: no se han encontrados datos OpenPGP válidos >>>>>>> >>>>>>> También intenté instalar el paquete manualmente como sugiere >>>>>>> Carlos Ortega pero tiene alrededor de 5 dependencias que a su >>>>>>> vez tiene otras más, es un proceso muy engorroso que comúnmente >>>>>>> da errores (Estoy instalándolo sobre Ubuntu Linux). >>>>>>> >>>>>>> Además reinstalé la PC como sugiere Luis Gutiérrez ahora tengo >>>>>>> la version 2.14 de R pero los problemas para instalar ggplot2 >>>>>>> son los mismos. >>>>>>> >>>>>>> Alguna otra ayuda????????? >>>>>>> >>>>>>> Un Saludos, >>>>>>> >>>>>>> Leonardo >>>>>> >>>>>> >>>>>> >>>>>> --- >>>>>> This message was processed by Kaspersky Mail Gateway 5.6.28/RELEASE running at host imx2.etecsa.cu >>>>>> Visit our web-site:<http://www.kaspersky.com>,<http://www.viruslist.com> >>>>>> >>>>>> >>>>>> _______________________________________________ >>>>>> R-help-es mailing list >>>>>> R-help-es en r-project.org >>>>>> https://stat.ethz.ch/mailman/listinfo/r-help-es >>>>> >>>>> >>>>> >>>>> >>>>> --- >>>>> This message was processed by Kaspersky Mail Gateway 5.6.28/RELEASE running at host imx2.etecsa.cu >>>>> Visit our web-site:<http://www.kaspersky.com>,<http://www.viruslist.com> >>>> >>>> >>>> >>>> >>>> --- >>>> This message was processed by Kaspersky Mail Gateway 5.6.28/RELEASE running at host imx2.etecsa.cu >>>> Visit our web-site:<http://www.kaspersky.com>,<http://www.viruslist.com> >>> >>> >>> >>> >>> --- >>> This message was processed by Kaspersky Mail Gateway 5.6.28/RELEASE running at host imx2.etecsa.cu >>> Visit our web-site:<http://www.kaspersky.com>,<http://www.viruslist.com> >> >> >> >> >> --- >> This message was processed by Kaspersky Mail Gateway 5.6.28/RELEASE running at host imx3.etecsa.cu >> Visit our web-site:<http://www.kaspersky.com>,<http://www.viruslist.com> > > > > > --- > This message was processed by Kaspersky Mail Gateway 5.6.28/RELEASE running at host imx2.etecsa.cu > Visit our web-site: <http://www.kaspersky.com>, <http://www.viruslist.com>------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20120704/8b10f126/attachment-0001.html> ------------ próxima parte ------------ An embedded and charset-unspecified text was scrubbed... Name: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20120704/8b10f126/attachment-0001.pl>