元测试代码错误

Question

我正在尝试对 R 执行 Yuen 的稳健测试，但我遇到以下错误代码的问题：

Error in formula.default(object, env = baseenv()) : invalid formula

我使用了以下包和代码：

library(WRS2)
library(car)
library(akima)
library(MASS)
library(robustbase)
library(rrcov)
library(scatterplot3d)
library(trimcluster)
library(parallel)
PremS <- c(0.06,    0.02,   0.04,   0.03,   0.03,   0.01,   0.04,   0.02,   0.07,   0.05,   0,  0.11,   0.04,   0.02,   0.03,   0.03,   0.04,   0.01,   0.06,   0.03,   0.02,   0.03,   0.03,   0.02,   0.02,   0.03,   0.02,   0.02,   0.04,   0.03,   0.03,   0.03,   0.02,   0.04,   0.06,   0.03)
PremR <- c(0.04,    0.02,   0.05,   0.02,   0.04,   0.02,   0.01,   0.01,   0.02,   0.04,   0,  0.02,   0.01,   0.01,   0.01,   0.03,   0.02,   0.03,   0.02,   0.02,   0.02,   0.01,   0.01,   0.04,   0.01,   0.03,   0.03,   0.03,   0.02)
yuen(PremR,PremS,tr=0.2)

如果有人能提供帮助，我将非常感激！

Answer 1

首先，为了调用函数yuen，你只需要加载包WRS2。其他的有点过分了。

现在回答问题。您需要使用类formula 的对象调用该函数。像这样：

yuen(PremR ~ PremS,tr=0.2)
Error in model.frame.default(formula) : 
  variable lengths differ (found for 'PremS')

现在抛出的错误不同了。 所以我检查了两个变量的长度。

length(PremS)
[1] 36
length(PremR)
[1] 29

您现在需要更正您的数据。