我试图根据线性回归和Theil回归(mblm函数w /重复= FALSE)来引导我的数据以获得2000次测量。
我的引导程序R代码完全适用于正常回归(从我所知道的),如下所示:
> fitfunc <- function(formula, data, index) {
+ d<- data[index,]
+ f<- lm(formula,data=d)
+ return(coef(f))
+ }
boot(dataframe, fitfunc, R=2000, formula=`Index A`~`Measurement B`)
但是在尝试Theil估算器引导程序时出现错误:
> fitfuncTheil <- function(formula,data,index) {
+ d<- data[index,]
+ f<- mblm(formula, data=d, repeated=FALSE)
+ return(coef(f))
+ }
> boot(dataframe, fitfuncTheil, R=2000, formula=`Index A`~`Measurement B`)
Error in order(x) : argument 1 is not a vector
In addition: Warning message:
In is.na(x) :
错误信息似乎是基本的,但我无法弄清楚为什么这会在一个案例中起作用而在另一个案例中不起作用。
答案 0 :(得分:0)
一旦我从列名中删除了空格(在公式字段中引用),问题就解决了。