Formato di lm di R() Formula con una Trasformazione

Question

Non riesco a capire come effettuare le seguenti operazioni in una sola riga:

data(attenu) 
x_temp = attenu$accel^(1/4) 
y_temp = log(attenu$dist) 
best_line = lm(y_temp ~ x_temp) 

Dal momento che i lavori di cui sopra, ho pensato che avrei potuto fare quanto segue:

data(attenu) 
best_line = lm(log(attenu$dist) ~ (attenu$accel^(1/4))) 

Ma questo dà l'errore:

Error in terms.formula(formula, data = data) : invalid power in formula 

C'è ovviamente qualcosa che mi manca quando si usa trasformato vari abilita nel formato formula R. Perché non funziona?

Answer 1

Si sta cercando la funzione I in modo che l'operatore ^ venga trattato come aritmetico nella formula, ad es.

x <- runif(1:100) 
y <- x + rnorm(100,0, 3) 
lm(log(y) ~ I(x^(1/4))