使用函数计算R中的子回归

Question

我目前正在使用内置数据集“ iris”。我已经为这三种植物类型分别计算了回归，并且效果很好：

reg1 = with(iris[iris$Species=="setosa",], lm(Sepal.Width~Sepal.Length))

reg2 = with(iris[iris$Species=="versicolor",], lm(Sepal.Width~Sepal.Length))

reg3 = with(iris[iris$Species=="virginica",], lm(Sepal.Width~Sepal.Length))

现在，我想知道是否还可以使用R中的by函数来计算三个子回归，并获得与上述相似的输出。

我已经这样尝试过，但是我只收到错误消息：

b <- by(iris, iris$Species, function(x){
  regby <- lm(Sepal.Width~Sepal.Length)
}

有人对我有什么提示吗？