[R-es] glm com etiquetas en las variables
Hola: Si aún hay alguien que no está de vacaciones, igual me pueden ayudar. Quiero ajustar unos modelos: REG_LOG - glm (low ~ X, family = binomial, data = DATOS) Ejemplo: library(MASS) data(birthwt, package=MASS) birthwt$low - factor(birthwt$low) birthwt$race - factor(birthwt$smoke) REG_LOG - glm (low ~ smoke, family = binomial, data = birthwt) summary (REG_LOG) Se pueden colocar etiquetas en las variable de la fórmula de manera que en los resultados salgan las etiquetas y no el nombre de la variable? Seria una cosa como: REG_LOG - glm (low = Bajo peso recién nacido ~ smoke = Consumo tabaco embarazo, family = binomial, data = birthwt) que, evidentemente no funciona. Muchas gracias y saludos ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es
Re: [R-es] glm com etiquetas en las variables: solucionado
Hola: Dos soluciones: On Tue, 21 Jul 2015 11:26:32 +0200 Carlos Ortega c...@qualityexcellence.es wrote: Esta es una alternativa sólo válida para el modelo de ejemplo, pero te da una idea de cómo se podría hacer... Se puede generalizar para que tengas un vector con las relaciones nombres de las variables y su descripción y al hacer la asignación (mySummary$call) las cambie programáticamente. #--- mySummary - summary (REG_LOG) mySummary$call - c(glm(formula = low = 'Bajo peso...' ~ smoke = 'Consumo tabaco..', family = 'binomial', data = birthwt)) mySummary #--- - El resultado de summary(), lo guardo en una variable. - El resultado de summary() es un objeto que tiene diferentes elementos. Lo que devuelve cada función, aparece en la ayuda de R bajo el epígrafe Value. Mira lo que devuelve summary.glm.. - Uno de los elementos que devuelve summary.glm es call que es la llamada, la fórmula del modelo. - Lo que hago es coger esa parte mySummary$call y modificarlo con el literal que yo quiero. - Una vez modificado, al ejecutar mySummary aparece la modificación, junto con el resto de valores que no he tocado. On Tue, 21 Jul 2015 11:31:44 +0200 Jose Luis Cañadas Reche canadasre...@gmail.com wrote: Otra idea es utilizar el paquete memisc library(MASS) data(birthwt, package=MASS) birthwt$low - factor(birthwt$low) birthwt$race - factor(birthwt$smoke) REG_LOG - glm (low ~ smoke, family = binomial, data = birthwt) summary (REG_LOG) library(memisc) mtable1 - mtable(REG_LOG) mtable1 - relabel(mtable1, (Intercept) = Constante, smoke = Consumo tabaco embarazo ) mtable1 Calls: REG_LOG: glm(formula = low ~ smoke, family = binomial, data = birthwt) == Constante-1.087*** (0.215) Consumo tabaco embarazo 0.704* (0.320) -- Aldrich-Nelson R-sq.0.025 McFadden R-sq. 0.021 Cox-Snell R-sq. 0.025 Nagelkerke R-sq.0.036 phi 1.000 Likelihood-ratio4.867 p 0.027 Log-likelihood -114.902 Deviance 229.805 AIC 233.805 BIC 240.288 N 189 == # o un fichero separado por comas write.mtable(mtable1, file = resumen.csv, colsep=,) Las dos funcionan sin problema. La Jose Luis Cañadas es mucho más sencilla, pero utiliza el enesimo paquete memisc y, si puedo, los evito: ¿que pasará si el que lo mantiene se jubila antes que yo? Implemento por tanto la primera. Como el código va dentro de una función que solo realiza reg logísticas, no me importa que tenga que escribir un poca más de código una vez. Quedaría así: library(MASS) data(birthwt, package=MASS) birthwt$low - factor(birthwt$low) birthwt$race - factor(birthwt$smoke) REG_LOG - glm (low ~ smoke, family = binomial, data = birthwt) mySummary - summary (REG_LOG) mySummary$call - c(glm(formula = low = 'Bajo peso...' ~ smoke = 'Consumo tabaco..', family = 'binomial', data = birthwt)) mySummary Call: glm(formula = low = 'Bajo peso...' ~ smoke = 'Consumo tabaco..', family = 'binomial', data = birthwt) Deviance Residuals: Min 1Q Median 3Q Max -1,0197 -0,7623 -0,7623 1,3438 1,6599 Coefficients: Estimate Std. Error z valuePr(|z|) (Intercept) -1,0871 0,2147 -5,062 0,00414 *** smoke 0,7041 0,3196 2,203 0,0276 * --- [borrado] Gracias y saludos. El 21 de julio de 2015, 10:57, Griera gri...@yandex.com escribió: Hola: Si aún hay alguien que no está de vacaciones, igual me pueden ayudar. Quiero ajustar unos modelos: REG_LOG - glm (low ~ X, family = binomial, data = DATOS) Ejemplo: library(MASS) data(birthwt, package=MASS) birthwt$low - factor(birthwt$low) birthwt$race - factor(birthwt$smoke) REG_LOG - glm (low ~ smoke, family = binomial, data = birthwt) summary (REG_LOG) Se pueden colocar etiquetas en las variable de la fórmula de manera que en los resultados salgan las etiquetas y no el nombre de la variable? Seria una cosa como: REG_LOG - glm (low = Bajo peso recién nacido ~ smoke = Consumo tabaco embarazo, family = binomial, data = birthwt) que, evidentemente no funciona. Muchas gracias y saludos ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es ___ R-help-es mailing list
Re: [R-es] Duda R
Hola Marcos. Utilicé hace mucho tiempo la función krige, pero en este caso me temo que el error está en los datos que estás suministrando. Comprueba que en plantilla.sp está la columna llamada RadPcp, porque el error lo que dice es que no encuentra ese objeto (primero mira a ver si es una variable de los datos que suministras, y al no encontrarla, lo intenta como objeto y tampoco lo encuentra). Creo, aunque sin unos datos para comprobarlo no estoy seguro, que el problema viene en este fragmento de código: plantilla.sp - as.data.frame(radar.mx)[,1:2] coordinates(plantilla.sp) - ~x+y class(plantilla.sp) donde creo que te estás dejando la columna RadPcp fuera. Espero que te sirva. ¡Un saludo! Víctor Granda García PhD. Plant Biology Dpto BOS, Área de Fisiología Vegetal Universidad de Oviedo El mar., 21 jul. 2015 a las 11:31, Marcos Bermejo ( markbermej...@hotmail.com) escribió: Hola, Estoy intentando hacer un kriging con deriva externa con la función krige(), pero me da un error: Error in eval(expr, envir, enclos) : object 'radPcp' not found El código principal es: radar.mx - rasterToPoints(radar.cu) colnames(radar.mx) - c(x,y,radPcp) radar.spdf - as.data.frame(radar.mx) lluvia.pluv - as.data.frame(t(rbind(pluvio.h[1,log.NA], coord.aux[,log.NA]))) colnames(lluvia.pluv) - c(pluvPcp, x, y) lluvia.pluv.spdf - lluvia.pluv coordinates(lluvia.pluv.spdf) - ~x+y #radarpx - as(radar.spdf, SpatialPixelsDataFrame) plantilla.sp - as.data.frame(radar.mx)[,1:2] coordinates(plantilla.sp) - ~x+y class(plantilla.sp) for (v in 1:length(pluvio.h[1,log.NA])) { lluvia.pluv.df - data.frame(lluvia.pluv.spdf) [-v,] lluvia.pluv.spdf - lluvia.pluv.df coordinates(lluvia.pluv.spdf) - ~x+y # Hacemos el KED. Ver función krige(): KED.rad - krige( formula=pluvPcp~radPcp, # covariable - radar locations=lluvia.pluv.spdf, newdata=plantilla.sp, # podría ser cualquier objeto Spatial model=v.fit, # modelo de semivariograma. maxdist=5000 ) ¿Me podríais ayudar? Un saludo, [[alternative HTML version deleted]] ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es [[alternative HTML version deleted]] ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es
Re: [R-es] glm com etiquetas en las variables
Hola, Esta es una alternativa sólo válida para el modelo de ejemplo, pero te da una idea de cómo se podría hacer... Se puede generalizar para que tengas un vector con las relaciones nombres de las variables y su descripción y al hacer la asignación (mySummary$call) las cambie programáticamente. #--- mySummary - summary (REG_LOG) mySummary$call - c(glm(formula = low = 'Bajo peso...' ~ smoke = 'Consumo tabaco..', family = 'binomial', data = birthwt)) mySummary #--- Saludos, Carlos Ortega www.qualityexcellence.es El 21 de julio de 2015, 10:57, Griera gri...@yandex.com escribió: Hola: Si aún hay alguien que no está de vacaciones, igual me pueden ayudar. Quiero ajustar unos modelos: REG_LOG - glm (low ~ X, family = binomial, data = DATOS) Ejemplo: library(MASS) data(birthwt, package=MASS) birthwt$low - factor(birthwt$low) birthwt$race - factor(birthwt$smoke) REG_LOG - glm (low ~ smoke, family = binomial, data = birthwt) summary (REG_LOG) Se pueden colocar etiquetas en las variable de la fórmula de manera que en los resultados salgan las etiquetas y no el nombre de la variable? Seria una cosa como: REG_LOG - glm (low = Bajo peso recién nacido ~ smoke = Consumo tabaco embarazo, family = binomial, data = birthwt) que, evidentemente no funciona. Muchas gracias y saludos ___ 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]] ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es
[R-es] Duda R
Hola, Estoy intentando hacer un kriging con deriva externa con la funci�n krige(), pero me da un error: Error in eval(expr, envir, enclos) : object 'radPcp' not found El c�digo principal es: radar.mx - rasterToPoints(radar.cu) colnames(radar.mx) - c(x,y,radPcp) radar.spdf - as.data.frame(radar.mx) lluvia.pluv - as.data.frame(t(rbind(pluvio.h[1,log.NA], coord.aux[,log.NA]))) colnames(lluvia.pluv) - c(pluvPcp, x, y) lluvia.pluv.spdf - lluvia.pluv coordinates(lluvia.pluv.spdf) - ~x+y #radarpx - as(radar.spdf, SpatialPixelsDataFrame) plantilla.sp - as.data.frame(radar.mx)[,1:2] coordinates(plantilla.sp) - ~x+y class(plantilla.sp) for (v in 1:length(pluvio.h[1,log.NA])) { lluvia.pluv.df - data.frame(lluvia.pluv.spdf) [-v,] lluvia.pluv.spdf - lluvia.pluv.df coordinates(lluvia.pluv.spdf) - ~x+y # Hacemos el KED. Ver funci�n krige(): KED.rad - krige( formula=pluvPcp~radPcp, # covariable - radar locations=lluvia.pluv.spdf, newdata=plantilla.sp, # podr�a ser cualquier objeto Spatial model=v.fit, # modelo de semivariograma. maxdist=5000 ) �Me podr�ais ayudar? Un saludo, [[alternative HTML version deleted]] ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es
Re: [R-es] glm com etiquetas en las variables
Hola. Otra idea es utilizar el paquete memisc library(MASS) data(birthwt, package=MASS) birthwt$low - factor(birthwt$low) birthwt$race - factor(birthwt$smoke) REG_LOG - glm (low ~ smoke, family = binomial, data = birthwt) summary (REG_LOG) library(memisc) mtable1 - mtable(REG_LOG) mtable1 - relabel(mtable1, (Intercept) = Constante, smoke = Consumo tabaco embarazo ) mtable1 Calls: REG_LOG: glm(formula = low ~ smoke, family = binomial, data = birthwt) == Constante-1.087*** (0.215) Consumo tabaco embarazo 0.704* (0.320) -- Aldrich-Nelson R-sq.0.025 McFadden R-sq. 0.021 Cox-Snell R-sq. 0.025 Nagelkerke R-sq.0.036 phi 1.000 Likelihood-ratio4.867 p 0.027 Log-likelihood -114.902 Deviance 229.805 AIC 233.805 BIC 240.288 N 189 == # o un fichero separado por comas write.mtable(mtable1, file = resumen.csv, colsep=,) El 21/07/15 a las 11:26, Carlos Ortega escribió: Hola, Esta es una alternativa sólo válida para el modelo de ejemplo, pero te da una idea de cómo se podría hacer... Se puede generalizar para que tengas un vector con las relaciones nombres de las variables y su descripción y al hacer la asignación (mySummary$call) las cambie programáticamente. #--- mySummary - summary (REG_LOG) mySummary$call - c(glm(formula = low = 'Bajo peso...' ~ smoke = 'Consumo tabaco..', family = 'binomial', data = birthwt)) mySummary #--- Saludos, Carlos Ortega www.qualityexcellence.es El 21 de julio de 2015, 10:57, Griera gri...@yandex.com escribió: Hola: Si aún hay alguien que no está de vacaciones, igual me pueden ayudar. Quiero ajustar unos modelos: REG_LOG - glm (low ~ X, family = binomial, data = DATOS) Ejemplo: library(MASS) data(birthwt, package=MASS) birthwt$low - factor(birthwt$low) birthwt$race - factor(birthwt$smoke) REG_LOG - glm (low ~ smoke, family = binomial, data = birthwt) summary (REG_LOG) Se pueden colocar etiquetas en las variable de la fórmula de manera que en los resultados salgan las etiquetas y no el nombre de la variable? Seria una cosa como: REG_LOG - glm (low = Bajo peso recién nacido ~ smoke = Consumo tabaco embarazo, family = binomial, data = birthwt) que, evidentemente no funciona. Muchas gracias y saludos ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es