基于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)

我想执行回归，例如使用lm on rpm对x：

进行线性回归

    lm(x ~ rpm)

因此，在这个例子中，我将有两个回归，一个用于实验35，一个用于实验36.我可以对数据进行分组，但我的实际数据将有超过200个实验，我想知道我是否可以只用一个几行代码？

Answer 1

有plyr：

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