Question

我的数据框是

library(car)
DF2 <- data.frame(WindVel=c(2.45,2.7,2.9,3.05,3.4,3.6,3.95,4.1,4.6,5),
                  DCOutput=c(0.123,0.5,0.653,0.558,1.057,1.137,1.144,NA,NA,NA))

后来我想通过排除NA来对数据帧执行Box-Tidwell转换。因此我使用了：

boxTidwell(DF2$DCOutput, DF2$WindVel, options(na.action="na.exclude"))

这会出错：

Error in model.frame.default(formula = y ~ cbind(x1, x2),  drop.unused.levels = TRUE) : 
invalid type (list) for variable 'cbind(x1, x2)'

我试过了：

boxTidwell(DF2$DCOutput, DF2$WindVel, na.action=na.exclude)

然而，它不起作用。请告诉我如何在使用boxTidwell时排除NA。我对用于NA的用法的car包的用户手册无法理解。

Answer 1

或者您可以使用pam_unix.so的公式方法。来自boxTidwell：

?boxTidwell

由于## S3 method for class 'formula' boxTidwell(formula, other.x=NULL, data=NULL, subset, na.action=getOption("na.action"), verbose=FALSE, tol=0.001, max.iter=25, ...) na.action: a function that indicates what should happen when the data contain NAs. The default is set by the na.action setting of options.选项的开箱即用价值为na.action，您将获得完整案例分析，就像您从na.omit那样。

在您的情况下，您可以使用

na.exclude

这种参数模式对于使用公式接口的统计建模函数很常见。

约翰

boxTidwell - 'NA'处理

1 个答案: