rrvglm: ошибка в lm.fit (X.vlm, y = z.vlm, ...): NA / NaN / Inf в «x»

Question

Я хочу подогнать модель rr-vglm, используя функцию rrvglm в package-VGAM.Я снабдил его простыми данными моделирования, но я получаю ошибку

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