r - 如何在lm()中子集一个值范围

Question

lm()的帮助文件不包含在subset参数的语法中。我不确定如何找到最适合我的一部分数据集的线。这个问题是相似的，但是我无法使用它解决我的特定问题。 How does the subset argument work in the lm() function?

这是我的代码:

    with(dat[dat$SIZE <7 & dat$SIZE > 0.8 ,], plot(SP.RICH~SIZE, log="x",
      xlim=c(1,9), ylim=c(60,180), ylab="plant species richness", 
      xlab="log area (ha)", type="n"))
   with(dat[dat$SIZE <7 & dat$SIZE > 0.8 ,], points(SP.RICH~SIZE, pch=20, cex=1))
   fit=lm(SP.RICH~SIZE, subset=c(1:7))

我想确保仅对我在plot()和points()命令中上面子集的值绘制回归线。

Answer 1

lm()和其他模型拟合函数中的子集参数将数据帧长度的逻辑矢量作为其参数，该长度是在数据帧的环境中评估的。因此，如果我对您的理解正确，我将使用以下内容:

fit <- lm(SP.RICH~SIZE, data=dat, subset=(SIZE>0.8 & SIZE<7))