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