R help es - Apr 2012 - Coeficientes Beta en un modelo de regresión lineal aplicando stepwise regression: ¿cómo los obtengo en R?

Hola a todos, el motivo de esta pregunta es que he estado trabajando con modelos
de regresión lineal aplicando stepwise regression. He consultado todas las
fuentes que he encontrado en internet pero no veo que el R te de los
coeficientes Beta directamente. Es decir, que ya estén programados por alguien y
puedas sacarlos haciendo algo así como:
library(stats)
resultado<-step(modelo,direction="both",trace=1,steps=5000,k=2)
resultado$betas
Sería una pena que no esté implementado, dado que en muchos artículos (p. ej. en
temas médicos) lo emplean, probablemente porque el SPSS te lo da directamente.
 
saludos,
 
Fernando
 
Os pongo aquí qué son eso de los coeficientes beta: 
Beta (standardised regression coefficients) The beta value
is a measure of how strongly each predictor variable influences the criterion
variable. The beta is measured in units of standard deviation. For example, a
beta value of 2.5 indicates that a change of one standard deviation in the
predictor variable will result in a change of 2.5 standard deviations in the
criterion variable. Thus, the higher the beta value the greater the impact of
the predictor variable on the criterion variable.
	[[alternative HTML version deleted]]

Fernando
Probá con

resultado$coefficients

Prof. Julio Di Rienzo
Estadística y Biometría
FCA- U.N. Córdoba
IBS-RARG President
http://sites.google.com/site/juliodirienzo
"Biometry, the active pursuit of biological
knowledge by quantitative methods."
(R.A. Fisher, 1948)



