*Buenas tardes*
*
*
*Los siguientes comando: *
library(coda)
library(spdep)
library(rgdal)
NY8 <- readOGR(".", "NY8_utm18")
NY_nb <- read.gal("NY_nb.gal", region.id = row.names(as(NY8,
"data.frame")))
nylm <- lm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data = NY8)
summary(nylm)
*Generan el siguiente resultado (en el archivo adjunto se puede visualizar
resultado y archivos para correr estos comandos ): *
*
*
Quisiera saber cómo puedo sacar los valores p. Al utilizar la opción $
únicamente me permite sacar los estimaciones de los coeficientes: 1.53e-02,
2.97e-08 y -1.92e-08 pero no los valores p.
Agradezco su ayuda
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20110630/cc26c704/attachment-0001.html>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: Consulta.docx
Type: application/vnd.openxmlformats-officedocument.wordprocessingml.document
Size: 76809 bytes
Desc: no disponible
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20110630/cc26c704/attachment-0001.bin>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: archivos.zip
Type: application/zip
Size: 285947 bytes
Desc: no disponible
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20110630/cc26c704/attachment-0001.zip>
Hola Ricardo, prueba a hacer: s <- summary(nylm) s$coefficients["(Intercept)", "Pr(>|z|)"] s$coefficients["PEXPOSURE", "Pr(>|z|)"] s$coefficients["PCTAGE65P", "Pr(>|z|)"] s$coefficients["PCTOWNHOME", "Pr(>|z|)"] Un saludo, Guillermo> *Buenas tardes* > > * > * > > *Los siguientes comando: * > > > library(coda) > > library(spdep) > > library(rgdal) > > NY8 <- readOGR(".", "NY8_utm18") > > NY_nb <- read.gal("NY_nb.gal", region.id = row.names(as(NY8,"data.frame")))> > nylm <- lm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data = NY8) > > summary(nylm) > > > > *Generan el siguiente resultado (en el archivo adjunto se puede visualizar > resultado y archivos para correr estos comandos ): * > > * > * > > Quisiera saber cómo puedo sacar los valores p. Al utilizar la opción $ > únicamente me permite sacar los estimaciones de los coeficientes:1.53e-02,> 2.97e-08 y -1.92e-08 pero no los valores p. > > > > Agradezco su ayuda >
Hola, Con los datos que has aportado, he podido reproducir completamente tu ejemplo. Y he podido acceder al p-value de la misma manera que ya te comentaba en el correo anterior. Adjunto el detalle con los datos de tu ejemplo: ######################> sCall: lm(formula = Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data = NY8) Residuals: Min 1Q Median 3Q Max -1.7417 -0.3957 -0.0326 0.3353 4.1398 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -0.51728 0.15856 -3.262 0.00124 ** PEXPOSURE 0.04884 0.03506 1.393 0.16480 PCTAGE65P 3.95089 0.60550 6.525 3.22e-10 *** PCTOWNHOME -0.56004 0.17031 -3.288 0.00114 ** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 0.6571 on 277 degrees of freedom Multiple R-squared: 0.1932, Adjusted R-squared: 0.1844 F-statistic: 22.1 on 3 and 277 DF, p-value: 7.306e-13> names(s) [1] "call" "terms" "residuals" "coefficients"[5] "aliased" "sigma" "df" "r.squared" [9] "adj.r.squared" "fstatistic" "cov.unscaled" > > cf.v<-s$fstatistic> cf.v value numdf dendf 22.10496 3.00000 277.00000 > > > pf(cf.v[1], cf.v[2], cf.v[3], lower.tail=F) value 7.306218e-13 > > p.value<-pf(cf.v[1], cf.v[2], cf.v[3], lower.tail=F)> p.value value 7.306218e-13 ###################### Saludos, Carlos Ortega www.qualityexcellence.es 2011/6/30 Ricardo Borda <ricardobordah@yahoo.es>> *Buenas tardes* > > * > * > > *Los siguientes comando: * > > > library(coda) > > library(spdep) > > library(rgdal) > > NY8 <- readOGR(".", "NY8_utm18") > > NY_nb <- read.gal("NY_nb.gal", region.id = row.names(as(NY8, > "data.frame"))) > > nylm <- lm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data = NY8) > > summary(nylm) > > > > *Generan el siguiente resultado (en el archivo adjunto se puede visualizar > resultado y archivos para correr estos comandos ): * > > * > * > > Quisiera saber cómo puedo sacar los valores p. Al utilizar la opción $ > únicamente me permite sacar los estimaciones de los coeficientes: > 1.53e-02, 2.97e-08 y -1.92e-08 pero no los valores p. > > > > Agradezco su ayuda > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > >[[alternative HTML version deleted]]
Hola,
cuando usas summary.lm (el resultado de aplicar summary() sobre el objeto
resultante al estimar un modelo con lm) el resultado es una lista con 11
elementos. El cuarto elemento de esa lista (de nombre coefficients) es una
matriz de dimensión número de coeficientes x 4 (4x4 en este caso)
donde las columnas son:
1) los coeficientes
2) sus desv. típicas estimadas
3) el estadístico t de la Ho: coef=0
4) y el p-value
A ti te interesa extraer esta cuarta columna:
summary(nylm)$coef[,4]
Jorge
PD: te adjunto el código con breves explicaciones:
library(coda)
library(spdep)
library(rgdal)
NY8 <- readOGR(".", "NY8_utm18")
NY_nb <- read.gal("NY_nb.gal", region.id = row.names(as(NY8,
"data.frame")))
nylm <- lm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data = NY8)
nylm.summary <- summary(nylm)
# Lista de 11 elementos
str(nylm.summary)
names(nylm.summary)
# El cuarto elemento denominado coefficients
# es una matriz de número de coeficientes x 4 (4x4 en este caso)
# donde las columnas son: los coeficientes, sus desv. típicas estimadas,
# el estadístico t de la Ho: coef=0 y el p-value
is.matrix(nylm.summary$coef)
dim(nylm.summary$coef)
nylm.summary$coef
# Los valore p, por tanto, son la cuarta columna
nylm.summary$coef[,4]
# O por el nombre
colnames(nylm.summary$coef)
nylm.summary$coef[,"Pr(>|t|)"]
El Thursday 30 June 2011 20:38:37 escribió:> *Buenas tardes*
>
> *
> *
>
> *Los siguientes comando: *
>
>
> library(coda)
>
> library(spdep)
>
> library(rgdal)
>
> NY8 <- readOGR(".", "NY8_utm18")
>
> NY_nb <- read.gal("NY_nb.gal", region.id = row.names(as(NY8,
> "data.frame")))
>
> nylm <- lm(Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME, data = NY8)
>
> summary(nylm)
>
>
>
> *Generan el siguiente resultado (en el archivo adjunto se puede visualizar
> resultado y archivos para correr estos comandos ): *
>
> *
> *
>
> Quisiera saber cómo puedo sacar los valores p. Al utilizar la opción $
> únicamente me permite sacar los estimaciones de los coeficientes:
> 1.53e-02, 2.97e-08 y -1.92e-08 pero no los valores p.
>
>
>
> Agradezco su ayuda