Eva Prieto Castro
2012-Jul-15 10:44 UTC
[R-es] Summary(lm) no devuelve intervalo de confianza
Buenos días: Tengo un modelo de las siguientes características: res <- lm(Y~ t + I(t^2) + I(t^3)) Resulta que el summary del res me devuelve ésto:> summary(res)Call: lm(formula = as.formula(modelo)) Residuals: Min 1Q Median 3Q Max -20.911 -8.699 -3.224 3.802 51.739 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 11.49248 9.93797 1.156 0.252 t -1.49054 3.99867 -0.373 0.711 I(t^2) 0.38444 0.43685 0.880 0.383 I(t^3) -0.01109 0.01370 -0.809 0.422 Residual standard error: 15.76 on 56 degrees of freedom Multiple R-squared: 0.4272, Adjusted R-squared: 0.3965 F-statistic: 13.92 on 3 and 56 DF, p-value: 6.765e-07 Ahí puedo consultar el valor estimado, el standard error, el t y el p. Pero, ¿cómo calculo el intervalo de confianza?. Gracias. Eva [[alternative HTML version deleted]]
J. Pelayo izquierdo
2012-Jul-15 10:48 UTC
[R-es] Summary(lm) no devuelve intervalo de confianza
Sí haces: summary(res)[-1] Te muestra todas las opciones del summary. Un saludo, Pelayo. El 15/07/2012 12:44, "Eva Prieto Castro" <evapcastro@yahoo.es> escribió:> > Buenos días: > > Tengo un modelo de las siguientes características: > > res <- lm(Y~ t + I(t^2) + I(t^3)) > > Resulta que el summary del res me devuelve ésto: > > > summary(res) > Call: > lm(formula = as.formula(modelo)) > Residuals: > Min 1Q Median 3Q Max > -20.911 -8.699 -3.224 3.802 51.739 > Coefficients: > Estimate Std. Error t value Pr(>|t|) > (Intercept) 11.49248 9.93797 1.156 0.252 > t -1.49054 3.99867 -0.373 0.711 > I(t^2) 0.38444 0.43685 0.880 0.383 > I(t^3) -0.01109 0.01370 -0.809 0.422 > Residual standard error: 15.76 on 56 degrees of freedom > Multiple R-squared: 0.4272, Adjusted R-squared: 0.3965 > F-statistic: 13.92 on 3 and 56 DF, p-value: 6.765e-07 > > Ahí puedo consultar el valor estimado, el standard error, el t y el p. > Pero, ¿cómo calculo el intervalo de confianza?. > > 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 > >[[alternative HTML version deleted]]
Hola, mira si esto te sirve: http://www.r-tutor.com/elementary-statistics/simple-linear-regression/confidence-interval-linear-regression Usan la función predict.lm, no sé si es lo que buscas. Saludos, Olmo. El Sun, 15 Jul 2012 11:44:27 +0100 (BST) Eva Prieto Castro <evapcastro en yahoo.es> escribió:> > Buenos días: > > Tengo un modelo de las siguientes características: > > res <- lm(Y~ t + I(t^2) + I(t^3)) > > Resulta que el summary del res me devuelve ésto: > > > summary(res) > Call: > lm(formula = as.formula(modelo)) > Residuals: > Min 1Q Median 3Q Max > -20.911 -8.699 -3.224 3.802 51.739 > Coefficients: > Estimate Std. Error t value Pr(>|t|) > (Intercept) 11.49248 9.93797 1.156 0.252 > t -1.49054 3.99867 -0.373 0.711 > I(t^2) 0.38444 0.43685 0.880 0.383 > I(t^3) -0.01109 0.01370 -0.809 0.422 > Residual standard error: 15.76 on 56 degrees of freedom > Multiple R-squared: 0.4272, Adjusted R-squared: 0.3965 > F-statistic: 13.92 on 3 and 56 DF, p-value: 6.765e-07 > > Ahí puedo consultar el valor estimado, el standard error, el t y el > p. Pero, ¿cómo calculo el intervalo de confianza?. > Gracias. > > Eva > [[alternative HTML version deleted]] >
Hola Eva, Revisa ?confint si necesitas calcular el intervalo de confianza para los coeficientes del modelo y ?predict.lm para el caso de las observaciones. Un saludo, Jorge.- 2012/7/15 Eva Prieto Castro <>> > Buenos días: > > Tengo un modelo de las siguientes características: > > res <- lm(Y~ t + I(t^2) + I(t^3)) > > Resulta que el summary del res me devuelve ésto: > > > summary(res) > Call: > lm(formula = as.formula(modelo)) > Residuals: > Min 1Q Median 3Q Max > -20.911 -8.699 -3.224 3.802 51.739 > Coefficients: > Estimate Std. Error t value Pr(>|t|) > (Intercept) 11.49248 9.93797 1.156 0.252 > t -1.49054 3.99867 -0.373 0.711 > I(t^2) 0.38444 0.43685 0.880 0.383 > I(t^3) -0.01109 0.01370 -0.809 0.422 > Residual standard error: 15.76 on 56 degrees of freedom > Multiple R-squared: 0.4272, Adjusted R-squared: 0.3965 > F-statistic: 13.92 on 3 and 56 DF, p-value: 6.765e-07 > > Ahí puedo consultar el valor estimado, el standard error, el t y el p. > Pero, ¿cómo calculo el intervalo de confianza?. > > 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 > >[[alternative HTML version deleted]]
Eva Prieto Castro
2012-Jul-15 14:07 UTC
[R-es] Summary(lm) no devuelve intervalo de confianza
Hola: Gracias a todos! He resuelto con la sugerencia de Jorge: confint! Un saludo. Eva --- El dom, 15/7/12, Jorge I Velez <jorgeivanvelez@gmail.com> escribió: De: Jorge I Velez <jorgeivanvelez@gmail.com> Asunto: Re: [R-es] Summary(lm) no devuelve intervalo de confianza Para: "Eva Prieto Castro" <evapcastro@yahoo.es> CC: r-help-es@r-project.org Fecha: domingo, 15 de julio, 2012 15:56 Hola Eva, Revisa ?confint si necesitas calcular el intervalo de confianza para los coeficientes del modelo y ?predict.lm para el caso de las observaciones. Un saludo, Jorge.- 2012/7/15 Eva Prieto Castro <> Buenos días: Tengo un modelo de las siguientes características: res <- lm(Y~ t + I(t^2) + I(t^3)) Resulta que el summary del res me devuelve ésto:> summary(res)Call: lm(formula = as.formula(modelo)) Residuals: Min 1Q Median 3Q Max -20.911 -8.699 -3.224 3.802 51.739 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 11.49248 9.93797 1.156 0.252 t -1.49054 3.99867 -0.373 0.711 I(t^2) 0.38444 0.43685 0.880 0.383 I(t^3) -0.01109 0.01370 -0.809 0.422 Residual standard error: 15.76 on 56 degrees of freedom Multiple R-squared: 0.4272, Adjusted R-squared: 0.3965 F-statistic: 13.92 on 3 and 56 DF, p-value: 6.765e-07 Ahí puedo consultar el valor estimado, el standard error, el t y el p. Pero, ¿cómo calculo el intervalo de confianza?. 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 [[alternative HTML version deleted]]