2012/4/18 Fernando Sanchez <fernandsanche@yahoo.es>
> Hola a todos, el motivo de esta pregunta es que he estado trabajando con
> modelos de regresión lineal aplicando stepwise regression. He consultado
> todas las fuentes que he encontrado en internet pero no veo que el R te de
> los coeficientes Beta directamente. Es decir, que ya estén programados por
> alguien y puedas sacarlos haciendo algo así como:
> library(stats)
>
resultado<-step(modelo,direction="both",trace=1,steps=5000,k=2)
> resultado$betas
> Sería una pena que no esté implementado, dado que en muchos artículos (p.
> ej. en temas médicos) lo emplean, probablemente porque el SPSS te lo da
> directamente.
>
> saludos,
>
> Fernando
>
> Os pongo aquí qué son eso de los coeficientes beta:
> Beta (standardised regression coefficients) The beta value
> is a measure of how strongly each predictor variable influences the
> criterion
> variable. The beta is measured in units of standard deviation. For
> example, a
> beta value of 2.5 indicates that a change of one standard deviation in the
> predictor variable will result in a change of 2.5 standard deviations in
> the
> criterion variable. Thus, the higher the beta value the greater the impact
> of
> the predictor variable on the criterion variable.
>        [[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 Julio, muchísimas gracias pero esos no son los coeficientes que necesito,
esos que salen si no me equivoco, son los de la propia ecuación. Voy a buscar un
ejemplo y os lo envío.
 
saludos,
 
Fernando
 
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
Fernando 
Probá con 
resultado$coefficients
Prof. Julio Di Rienzo
Estadística y Biometría
FCA- U.N. Córdoba
IBS-RARG President
http://sites.google.com/site/juliodirienzo
"Biometry, the active pursuit of biological
knowledge by quantitative methods." 
(R.A. Fisher, 1948)
 
 


Hola a todos, el motivo de esta pregunta es que he estado trabajando con modelos
de regresión lineal aplicando stepwise regression. He consultado todas las
fuentes que he encontrado en internet pero no veo que el R te de los
coeficientes Beta directamente. Es decir, que ya estén programados por alguien y
puedas sacarlos haciendo algo así como:
library(stats)
resultado<-step(modelo,direction="both",trace=1,steps=5000,k=2)
resultado$betas
Sería una pena que no esté implementado, dado que en muchos artículos (p. ej. en
temas médicos) lo emplean, probablemente porque el SPSS te lo da directamente.
 
saludos,
 
Fernando
 
Os pongo aquí qué son eso de los coeficientes beta:
Beta (standardised regression coefficients) The beta value
is a measure of how strongly each predictor variable influences the criterion
variable. The beta is measured in units of standard deviation. For example, a
beta value of 2.5 indicates that a change of one standard deviation in the
predictor variable will result in a change of 2.5 standard deviations in the
criterion variable. Thus, the higher the beta value the greater the impact of
the predictor variable on the criterion variable.
       [[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 Fernando.

Quizá esto te pueda ser útil:

http://sas-and-r.blogspot.com/2010/11/example-814-generating-standardized.html?utm_source=feedburner&utm_medium=feed&utm_campaign=Feed%3A+SASandR+%28SAS+and+R%29

Salud.

On Wed, Apr 18, 2012 at 16:57, Fernando Sanchez
<fernandsanche@yahoo.es>wrote:
> Hola Julio, muchísimas gracias pero esos no son los coeficientes que
> necesito, esos que salen si no me equivoco, son los de la propia ecuación.
> Voy a buscar un ejemplo y os lo envío.
>
> saludos,
>
> Fernando
>
>
>
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
>
> Fernando
> Probá con
> resultado$coefficients
> Prof. Julio Di Rienzo
> Estadística y Biometría
> FCA- U.N. Córdoba
> IBS-RARG President
> http://sites.google.com/site/juliodirienzo
> "Biometry, the active pursuit of biological
> knowledge by quantitative methods."
> (R.A. Fisher, 1948)
>
>
>
>
> Hola a todos, el motivo de esta pregunta es que he estado trabajando con
> modelos de regresión lineal aplicando stepwise regression. He consultado
> todas las fuentes que he encontrado en internet pero no veo que el R te de
> los coeficientes Beta directamente. Es decir, que ya estén programados por
> alguien y puedas sacarlos haciendo algo así como:
> library(stats)
>
resultado<-step(modelo,direction="both",trace=1,steps=5000,k=2)
> resultado$betas
> Sería una pena que no esté implementado, dado que en muchos artículos (p.
> ej. en temas médicos) lo emplean, probablemente porque el SPSS te lo da
> directamente.
>
> saludos,
>
> Fernando
>
> Os pongo aquí qué son eso de los coeficientes beta:
> Beta (standardised regression coefficients) The beta value
> is a measure of how strongly each predictor variable influences the
> criterion
> variable. The beta is measured in units of standard deviation. For
> example, a
> beta value of 2.5 indicates that a change of one standard deviation in the
> predictor variable will result in a change of 2.5 standard deviations in
> the
> criterion variable. Thus, the higher the beta value the greater the impact
> of
> the predictor variable on the criterion variable.
>        [[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
>
>

-- 
«But Gwindor answered: ''The doom lies in yourself, not in your
name.''»

JRR Tolkien

	[[alternative HTML version deleted]]

Buenas noches Fernando,

No responde directamente a tu pregunta, pero la discusion/resumen en [1]
podria ser de gran utilidad, sobretodo porque menciona los peligros de usar
procedimientos stepwise para la seleccion de modelos de regresion.

En cuanto a tu pregunta, creo que

betas <- coefficients(summary(resultado))
betas[, 1]/betas[,2]

es lo que necesitas.

Un saludo,
Jorge.-

[1] http://core.ecu.edu/psyc/wuenschk/stathelp/Stepwise-Voodoo.htm


2012/4/18 Fernando Sanchez <>
> Hola a todos, el motivo de esta pregunta es que he estado trabajando con
> modelos de regresión lineal aplicando stepwise regression. He consultado
> todas las fuentes que he encontrado en internet pero no veo que el R te de
> los coeficientes Beta directamente. Es decir, que ya estén programados por
> alguien y puedas sacarlos haciendo algo así como:
> library(stats)
>
resultado<-step(modelo,direction="both",trace=1,steps=5000,k=2)
> resultado$betas
> Sería una pena que no esté implementado, dado que en muchos artículos (p.
> ej. en temas médicos) lo emplean, probablemente porque el SPSS te lo da
> directamente.
>
> saludos,
>
> Fernando
>
> Os pongo aquí qué son eso de los coeficientes beta:
> Beta (standardised regression coefficients) The beta value
> is a measure of how strongly each predictor variable influences the
> criterion
> variable. The beta is measured in units of standard deviation. For
> example, a
> beta value of 2.5 indicates that a change of one standard deviation in the
> predictor variable will result in a change of 2.5 standard deviations in
> the
> criterion variable. Thus, the higher the beta value the greater the impact
> of
> the predictor variable on the criterion variable.
>        [[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]]

Gracias Freddy. Esta explicación es bastante mejor que la de Fernando ya 
que la de Fernando:

  The beta is measured in units of standard deviation

Al no especificar las sd de qué variable, hacía pensar simplemente en 
los valores tipificados de los coeficientes (y evidentemente no es el 
caso). Se trata de "... in standard deviations [of outcome], per 
standard deviation change in the predictors."

Gracias


El 19/04/12 03:09, Freddy López escribió:
> Hola Fernando.
>
> Quizá esto te pueda ser útil:
>
>
http://sas-and-r.blogspot.com/2010/11/example-814-generating-standardized.html?utm_source=feedburner&utm_medium=feed&utm_campaign=Feed%3A+SASandR+%28SAS+and+R%29
>
> Salud.
>
> On Wed, Apr 18, 2012 at 16:57, Fernando
Sanchez<fernandsanche@yahoo.es>wrote:
>
>    
>> Hola Julio, muchísimas gracias pero esos no son los coeficientes que
>> necesito, esos que salen si no me equivoco, son los de la propia
ecuación.
>> Voy a buscar un ejemplo y os lo envío.
>>
>> saludos,
>>
>> Fernando
>>
>>
>>
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
>>
>> Fernando
>> Probá con
>> resultado$coefficients
>> Prof. Julio Di Rienzo
>> Estadística y Biometría
>> FCA- U.N. Córdoba
>> IBS-RARG President
>> http://sites.google.com/site/juliodirienzo
>> "Biometry, the active pursuit of biological
>> knowledge by quantitative methods."
>> (R.A. Fisher, 1948)
>>
>>
>>
>>
>> Hola a todos, el motivo de esta pregunta es que he estado trabajando
con
>> modelos de regresión lineal aplicando stepwise regression. He
consultado
>> todas las fuentes que he encontrado en internet pero no veo que el R te
de
>> los coeficientes Beta directamente. Es decir, que ya estén programados
por
>> alguien y puedas sacarlos haciendo algo así como:
>> library(stats)
>>
resultado<-step(modelo,direction="both",trace=1,steps=5000,k=2)
>> resultado$betas
>> Sería una pena que no esté implementado, dado que en muchos artículos
(p.
>> ej. en temas médicos) lo emplean, probablemente porque el SPSS te lo da
>> directamente.
>>
>> saludos,
>>
>> Fernando
>>
>> Os pongo aquí qué son eso de los coeficientes beta:
>> Beta (standardised regression coefficients) The beta value
>> is a measure of how strongly each predictor variable influences the
>> criterion
>> variable. The beta is measured in units of standard deviation. For
>> example, a
>> beta value of 2.5 indicates that a change of one standard deviation in
the
>> predictor variable will result in a change of 2.5 standard deviations
in
>> the
>> criterion variable. Thus, the higher the beta value the greater the
impact
>> of
>> the predictor variable on the criterion variable.
>>         [[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
>>
>>
>>      
>
>    
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>    
-- 
_____---^---_____

Univ. de Extremadura
Dept. Matemáticas.
Despacho B29
Tf: + 34 924 289 300
Ext. 86823


	[[alternative HTML version deleted]]

Hola,

Aquí tienes una referencia y solución en R (segunda respuesta):

http://stats.stackexchange.com/questions/22346/standardized-beta-weights-for-a-multilevel-regression


Saludos,
Carlos Ortega
www.qualityexcellence.es

El 18 de abril de 2012 22:55, Fernando Sanchez
<fernandsanche@yahoo.es>escribió:
> Hola a todos, el motivo de esta pregunta es que he estado trabajando con
> modelos de regresión lineal aplicando stepwise regression. He consultado
> todas las fuentes que he encontrado en internet pero no veo que el R te de
> los coeficientes Beta directamente. Es decir, que ya estén programados por
> alguien y puedas sacarlos haciendo algo así como:
> library(stats)
>
resultado<-step(modelo,direction="both",trace=1,steps=5000,k=2)
> resultado$betas
> Sería una pena que no esté implementado, dado que en muchos artículos (p.
> ej. en temas médicos) lo emplean, probablemente porque el SPSS te lo da
> directamente.
>
> saludos,
>
> Fernando
>
> Os pongo aquí qué son eso de los coeficientes beta:
> Beta (standardised regression coefficients) The beta value
> is a measure of how strongly each predictor variable influences the
> criterion
> variable. The beta is measured in units of standard deviation. For
> example, a
> beta value of 2.5 indicates that a change of one standard deviation in the
> predictor variable will result in a change of 2.5 standard deviations in
> the
> criterion variable. Thus, the higher the beta value the greater the impact
> of
> the predictor variable on the criterion variable.
>        [[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 a todos, quiero agradeceros vuestros buenos consejos. Finalmente lo tengo
resuelto empleando el link de Carlos Ortega.Os adjunto un código que al
ejecutarse da los coeficientes Beta:
 
datos =
matrix(c(20,65,174,22,70,180,19,68,170,12,142,123,23,24,12,234,524,534,54,756),nrow=5,byrow=T)
colnames(datos)=c("a","b","c","d")
datos<-as.data.frame(datos)
lm1<-lm(d ~ .,data=datos)
library(QuantPsyc)
lm.beta(lm1)
 
Con el fin de garantizar que esto era lo que quería, he metido la misma
información en SPSS y he obtenido el mismo resultado así que doy el tema por
cerrado. De nuevo gracias!!!!!
 
Disculpad mi tardanza en responder pero el modelo lineal con mis datos me da el
siguiente error:
Mensajes de aviso perdidos
In b * sx :
  longitud de objeto mayor no es múltiplo de la longitud de uno menor
 
Y eso todavía no lo he resuelto pero creo que es algo al margen de estos
coeficientes y está relacionado con que algunos coeficientes el modelo lineal da
NA.
 
Saludos,
Fernando

________________________________________________________________________________________________________________________________________________________________
Hola,
Aquí tienes una referencia y solución en R (segunda respuesta):
http://stats.stackexchange.com/questions/22346/standardized-beta-weights-for-a-multilevel-regression
Saludos,
Carlos Ortega
www.qualityexcellence.es

ribió:
Hola a todos, el motivo de esta pregunta es que he estado trabajando con modelos
de regresión lineal aplicando stepwise regression. He consultado todas las
fuentes que he encontrado en internet pero no veo que el R te de los
coeficientes Beta directamente. Es decir, que ya estén programados por alguien y
puedas sacarlos haciendo algo así como:
library(stats)
resultado<-step(modelo,direction="both",trace=1,steps=5000,k=2)
resultado$betas
Sería una pena que no esté implementado, dado que en muchos artículos (p. ej. en
temas médicos) lo emplean, probablemente porque el SPSS te lo da directamente.
 
saludos,
 
Fernando
 
Os pongo aquí qué son eso de los coeficientes beta:
Beta (standardised regression coefficients) The beta value
is a measure of how strongly each predictor variable influences the criterion
variable. The beta is measured in units of standard deviation. For example, a
beta value of 2.5 indicates that a change of one standard deviation in the
predictor variable will result in a change of 2.5 standard deviations in the
criterion variable. Thus, the higher the beta value the greater the impact of
the predictor variable on the criterion variable.
       [[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
____________________________________________________________________________________________________________________________________________________
Gracias Freddy. Esta explicación es bastante mejor que la de Fernando ya que la
de Fernando:
  The beta is measured in units of standard deviation
Al no especificar las sd de qué variable, hacía pensar simplemente en los
valores tipificados de los coeficientes (y evidentemente no es el caso). Se
trata de "... in standard deviations [of outcome], per standard deviation
change in the predictors."
Gracias

El 19/04/12 03:09, Freddy López escribió:
> Hola Fernando.
>
> Quizá esto te pueda ser útil:
>
> http://sas-and-r.blogspot.com/2010/11/example-814-generating-standardi
>
zed.html?utm_source=feedburner&utm_medium=feed&utm_campaign=Feed%3A+SA
> SandR+%28SAS+and+R%29
>
> Salud.
>
ote:
>
>    
>> Hola Julio, muchísimas gracias pero esos no son los coeficientes que 
>> necesito, esos que salen si no me equivoco, son los de la propia
ecuación.
>> Voy a buscar un ejemplo y os lo envío.
>>
>> saludos,
>>
>> Fernando
>>
>>
>> ---------------------------------------------------------------------
>> ---------------------------------------------------------------------
>> ---------------------------------------------------------------------
>> ---------------------------------------------------------
>>
>> Fernando
>> Probá con
>> resultado$coefficients
>> Prof. Julio Di Rienzo
>> Estadística y Biometría
>> FCA- U.N. Córdoba
>> IBS-RARG President
>> http://sites.google.com/site/juliodirienzo
>> "Biometry, the active pursuit of biological knowledge by
quantitative
>> methods."
>> (R.A. Fisher, 1948)
>>
>>
>>
>>
>> Hola a todos, el motivo de esta pregunta es que he estado trabajando 
>> con modelos de regresión lineal aplicando stepwise regression. He 
>> consultado todas las fuentes que he encontrado en internet pero no 
>> veo que el R te de los coeficientes Beta directamente. Es decir, que 
>> ya estén programados por alguien y puedas sacarlos haciendo algo así
como:
>> library(stats)
>>
resultado<-step(modelo,direction="both",trace=1,steps=5000,k=2)
>> resultado$betas
>> Sería una pena que no esté implementado, dado que en muchos artículos
(p.
>> ej. en temas médicos) lo emplean, probablemente porque el SPSS te lo 
>> da directamente.
>>
>> saludos,
>>
>> Fernando
>>
>> Os pongo aquí qué son eso de los coeficientes beta:
>> Beta (standardised regression coefficients) The beta value is a 
>> measure of how strongly each predictor variable influences the 
>> criterion variable. The beta is measured in units of standard 
>> deviation. For example, a beta value of 2.5 indicates that a change 
>> of one standard deviation in the predictor variable will result in a 
>> change of 2.5 standard deviations in the criterion variable. Thus, 
>> the higher the beta value the greater the impact of the predictor 
>> variable on the criterion variable.
>>         [[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
>>
>>
>>      
>
>    
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>    
--
_____---^---_____
Univ. de Extremadura
Dept. Matemáticas.
Despacho B29
Tf: + 34 924 289 300
Ext. 86823

 [[alternative HTML version deleted]]
	[[alternative HTML version deleted]]