R函数错误，“没有适用于'预测'的方法应用于“NULL”类的对象

Question

我正在尝试使用一个函数来生成一系列 KS 测试的 p 值。就我而言，我希望从我的 aaa 数据集中选择不同的前 M[i] 行，并经过几个步骤生成预测 glm 回归值 dd，最后比较每两个数据集以获得 KS 测试的 p 值。但是我的代码效果不佳。我使用模拟数据在下面附加了我的代码：

set.seed(1)
aaa <- matrix(runif(8000*950), 8000, 950)
ph = sample(rep(0:1,l=950))
kst <- function(M){
  N = length(M)
  aa = list()
  bb = list()
  cc = list()
  dd = list()
  ee = list()
  p = 0
  for (i in 2:N){
    aa[[i]] = aaa[1:M[i],] ## class(aaa) is matrix 
    bb[[i]] = rbind(ph, aa[[i]])
    cc[[i]] = as.data.frame(t(bb[[i]]))
    dd[[i]] = glm(data = cc[[i]], formula = ph~., family = binomial)
    ee[[i]] = ks.test(predict(dd[[i-1]]), predict(dd[[i]]))
    p[i] = ee[[i]]$p.value
  }
  return(p)
}
kst(M = seq(50,950,by=50))

但是R说：

> kst(M = seq(50,950,by=50))
Error in UseMethod("predict") : 
  no applicable method for 'predict' applied to an object of class "NULL"

我认为我的 ks.test 循环也有错误，但我不知道如何修复它。 谁能帮我解决这个问题？提前致谢。

Answer 1

问题可能出在for (i in 1:N-1)

R 将1:N-1 解释为序列 1 到 N，然后从整个向量中减去 1。 试试1:5-1 看看我的意思。

试试for (i in 1:(N-1))

也在这一行

 ee[[i]] = ks.test(predict(dd[[i]]), predict(dd[[i+1]])) ## this line may also have problem.

dd[[i+1]] 尚未定义，因为您的循环仅将dd 填充到i。

不确定这是您想要做的，但也许：

  if(i>1L)ee[[i]] = ks.test(predict(dd[[i-1]]), predict(dd[[i]]))

这样，在您填写了dd 的第一个条目后，您就开始测试当前的dd 与之前的对比。