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