В настоящее время я работаю со встроенным набором данных «Ирис». Я рассчитал отдельные регрессии для трех типов растений, как это, и это сработало:
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))
Теперь мне интересно, можно ли также рассчитать три субрегрессии, используя функцию by в R, и получить результат, аналогичный приведенному выше.
Я пробовал вот так, но я получаю только сообщения об ошибках:
b <- by(iris, iris$Species, function(x){
regby <- lm(Sepal.Width~Sepal.Length)
}
У кого-нибудь есть советы для меня?