Я использую Bicopselect, используя приведенную ниже цитату с набором семейств, ограниченным от 1 до 7 и 9. Но мои результаты по-прежнему показывают копулы, которые находятся за пределами моих настроек, например, 17 (повернут BB1) для некоторых пар.
Может ли кто-нибудь сообщить мне, если что-то не так, когда я устанавливаю семейный набор в своих кодах.
Ваша помощь очень ценится.
myfunc<-function(rankresi) {
pairs<-combn(c(1:ncol(rankresi)),2,simplify=FALSE)
res<-sapply(pairs,function(x)
BiCopSelect(rankresi[,x[1]],rankresi[,x[2]], familyset=c(1:7,9)))
res
}