Displaying 1 result from an estimated 1 matches for "up_pi1".
Did you mean:
up_pi2
2017 Aug 06
1
Help with optim function in R, please?
...copula.=copula){
??? copula[[1]]@parameters <- par[1]
??? copula[[2]]@parameters <- par[2]
??? tau_1 <- par[3]*dCopula(data,copula[[1]],log = F)
??? tau_2 <- par[4]*dCopula(data,copula[[2]],log=F)
??? Tau_1 <- tau_1 / sum(tau_1 + tau_2)
??? Tau_2 <- tau_2 / sum(tau_1 + tau_2)
??? up_pi1 <- sum(Tau_1)/ 100
??? up_pi2 <- sum(Tau_2) / 100
??? # Tau <- c(Tau_1, Tau_2)
??? ll <- (-sum(Tau_1*dCopula(data,copula[[1]],log=T)))
??? #ll_1 <- -sum(Tau_2*dCopula(data,copula[[2]],log=T))
???
???
????if (is.finite(ll)) {
? ????return(ll)
??? } else {
????? if (is.na(ll)) {
?????...