Skip to main content
se añadieron 849 caracteres en el cuerpo
Origen Enlace
Patricio Moracho
  • 61.1k
  • 12
  • 42
  • 72

SiFinalmente, lo que buscas es agregar los coeficientes estándarizadossería generar un modelo de regresión simple por cada variable del vector lm.betaindependiente a, y el resultado dejarlo salvarlo como columna de tabla_b. Para esto necesitarás:

  • Un ciclo que itere por cada variable independiente y realice la regresión, la formula puede construirse directamente desde cadenas y trasformarse en formula mediante as.formula()
  • Una lista dónde salvar el coficiente que buscas, o en definitiva cualquier otro dato que necesites

Finalmente, puedes hacer un "matching" porcon el nombreciclo completo, tendremos una lista de lacoeficientes que simplemente habrá que transformar en un vector y agregarlo como columna, por ejemplo a table_b:

coeficienteslibrary(lm.beta)

set.seed(666)
tabla_a <- lmdata.betaframe(lmmatrix(formulasample(1:5,3000,replace = X1TRUE), ~ncol=30))

dependiente X2= +"X1"
independiente X3= +c("X2","X3","X12","X20", X12"X21")
coeficientes +<- X20list()
for +(var X21in independiente) {
  formula <- as.formula(paste(dependiente, "~ ", var))
  fit <- lm(data = tabla_a, formula = formula)
  coeficientes[[var]] <- lm.beta(fit)$standardized.coefficientscoefficients[2]
}
tabla_b$coeficientes
tabla_b <-= coeficientes[matchdata.frame(tabla_b$independientedependiente, namesindependiente)
tabla_b$coeficiente <- unlist(coeficientes))] 

tabla_b

  dependiente independiente coeficientescoeficiente
1          X1            X2  -0.0726077808760396
2          X1            X3   0.0336804904449120
3          X1           X12   0.1190467012002269
4          X1           X20  -0.1310141721159223
5          X1           X21  -0.2438325926675073

Si lo que buscas es agregar los coeficientes estándarizados de lm.beta a tabla_b, puedes hacer un "matching" por el nombre de la columna, por ejemplo:

coeficientes <- lm.beta(lm(formula = X1 ~ X2 + X3 + X12 + X20 + X21, data = tabla_a))$standardized.coefficients

tabla_b$coeficientes <- coeficientes[match(tabla_b$independiente, names(coeficientes))]

tabla_b

  dependiente independiente coeficientes
1          X1            X2  -0.07260778
2          X1            X3   0.03368049
3          X1           X12   0.11904670
4          X1           X20  -0.13101417
5          X1           X21  -0.24383259

Finalmente, lo que buscas sería generar un modelo de regresión simple por cada variable del vector independiente, y el resultado dejarlo salvarlo como columna de tabla_b. Para esto necesitarás:

  • Un ciclo que itere por cada variable independiente y realice la regresión, la formula puede construirse directamente desde cadenas y trasformarse en formula mediante as.formula()
  • Una lista dónde salvar el coficiente que buscas, o en definitiva cualquier otro dato que necesites

Finalmente, con el ciclo completo, tendremos una lista de coeficientes que simplemente habrá que transformar en un vector y agregarlo como columna a table_b:

library(lm.beta)

set.seed(666)
tabla_a <-data.frame(matrix(sample(1:5,3000,replace = TRUE), ncol=30))

dependiente = "X1"
independiente = c("X2","X3","X12","X20", "X21")
coeficientes <- list()
for (var in independiente) {
  formula <- as.formula(paste(dependiente, "~ ", var))
  fit <- lm(data = tabla_a, formula = formula)
  coeficientes[[var]] <- lm.beta(fit)$standardized.coefficients[2]
}

tabla_b = data.frame(dependiente, independiente)
tabla_b$coeficiente <- unlist(coeficientes) 

tabla_b

  dependiente independiente coeficiente
1          X1            X2 -0.08760396
2          X1            X3  0.04449120
3          X1           X12  0.12002269
4          X1           X20 -0.21159223
5          X1           X21 -0.26675073
Origen Enlace
Patricio Moracho
  • 61.1k
  • 12
  • 42
  • 72

Si lo que buscas es agregar los coeficientes estándarizados de lm.beta a tabla_b, puedes hacer un "matching" por el nombre de la columna, por ejemplo:

coeficientes <- lm.beta(lm(formula = X1 ~ X2 + X3 + X12 + X20 + X21, data = tabla_a))$standardized.coefficients

tabla_b$coeficientes <- coeficientes[match(tabla_b$independiente, names(coeficientes))]

tabla_b

  dependiente independiente coeficientes
1          X1            X2  -0.07260778
2          X1            X3   0.03368049
3          X1           X12   0.11904670
4          X1           X20  -0.13101417
5          X1           X21  -0.24383259