rrvglm：lm.fit（X.vlm，y = z.vlm，...）中的错误：NA / NaN / Inf in'x'

Question

我想使用VGAM包中的rrvglm函数来拟合rr-vglm模型。我用简单的模拟数据拟合了它，但出现错误

Error in lm.fit(X.vlm, y = z.vlm, ...) : NA/NaN/Inf in 'x'

以下是代码：

library("VGAM")
n <- 100
p <- 200
q <- 6
X <- matrix(ncol = p, nrow = n, rnorm(n * p))
Y <- matrix(0, n, q)
set <- sample(1:q, n, replace = T)
for(i in 1:n){
  Y[i, set[i]] <- 1
}
fit.vglm <- rrvglm(Y ~ X, multinomial(refLevel = 1), Rank = 2)

rrvglm是否可以在“ n