r-在一小段时间间隔内绘制lm

Question

我正在编写一个函数，该函数在变量的一个子集中绘制线性回归图。

FitWeibull <- function(data, xinf, xsup){
  xlim <- c(xinf, xsup)
  sub.data <- data
  my.lm <- lm(Y~X, data = sub.data, subset = xlim)
  return(my.lm)
 }

这个函数会让我在xinf和xsup之间绘制回归图，并保留图上的其他点吗？

Answer 1

您需要指定要设置的变量。最简单的方法是将sub.data <- data行替换为：sub.data <- data[(X>xinf & X<xsup),]，这样您就可以获得

FitWeibull <- function(data, xinf, xsup){
  sub.data <- data[(data$X > xinf & data$X < xsup),]
  my.lm <- lm(Y~X, data = sub.data)
  return(my.lm)
}

然后，您可以对返回的对象执行任何打印/汇总操作。