在一个小间隔中绘制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)
}

然后，您可以使用返回的对象执行所需的任何绘图/摘要。