tcopula - не может сгенерировать действительные случайные числа из многомерного распределения - PullRequest
0 голосов
/ 31 августа 2018

Я использую пакет copula, чтобы оценить многомерное распределение, я использую tcopula, чтобы связать маргиналы. Однако генератор случайных чисел из многомерного распределения возвращает бессмыслицу.

library("forecast")
library("copula")
library("MASS")
library("LaplacesDemon")
lpm <- list()
for (i in 1:11) {
  pm <- fitdistr(na.omit(RVV[,i]), densfun = "t")$estimate
  lpm[[i]] <- list(mu = pm[1], tau = pm[2], nu = pm[3])
}
t.cop <- tCopula(param = cor(RVV)[,1][-1], dim = 11, dispstr = "toep",
                 df = 4, df.fixed = TRUE)
t.cop2 <- fitCopula(t.cop, pobs(RVV))@copula
mvdc1 <- mvdc(t.cop2, margins = c("stp","stp","stp","stp","stp","stp","stp","stp","stp","stp","stp"),
              paramMargins = lpm )
rvv <- rMvdc(252, mvdc1)

rvv содержит сгенерированное случайное число.

использованные данные, "RVV"

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...