我不太清楚如何在一行中执行以下操作:
data(attenu)
x_temp = attenu$accel^(1/4)
y_temp = log(attenu$dist)
best_line = lm(y_temp ~ x_temp)
由于上述工作,我想我可以做到以下几点:
data(attenu)
best_line = lm( log(attenu$dist) ~ (attenu$accel^(1/4)) )
但这给出了错误:
Error in terms.formula(formula, data = data) : invalid power in formula
当以R的公式格式使用转换后的变量时,显然缺少某些东西。为什么不起作用?
最佳答案
您正在寻找函数I
,以便^
运算符在公式(即)中被视为算术。
x <- runif(1:100)
y <- x + rnorm(100,0, 3)
lm(log(y) ~ I(x^(1/4))
关于r - R的lm()公式的格式与转换,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/30929115/