r - lm.fit(x,y,offset = offset, single.ok,...) 中的错误 0 个使用 boxcox 公式的非 NA 情况

Question

我正在尝试使用以下代码运行 boxcox 转换:

urban1 <- subset(ski,urban <= 4,na.rm=TRUE)
ski$gender <- as.numeric((as.character(ski$gender)),na.rm=TRUE)
urban1 <- as.numeric((as.character(urban1)))
x <- (ski$gender*urban1)
y <- ski$EPSI.
bc <- boxcox(y ~ x) 
(trans <- bc$x[which.max(bc$y)]) 
model3 <- lm(y ~ x) 
model3new <- lm(y^trans ~ x)
ski$EPSI. <- ski$EPSI. + 1

但我不断收到此错误:

lm.fit(x,y,offset = offset, Single.ok = Single.ok, ...) 中的错误:0(非 NA)情况调用:... eval -> eval -> boxcar -> boxcar.formula -> lm -> lm.fit 执行已停止

提前致谢!

Answer 1

错误信息

lm.fit(x,y,offset = offset, singular.ok = singular.ok, ...) : 0 (non-NA) cases

是当变量x或y(或两者)仅具有NA时由lm(y ~ x)命令生成的。
这是一个例子:

n <- 10
x <- rnorm(n,1)
y <- rep(NA,n)
lm(y ~ x)

Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : 
  0 (non-NA) cases

在您的代码中，我建议测试(就在您的 lm 命令之前)您的变量之一是否具有所有 NA 使用:

all(is.na(x))
all(is.na(y))
all(is.na(y^trans))

在我的例子中:

all(is.na(y))
[1] TRUE