0
두 개의 매개 변수를 포함하는 하나의 관측 세트가 있습니다.두 개의 매개 변수 관측을 코 콥터로 연결
copula (합치기 및 여백 함수의 매개 변수 추정)에 맞게하려면 어떻게해야합니까?
마진 분포가 로그 정규 분포이고, 접합부가 Gumbel copula라고 가정 해 봅시다.
데이터는 아래와 같다 :
1 974.0304 1010
2 6094.2672 1150
3 3103.2720 1490
4 1746.1872 1210
5 6683.7744 3060
6 6299.6832 3330
7 4784.0112 1550
8 1472.4288 607
9 3758.5728 1970
10 4381.2144 1350
Library(copula)
gumbel.cop <- gumbelCopula(dim=2)
myMvd <- mvdc(gumbel.cop, c("lnorm","lnorm"), list(list(meanlog = 7.1445391,sdlog=0.4568783), list(meanlog = 7.957392,sdlog=0.559831)))
x <- rmvdc(myMvd, 1000)
fit <- fitMvdc(x, myMvd, c(7.1445391,0.4568783,7.957392,0.559831))
meanlog 및 sdlog 값은 데이터 세트로부터 도출된다. 오류 메시지 :
"Error in if (alpha - 1 < .Machine$double.eps^(1/3)) return(rCopula(n, :
missing value where TRUE/FALSE needed"
주어진 데이터 및 데이터 세트에서 파생 된 여백 분포를 사용하여 copula 매개 변수를 선택하는 방법은 무엇입니까?
당신은 무엇을 시도 했습니까? 어떤 문제에 직면 해 있습니까? – Llopis
답장을 보내 주셔서 대단히 감사합니다. 나는 두 매개 변수에 대한 단 변수 분석을 마쳤으며, 분포가 "로그 정규 분포"로 추정되었다. 그리고 "copula"패키지를 설치합니다. "fitCopula"함수는 코 펄라에 대한 매개 변수를 필요로합니다. 이것은 내가 예상 할 수있는 것이기 때문에 이해하지 못했습니다. 또한, 코 플레 분석을 위해 분포 (log-normal)를 균일 한 한계 분포로 변환하는 방법은 무엇입니까? – simonyy
이 질문과 관련하여 더 많은 질문이 있으면 질문을 편집하여 추가하십시오. 따라서이 매개 변수없이이 함수를 사용하면 정확한 오류는 무엇입니까? 복사 해 주시겠습니까? 이 매개 변수에 대한 패키지 도움말을 읽었습니까? 보여준 데이터를 줄이거 나 (단 5 줄 또는 10 줄로 할 것인가?) 당신이 그것을 분석하는데 사용하는 코드를 복사 할 수 있습니까? 이렇게하면 문제를보기 위해 코드를 직접 복사하여 붙여 넣을 수 있습니다. – Llopis