基於r中分類變量的線性回歸

Question

我看過各種資源，但找不到解決我問題的明顯方法。 如果在其他地方回答此問題，則指向正確方向的指針將很有用。 使用以下數據集：

    tmp <- c(35,35,35,35,35,36,36,36,36,36)
    experiment <- factor(tmp)
    rpm <- c(10,20,40,80,120,10,20,40,80,120)
    x <- c(678,1889,3416,8916,17917,665,1385,3377,8551,16793)
    test <- data.frame(experiment,rpm,x)

我想執行回歸，例如對x使用在rpm上使用lm的線性回歸：

    lm(x ~ rpm)

因此，在此示例中，我將進行兩個回歸，一個用於實驗35，一個用於實驗36。我可以對數據進行子集化，但是我的實際數據將具有200多個實驗，我想知道是否可以僅使用幾行碼？

Answer 1

有plyr ：

library(plyr)
dlply(test, .(experiment), function(df) with(df, lm(x ~ rpm)))