Nothing
R/claytonExpr.R
In copula: Multivariate Dependence with Copulas
## This file is not to be edited. It is generated by ../inst/docs/mathnb/exprPrep.R
`claytonCopula.cdf.expr` <-
expression((u1^(-alpha))^(-1/alpha), (-1 + u1^(-alpha) + u2^(-alpha))^(-1/alpha),
(-2 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha))^(-1/alpha),
(-3 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha))^(-1/alpha),
(-4 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha))^(-1/alpha), (-5 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha))^(-1/alpha),
(-6 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha))^(-1/alpha),
(-7 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^(-1/alpha),
(-8 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha) +
u9^(-alpha))^(-1/alpha), (-9 + u1^(-alpha) + u10^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-1/alpha))
`claytonCopula.cdf.algr` <-
expression({
.value <- (u1^-alpha)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr1 <- -1
.expr2 <- -alpha
.value <- (.expr1 + u1^.expr2 + u2^.expr2)^(.expr1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-2 + u1^.expr2 + u2^.expr2 + u3^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-3 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-4 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-5 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-6 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2 + u7^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-7 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2 + u7^.expr2 + u8^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-8 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2 + u7^.expr2 + u8^.expr2 + u9^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-9 + u1^.expr2 + u10^.expr2 + u2^.expr2 + u3^.expr2 +
u4^.expr2 + u5^.expr2 + u6^.expr2 + u7^.expr2 + u8^.expr2 +
u9^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
})
`claytonCopula.cdfDerWrtArg.expr` <-
expression((u1^(-alpha))^((-1/alpha) - 1) * ((-1/alpha) * (u1^((-alpha) -
1) * (-alpha))), (-1 + u1^(-alpha) + u2^(-alpha))^((-1/alpha) -
1) * ((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))), (-2 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha))^((-1/alpha) - 1) *
((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))), (-3 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha))^((-1/alpha) - 1) *
((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))), (-4 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha))^((-1/alpha) -
1) * ((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))), (-5 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha))^((-1/alpha) - 1) * ((-1/alpha) * (u1^((-alpha) -
1) * (-alpha))), (-6 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha))^((-1/alpha) -
1) * ((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))), (-7 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^((-1/alpha) - 1) *
((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))), (-8 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^((-1/alpha) - 1) *
((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))), (-9 + u1^(-alpha) +
u10^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^((-1/alpha) -
1) * ((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))))
`claytonCopula.cdfDerWrtArg.algr` <-
expression({
.expr1 <- -alpha
.expr4 <- -1/alpha
.value <- (u1^.expr1)^(.expr4 - 1) * (.expr4 * (u1^(.expr1 -
1) * .expr1))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr1 <- -1
.expr2 <- -alpha
.expr7 <- .expr1/alpha
.value <- (.expr1 + u1^.expr2 + u2^.expr2)^(.expr7 - 1) *
(.expr7 * (u1^(.expr2 - 1) * .expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr10 <- -1/alpha
.value <- (-2 + u1^.expr2 + u2^.expr2 + u3^.expr2)^(.expr10 -
1) * (.expr10 * (u1^(.expr2 - 1) * .expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr12 <- -1/alpha
.value <- (-3 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2)^(.expr12 -
1) * (.expr12 * (u1^(.expr2 - 1) * .expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr14 <- -1/alpha
.value <- (-4 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2)^(.expr14 - 1) * (.expr14 * (u1^(.expr2 - 1) *
.expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr16 <- -1/alpha
.value <- (-5 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2)^(.expr16 - 1) * (.expr16 * (u1^(.expr2 -
1) * .expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr18 <- -1/alpha
.value <- (-6 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2 + u7^.expr2)^(.expr18 - 1) * (.expr18 *
(u1^(.expr2 - 1) * .expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr20 <- -1/alpha
.value <- (-7 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2 + u7^.expr2 + u8^.expr2)^(.expr20 -
1) * (.expr20 * (u1^(.expr2 - 1) * .expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr22 <- -1/alpha
.value <- (-8 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2 + u7^.expr2 + u8^.expr2 + u9^.expr2)^(.expr22 -
1) * (.expr22 * (u1^(.expr2 - 1) * .expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr24 <- -1/alpha
.value <- (-9 + u1^.expr2 + u10^.expr2 + u2^.expr2 + u3^.expr2 +
u4^.expr2 + u5^.expr2 + u6^.expr2 + u7^.expr2 + u8^.expr2 +
u9^.expr2)^(.expr24 - 1) * (.expr24 * (u1^(.expr2 - 1) *
.expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
})
`claytonCopula.cdfDerWrtPar.expr` <-
expression((u1^(-alpha))^(-1/alpha) * (log((u1^(-alpha))) * (1/alpha^2)) -
(u1^(-alpha))^((-1/alpha) - 1) * ((-1/alpha) * (u1^(-alpha) *
log(u1))), (-1 + u1^(-alpha) + u2^(-alpha))^(-1/alpha) *
(log((-1 + u1^(-alpha) + u2^(-alpha))) * (1/alpha^2)) - (-1 +
u1^(-alpha) + u2^(-alpha))^((-1/alpha) - 1) * ((-1/alpha) *
(u2^(-alpha) * log(u2) + u1^(-alpha) * log(u1))), (-2 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha))^(-1/alpha) * (log((-2 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha))) * (1/alpha^2)) - (-2 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha))^((-1/alpha) - 1) * ((-1/alpha) *
(u3^(-alpha) * log(u3) + (u2^(-alpha) * log(u2) + u1^(-alpha) *
log(u1)))), (-3 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha))^(-1/alpha) * (log((-3 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha))) * (1/alpha^2)) - (-3 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha))^((-1/alpha) - 1) *
((-1/alpha) * (u4^(-alpha) * log(u4) + (u3^(-alpha) * log(u3) +
(u2^(-alpha) * log(u2) + u1^(-alpha) * log(u1))))), (-4 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha))^(-1/alpha) *
(log((-4 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha))) * (1/alpha^2)) - (-4 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha))^((-1/alpha) - 1) *
((-1/alpha) * (u5^(-alpha) * log(u5) + (u4^(-alpha) * log(u4) +
(u3^(-alpha) * log(u3) + (u2^(-alpha) * log(u2) + u1^(-alpha) *
log(u1)))))), (-5 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha))^(-1/alpha) * (log((-5 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha))) * (1/alpha^2)) - (-5 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha))^((-1/alpha) -
1) * ((-1/alpha) * (u6^(-alpha) * log(u6) + (u5^(-alpha) *
log(u5) + (u4^(-alpha) * log(u4) + (u3^(-alpha) * log(u3) +
(u2^(-alpha) * log(u2) + u1^(-alpha) * log(u1))))))), (-6 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha))^(-1/alpha) * (log((-6 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha))) * (1/alpha^2)) - (-6 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha))^((-1/alpha) -
1) * ((-1/alpha) * (u7^(-alpha) * log(u7) + (u6^(-alpha) *
log(u6) + (u5^(-alpha) * log(u5) + (u4^(-alpha) * log(u4) +
(u3^(-alpha) * log(u3) + (u2^(-alpha) * log(u2) + u1^(-alpha) *
log(u1)))))))), (-7 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^(-1/alpha) *
(log((-7 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha))) *
(1/alpha^2)) - (-7 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^((-1/alpha) -
1) * ((-1/alpha) * (u8^(-alpha) * log(u8) + (u7^(-alpha) *
log(u7) + (u6^(-alpha) * log(u6) + (u5^(-alpha) * log(u5) +
(u4^(-alpha) * log(u4) + (u3^(-alpha) * log(u3) + (u2^(-alpha) *
log(u2) + u1^(-alpha) * log(u1))))))))), (-8 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-1/alpha) * (log((-8 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))) *
(1/alpha^2)) - (-8 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha) +
u9^(-alpha))^((-1/alpha) - 1) * ((-1/alpha) * (u9^(-alpha) *
log(u9) + (u8^(-alpha) * log(u8) + (u7^(-alpha) * log(u7) +
(u6^(-alpha) * log(u6) + (u5^(-alpha) * log(u5) + (u4^(-alpha) *
log(u4) + (u3^(-alpha) * log(u3) + (u2^(-alpha) * log(u2) +
u1^(-alpha) * log(u1)))))))))), (-9 + u1^(-alpha) + u10^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-1/alpha) * (log((-9 +
u1^(-alpha) + u10^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha) +
u9^(-alpha))) * (1/alpha^2)) - (-9 + u1^(-alpha) + u10^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^((-1/alpha) - 1) *
((-1/alpha) * (u9^(-alpha) * log(u9) + (u8^(-alpha) * log(u8) +
(u7^(-alpha) * log(u7) + (u6^(-alpha) * log(u6) + (u5^(-alpha) *
log(u5) + (u4^(-alpha) * log(u4) + (u3^(-alpha) *
log(u3) + (u2^(-alpha) * log(u2) + (u10^(-alpha) *
log(u10) + u1^(-alpha) * log(u1))))))))))))
`claytonCopula.cdfDerWrtPar.algr` <-
expression({
.expr2 <- u1^-alpha
.expr4 <- -1/alpha
.value <- .expr2^.expr4 * (log(.expr2) * (1/alpha^2)) - .expr2^(.expr4 -
1) * (.expr4 * (.expr2 * log(u1)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr1 <- -1
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr6 <- .expr1 + .expr3 + .expr5
.expr7 <- .expr1/alpha
.value <- .expr6^.expr7 * (log(.expr6) * (1/alpha^2)) - .expr6^(.expr7 -
1) * (.expr7 * (.expr5 * log(u2) + .expr3 * log(u1)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr7 <- u3^.expr2
.expr8 <- -2 + .expr3 + .expr5 + .expr7
.expr10 <- -1/alpha
.value <- .expr8^.expr10 * (log(.expr8) * (1/alpha^2)) -
.expr8^(.expr10 - 1) * (.expr10 * (.expr7 * log(u3) +
(.expr5 * log(u2) + .expr3 * log(u1))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr7 <- u3^.expr2
.expr9 <- u4^.expr2
.expr10 <- -3 + .expr3 + .expr5 + .expr7 + .expr9
.expr12 <- -1/alpha
.value <- .expr10^.expr12 * (log(.expr10) * (1/alpha^2)) -
.expr10^(.expr12 - 1) * (.expr12 * (.expr9 * log(u4) +
(.expr7 * log(u3) + (.expr5 * log(u2) + .expr3 *
log(u1)))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr7 <- u3^.expr2
.expr9 <- u4^.expr2
.expr11 <- u5^.expr2
.expr12 <- -4 + .expr3 + .expr5 + .expr7 + .expr9 + .expr11
.expr14 <- -1/alpha
.value <- .expr12^.expr14 * (log(.expr12) * (1/alpha^2)) -
.expr12^(.expr14 - 1) * (.expr14 * (.expr11 * log(u5) +
(.expr9 * log(u4) + (.expr7 * log(u3) + (.expr5 *
log(u2) + .expr3 * log(u1))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr7 <- u3^.expr2
.expr9 <- u4^.expr2
.expr11 <- u5^.expr2
.expr13 <- u6^.expr2
.expr14 <- -5 + .expr3 + .expr5 + .expr7 + .expr9 + .expr11 +
.expr13
.expr16 <- -1/alpha
.value <- .expr14^.expr16 * (log(.expr14) * (1/alpha^2)) -
.expr14^(.expr16 - 1) * (.expr16 * (.expr13 * log(u6) +
(.expr11 * log(u5) + (.expr9 * log(u4) + (.expr7 *
log(u3) + (.expr5 * log(u2) + .expr3 * log(u1)))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr7 <- u3^.expr2
.expr9 <- u4^.expr2
.expr11 <- u5^.expr2
.expr13 <- u6^.expr2
.expr15 <- u7^.expr2
.expr16 <- -6 + .expr3 + .expr5 + .expr7 + .expr9 + .expr11 +
.expr13 + .expr15
.expr18 <- -1/alpha
.value <- .expr16^.expr18 * (log(.expr16) * (1/alpha^2)) -
.expr16^(.expr18 - 1) * (.expr18 * (.expr15 * log(u7) +
(.expr13 * log(u6) + (.expr11 * log(u5) + (.expr9 *
log(u4) + (.expr7 * log(u3) + (.expr5 * log(u2) +
.expr3 * log(u1))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr7 <- u3^.expr2
.expr9 <- u4^.expr2
.expr11 <- u5^.expr2
.expr13 <- u6^.expr2
.expr15 <- u7^.expr2
.expr17 <- u8^.expr2
.expr18 <- -7 + .expr3 + .expr5 + .expr7 + .expr9 + .expr11 +
.expr13 + .expr15 + .expr17
.expr20 <- -1/alpha
.value <- .expr18^.expr20 * (log(.expr18) * (1/alpha^2)) -
.expr18^(.expr20 - 1) * (.expr20 * (.expr17 * log(u8) +
(.expr15 * log(u7) + (.expr13 * log(u6) + (.expr11 *
log(u5) + (.expr9 * log(u4) + (.expr7 * log(u3) +
(.expr5 * log(u2) + .expr3 * log(u1)))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr7 <- u3^.expr2
.expr9 <- u4^.expr2
.expr11 <- u5^.expr2
.expr13 <- u6^.expr2
.expr15 <- u7^.expr2
.expr17 <- u8^.expr2
.expr19 <- u9^.expr2
.expr20 <- -8 + .expr3 + .expr5 + .expr7 + .expr9 + .expr11 +
.expr13 + .expr15 + .expr17 + .expr19
.expr22 <- -1/alpha
.value <- .expr20^.expr22 * (log(.expr20) * (1/alpha^2)) -
.expr20^(.expr22 - 1) * (.expr22 * (.expr19 * log(u9) +
(.expr17 * log(u8) + (.expr15 * log(u7) + (.expr13 *
log(u6) + (.expr11 * log(u5) + (.expr9 * log(u4) +
(.expr7 * log(u3) + (.expr5 * log(u2) + .expr3 *
log(u1))))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u10^.expr2
.expr7 <- u2^.expr2
.expr9 <- u3^.expr2
.expr11 <- u4^.expr2
.expr13 <- u5^.expr2
.expr15 <- u6^.expr2
.expr17 <- u7^.expr2
.expr19 <- u8^.expr2
.expr21 <- u9^.expr2
.expr22 <- -9 + .expr3 + .expr5 + .expr7 + .expr9 + .expr11 +
.expr13 + .expr15 + .expr17 + .expr19 + .expr21
.expr24 <- -1/alpha
.value <- .expr22^.expr24 * (log(.expr22) * (1/alpha^2)) -
.expr22^(.expr24 - 1) * (.expr24 * (.expr21 * log(u9) +
(.expr19 * log(u8) + (.expr17 * log(u7) + (.expr15 *
log(u6) + (.expr13 * log(u5) + (.expr11 * log(u4) +
(.expr9 * log(u3) + (.expr7 * log(u2) + (.expr5 *
log(u10) + .expr3 * log(u1)))))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
})
`claytonCopula.pdf.expr` <-
expression(1/(u1 * (u1^(-alpha))^(1/alpha)), (1 + alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * (-1 + u1^(-alpha) + u2^(-alpha))^(-2 -
1/alpha), (1 + alpha) * (1 + 2 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * (-2 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha))^(-3 - 1/alpha), (1 + alpha) * (1 + 2 * alpha) *
(1 + 3 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 -
alpha) * u4^(-1 - alpha) * (-3 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha))^(-4 - 1/alpha), (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * (-4 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha))^(-5 - 1/alpha), (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 +
5 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 -
alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) *
(-5 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha))^(-6 - 1/alpha), (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 +
5 * alpha) * (1 + 6 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) *
u6^(-1 - alpha) * u7^(-1 - alpha) * (-6 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha))^(-7 -
1/alpha), (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * (1 +
7 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 -
alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) *
u7^(-1 - alpha) * u8^(-1 - alpha) * (-7 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) +
u8^(-alpha))^(-8 - 1/alpha), (1 + alpha) * (1 + 2 * alpha) *
(1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * (1 + 7 * alpha) * (1 + 8 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) * u8^(-1 -
alpha) * u9^(-1 - alpha) * (-8 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) +
u8^(-alpha) + u9^(-alpha))^(-9 - 1/alpha), (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 +
5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) * (1 + 8 *
alpha) * (1 + 9 * alpha) * u1^(-1 - alpha) * u10^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 -
alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) * u8^(-1 - alpha) *
u9^(-1 - alpha) * (-9 + u1^(-alpha) + u10^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) +
u8^(-alpha) + u9^(-alpha))^(-10 - 1/alpha))
`claytonCopula.pdf.algr` <-
expression({
.value <- 1/(u1 * (u1^-alpha)^(1/alpha))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -1
.expr3 <- .expr2 - alpha
.expr8 <- -alpha
.value <- (1 + alpha) * u1^.expr3 * u2^.expr3 * (.expr2 +
u1^.expr8 + u2^.expr8)^(-2 - 1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr6 <- -1 - alpha
.expr14 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * u1^.expr6 * u2^.expr6 *
u3^.expr6 * (-2 + u1^.expr14 + u2^.expr14 + u3^.expr14)^(-3 -
1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr9 <- -1 - alpha
.expr19 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
u1^.expr9 * u2^.expr9 * u3^.expr9 * u4^.expr9 * (-3 +
u1^.expr19 + u2^.expr19 + u3^.expr19 + u4^.expr19)^(-4 -
1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr12 <- -1 - alpha
.expr24 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * u1^.expr12 * u2^.expr12 * u3^.expr12 *
u4^.expr12 * u5^.expr12 * (-4 + u1^.expr24 + u2^.expr24 +
u3^.expr24 + u4^.expr24 + u5^.expr24)^(-5 - 1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr15 <- -1 - alpha
.expr29 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * u1^.expr15 * u2^.expr15 *
u3^.expr15 * u4^.expr15 * u5^.expr15 * u6^.expr15 * (-5 +
u1^.expr29 + u2^.expr29 + u3^.expr29 + u4^.expr29 + u5^.expr29 +
u6^.expr29)^(-6 - 1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr18 <- -1 - alpha
.expr34 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
u1^.expr18 * u2^.expr18 * u3^.expr18 * u4^.expr18 * u5^.expr18 *
u6^.expr18 * u7^.expr18 * (-6 + u1^.expr34 + u2^.expr34 +
u3^.expr34 + u4^.expr34 + u5^.expr34 + u6^.expr34 + u7^.expr34)^(-7 -
1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr21 <- -1 - alpha
.expr39 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * u1^.expr21 * u2^.expr21 * u3^.expr21 *
u4^.expr21 * u5^.expr21 * u6^.expr21 * u7^.expr21 * u8^.expr21 *
(-7 + u1^.expr39 + u2^.expr39 + u3^.expr39 + u4^.expr39 +
u5^.expr39 + u6^.expr39 + u7^.expr39 + u8^.expr39)^(-8 -
1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr24 <- -1 - alpha
.expr44 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * u1^.expr24 * u2^.expr24 *
u3^.expr24 * u4^.expr24 * u5^.expr24 * u6^.expr24 * u7^.expr24 *
u8^.expr24 * u9^.expr24 * (-8 + u1^.expr44 + u2^.expr44 +
u3^.expr44 + u4^.expr44 + u5^.expr44 + u6^.expr44 + u7^.expr44 +
u8^.expr44 + u9^.expr44)^(-9 - 1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr27 <- -1 - alpha
.expr49 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * (1 + 9 * alpha) *
u1^.expr27 * u10^.expr27 * u2^.expr27 * u3^.expr27 *
u4^.expr27 * u5^.expr27 * u6^.expr27 * u7^.expr27 * u8^.expr27 *
u9^.expr27 * (-9 + u1^.expr49 + u10^.expr49 + u2^.expr49 +
u3^.expr49 + u4^.expr49 + u5^.expr49 + u6^.expr49 + u7^.expr49 +
u8^.expr49 + u9^.expr49)^(-10 - 1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
})
`claytonCopula.pdfDerWrtArg.expr` <-
expression(-(((u1^(-alpha))^(1/alpha) + u1 * ((u1^(-alpha))^((1/alpha) -
1) * ((1/alpha) * (u1^((-alpha) - 1) * (-alpha)))))/(u1 *
(u1^(-alpha))^(1/alpha))^2), (1 + alpha) * (u1^((-1 - alpha) -
1) * (-1 - alpha)) * u2^(-1 - alpha) * (-1 + u1^(-alpha) +
u2^(-alpha))^(-2 - 1/alpha) + (1 + alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * ((-1 + u1^(-alpha) + u2^(-alpha))^((-2 -
1/alpha) - 1) * ((-2 - 1/alpha) * (u1^((-alpha) - 1) * (-alpha)))),
(1 + alpha) * (1 + 2 * alpha) * (u1^((-1 - alpha) - 1) *
(-1 - alpha)) * u2^(-1 - alpha) * u3^(-1 - alpha) * (-2 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha))^(-3 - 1/alpha) +
(1 + alpha) * (1 + 2 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * ((-2 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha))^((-3 - 1/alpha) - 1) * ((-3 - 1/alpha) *
(u1^((-alpha) - 1) * (-alpha)))), (1 + alpha) * (1 +
2 * alpha) * (1 + 3 * alpha) * (u1^((-1 - alpha) - 1) *
(-1 - alpha)) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 -
alpha) * (-3 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha))^(-4 - 1/alpha) + (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * ((-3 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha))^((-4 - 1/alpha) -
1) * ((-4 - 1/alpha) * (u1^((-alpha) - 1) * (-alpha)))),
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (u1^((-1 - alpha) - 1) * (-1 - alpha)) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 -
alpha) * (-4 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha))^(-5 - 1/alpha) + (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) *
u4^(-1 - alpha) * u5^(-1 - alpha) * ((-4 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha))^((-5 -
1/alpha) - 1) * ((-5 - 1/alpha) * (u1^((-alpha) - 1) *
(-alpha)))), (1 + alpha) * (1 + 2 * alpha) * (1 + 3 *
alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (u1^((-1 -
alpha) - 1) * (-1 - alpha)) * u2^(-1 - alpha) * u3^(-1 -
alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 -
alpha) * (-5 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha))^(-6 - 1/alpha) +
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 +
4 * alpha) * (1 + 5 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * u6^(-1 - alpha) * ((-5 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha))^((-6 - 1/alpha) - 1) * ((-6 - 1/alpha) *
(u1^((-alpha) - 1) * (-alpha)))), (1 + alpha) * (1 +
2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 +
5 * alpha) * (1 + 6 * alpha) * (u1^((-1 - alpha) - 1) *
(-1 - alpha)) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 -
alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 -
alpha) * (-6 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha))^(-7 -
1/alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) *
u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) *
u7^(-1 - alpha) * ((-6 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha))^((-7 - 1/alpha) - 1) * ((-7 - 1/alpha) *
(u1^((-alpha) - 1) * (-alpha)))), (1 + alpha) * (1 +
2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 +
5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) * (u1^((-1 -
alpha) - 1) * (-1 - alpha)) * u2^(-1 - alpha) * u3^(-1 -
alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 -
alpha) * u7^(-1 - alpha) * u8^(-1 - alpha) * (-7 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^(-8 - 1/alpha) +
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 +
4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) *
u6^(-1 - alpha) * u7^(-1 - alpha) * u8^(-1 - alpha) *
((-7 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) +
u8^(-alpha))^((-8 - 1/alpha) - 1) * ((-8 - 1/alpha) *
(u1^((-alpha) - 1) * (-alpha)))), (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
(1 + 8 * alpha) * (u1^((-1 - alpha) - 1) * (-1 - alpha)) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) *
u8^(-1 - alpha) * u9^(-1 - alpha) * (-8 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-9 -
1/alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) *
u8^(-1 - alpha) * u9^(-1 - alpha) * ((-8 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^((-9 -
1/alpha) - 1) * ((-9 - 1/alpha) * (u1^((-alpha) - 1) *
(-alpha)))), (1 + alpha) * (1 + 2 * alpha) * (1 + 3 *
alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 *
alpha) * (1 + 7 * alpha) * (1 + 8 * alpha) * (1 + 9 *
alpha) * (u1^((-1 - alpha) - 1) * (-1 - alpha)) * u10^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 -
alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 -
alpha) * u8^(-1 - alpha) * u9^(-1 - alpha) * (-9 + u1^(-alpha) +
u10^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha) +
u9^(-alpha))^(-10 - 1/alpha) + (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) * (1 + 8 *
alpha) * (1 + 9 * alpha) * u1^(-1 - alpha) * u10^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 -
alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 -
alpha) * u8^(-1 - alpha) * u9^(-1 - alpha) * ((-9 + u1^(-alpha) +
u10^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha) +
u9^(-alpha))^((-10 - 1/alpha) - 1) * ((-10 - 1/alpha) *
(u1^((-alpha) - 1) * (-alpha)))))
`claytonCopula.pdfDerWrtArg.algr` <-
expression({
.expr1 <- -alpha
.expr2 <- u1^.expr1
.expr3 <- 1/alpha
.expr4 <- .expr2^.expr3
.value <- -((.expr4 + u1 * (.expr2^(.expr3 - 1) * (.expr3 *
(u1^(.expr1 - 1) * .expr1))))/(u1 * .expr4)^2)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr1 <- 1 + alpha
.expr2 <- -1
.expr3 <- .expr2 - alpha
.expr8 <- u2^.expr3
.expr10 <- -alpha
.expr14 <- .expr2 + u1^.expr10 + u2^.expr10
.expr17 <- -2 - 1/alpha
.value <- .expr1 * (u1^(.expr3 - 1) * .expr3) * .expr8 *
.expr14^.expr17 + .expr1 * u1^.expr3 * .expr8 * (.expr14^(.expr17 -
1) * (.expr17 * (u1^(.expr10 - 1) * .expr10)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr4 <- (1 + alpha) * (1 + 2 * alpha)
.expr6 <- -1 - alpha
.expr11 <- u2^.expr6
.expr13 <- u3^.expr6
.expr16 <- -alpha
.expr22 <- -2 + u1^.expr16 + u2^.expr16 + u3^.expr16
.expr25 <- -3 - 1/alpha
.value <- .expr4 * (u1^(.expr6 - 1) * .expr6) * .expr11 *
.expr13 * .expr22^.expr25 + .expr4 * u1^.expr6 * .expr11 *
.expr13 * (.expr22^(.expr25 - 1) * (.expr25 * (u1^(.expr16 -
1) * .expr16)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr7 <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha)
.expr9 <- -1 - alpha
.expr14 <- u2^.expr9
.expr16 <- u3^.expr9
.expr18 <- u4^.expr9
.expr21 <- -alpha
.expr29 <- -3 + u1^.expr21 + u2^.expr21 + u3^.expr21 + u4^.expr21
.expr32 <- -4 - 1/alpha
.value <- .expr7 * (u1^(.expr9 - 1) * .expr9) * .expr14 *
.expr16 * .expr18 * .expr29^.expr32 + .expr7 * u1^.expr9 *
.expr14 * .expr16 * .expr18 * (.expr29^(.expr32 - 1) *
(.expr32 * (u1^(.expr21 - 1) * .expr21)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr10 <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha)
.expr12 <- -1 - alpha
.expr17 <- u2^.expr12
.expr19 <- u3^.expr12
.expr21 <- u4^.expr12
.expr23 <- u5^.expr12
.expr26 <- -alpha
.expr36 <- -4 + u1^.expr26 + u2^.expr26 + u3^.expr26 + u4^.expr26 +
u5^.expr26
.expr39 <- -5 - 1/alpha
.value <- .expr10 * (u1^(.expr12 - 1) * .expr12) * .expr17 *
.expr19 * .expr21 * .expr23 * .expr36^.expr39 + .expr10 *
u1^.expr12 * .expr17 * .expr19 * .expr21 * .expr23 *
(.expr36^(.expr39 - 1) * (.expr39 * (u1^(.expr26 - 1) *
.expr26)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr13 <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha)
.expr15 <- -1 - alpha
.expr20 <- u2^.expr15
.expr22 <- u3^.expr15
.expr24 <- u4^.expr15
.expr26 <- u5^.expr15
.expr28 <- u6^.expr15
.expr31 <- -alpha
.expr43 <- -5 + u1^.expr31 + u2^.expr31 + u3^.expr31 + u4^.expr31 +
u5^.expr31 + u6^.expr31
.expr46 <- -6 - 1/alpha
.value <- .expr13 * (u1^(.expr15 - 1) * .expr15) * .expr20 *
.expr22 * .expr24 * .expr26 * .expr28 * .expr43^.expr46 +
.expr13 * u1^.expr15 * .expr20 * .expr22 * .expr24 *
.expr26 * .expr28 * (.expr43^(.expr46 - 1) * (.expr46 *
(u1^(.expr31 - 1) * .expr31)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr16 <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha)
.expr18 <- -1 - alpha
.expr23 <- u2^.expr18
.expr25 <- u3^.expr18
.expr27 <- u4^.expr18
.expr29 <- u5^.expr18
.expr31 <- u6^.expr18
.expr33 <- u7^.expr18
.expr36 <- -alpha
.expr50 <- -6 + u1^.expr36 + u2^.expr36 + u3^.expr36 + u4^.expr36 +
u5^.expr36 + u6^.expr36 + u7^.expr36
.expr53 <- -7 - 1/alpha
.value <- .expr16 * (u1^(.expr18 - 1) * .expr18) * .expr23 *
.expr25 * .expr27 * .expr29 * .expr31 * .expr33 * .expr50^.expr53 +
.expr16 * u1^.expr18 * .expr23 * .expr25 * .expr27 *
.expr29 * .expr31 * .expr33 * (.expr50^(.expr53 -
1) * (.expr53 * (u1^(.expr36 - 1) * .expr36)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr19 <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha)
.expr21 <- -1 - alpha
.expr26 <- u2^.expr21
.expr28 <- u3^.expr21
.expr30 <- u4^.expr21
.expr32 <- u5^.expr21
.expr34 <- u6^.expr21
.expr36 <- u7^.expr21
.expr38 <- u8^.expr21
.expr41 <- -alpha
.expr57 <- -7 + u1^.expr41 + u2^.expr41 + u3^.expr41 + u4^.expr41 +
u5^.expr41 + u6^.expr41 + u7^.expr41 + u8^.expr41
.expr60 <- -8 - 1/alpha
.value <- .expr19 * (u1^(.expr21 - 1) * .expr21) * .expr26 *
.expr28 * .expr30 * .expr32 * .expr34 * .expr36 * .expr38 *
.expr57^.expr60 + .expr19 * u1^.expr21 * .expr26 * .expr28 *
.expr30 * .expr32 * .expr34 * .expr36 * .expr38 * (.expr57^(.expr60 -
1) * (.expr60 * (u1^(.expr41 - 1) * .expr41)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr22 <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha)
.expr24 <- -1 - alpha
.expr29 <- u2^.expr24
.expr31 <- u3^.expr24
.expr33 <- u4^.expr24
.expr35 <- u5^.expr24
.expr37 <- u6^.expr24
.expr39 <- u7^.expr24
.expr41 <- u8^.expr24
.expr43 <- u9^.expr24
.expr46 <- -alpha
.expr64 <- -8 + u1^.expr46 + u2^.expr46 + u3^.expr46 + u4^.expr46 +
u5^.expr46 + u6^.expr46 + u7^.expr46 + u8^.expr46 + u9^.expr46
.expr67 <- -9 - 1/alpha
.value <- .expr22 * (u1^(.expr24 - 1) * .expr24) * .expr29 *
.expr31 * .expr33 * .expr35 * .expr37 * .expr39 * .expr41 *
.expr43 * .expr64^.expr67 + .expr22 * u1^.expr24 * .expr29 *
.expr31 * .expr33 * .expr35 * .expr37 * .expr39 * .expr41 *
.expr43 * (.expr64^(.expr67 - 1) * (.expr67 * (u1^(.expr46 -
1) * .expr46)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr25 <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * (1 + 9 * alpha)
.expr27 <- -1 - alpha
.expr32 <- u10^.expr27
.expr34 <- u2^.expr27
.expr36 <- u3^.expr27
.expr38 <- u4^.expr27
.expr40 <- u5^.expr27
.expr42 <- u6^.expr27
.expr44 <- u7^.expr27
.expr46 <- u8^.expr27
.expr48 <- u9^.expr27
.expr51 <- -alpha
.expr71 <- -9 + u1^.expr51 + u10^.expr51 + u2^.expr51 + u3^.expr51 +
u4^.expr51 + u5^.expr51 + u6^.expr51 + u7^.expr51 + u8^.expr51 +
u9^.expr51
.expr74 <- -10 - 1/alpha
.value <- .expr25 * (u1^(.expr27 - 1) * .expr27) * .expr32 *
.expr34 * .expr36 * .expr38 * .expr40 * .expr42 * .expr44 *
.expr46 * .expr48 * .expr71^.expr74 + .expr25 * u1^.expr27 *
.expr32 * .expr34 * .expr36 * .expr38 * .expr40 * .expr42 *
.expr44 * .expr46 * .expr48 * (.expr71^(.expr74 - 1) *
(.expr74 * (u1^(.expr51 - 1) * .expr51)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
})
`claytonCopula.pdfDerWrtPar.expr` <-
expression(u1 * ((u1^(-alpha))^(1/alpha) * (log((u1^(-alpha))) *
(1/alpha^2)) + (u1^(-alpha))^((1/alpha) - 1) * ((1/alpha) *
(u1^(-alpha) * log(u1))))/(u1 * (u1^(-alpha))^(1/alpha))^2,
((u1^(-1 - alpha) - (1 + alpha) * (u1^(-1 - alpha) * log(u1))) *
u2^(-1 - alpha) - (1 + alpha) * u1^(-1 - alpha) * (u2^(-1 -
alpha) * log(u2))) * (-1 + u1^(-alpha) + u2^(-alpha))^(-2 -
1/alpha) + (1 + alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
((-1 + u1^(-alpha) + u2^(-alpha))^(-2 - 1/alpha) * (log((-1 +
u1^(-alpha) + u2^(-alpha))) * (1/alpha^2)) - (-1 +
u1^(-alpha) + u2^(-alpha))^((-2 - 1/alpha) - 1) *
((-2 - 1/alpha) * (u2^(-alpha) * log(u2) + u1^(-alpha) *
log(u1)))), (((((1 + 2 * alpha) + (1 + alpha) *
2) * u1^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) *
(u1^(-1 - alpha) * log(u1))) * u2^(-1 - alpha) - (1 +
alpha) * (1 + 2 * alpha) * u1^(-1 - alpha) * (u2^(-1 -
alpha) * log(u2))) * u3^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
(u3^(-1 - alpha) * log(u3))) * (-2 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha))^(-3 - 1/alpha) + (1 + alpha) * (1 + 2 *
alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 -
alpha) * ((-2 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha))^(-3 -
1/alpha) * (log((-2 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha))) *
(1/alpha^2)) - (-2 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha))^((-3 -
1/alpha) - 1) * ((-3 - 1/alpha) * (u3^(-alpha) * log(u3) +
(u2^(-alpha) * log(u2) + u1^(-alpha) * log(u1))))), (((((((1 +
2 * alpha) + (1 + alpha) * 2) * (1 + 3 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * 3) * u1^(-1 - alpha) - (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (u1^(-1 -
alpha) * log(u1))) * u2^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * u1^(-1 - alpha) *
(u2^(-1 - alpha) * log(u2))) * u3^(-1 - alpha) - (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * (u3^(-1 - alpha) * log(u3))) *
u4^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 -
alpha) * (u4^(-1 - alpha) * log(u4))) * (-3 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha))^(-4 - 1/alpha) +
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 -
alpha) * ((-3 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha))^(-4 - 1/alpha) * (log((-3 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha))) * (1/alpha^2)) -
(-3 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha))^((-4 -
1/alpha) - 1) * ((-4 - 1/alpha) * (u4^(-alpha) *
log(u4) + (u3^(-alpha) * log(u3) + (u2^(-alpha) *
log(u2) + u1^(-alpha) * log(u1)))))), (((((((((1 +
2 * alpha) + (1 + alpha) * 2) * (1 + 3 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * 3) * (1 + 4 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * 4) * u1^(-1 -
alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (u1^(-1 - alpha) * log(u1))) * u2^(-1 -
alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * u1^(-1 - alpha) * (u2^(-1 - alpha) *
log(u2))) * u3^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * (u3^(-1 - alpha) * log(u3))) *
u4^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * (u4^(-1 - alpha) * log(u4))) *
u5^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * (u5^(-1 -
alpha) * log(u5))) * (-4 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha))^(-5 - 1/alpha) +
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 +
4 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) *
((-4 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha))^(-5 - 1/alpha) * (log((-4 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha))) * (1/alpha^2)) - (-4 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha))^((-5 -
1/alpha) - 1) * ((-5 - 1/alpha) * (u5^(-alpha) *
log(u5) + (u4^(-alpha) * log(u4) + (u3^(-alpha) *
log(u3) + (u2^(-alpha) * log(u2) + u1^(-alpha) *
log(u1))))))), (((((((((((1 + 2 * alpha) + (1 +
alpha) * 2) * (1 + 3 * alpha) + (1 + alpha) * (1 + 2 *
alpha) * 3) * (1 + 4 * alpha) + (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * 4) * (1 + 5 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * 5) * u1^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (u1^(-1 - alpha) * log(u1))) * u2^(-1 - alpha) -
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 +
4 * alpha) * (1 + 5 * alpha) * u1^(-1 - alpha) *
(u2^(-1 - alpha) * log(u2))) * u3^(-1 - alpha) -
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 +
4 * alpha) * (1 + 5 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * (u3^(-1 - alpha) * log(u3))) *
u4^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * (u4^(-1 -
alpha) * log(u4))) * u5^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * (u5^(-1 - alpha) *
log(u5))) * u6^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 -
alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) * (u6^(-1 -
alpha) * log(u6))) * (-5 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha))^(-6 -
1/alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * u6^(-1 - alpha) * ((-5 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha))^(-6 - 1/alpha) * (log((-5 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha))) * (1/alpha^2)) - (-5 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha))^((-6 -
1/alpha) - 1) * ((-6 - 1/alpha) * (u6^(-alpha) * log(u6) +
(u5^(-alpha) * log(u5) + (u4^(-alpha) * log(u4) + (u3^(-alpha) *
log(u3) + (u2^(-alpha) * log(u2) + u1^(-alpha) *
log(u1)))))))), (((((((((((((1 + 2 * alpha) + (1 +
alpha) * 2) * (1 + 3 * alpha) + (1 + alpha) * (1 + 2 *
alpha) * 3) * (1 + 4 * alpha) + (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * 4) * (1 + 5 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * 5) * (1 + 6 * alpha) + (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * 6) * u1^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * (u1^(-1 - alpha) * log(u1))) *
u2^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * u1^(-1 - alpha) * (u2^(-1 - alpha) * log(u2))) *
u3^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * (u3^(-1 -
alpha) * log(u3))) * u4^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * (u4^(-1 - alpha) *
log(u4))) * u5^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * (u5^(-1 -
alpha) * log(u5))) * u6^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * (u6^(-1 - alpha) * log(u6))) * u7^(-1 -
alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) *
u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) *
(u7^(-1 - alpha) * log(u7))) * (-6 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha))^(-7 - 1/alpha) + (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 -
alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) * ((-6 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha))^(-7 - 1/alpha) * (log((-6 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha))) * (1/alpha^2)) -
(-6 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha))^((-7 - 1/alpha) -
1) * ((-7 - 1/alpha) * (u7^(-alpha) * log(u7) + (u6^(-alpha) *
log(u6) + (u5^(-alpha) * log(u5) + (u4^(-alpha) *
log(u4) + (u3^(-alpha) * log(u3) + (u2^(-alpha) *
log(u2) + u1^(-alpha) * log(u1))))))))), (((((((((((((((1 +
2 * alpha) + (1 + alpha) * 2) * (1 + 3 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * 3) * (1 + 4 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * 4) * (1 +
5 * alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 *
alpha) * (1 + 4 * alpha) * 5) * (1 + 6 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * 6) * (1 + 7 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * 7) * u1^(-1 -
alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (u1^(-1 - alpha) * log(u1))) * u2^(-1 -
alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * u1^(-1 - alpha) * (u2^(-1 - alpha) *
log(u2))) * u3^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * (u3^(-1 - alpha) * log(u3))) *
u4^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * (1 + 7 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * (u4^(-1 - alpha) * log(u4))) *
u5^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * (1 + 7 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * (u5^(-1 -
alpha) * log(u5))) * u6^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) *
u4^(-1 - alpha) * u5^(-1 - alpha) * (u6^(-1 - alpha) *
log(u6))) * u7^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 -
alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) * (u7^(-1 -
alpha) * log(u7))) * u8^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) *
u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) *
u7^(-1 - alpha) * (u8^(-1 - alpha) * log(u8))) * (-7 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^(-8 -
1/alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) *
u6^(-1 - alpha) * u7^(-1 - alpha) * u8^(-1 - alpha) *
((-7 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^(-8 -
1/alpha) * (log((-7 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha) + u8^(-alpha))) * (1/alpha^2)) - (-7 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^((-8 -
1/alpha) - 1) * ((-8 - 1/alpha) * (u8^(-alpha) *
log(u8) + (u7^(-alpha) * log(u7) + (u6^(-alpha) *
log(u6) + (u5^(-alpha) * log(u5) + (u4^(-alpha) *
log(u4) + (u3^(-alpha) * log(u3) + (u2^(-alpha) *
log(u2) + u1^(-alpha) * log(u1)))))))))), (((((((((((((((((1 +
2 * alpha) + (1 + alpha) * 2) * (1 + 3 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * 3) * (1 + 4 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * 4) * (1 +
5 * alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 *
alpha) * (1 + 4 * alpha) * 5) * (1 + 6 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * 6) * (1 + 7 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * 7) * (1 +
8 * alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 *
alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 *
alpha) * (1 + 7 * alpha) * 8) * u1^(-1 - alpha) - (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 *
alpha) * (1 + 8 * alpha) * (u1^(-1 - alpha) * log(u1))) *
u2^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * (1 + 7 * alpha) * (1 + 8 * alpha) * u1^(-1 -
alpha) * (u2^(-1 - alpha) * log(u2))) * u3^(-1 - alpha) -
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 +
4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * (u3^(-1 - alpha) * log(u3))) *
u4^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * (1 + 7 * alpha) * (1 + 8 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * (u4^(-1 -
alpha) * log(u4))) * u5^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
(1 + 8 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * (u5^(-1 - alpha) *
log(u5))) * u6^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) * (1 + 8 *
alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 -
alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) * (u6^(-1 -
alpha) * log(u6))) * u7^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
(1 + 8 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) *
u6^(-1 - alpha) * (u7^(-1 - alpha) * log(u7))) * u8^(-1 -
alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) *
(u8^(-1 - alpha) * log(u8))) * u9^(-1 - alpha) - (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 *
alpha) * (1 + 8 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 -
alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) * u8^(-1 -
alpha) * (u9^(-1 - alpha) * log(u9))) * (-8 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-9 -
1/alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) *
u8^(-1 - alpha) * u9^(-1 - alpha) * ((-8 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-9 -
1/alpha) * (log((-8 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) +
u8^(-alpha) + u9^(-alpha))) * (1/alpha^2)) - (-8 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^((-9 -
1/alpha) - 1) * ((-9 - 1/alpha) * (u9^(-alpha) * log(u9) +
(u8^(-alpha) * log(u8) + (u7^(-alpha) * log(u7) + (u6^(-alpha) *
log(u6) + (u5^(-alpha) * log(u5) + (u4^(-alpha) *
log(u4) + (u3^(-alpha) * log(u3) + (u2^(-alpha) *
log(u2) + u1^(-alpha) * log(u1))))))))))), (((((((((((((((((((1 +
2 * alpha) + (1 + alpha) * 2) * (1 + 3 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * 3) * (1 + 4 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * 4) * (1 +
5 * alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 *
alpha) * (1 + 4 * alpha) * 5) * (1 + 6 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * 6) * (1 + 7 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * 7) * (1 +
8 * alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 *
alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 *
alpha) * (1 + 7 * alpha) * 8) * (1 + 9 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 *
alpha) * (1 + 8 * alpha) * 9) * u1^(-1 - alpha) - (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 *
alpha) * (1 + 8 * alpha) * (1 + 9 * alpha) * (u1^(-1 -
alpha) * log(u1))) * u10^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
(1 + 8 * alpha) * (1 + 9 * alpha) * u1^(-1 - alpha) *
(u10^(-1 - alpha) * log(u10))) * u2^(-1 - alpha) - (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 *
alpha) * (1 + 8 * alpha) * (1 + 9 * alpha) * u1^(-1 -
alpha) * u10^(-1 - alpha) * (u2^(-1 - alpha) * log(u2))) *
u3^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * (1 + 7 * alpha) * (1 + 8 * alpha) * (1 +
9 * alpha) * u1^(-1 - alpha) * u10^(-1 - alpha) * u2^(-1 -
alpha) * (u3^(-1 - alpha) * log(u3))) * u4^(-1 - alpha) -
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 +
4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * (1 + 9 * alpha) *
u1^(-1 - alpha) * u10^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * (u4^(-1 - alpha) * log(u4))) *
u5^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * (1 + 7 * alpha) * (1 + 8 * alpha) * (1 +
9 * alpha) * u1^(-1 - alpha) * u10^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * (u5^(-1 -
alpha) * log(u5))) * u6^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
(1 + 8 * alpha) * (1 + 9 * alpha) * u1^(-1 - alpha) *
u10^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) *
u4^(-1 - alpha) * u5^(-1 - alpha) * (u6^(-1 - alpha) *
log(u6))) * u7^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) * (1 + 8 *
alpha) * (1 + 9 * alpha) * u1^(-1 - alpha) * u10^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 -
alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) * (u7^(-1 -
alpha) * log(u7))) * u8^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
(1 + 8 * alpha) * (1 + 9 * alpha) * u1^(-1 - alpha) *
u10^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) *
u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) *
u7^(-1 - alpha) * (u8^(-1 - alpha) * log(u8))) * u9^(-1 -
alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * (1 + 9 * alpha) *
u1^(-1 - alpha) * u10^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) *
u6^(-1 - alpha) * u7^(-1 - alpha) * u8^(-1 - alpha) *
(u9^(-1 - alpha) * log(u9))) * (-9 + u1^(-alpha) + u10^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-10 -
1/alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * (1 + 9 * alpha) *
u1^(-1 - alpha) * u10^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) *
u6^(-1 - alpha) * u7^(-1 - alpha) * u8^(-1 - alpha) *
u9^(-1 - alpha) * ((-9 + u1^(-alpha) + u10^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-10 -
1/alpha) * (log((-9 + u1^(-alpha) + u10^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha) + u8^(-alpha) + u9^(-alpha))) * (1/alpha^2)) -
(-9 + u1^(-alpha) + u10^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) +
u8^(-alpha) + u9^(-alpha))^((-10 - 1/alpha) - 1) *
((-10 - 1/alpha) * (u9^(-alpha) * log(u9) + (u8^(-alpha) *
log(u8) + (u7^(-alpha) * log(u7) + (u6^(-alpha) *
log(u6) + (u5^(-alpha) * log(u5) + (u4^(-alpha) *
log(u4) + (u3^(-alpha) * log(u3) + (u2^(-alpha) *
log(u2) + (u10^(-alpha) * log(u10) + u1^(-alpha) *
log(u1)))))))))))))
`claytonCopula.pdfDerWrtPar.algr` <-
expression({
.expr2 <- u1^-alpha
.expr3 <- 1/alpha
.expr4 <- .expr2^.expr3
.value <- u1 * (.expr4 * (log(.expr2) * (1/alpha^2)) + .expr2^(.expr3 -
1) * (.expr3 * (.expr2 * log(u1))))/(u1 * .expr4)^2
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr1 <- -1
.expr2 <- .expr1 - alpha
.expr3 <- u1^.expr2
.expr4 <- 1 + alpha
.expr5 <- log(u1)
.expr9 <- u2^.expr2
.expr11 <- .expr4 * .expr3
.expr12 <- log(u2)
.expr16 <- -alpha
.expr17 <- u1^.expr16
.expr19 <- u2^.expr16
.expr20 <- .expr1 + .expr17 + .expr19
.expr23 <- -2 - 1/alpha
.expr24 <- .expr20^.expr23
.value <- ((.expr3 - .expr4 * (.expr3 * .expr5)) * .expr9 -
.expr11 * (.expr9 * .expr12)) * .expr24 + .expr11 * .expr9 *
(.expr24 * (log(.expr20) * (1/alpha^2)) - .expr20^(.expr23 -
1) * (.expr23 * (.expr19 * .expr12 + .expr17 * .expr5)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- -1 - alpha
.expr8 <- u1^.expr7
.expr10 <- .expr3 * .expr2
.expr11 <- log(u1)
.expr15 <- u2^.expr7
.expr17 <- .expr10 * .expr8
.expr18 <- log(u2)
.expr22 <- u3^.expr7
.expr24 <- .expr17 * .expr15
.expr25 <- log(u3)
.expr30 <- -alpha
.expr31 <- u1^.expr30
.expr33 <- u2^.expr30
.expr35 <- u3^.expr30
.expr36 <- -2 + .expr31 + .expr33 + .expr35
.expr39 <- -3 - 1/alpha
.expr40 <- .expr36^.expr39
.value <- ((((.expr2 + .expr3 * 2) * .expr8 - .expr10 * (.expr8 *
.expr11)) * .expr15 - .expr17 * (.expr15 * .expr18)) *
.expr22 - .expr24 * (.expr22 * .expr25)) * .expr40 +
.expr24 * .expr22 * (.expr40 * (log(.expr36) * (1/alpha^2)) -
.expr36^(.expr39 - 1) * (.expr39 * (.expr35 * .expr25 +
(.expr33 * .expr18 + .expr31 * .expr11))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- 1 + 3 * alpha
.expr9 <- .expr3 * .expr2
.expr13 <- -1 - alpha
.expr14 <- u1^.expr13
.expr16 <- .expr9 * .expr7
.expr17 <- log(u1)
.expr21 <- u2^.expr13
.expr23 <- .expr16 * .expr14
.expr24 <- log(u2)
.expr28 <- u3^.expr13
.expr30 <- .expr23 * .expr21
.expr31 <- log(u3)
.expr35 <- u4^.expr13
.expr37 <- .expr30 * .expr28
.expr38 <- log(u4)
.expr43 <- -alpha
.expr44 <- u1^.expr43
.expr46 <- u2^.expr43
.expr48 <- u3^.expr43
.expr50 <- u4^.expr43
.expr51 <- -3 + .expr44 + .expr46 + .expr48 + .expr50
.expr54 <- -4 - 1/alpha
.expr55 <- .expr51^.expr54
.value <- ((((((.expr2 + .expr3 * 2) * .expr7 + .expr9 *
3) * .expr14 - .expr16 * (.expr14 * .expr17)) * .expr21 -
.expr23 * (.expr21 * .expr24)) * .expr28 - .expr30 *
(.expr28 * .expr31)) * .expr35 - .expr37 * (.expr35 *
.expr38)) * .expr55 + .expr37 * .expr35 * (.expr55 *
(log(.expr51) * (1/alpha^2)) - .expr51^(.expr54 - 1) *
(.expr54 * (.expr50 * .expr38 + (.expr48 * .expr31 +
(.expr46 * .expr24 + .expr44 * .expr17)))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- 1 + 3 * alpha
.expr9 <- .expr3 * .expr2
.expr13 <- 1 + 4 * alpha
.expr15 <- .expr9 * .expr7
.expr19 <- -1 - alpha
.expr20 <- u1^.expr19
.expr22 <- .expr15 * .expr13
.expr23 <- log(u1)
.expr27 <- u2^.expr19
.expr29 <- .expr22 * .expr20
.expr30 <- log(u2)
.expr34 <- u3^.expr19
.expr36 <- .expr29 * .expr27
.expr37 <- log(u3)
.expr41 <- u4^.expr19
.expr43 <- .expr36 * .expr34
.expr44 <- log(u4)
.expr48 <- u5^.expr19
.expr50 <- .expr43 * .expr41
.expr51 <- log(u5)
.expr56 <- -alpha
.expr57 <- u1^.expr56
.expr59 <- u2^.expr56
.expr61 <- u3^.expr56
.expr63 <- u4^.expr56
.expr65 <- u5^.expr56
.expr66 <- -4 + .expr57 + .expr59 + .expr61 + .expr63 + .expr65
.expr69 <- -5 - 1/alpha
.expr70 <- .expr66^.expr69
.value <- ((((((((.expr2 + .expr3 * 2) * .expr7 + .expr9 *
3) * .expr13 + .expr15 * 4) * .expr20 - .expr22 * (.expr20 *
.expr23)) * .expr27 - .expr29 * (.expr27 * .expr30)) *
.expr34 - .expr36 * (.expr34 * .expr37)) * .expr41 -
.expr43 * (.expr41 * .expr44)) * .expr48 - .expr50 *
(.expr48 * .expr51)) * .expr70 + .expr50 * .expr48 *
(.expr70 * (log(.expr66) * (1/alpha^2)) - .expr66^(.expr69 -
1) * (.expr69 * (.expr65 * .expr51 + (.expr63 * .expr44 +
(.expr61 * .expr37 + (.expr59 * .expr30 + .expr57 *
.expr23))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- 1 + 3 * alpha
.expr9 <- .expr3 * .expr2
.expr13 <- 1 + 4 * alpha
.expr15 <- .expr9 * .expr7
.expr19 <- 1 + 5 * alpha
.expr21 <- .expr15 * .expr13
.expr25 <- -1 - alpha
.expr26 <- u1^.expr25
.expr28 <- .expr21 * .expr19
.expr29 <- log(u1)
.expr33 <- u2^.expr25
.expr35 <- .expr28 * .expr26
.expr36 <- log(u2)
.expr40 <- u3^.expr25
.expr42 <- .expr35 * .expr33
.expr43 <- log(u3)
.expr47 <- u4^.expr25
.expr49 <- .expr42 * .expr40
.expr50 <- log(u4)
.expr54 <- u5^.expr25
.expr56 <- .expr49 * .expr47
.expr57 <- log(u5)
.expr61 <- u6^.expr25
.expr63 <- .expr56 * .expr54
.expr64 <- log(u6)
.expr69 <- -alpha
.expr70 <- u1^.expr69
.expr72 <- u2^.expr69
.expr74 <- u3^.expr69
.expr76 <- u4^.expr69
.expr78 <- u5^.expr69
.expr80 <- u6^.expr69
.expr81 <- -5 + .expr70 + .expr72 + .expr74 + .expr76 + .expr78 +
.expr80
.expr84 <- -6 - 1/alpha
.expr85 <- .expr81^.expr84
.value <- ((((((((((.expr2 + .expr3 * 2) * .expr7 + .expr9 *
3) * .expr13 + .expr15 * 4) * .expr19 + .expr21 * 5) *
.expr26 - .expr28 * (.expr26 * .expr29)) * .expr33 -
.expr35 * (.expr33 * .expr36)) * .expr40 - .expr42 *
(.expr40 * .expr43)) * .expr47 - .expr49 * (.expr47 *
.expr50)) * .expr54 - .expr56 * (.expr54 * .expr57)) *
.expr61 - .expr63 * (.expr61 * .expr64)) * .expr85 +
.expr63 * .expr61 * (.expr85 * (log(.expr81) * (1/alpha^2)) -
.expr81^(.expr84 - 1) * (.expr84 * (.expr80 * .expr64 +
(.expr78 * .expr57 + (.expr76 * .expr50 + (.expr74 *
.expr43 + (.expr72 * .expr36 + .expr70 * .expr29)))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- 1 + 3 * alpha
.expr9 <- .expr3 * .expr2
.expr13 <- 1 + 4 * alpha
.expr15 <- .expr9 * .expr7
.expr19 <- 1 + 5 * alpha
.expr21 <- .expr15 * .expr13
.expr25 <- 1 + 6 * alpha
.expr27 <- .expr21 * .expr19
.expr31 <- -1 - alpha
.expr32 <- u1^.expr31
.expr34 <- .expr27 * .expr25
.expr35 <- log(u1)
.expr39 <- u2^.expr31
.expr41 <- .expr34 * .expr32
.expr42 <- log(u2)
.expr46 <- u3^.expr31
.expr48 <- .expr41 * .expr39
.expr49 <- log(u3)
.expr53 <- u4^.expr31
.expr55 <- .expr48 * .expr46
.expr56 <- log(u4)
.expr60 <- u5^.expr31
.expr62 <- .expr55 * .expr53
.expr63 <- log(u5)
.expr67 <- u6^.expr31
.expr69 <- .expr62 * .expr60
.expr70 <- log(u6)
.expr74 <- u7^.expr31
.expr76 <- .expr69 * .expr67
.expr77 <- log(u7)
.expr82 <- -alpha
.expr83 <- u1^.expr82
.expr85 <- u2^.expr82
.expr87 <- u3^.expr82
.expr89 <- u4^.expr82
.expr91 <- u5^.expr82
.expr93 <- u6^.expr82
.expr95 <- u7^.expr82
.expr96 <- -6 + .expr83 + .expr85 + .expr87 + .expr89 + .expr91 +
.expr93 + .expr95
.expr99 <- -7 - 1/alpha
.expr100 <- .expr96^.expr99
.value <- ((((((((((((.expr2 + .expr3 * 2) * .expr7 + .expr9 *
3) * .expr13 + .expr15 * 4) * .expr19 + .expr21 * 5) *
.expr25 + .expr27 * 6) * .expr32 - .expr34 * (.expr32 *
.expr35)) * .expr39 - .expr41 * (.expr39 * .expr42)) *
.expr46 - .expr48 * (.expr46 * .expr49)) * .expr53 -
.expr55 * (.expr53 * .expr56)) * .expr60 - .expr62 *
(.expr60 * .expr63)) * .expr67 - .expr69 * (.expr67 *
.expr70)) * .expr74 - .expr76 * (.expr74 * .expr77)) *
.expr100 + .expr76 * .expr74 * (.expr100 * (log(.expr96) *
(1/alpha^2)) - .expr96^(.expr99 - 1) * (.expr99 * (.expr95 *
.expr77 + (.expr93 * .expr70 + (.expr91 * .expr63 + (.expr89 *
.expr56 + (.expr87 * .expr49 + (.expr85 * .expr42 + .expr83 *
.expr35))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- 1 + 3 * alpha
.expr9 <- .expr3 * .expr2
.expr13 <- 1 + 4 * alpha
.expr15 <- .expr9 * .expr7
.expr19 <- 1 + 5 * alpha
.expr21 <- .expr15 * .expr13
.expr25 <- 1 + 6 * alpha
.expr27 <- .expr21 * .expr19
.expr31 <- 1 + 7 * alpha
.expr33 <- .expr27 * .expr25
.expr37 <- -1 - alpha
.expr38 <- u1^.expr37
.expr40 <- .expr33 * .expr31
.expr41 <- log(u1)
.expr45 <- u2^.expr37
.expr47 <- .expr40 * .expr38
.expr48 <- log(u2)
.expr52 <- u3^.expr37
.expr54 <- .expr47 * .expr45
.expr55 <- log(u3)
.expr59 <- u4^.expr37
.expr61 <- .expr54 * .expr52
.expr62 <- log(u4)
.expr66 <- u5^.expr37
.expr68 <- .expr61 * .expr59
.expr69 <- log(u5)
.expr73 <- u6^.expr37
.expr75 <- .expr68 * .expr66
.expr76 <- log(u6)
.expr80 <- u7^.expr37
.expr82 <- .expr75 * .expr73
.expr83 <- log(u7)
.expr87 <- u8^.expr37
.expr89 <- .expr82 * .expr80
.expr90 <- log(u8)
.expr95 <- -alpha
.expr96 <- u1^.expr95
.expr98 <- u2^.expr95
.expr100 <- u3^.expr95
.expr102 <- u4^.expr95
.expr104 <- u5^.expr95
.expr106 <- u6^.expr95
.expr108 <- u7^.expr95
.expr110 <- u8^.expr95
.expr111 <- -7 + .expr96 + .expr98 + .expr100 + .expr102 +
.expr104 + .expr106 + .expr108 + .expr110
.expr114 <- -8 - 1/alpha
.expr115 <- .expr111^.expr114
.value <- ((((((((((((((.expr2 + .expr3 * 2) * .expr7 + .expr9 *
3) * .expr13 + .expr15 * 4) * .expr19 + .expr21 * 5) *
.expr25 + .expr27 * 6) * .expr31 + .expr33 * 7) * .expr38 -
.expr40 * (.expr38 * .expr41)) * .expr45 - .expr47 *
(.expr45 * .expr48)) * .expr52 - .expr54 * (.expr52 *
.expr55)) * .expr59 - .expr61 * (.expr59 * .expr62)) *
.expr66 - .expr68 * (.expr66 * .expr69)) * .expr73 -
.expr75 * (.expr73 * .expr76)) * .expr80 - .expr82 *
(.expr80 * .expr83)) * .expr87 - .expr89 * (.expr87 *
.expr90)) * .expr115 + .expr89 * .expr87 * (.expr115 *
(log(.expr111) * (1/alpha^2)) - .expr111^(.expr114 -
1) * (.expr114 * (.expr110 * .expr90 + (.expr108 * .expr83 +
(.expr106 * .expr76 + (.expr104 * .expr69 + (.expr102 *
.expr62 + (.expr100 * .expr55 + (.expr98 * .expr48 +
.expr96 * .expr41)))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- 1 + 3 * alpha
.expr9 <- .expr3 * .expr2
.expr13 <- 1 + 4 * alpha
.expr15 <- .expr9 * .expr7
.expr19 <- 1 + 5 * alpha
.expr21 <- .expr15 * .expr13
.expr25 <- 1 + 6 * alpha
.expr27 <- .expr21 * .expr19
.expr31 <- 1 + 7 * alpha
.expr33 <- .expr27 * .expr25
.expr37 <- 1 + 8 * alpha
.expr39 <- .expr33 * .expr31
.expr43 <- -1 - alpha
.expr44 <- u1^.expr43
.expr46 <- .expr39 * .expr37
.expr47 <- log(u1)
.expr51 <- u2^.expr43
.expr53 <- .expr46 * .expr44
.expr54 <- log(u2)
.expr58 <- u3^.expr43
.expr60 <- .expr53 * .expr51
.expr61 <- log(u3)
.expr65 <- u4^.expr43
.expr67 <- .expr60 * .expr58
.expr68 <- log(u4)
.expr72 <- u5^.expr43
.expr74 <- .expr67 * .expr65
.expr75 <- log(u5)
.expr79 <- u6^.expr43
.expr81 <- .expr74 * .expr72
.expr82 <- log(u6)
.expr86 <- u7^.expr43
.expr88 <- .expr81 * .expr79
.expr89 <- log(u7)
.expr93 <- u8^.expr43
.expr95 <- .expr88 * .expr86
.expr96 <- log(u8)
.expr100 <- u9^.expr43
.expr102 <- .expr95 * .expr93
.expr103 <- log(u9)
.expr108 <- -alpha
.expr109 <- u1^.expr108
.expr111 <- u2^.expr108
.expr113 <- u3^.expr108
.expr115 <- u4^.expr108
.expr117 <- u5^.expr108
.expr119 <- u6^.expr108
.expr121 <- u7^.expr108
.expr123 <- u8^.expr108
.expr125 <- u9^.expr108
.expr126 <- -8 + .expr109 + .expr111 + .expr113 + .expr115 +
.expr117 + .expr119 + .expr121 + .expr123 + .expr125
.expr129 <- -9 - 1/alpha
.expr130 <- .expr126^.expr129
.value <- ((((((((((((((((.expr2 + .expr3 * 2) * .expr7 +
.expr9 * 3) * .expr13 + .expr15 * 4) * .expr19 + .expr21 *
5) * .expr25 + .expr27 * 6) * .expr31 + .expr33 * 7) *
.expr37 + .expr39 * 8) * .expr44 - .expr46 * (.expr44 *
.expr47)) * .expr51 - .expr53 * (.expr51 * .expr54)) *
.expr58 - .expr60 * (.expr58 * .expr61)) * .expr65 -
.expr67 * (.expr65 * .expr68)) * .expr72 - .expr74 *
(.expr72 * .expr75)) * .expr79 - .expr81 * (.expr79 *
.expr82)) * .expr86 - .expr88 * (.expr86 * .expr89)) *
.expr93 - .expr95 * (.expr93 * .expr96)) * .expr100 -
.expr102 * (.expr100 * .expr103)) * .expr130 + .expr102 *
.expr100 * (.expr130 * (log(.expr126) * (1/alpha^2)) -
.expr126^(.expr129 - 1) * (.expr129 * (.expr125 * .expr103 +
(.expr123 * .expr96 + (.expr121 * .expr89 + (.expr119 *
.expr82 + (.expr117 * .expr75 + (.expr115 * .expr68 +
(.expr113 * .expr61 + (.expr111 * .expr54 + .expr109 *
.expr47))))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- 1 + 3 * alpha
.expr9 <- .expr3 * .expr2
.expr13 <- 1 + 4 * alpha
.expr15 <- .expr9 * .expr7
.expr19 <- 1 + 5 * alpha
.expr21 <- .expr15 * .expr13
.expr25 <- 1 + 6 * alpha
.expr27 <- .expr21 * .expr19
.expr31 <- 1 + 7 * alpha
.expr33 <- .expr27 * .expr25
.expr37 <- 1 + 8 * alpha
.expr39 <- .expr33 * .expr31
.expr43 <- 1 + 9 * alpha
.expr45 <- .expr39 * .expr37
.expr49 <- -1 - alpha
.expr50 <- u1^.expr49
.expr52 <- .expr45 * .expr43
.expr53 <- log(u1)
.expr57 <- u10^.expr49
.expr59 <- .expr52 * .expr50
.expr60 <- log(u10)
.expr64 <- u2^.expr49
.expr66 <- .expr59 * .expr57
.expr67 <- log(u2)
.expr71 <- u3^.expr49
.expr73 <- .expr66 * .expr64
.expr74 <- log(u3)
.expr78 <- u4^.expr49
.expr80 <- .expr73 * .expr71
.expr81 <- log(u4)
.expr85 <- u5^.expr49
.expr87 <- .expr80 * .expr78
.expr88 <- log(u5)
.expr92 <- u6^.expr49
.expr94 <- .expr87 * .expr85
.expr95 <- log(u6)
.expr99 <- u7^.expr49
.expr101 <- .expr94 * .expr92
.expr102 <- log(u7)
.expr106 <- u8^.expr49
.expr108 <- .expr101 * .expr99
.expr109 <- log(u8)
.expr113 <- u9^.expr49
.expr115 <- .expr108 * .expr106
.expr116 <- log(u9)
.expr121 <- -alpha
.expr122 <- u1^.expr121
.expr124 <- u10^.expr121
.expr126 <- u2^.expr121
.expr128 <- u3^.expr121
.expr130 <- u4^.expr121
.expr132 <- u5^.expr121
.expr134 <- u6^.expr121
.expr136 <- u7^.expr121
.expr138 <- u8^.expr121
.expr140 <- u9^.expr121
.expr141 <- -9 + .expr122 + .expr124 + .expr126 + .expr128 +
.expr130 + .expr132 + .expr134 + .expr136 + .expr138 +
.expr140
.expr144 <- -10 - 1/alpha
.expr145 <- .expr141^.expr144
.value <- ((((((((((((((((((.expr2 + .expr3 * 2) * .expr7 +
.expr9 * 3) * .expr13 + .expr15 * 4) * .expr19 + .expr21 *
5) * .expr25 + .expr27 * 6) * .expr31 + .expr33 * 7) *
.expr37 + .expr39 * 8) * .expr43 + .expr45 * 9) * .expr50 -
.expr52 * (.expr50 * .expr53)) * .expr57 - .expr59 *
(.expr57 * .expr60)) * .expr64 - .expr66 * (.expr64 *
.expr67)) * .expr71 - .expr73 * (.expr71 * .expr74)) *
.expr78 - .expr80 * (.expr78 * .expr81)) * .expr85 -
.expr87 * (.expr85 * .expr88)) * .expr92 - .expr94 *
(.expr92 * .expr95)) * .expr99 - .expr101 * (.expr99 *
.expr102)) * .expr106 - .expr108 * (.expr106 * .expr109)) *
.expr113 - .expr115 * (.expr113 * .expr116)) * .expr145 +
.expr115 * .expr113 * (.expr145 * (log(.expr141) * (1/alpha^2)) -
.expr141^(.expr144 - 1) * (.expr144 * (.expr140 *
.expr116 + (.expr138 * .expr109 + (.expr136 *
.expr102 + (.expr134 * .expr95 + (.expr132 *
.expr88 + (.expr130 * .expr81 + (.expr128 * .expr74 +
(.expr126 * .expr67 + (.expr124 * .expr60 + .expr122 *
.expr53)))))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
})
`claytonCopula.genfunDer.expr` <-
expression(-(alpha * u^(-1 - alpha)), alpha * (1 + alpha) * u^(-2 -
alpha))
`claytonCopula.genfunDer.algr` <-
expression({
.value <- -(alpha * u^(-1 - alpha))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.value <- alpha * (1 + alpha) * u^(-2 - alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
})
Try the copula package in your browser
Any scripts or data that you put into this service are public.
copula documentation built on Sept. 11, 2024, 7:48 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
## This file is not to be edited. It is generated by ../inst/docs/mathnb/exprPrep.R
`claytonCopula.cdf.expr` <-
expression((u1^(-alpha))^(-1/alpha), (-1 + u1^(-alpha) + u2^(-alpha))^(-1/alpha),
(-2 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha))^(-1/alpha),
(-3 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha))^(-1/alpha),
(-4 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha))^(-1/alpha), (-5 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha))^(-1/alpha),
(-6 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha))^(-1/alpha),
(-7 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^(-1/alpha),
(-8 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha) +
u9^(-alpha))^(-1/alpha), (-9 + u1^(-alpha) + u10^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-1/alpha))
`claytonCopula.cdf.algr` <-
expression({
.value <- (u1^-alpha)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr1 <- -1
.expr2 <- -alpha
.value <- (.expr1 + u1^.expr2 + u2^.expr2)^(.expr1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-2 + u1^.expr2 + u2^.expr2 + u3^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-3 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-4 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-5 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-6 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2 + u7^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-7 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2 + u7^.expr2 + u8^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-8 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2 + u7^.expr2 + u8^.expr2 + u9^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.value <- (-9 + u1^.expr2 + u10^.expr2 + u2^.expr2 + u3^.expr2 +
u4^.expr2 + u5^.expr2 + u6^.expr2 + u7^.expr2 + u8^.expr2 +
u9^.expr2)^(-1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
})
`claytonCopula.cdfDerWrtArg.expr` <-
expression((u1^(-alpha))^((-1/alpha) - 1) * ((-1/alpha) * (u1^((-alpha) -
1) * (-alpha))), (-1 + u1^(-alpha) + u2^(-alpha))^((-1/alpha) -
1) * ((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))), (-2 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha))^((-1/alpha) - 1) *
((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))), (-3 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha))^((-1/alpha) - 1) *
((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))), (-4 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha))^((-1/alpha) -
1) * ((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))), (-5 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha))^((-1/alpha) - 1) * ((-1/alpha) * (u1^((-alpha) -
1) * (-alpha))), (-6 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha))^((-1/alpha) -
1) * ((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))), (-7 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^((-1/alpha) - 1) *
((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))), (-8 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^((-1/alpha) - 1) *
((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))), (-9 + u1^(-alpha) +
u10^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^((-1/alpha) -
1) * ((-1/alpha) * (u1^((-alpha) - 1) * (-alpha))))
`claytonCopula.cdfDerWrtArg.algr` <-
expression({
.expr1 <- -alpha
.expr4 <- -1/alpha
.value <- (u1^.expr1)^(.expr4 - 1) * (.expr4 * (u1^(.expr1 -
1) * .expr1))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr1 <- -1
.expr2 <- -alpha
.expr7 <- .expr1/alpha
.value <- (.expr1 + u1^.expr2 + u2^.expr2)^(.expr7 - 1) *
(.expr7 * (u1^(.expr2 - 1) * .expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr10 <- -1/alpha
.value <- (-2 + u1^.expr2 + u2^.expr2 + u3^.expr2)^(.expr10 -
1) * (.expr10 * (u1^(.expr2 - 1) * .expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr12 <- -1/alpha
.value <- (-3 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2)^(.expr12 -
1) * (.expr12 * (u1^(.expr2 - 1) * .expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr14 <- -1/alpha
.value <- (-4 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2)^(.expr14 - 1) * (.expr14 * (u1^(.expr2 - 1) *
.expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr16 <- -1/alpha
.value <- (-5 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2)^(.expr16 - 1) * (.expr16 * (u1^(.expr2 -
1) * .expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr18 <- -1/alpha
.value <- (-6 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2 + u7^.expr2)^(.expr18 - 1) * (.expr18 *
(u1^(.expr2 - 1) * .expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr20 <- -1/alpha
.value <- (-7 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2 + u7^.expr2 + u8^.expr2)^(.expr20 -
1) * (.expr20 * (u1^(.expr2 - 1) * .expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr22 <- -1/alpha
.value <- (-8 + u1^.expr2 + u2^.expr2 + u3^.expr2 + u4^.expr2 +
u5^.expr2 + u6^.expr2 + u7^.expr2 + u8^.expr2 + u9^.expr2)^(.expr22 -
1) * (.expr22 * (u1^(.expr2 - 1) * .expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr24 <- -1/alpha
.value <- (-9 + u1^.expr2 + u10^.expr2 + u2^.expr2 + u3^.expr2 +
u4^.expr2 + u5^.expr2 + u6^.expr2 + u7^.expr2 + u8^.expr2 +
u9^.expr2)^(.expr24 - 1) * (.expr24 * (u1^(.expr2 - 1) *
.expr2))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
})
`claytonCopula.cdfDerWrtPar.expr` <-
expression((u1^(-alpha))^(-1/alpha) * (log((u1^(-alpha))) * (1/alpha^2)) -
(u1^(-alpha))^((-1/alpha) - 1) * ((-1/alpha) * (u1^(-alpha) *
log(u1))), (-1 + u1^(-alpha) + u2^(-alpha))^(-1/alpha) *
(log((-1 + u1^(-alpha) + u2^(-alpha))) * (1/alpha^2)) - (-1 +
u1^(-alpha) + u2^(-alpha))^((-1/alpha) - 1) * ((-1/alpha) *
(u2^(-alpha) * log(u2) + u1^(-alpha) * log(u1))), (-2 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha))^(-1/alpha) * (log((-2 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha))) * (1/alpha^2)) - (-2 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha))^((-1/alpha) - 1) * ((-1/alpha) *
(u3^(-alpha) * log(u3) + (u2^(-alpha) * log(u2) + u1^(-alpha) *
log(u1)))), (-3 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha))^(-1/alpha) * (log((-3 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha))) * (1/alpha^2)) - (-3 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha))^((-1/alpha) - 1) *
((-1/alpha) * (u4^(-alpha) * log(u4) + (u3^(-alpha) * log(u3) +
(u2^(-alpha) * log(u2) + u1^(-alpha) * log(u1))))), (-4 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha))^(-1/alpha) *
(log((-4 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha))) * (1/alpha^2)) - (-4 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha))^((-1/alpha) - 1) *
((-1/alpha) * (u5^(-alpha) * log(u5) + (u4^(-alpha) * log(u4) +
(u3^(-alpha) * log(u3) + (u2^(-alpha) * log(u2) + u1^(-alpha) *
log(u1)))))), (-5 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha))^(-1/alpha) * (log((-5 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha))) * (1/alpha^2)) - (-5 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha))^((-1/alpha) -
1) * ((-1/alpha) * (u6^(-alpha) * log(u6) + (u5^(-alpha) *
log(u5) + (u4^(-alpha) * log(u4) + (u3^(-alpha) * log(u3) +
(u2^(-alpha) * log(u2) + u1^(-alpha) * log(u1))))))), (-6 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha))^(-1/alpha) * (log((-6 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha))) * (1/alpha^2)) - (-6 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha))^((-1/alpha) -
1) * ((-1/alpha) * (u7^(-alpha) * log(u7) + (u6^(-alpha) *
log(u6) + (u5^(-alpha) * log(u5) + (u4^(-alpha) * log(u4) +
(u3^(-alpha) * log(u3) + (u2^(-alpha) * log(u2) + u1^(-alpha) *
log(u1)))))))), (-7 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^(-1/alpha) *
(log((-7 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha))) *
(1/alpha^2)) - (-7 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^((-1/alpha) -
1) * ((-1/alpha) * (u8^(-alpha) * log(u8) + (u7^(-alpha) *
log(u7) + (u6^(-alpha) * log(u6) + (u5^(-alpha) * log(u5) +
(u4^(-alpha) * log(u4) + (u3^(-alpha) * log(u3) + (u2^(-alpha) *
log(u2) + u1^(-alpha) * log(u1))))))))), (-8 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-1/alpha) * (log((-8 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))) *
(1/alpha^2)) - (-8 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha) +
u9^(-alpha))^((-1/alpha) - 1) * ((-1/alpha) * (u9^(-alpha) *
log(u9) + (u8^(-alpha) * log(u8) + (u7^(-alpha) * log(u7) +
(u6^(-alpha) * log(u6) + (u5^(-alpha) * log(u5) + (u4^(-alpha) *
log(u4) + (u3^(-alpha) * log(u3) + (u2^(-alpha) * log(u2) +
u1^(-alpha) * log(u1)))))))))), (-9 + u1^(-alpha) + u10^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-1/alpha) * (log((-9 +
u1^(-alpha) + u10^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha) +
u9^(-alpha))) * (1/alpha^2)) - (-9 + u1^(-alpha) + u10^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^((-1/alpha) - 1) *
((-1/alpha) * (u9^(-alpha) * log(u9) + (u8^(-alpha) * log(u8) +
(u7^(-alpha) * log(u7) + (u6^(-alpha) * log(u6) + (u5^(-alpha) *
log(u5) + (u4^(-alpha) * log(u4) + (u3^(-alpha) *
log(u3) + (u2^(-alpha) * log(u2) + (u10^(-alpha) *
log(u10) + u1^(-alpha) * log(u1))))))))))))
`claytonCopula.cdfDerWrtPar.algr` <-
expression({
.expr2 <- u1^-alpha
.expr4 <- -1/alpha
.value <- .expr2^.expr4 * (log(.expr2) * (1/alpha^2)) - .expr2^(.expr4 -
1) * (.expr4 * (.expr2 * log(u1)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr1 <- -1
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr6 <- .expr1 + .expr3 + .expr5
.expr7 <- .expr1/alpha
.value <- .expr6^.expr7 * (log(.expr6) * (1/alpha^2)) - .expr6^(.expr7 -
1) * (.expr7 * (.expr5 * log(u2) + .expr3 * log(u1)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr7 <- u3^.expr2
.expr8 <- -2 + .expr3 + .expr5 + .expr7
.expr10 <- -1/alpha
.value <- .expr8^.expr10 * (log(.expr8) * (1/alpha^2)) -
.expr8^(.expr10 - 1) * (.expr10 * (.expr7 * log(u3) +
(.expr5 * log(u2) + .expr3 * log(u1))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr7 <- u3^.expr2
.expr9 <- u4^.expr2
.expr10 <- -3 + .expr3 + .expr5 + .expr7 + .expr9
.expr12 <- -1/alpha
.value <- .expr10^.expr12 * (log(.expr10) * (1/alpha^2)) -
.expr10^(.expr12 - 1) * (.expr12 * (.expr9 * log(u4) +
(.expr7 * log(u3) + (.expr5 * log(u2) + .expr3 *
log(u1)))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr7 <- u3^.expr2
.expr9 <- u4^.expr2
.expr11 <- u5^.expr2
.expr12 <- -4 + .expr3 + .expr5 + .expr7 + .expr9 + .expr11
.expr14 <- -1/alpha
.value <- .expr12^.expr14 * (log(.expr12) * (1/alpha^2)) -
.expr12^(.expr14 - 1) * (.expr14 * (.expr11 * log(u5) +
(.expr9 * log(u4) + (.expr7 * log(u3) + (.expr5 *
log(u2) + .expr3 * log(u1))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr7 <- u3^.expr2
.expr9 <- u4^.expr2
.expr11 <- u5^.expr2
.expr13 <- u6^.expr2
.expr14 <- -5 + .expr3 + .expr5 + .expr7 + .expr9 + .expr11 +
.expr13
.expr16 <- -1/alpha
.value <- .expr14^.expr16 * (log(.expr14) * (1/alpha^2)) -
.expr14^(.expr16 - 1) * (.expr16 * (.expr13 * log(u6) +
(.expr11 * log(u5) + (.expr9 * log(u4) + (.expr7 *
log(u3) + (.expr5 * log(u2) + .expr3 * log(u1)))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr7 <- u3^.expr2
.expr9 <- u4^.expr2
.expr11 <- u5^.expr2
.expr13 <- u6^.expr2
.expr15 <- u7^.expr2
.expr16 <- -6 + .expr3 + .expr5 + .expr7 + .expr9 + .expr11 +
.expr13 + .expr15
.expr18 <- -1/alpha
.value <- .expr16^.expr18 * (log(.expr16) * (1/alpha^2)) -
.expr16^(.expr18 - 1) * (.expr18 * (.expr15 * log(u7) +
(.expr13 * log(u6) + (.expr11 * log(u5) + (.expr9 *
log(u4) + (.expr7 * log(u3) + (.expr5 * log(u2) +
.expr3 * log(u1))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr7 <- u3^.expr2
.expr9 <- u4^.expr2
.expr11 <- u5^.expr2
.expr13 <- u6^.expr2
.expr15 <- u7^.expr2
.expr17 <- u8^.expr2
.expr18 <- -7 + .expr3 + .expr5 + .expr7 + .expr9 + .expr11 +
.expr13 + .expr15 + .expr17
.expr20 <- -1/alpha
.value <- .expr18^.expr20 * (log(.expr18) * (1/alpha^2)) -
.expr18^(.expr20 - 1) * (.expr20 * (.expr17 * log(u8) +
(.expr15 * log(u7) + (.expr13 * log(u6) + (.expr11 *
log(u5) + (.expr9 * log(u4) + (.expr7 * log(u3) +
(.expr5 * log(u2) + .expr3 * log(u1)))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u2^.expr2
.expr7 <- u3^.expr2
.expr9 <- u4^.expr2
.expr11 <- u5^.expr2
.expr13 <- u6^.expr2
.expr15 <- u7^.expr2
.expr17 <- u8^.expr2
.expr19 <- u9^.expr2
.expr20 <- -8 + .expr3 + .expr5 + .expr7 + .expr9 + .expr11 +
.expr13 + .expr15 + .expr17 + .expr19
.expr22 <- -1/alpha
.value <- .expr20^.expr22 * (log(.expr20) * (1/alpha^2)) -
.expr20^(.expr22 - 1) * (.expr22 * (.expr19 * log(u9) +
(.expr17 * log(u8) + (.expr15 * log(u7) + (.expr13 *
log(u6) + (.expr11 * log(u5) + (.expr9 * log(u4) +
(.expr7 * log(u3) + (.expr5 * log(u2) + .expr3 *
log(u1))))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -alpha
.expr3 <- u1^.expr2
.expr5 <- u10^.expr2
.expr7 <- u2^.expr2
.expr9 <- u3^.expr2
.expr11 <- u4^.expr2
.expr13 <- u5^.expr2
.expr15 <- u6^.expr2
.expr17 <- u7^.expr2
.expr19 <- u8^.expr2
.expr21 <- u9^.expr2
.expr22 <- -9 + .expr3 + .expr5 + .expr7 + .expr9 + .expr11 +
.expr13 + .expr15 + .expr17 + .expr19 + .expr21
.expr24 <- -1/alpha
.value <- .expr22^.expr24 * (log(.expr22) * (1/alpha^2)) -
.expr22^(.expr24 - 1) * (.expr24 * (.expr21 * log(u9) +
(.expr19 * log(u8) + (.expr17 * log(u7) + (.expr15 *
log(u6) + (.expr13 * log(u5) + (.expr11 * log(u4) +
(.expr9 * log(u3) + (.expr7 * log(u2) + (.expr5 *
log(u10) + .expr3 * log(u1)))))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
})
`claytonCopula.pdf.expr` <-
expression(1/(u1 * (u1^(-alpha))^(1/alpha)), (1 + alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * (-1 + u1^(-alpha) + u2^(-alpha))^(-2 -
1/alpha), (1 + alpha) * (1 + 2 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * (-2 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha))^(-3 - 1/alpha), (1 + alpha) * (1 + 2 * alpha) *
(1 + 3 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 -
alpha) * u4^(-1 - alpha) * (-3 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha))^(-4 - 1/alpha), (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * (-4 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha))^(-5 - 1/alpha), (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 +
5 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 -
alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) *
(-5 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha))^(-6 - 1/alpha), (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 +
5 * alpha) * (1 + 6 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) *
u6^(-1 - alpha) * u7^(-1 - alpha) * (-6 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha))^(-7 -
1/alpha), (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * (1 +
7 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 -
alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) *
u7^(-1 - alpha) * u8^(-1 - alpha) * (-7 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) +
u8^(-alpha))^(-8 - 1/alpha), (1 + alpha) * (1 + 2 * alpha) *
(1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * (1 + 7 * alpha) * (1 + 8 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) * u8^(-1 -
alpha) * u9^(-1 - alpha) * (-8 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) +
u8^(-alpha) + u9^(-alpha))^(-9 - 1/alpha), (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 +
5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) * (1 + 8 *
alpha) * (1 + 9 * alpha) * u1^(-1 - alpha) * u10^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 -
alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) * u8^(-1 - alpha) *
u9^(-1 - alpha) * (-9 + u1^(-alpha) + u10^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) +
u8^(-alpha) + u9^(-alpha))^(-10 - 1/alpha))
`claytonCopula.pdf.algr` <-
expression({
.value <- 1/(u1 * (u1^-alpha)^(1/alpha))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- -1
.expr3 <- .expr2 - alpha
.expr8 <- -alpha
.value <- (1 + alpha) * u1^.expr3 * u2^.expr3 * (.expr2 +
u1^.expr8 + u2^.expr8)^(-2 - 1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr6 <- -1 - alpha
.expr14 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * u1^.expr6 * u2^.expr6 *
u3^.expr6 * (-2 + u1^.expr14 + u2^.expr14 + u3^.expr14)^(-3 -
1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr9 <- -1 - alpha
.expr19 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
u1^.expr9 * u2^.expr9 * u3^.expr9 * u4^.expr9 * (-3 +
u1^.expr19 + u2^.expr19 + u3^.expr19 + u4^.expr19)^(-4 -
1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr12 <- -1 - alpha
.expr24 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * u1^.expr12 * u2^.expr12 * u3^.expr12 *
u4^.expr12 * u5^.expr12 * (-4 + u1^.expr24 + u2^.expr24 +
u3^.expr24 + u4^.expr24 + u5^.expr24)^(-5 - 1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr15 <- -1 - alpha
.expr29 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * u1^.expr15 * u2^.expr15 *
u3^.expr15 * u4^.expr15 * u5^.expr15 * u6^.expr15 * (-5 +
u1^.expr29 + u2^.expr29 + u3^.expr29 + u4^.expr29 + u5^.expr29 +
u6^.expr29)^(-6 - 1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr18 <- -1 - alpha
.expr34 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
u1^.expr18 * u2^.expr18 * u3^.expr18 * u4^.expr18 * u5^.expr18 *
u6^.expr18 * u7^.expr18 * (-6 + u1^.expr34 + u2^.expr34 +
u3^.expr34 + u4^.expr34 + u5^.expr34 + u6^.expr34 + u7^.expr34)^(-7 -
1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr21 <- -1 - alpha
.expr39 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * u1^.expr21 * u2^.expr21 * u3^.expr21 *
u4^.expr21 * u5^.expr21 * u6^.expr21 * u7^.expr21 * u8^.expr21 *
(-7 + u1^.expr39 + u2^.expr39 + u3^.expr39 + u4^.expr39 +
u5^.expr39 + u6^.expr39 + u7^.expr39 + u8^.expr39)^(-8 -
1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr24 <- -1 - alpha
.expr44 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * u1^.expr24 * u2^.expr24 *
u3^.expr24 * u4^.expr24 * u5^.expr24 * u6^.expr24 * u7^.expr24 *
u8^.expr24 * u9^.expr24 * (-8 + u1^.expr44 + u2^.expr44 +
u3^.expr44 + u4^.expr44 + u5^.expr44 + u6^.expr44 + u7^.expr44 +
u8^.expr44 + u9^.expr44)^(-9 - 1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr27 <- -1 - alpha
.expr49 <- -alpha
.value <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * (1 + 9 * alpha) *
u1^.expr27 * u10^.expr27 * u2^.expr27 * u3^.expr27 *
u4^.expr27 * u5^.expr27 * u6^.expr27 * u7^.expr27 * u8^.expr27 *
u9^.expr27 * (-9 + u1^.expr49 + u10^.expr49 + u2^.expr49 +
u3^.expr49 + u4^.expr49 + u5^.expr49 + u6^.expr49 + u7^.expr49 +
u8^.expr49 + u9^.expr49)^(-10 - 1/alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
})
`claytonCopula.pdfDerWrtArg.expr` <-
expression(-(((u1^(-alpha))^(1/alpha) + u1 * ((u1^(-alpha))^((1/alpha) -
1) * ((1/alpha) * (u1^((-alpha) - 1) * (-alpha)))))/(u1 *
(u1^(-alpha))^(1/alpha))^2), (1 + alpha) * (u1^((-1 - alpha) -
1) * (-1 - alpha)) * u2^(-1 - alpha) * (-1 + u1^(-alpha) +
u2^(-alpha))^(-2 - 1/alpha) + (1 + alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * ((-1 + u1^(-alpha) + u2^(-alpha))^((-2 -
1/alpha) - 1) * ((-2 - 1/alpha) * (u1^((-alpha) - 1) * (-alpha)))),
(1 + alpha) * (1 + 2 * alpha) * (u1^((-1 - alpha) - 1) *
(-1 - alpha)) * u2^(-1 - alpha) * u3^(-1 - alpha) * (-2 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha))^(-3 - 1/alpha) +
(1 + alpha) * (1 + 2 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * ((-2 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha))^((-3 - 1/alpha) - 1) * ((-3 - 1/alpha) *
(u1^((-alpha) - 1) * (-alpha)))), (1 + alpha) * (1 +
2 * alpha) * (1 + 3 * alpha) * (u1^((-1 - alpha) - 1) *
(-1 - alpha)) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 -
alpha) * (-3 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha))^(-4 - 1/alpha) + (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * ((-3 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha))^((-4 - 1/alpha) -
1) * ((-4 - 1/alpha) * (u1^((-alpha) - 1) * (-alpha)))),
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (u1^((-1 - alpha) - 1) * (-1 - alpha)) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 -
alpha) * (-4 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha))^(-5 - 1/alpha) + (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) *
u4^(-1 - alpha) * u5^(-1 - alpha) * ((-4 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha))^((-5 -
1/alpha) - 1) * ((-5 - 1/alpha) * (u1^((-alpha) - 1) *
(-alpha)))), (1 + alpha) * (1 + 2 * alpha) * (1 + 3 *
alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (u1^((-1 -
alpha) - 1) * (-1 - alpha)) * u2^(-1 - alpha) * u3^(-1 -
alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 -
alpha) * (-5 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha))^(-6 - 1/alpha) +
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 +
4 * alpha) * (1 + 5 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * u6^(-1 - alpha) * ((-5 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha))^((-6 - 1/alpha) - 1) * ((-6 - 1/alpha) *
(u1^((-alpha) - 1) * (-alpha)))), (1 + alpha) * (1 +
2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 +
5 * alpha) * (1 + 6 * alpha) * (u1^((-1 - alpha) - 1) *
(-1 - alpha)) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 -
alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 -
alpha) * (-6 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha))^(-7 -
1/alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) *
u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) *
u7^(-1 - alpha) * ((-6 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha))^((-7 - 1/alpha) - 1) * ((-7 - 1/alpha) *
(u1^((-alpha) - 1) * (-alpha)))), (1 + alpha) * (1 +
2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 +
5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) * (u1^((-1 -
alpha) - 1) * (-1 - alpha)) * u2^(-1 - alpha) * u3^(-1 -
alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 -
alpha) * u7^(-1 - alpha) * u8^(-1 - alpha) * (-7 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^(-8 - 1/alpha) +
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 +
4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) *
u6^(-1 - alpha) * u7^(-1 - alpha) * u8^(-1 - alpha) *
((-7 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) +
u8^(-alpha))^((-8 - 1/alpha) - 1) * ((-8 - 1/alpha) *
(u1^((-alpha) - 1) * (-alpha)))), (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
(1 + 8 * alpha) * (u1^((-1 - alpha) - 1) * (-1 - alpha)) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) *
u8^(-1 - alpha) * u9^(-1 - alpha) * (-8 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-9 -
1/alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) *
u8^(-1 - alpha) * u9^(-1 - alpha) * ((-8 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^((-9 -
1/alpha) - 1) * ((-9 - 1/alpha) * (u1^((-alpha) - 1) *
(-alpha)))), (1 + alpha) * (1 + 2 * alpha) * (1 + 3 *
alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 *
alpha) * (1 + 7 * alpha) * (1 + 8 * alpha) * (1 + 9 *
alpha) * (u1^((-1 - alpha) - 1) * (-1 - alpha)) * u10^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 -
alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 -
alpha) * u8^(-1 - alpha) * u9^(-1 - alpha) * (-9 + u1^(-alpha) +
u10^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha) +
u9^(-alpha))^(-10 - 1/alpha) + (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) * (1 + 8 *
alpha) * (1 + 9 * alpha) * u1^(-1 - alpha) * u10^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 -
alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 -
alpha) * u8^(-1 - alpha) * u9^(-1 - alpha) * ((-9 + u1^(-alpha) +
u10^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha) +
u9^(-alpha))^((-10 - 1/alpha) - 1) * ((-10 - 1/alpha) *
(u1^((-alpha) - 1) * (-alpha)))))
`claytonCopula.pdfDerWrtArg.algr` <-
expression({
.expr1 <- -alpha
.expr2 <- u1^.expr1
.expr3 <- 1/alpha
.expr4 <- .expr2^.expr3
.value <- -((.expr4 + u1 * (.expr2^(.expr3 - 1) * (.expr3 *
(u1^(.expr1 - 1) * .expr1))))/(u1 * .expr4)^2)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr1 <- 1 + alpha
.expr2 <- -1
.expr3 <- .expr2 - alpha
.expr8 <- u2^.expr3
.expr10 <- -alpha
.expr14 <- .expr2 + u1^.expr10 + u2^.expr10
.expr17 <- -2 - 1/alpha
.value <- .expr1 * (u1^(.expr3 - 1) * .expr3) * .expr8 *
.expr14^.expr17 + .expr1 * u1^.expr3 * .expr8 * (.expr14^(.expr17 -
1) * (.expr17 * (u1^(.expr10 - 1) * .expr10)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr4 <- (1 + alpha) * (1 + 2 * alpha)
.expr6 <- -1 - alpha
.expr11 <- u2^.expr6
.expr13 <- u3^.expr6
.expr16 <- -alpha
.expr22 <- -2 + u1^.expr16 + u2^.expr16 + u3^.expr16
.expr25 <- -3 - 1/alpha
.value <- .expr4 * (u1^(.expr6 - 1) * .expr6) * .expr11 *
.expr13 * .expr22^.expr25 + .expr4 * u1^.expr6 * .expr11 *
.expr13 * (.expr22^(.expr25 - 1) * (.expr25 * (u1^(.expr16 -
1) * .expr16)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr7 <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha)
.expr9 <- -1 - alpha
.expr14 <- u2^.expr9
.expr16 <- u3^.expr9
.expr18 <- u4^.expr9
.expr21 <- -alpha
.expr29 <- -3 + u1^.expr21 + u2^.expr21 + u3^.expr21 + u4^.expr21
.expr32 <- -4 - 1/alpha
.value <- .expr7 * (u1^(.expr9 - 1) * .expr9) * .expr14 *
.expr16 * .expr18 * .expr29^.expr32 + .expr7 * u1^.expr9 *
.expr14 * .expr16 * .expr18 * (.expr29^(.expr32 - 1) *
(.expr32 * (u1^(.expr21 - 1) * .expr21)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr10 <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha)
.expr12 <- -1 - alpha
.expr17 <- u2^.expr12
.expr19 <- u3^.expr12
.expr21 <- u4^.expr12
.expr23 <- u5^.expr12
.expr26 <- -alpha
.expr36 <- -4 + u1^.expr26 + u2^.expr26 + u3^.expr26 + u4^.expr26 +
u5^.expr26
.expr39 <- -5 - 1/alpha
.value <- .expr10 * (u1^(.expr12 - 1) * .expr12) * .expr17 *
.expr19 * .expr21 * .expr23 * .expr36^.expr39 + .expr10 *
u1^.expr12 * .expr17 * .expr19 * .expr21 * .expr23 *
(.expr36^(.expr39 - 1) * (.expr39 * (u1^(.expr26 - 1) *
.expr26)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr13 <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha)
.expr15 <- -1 - alpha
.expr20 <- u2^.expr15
.expr22 <- u3^.expr15
.expr24 <- u4^.expr15
.expr26 <- u5^.expr15
.expr28 <- u6^.expr15
.expr31 <- -alpha
.expr43 <- -5 + u1^.expr31 + u2^.expr31 + u3^.expr31 + u4^.expr31 +
u5^.expr31 + u6^.expr31
.expr46 <- -6 - 1/alpha
.value <- .expr13 * (u1^(.expr15 - 1) * .expr15) * .expr20 *
.expr22 * .expr24 * .expr26 * .expr28 * .expr43^.expr46 +
.expr13 * u1^.expr15 * .expr20 * .expr22 * .expr24 *
.expr26 * .expr28 * (.expr43^(.expr46 - 1) * (.expr46 *
(u1^(.expr31 - 1) * .expr31)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr16 <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha)
.expr18 <- -1 - alpha
.expr23 <- u2^.expr18
.expr25 <- u3^.expr18
.expr27 <- u4^.expr18
.expr29 <- u5^.expr18
.expr31 <- u6^.expr18
.expr33 <- u7^.expr18
.expr36 <- -alpha
.expr50 <- -6 + u1^.expr36 + u2^.expr36 + u3^.expr36 + u4^.expr36 +
u5^.expr36 + u6^.expr36 + u7^.expr36
.expr53 <- -7 - 1/alpha
.value <- .expr16 * (u1^(.expr18 - 1) * .expr18) * .expr23 *
.expr25 * .expr27 * .expr29 * .expr31 * .expr33 * .expr50^.expr53 +
.expr16 * u1^.expr18 * .expr23 * .expr25 * .expr27 *
.expr29 * .expr31 * .expr33 * (.expr50^(.expr53 -
1) * (.expr53 * (u1^(.expr36 - 1) * .expr36)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr19 <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha)
.expr21 <- -1 - alpha
.expr26 <- u2^.expr21
.expr28 <- u3^.expr21
.expr30 <- u4^.expr21
.expr32 <- u5^.expr21
.expr34 <- u6^.expr21
.expr36 <- u7^.expr21
.expr38 <- u8^.expr21
.expr41 <- -alpha
.expr57 <- -7 + u1^.expr41 + u2^.expr41 + u3^.expr41 + u4^.expr41 +
u5^.expr41 + u6^.expr41 + u7^.expr41 + u8^.expr41
.expr60 <- -8 - 1/alpha
.value <- .expr19 * (u1^(.expr21 - 1) * .expr21) * .expr26 *
.expr28 * .expr30 * .expr32 * .expr34 * .expr36 * .expr38 *
.expr57^.expr60 + .expr19 * u1^.expr21 * .expr26 * .expr28 *
.expr30 * .expr32 * .expr34 * .expr36 * .expr38 * (.expr57^(.expr60 -
1) * (.expr60 * (u1^(.expr41 - 1) * .expr41)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr22 <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha)
.expr24 <- -1 - alpha
.expr29 <- u2^.expr24
.expr31 <- u3^.expr24
.expr33 <- u4^.expr24
.expr35 <- u5^.expr24
.expr37 <- u6^.expr24
.expr39 <- u7^.expr24
.expr41 <- u8^.expr24
.expr43 <- u9^.expr24
.expr46 <- -alpha
.expr64 <- -8 + u1^.expr46 + u2^.expr46 + u3^.expr46 + u4^.expr46 +
u5^.expr46 + u6^.expr46 + u7^.expr46 + u8^.expr46 + u9^.expr46
.expr67 <- -9 - 1/alpha
.value <- .expr22 * (u1^(.expr24 - 1) * .expr24) * .expr29 *
.expr31 * .expr33 * .expr35 * .expr37 * .expr39 * .expr41 *
.expr43 * .expr64^.expr67 + .expr22 * u1^.expr24 * .expr29 *
.expr31 * .expr33 * .expr35 * .expr37 * .expr39 * .expr41 *
.expr43 * (.expr64^(.expr67 - 1) * (.expr67 * (u1^(.expr46 -
1) * .expr46)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr25 <- (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * (1 + 9 * alpha)
.expr27 <- -1 - alpha
.expr32 <- u10^.expr27
.expr34 <- u2^.expr27
.expr36 <- u3^.expr27
.expr38 <- u4^.expr27
.expr40 <- u5^.expr27
.expr42 <- u6^.expr27
.expr44 <- u7^.expr27
.expr46 <- u8^.expr27
.expr48 <- u9^.expr27
.expr51 <- -alpha
.expr71 <- -9 + u1^.expr51 + u10^.expr51 + u2^.expr51 + u3^.expr51 +
u4^.expr51 + u5^.expr51 + u6^.expr51 + u7^.expr51 + u8^.expr51 +
u9^.expr51
.expr74 <- -10 - 1/alpha
.value <- .expr25 * (u1^(.expr27 - 1) * .expr27) * .expr32 *
.expr34 * .expr36 * .expr38 * .expr40 * .expr42 * .expr44 *
.expr46 * .expr48 * .expr71^.expr74 + .expr25 * u1^.expr27 *
.expr32 * .expr34 * .expr36 * .expr38 * .expr40 * .expr42 *
.expr44 * .expr46 * .expr48 * (.expr71^(.expr74 - 1) *
(.expr74 * (u1^(.expr51 - 1) * .expr51)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
})
`claytonCopula.pdfDerWrtPar.expr` <-
expression(u1 * ((u1^(-alpha))^(1/alpha) * (log((u1^(-alpha))) *
(1/alpha^2)) + (u1^(-alpha))^((1/alpha) - 1) * ((1/alpha) *
(u1^(-alpha) * log(u1))))/(u1 * (u1^(-alpha))^(1/alpha))^2,
((u1^(-1 - alpha) - (1 + alpha) * (u1^(-1 - alpha) * log(u1))) *
u2^(-1 - alpha) - (1 + alpha) * u1^(-1 - alpha) * (u2^(-1 -
alpha) * log(u2))) * (-1 + u1^(-alpha) + u2^(-alpha))^(-2 -
1/alpha) + (1 + alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
((-1 + u1^(-alpha) + u2^(-alpha))^(-2 - 1/alpha) * (log((-1 +
u1^(-alpha) + u2^(-alpha))) * (1/alpha^2)) - (-1 +
u1^(-alpha) + u2^(-alpha))^((-2 - 1/alpha) - 1) *
((-2 - 1/alpha) * (u2^(-alpha) * log(u2) + u1^(-alpha) *
log(u1)))), (((((1 + 2 * alpha) + (1 + alpha) *
2) * u1^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) *
(u1^(-1 - alpha) * log(u1))) * u2^(-1 - alpha) - (1 +
alpha) * (1 + 2 * alpha) * u1^(-1 - alpha) * (u2^(-1 -
alpha) * log(u2))) * u3^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
(u3^(-1 - alpha) * log(u3))) * (-2 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha))^(-3 - 1/alpha) + (1 + alpha) * (1 + 2 *
alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 -
alpha) * ((-2 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha))^(-3 -
1/alpha) * (log((-2 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha))) *
(1/alpha^2)) - (-2 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha))^((-3 -
1/alpha) - 1) * ((-3 - 1/alpha) * (u3^(-alpha) * log(u3) +
(u2^(-alpha) * log(u2) + u1^(-alpha) * log(u1))))), (((((((1 +
2 * alpha) + (1 + alpha) * 2) * (1 + 3 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * 3) * u1^(-1 - alpha) - (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (u1^(-1 -
alpha) * log(u1))) * u2^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * u1^(-1 - alpha) *
(u2^(-1 - alpha) * log(u2))) * u3^(-1 - alpha) - (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * (u3^(-1 - alpha) * log(u3))) *
u4^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 -
alpha) * (u4^(-1 - alpha) * log(u4))) * (-3 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha))^(-4 - 1/alpha) +
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 -
alpha) * ((-3 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha))^(-4 - 1/alpha) * (log((-3 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha))) * (1/alpha^2)) -
(-3 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha))^((-4 -
1/alpha) - 1) * ((-4 - 1/alpha) * (u4^(-alpha) *
log(u4) + (u3^(-alpha) * log(u3) + (u2^(-alpha) *
log(u2) + u1^(-alpha) * log(u1)))))), (((((((((1 +
2 * alpha) + (1 + alpha) * 2) * (1 + 3 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * 3) * (1 + 4 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * 4) * u1^(-1 -
alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (u1^(-1 - alpha) * log(u1))) * u2^(-1 -
alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * u1^(-1 - alpha) * (u2^(-1 - alpha) *
log(u2))) * u3^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * (u3^(-1 - alpha) * log(u3))) *
u4^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * (u4^(-1 - alpha) * log(u4))) *
u5^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * (u5^(-1 -
alpha) * log(u5))) * (-4 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha))^(-5 - 1/alpha) +
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 +
4 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) *
((-4 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha))^(-5 - 1/alpha) * (log((-4 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha))) * (1/alpha^2)) - (-4 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha))^((-5 -
1/alpha) - 1) * ((-5 - 1/alpha) * (u5^(-alpha) *
log(u5) + (u4^(-alpha) * log(u4) + (u3^(-alpha) *
log(u3) + (u2^(-alpha) * log(u2) + u1^(-alpha) *
log(u1))))))), (((((((((((1 + 2 * alpha) + (1 +
alpha) * 2) * (1 + 3 * alpha) + (1 + alpha) * (1 + 2 *
alpha) * 3) * (1 + 4 * alpha) + (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * 4) * (1 + 5 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * 5) * u1^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (u1^(-1 - alpha) * log(u1))) * u2^(-1 - alpha) -
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 +
4 * alpha) * (1 + 5 * alpha) * u1^(-1 - alpha) *
(u2^(-1 - alpha) * log(u2))) * u3^(-1 - alpha) -
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 +
4 * alpha) * (1 + 5 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * (u3^(-1 - alpha) * log(u3))) *
u4^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * (u4^(-1 -
alpha) * log(u4))) * u5^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * (u5^(-1 - alpha) *
log(u5))) * u6^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 -
alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) * (u6^(-1 -
alpha) * log(u6))) * (-5 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha))^(-6 -
1/alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * u6^(-1 - alpha) * ((-5 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha))^(-6 - 1/alpha) * (log((-5 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha))) * (1/alpha^2)) - (-5 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha))^((-6 -
1/alpha) - 1) * ((-6 - 1/alpha) * (u6^(-alpha) * log(u6) +
(u5^(-alpha) * log(u5) + (u4^(-alpha) * log(u4) + (u3^(-alpha) *
log(u3) + (u2^(-alpha) * log(u2) + u1^(-alpha) *
log(u1)))))))), (((((((((((((1 + 2 * alpha) + (1 +
alpha) * 2) * (1 + 3 * alpha) + (1 + alpha) * (1 + 2 *
alpha) * 3) * (1 + 4 * alpha) + (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * 4) * (1 + 5 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * 5) * (1 + 6 * alpha) + (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * 6) * u1^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * (u1^(-1 - alpha) * log(u1))) *
u2^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * u1^(-1 - alpha) * (u2^(-1 - alpha) * log(u2))) *
u3^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * (u3^(-1 -
alpha) * log(u3))) * u4^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * (u4^(-1 - alpha) *
log(u4))) * u5^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * (u5^(-1 -
alpha) * log(u5))) * u6^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * (u6^(-1 - alpha) * log(u6))) * u7^(-1 -
alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) *
u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) *
(u7^(-1 - alpha) * log(u7))) * (-6 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha))^(-7 - 1/alpha) + (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 -
alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) * ((-6 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha))^(-7 - 1/alpha) * (log((-6 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha))) * (1/alpha^2)) -
(-6 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha))^((-7 - 1/alpha) -
1) * ((-7 - 1/alpha) * (u7^(-alpha) * log(u7) + (u6^(-alpha) *
log(u6) + (u5^(-alpha) * log(u5) + (u4^(-alpha) *
log(u4) + (u3^(-alpha) * log(u3) + (u2^(-alpha) *
log(u2) + u1^(-alpha) * log(u1))))))))), (((((((((((((((1 +
2 * alpha) + (1 + alpha) * 2) * (1 + 3 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * 3) * (1 + 4 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * 4) * (1 +
5 * alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 *
alpha) * (1 + 4 * alpha) * 5) * (1 + 6 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * 6) * (1 + 7 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * 7) * u1^(-1 -
alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (u1^(-1 - alpha) * log(u1))) * u2^(-1 -
alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * u1^(-1 - alpha) * (u2^(-1 - alpha) *
log(u2))) * u3^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * (u3^(-1 - alpha) * log(u3))) *
u4^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * (1 + 7 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * (u4^(-1 - alpha) * log(u4))) *
u5^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * (1 + 7 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * (u5^(-1 -
alpha) * log(u5))) * u6^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) *
u4^(-1 - alpha) * u5^(-1 - alpha) * (u6^(-1 - alpha) *
log(u6))) * u7^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 -
alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) * (u7^(-1 -
alpha) * log(u7))) * u8^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) *
u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) *
u7^(-1 - alpha) * (u8^(-1 - alpha) * log(u8))) * (-7 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^(-8 -
1/alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) *
u6^(-1 - alpha) * u7^(-1 - alpha) * u8^(-1 - alpha) *
((-7 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^(-8 -
1/alpha) * (log((-7 + u1^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha) + u8^(-alpha))) * (1/alpha^2)) - (-7 +
u1^(-alpha) + u2^(-alpha) + u3^(-alpha) + u4^(-alpha) +
u5^(-alpha) + u6^(-alpha) + u7^(-alpha) + u8^(-alpha))^((-8 -
1/alpha) - 1) * ((-8 - 1/alpha) * (u8^(-alpha) *
log(u8) + (u7^(-alpha) * log(u7) + (u6^(-alpha) *
log(u6) + (u5^(-alpha) * log(u5) + (u4^(-alpha) *
log(u4) + (u3^(-alpha) * log(u3) + (u2^(-alpha) *
log(u2) + u1^(-alpha) * log(u1)))))))))), (((((((((((((((((1 +
2 * alpha) + (1 + alpha) * 2) * (1 + 3 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * 3) * (1 + 4 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * 4) * (1 +
5 * alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 *
alpha) * (1 + 4 * alpha) * 5) * (1 + 6 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * 6) * (1 + 7 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * 7) * (1 +
8 * alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 *
alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 *
alpha) * (1 + 7 * alpha) * 8) * u1^(-1 - alpha) - (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 *
alpha) * (1 + 8 * alpha) * (u1^(-1 - alpha) * log(u1))) *
u2^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * (1 + 7 * alpha) * (1 + 8 * alpha) * u1^(-1 -
alpha) * (u2^(-1 - alpha) * log(u2))) * u3^(-1 - alpha) -
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 +
4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * (u3^(-1 - alpha) * log(u3))) *
u4^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * (1 + 7 * alpha) * (1 + 8 * alpha) * u1^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * (u4^(-1 -
alpha) * log(u4))) * u5^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
(1 + 8 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * (u5^(-1 - alpha) *
log(u5))) * u6^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) * (1 + 8 *
alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 -
alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) * (u6^(-1 -
alpha) * log(u6))) * u7^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
(1 + 8 * alpha) * u1^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) *
u6^(-1 - alpha) * (u7^(-1 - alpha) * log(u7))) * u8^(-1 -
alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) *
(u8^(-1 - alpha) * log(u8))) * u9^(-1 - alpha) - (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 *
alpha) * (1 + 8 * alpha) * u1^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 -
alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) * u8^(-1 -
alpha) * (u9^(-1 - alpha) * log(u9))) * (-8 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-9 -
1/alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * u1^(-1 - alpha) *
u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) *
u5^(-1 - alpha) * u6^(-1 - alpha) * u7^(-1 - alpha) *
u8^(-1 - alpha) * u9^(-1 - alpha) * ((-8 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-9 -
1/alpha) * (log((-8 + u1^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) +
u8^(-alpha) + u9^(-alpha))) * (1/alpha^2)) - (-8 + u1^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^((-9 -
1/alpha) - 1) * ((-9 - 1/alpha) * (u9^(-alpha) * log(u9) +
(u8^(-alpha) * log(u8) + (u7^(-alpha) * log(u7) + (u6^(-alpha) *
log(u6) + (u5^(-alpha) * log(u5) + (u4^(-alpha) *
log(u4) + (u3^(-alpha) * log(u3) + (u2^(-alpha) *
log(u2) + u1^(-alpha) * log(u1))))))))))), (((((((((((((((((((1 +
2 * alpha) + (1 + alpha) * 2) * (1 + 3 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * 3) * (1 + 4 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * 4) * (1 +
5 * alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 *
alpha) * (1 + 4 * alpha) * 5) * (1 + 6 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * 6) * (1 + 7 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * 7) * (1 +
8 * alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 *
alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 *
alpha) * (1 + 7 * alpha) * 8) * (1 + 9 * alpha) + (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 *
alpha) * (1 + 8 * alpha) * 9) * u1^(-1 - alpha) - (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 *
alpha) * (1 + 8 * alpha) * (1 + 9 * alpha) * (u1^(-1 -
alpha) * log(u1))) * u10^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
(1 + 8 * alpha) * (1 + 9 * alpha) * u1^(-1 - alpha) *
(u10^(-1 - alpha) * log(u10))) * u2^(-1 - alpha) - (1 +
alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 *
alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 *
alpha) * (1 + 8 * alpha) * (1 + 9 * alpha) * u1^(-1 -
alpha) * u10^(-1 - alpha) * (u2^(-1 - alpha) * log(u2))) *
u3^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * (1 + 7 * alpha) * (1 + 8 * alpha) * (1 +
9 * alpha) * u1^(-1 - alpha) * u10^(-1 - alpha) * u2^(-1 -
alpha) * (u3^(-1 - alpha) * log(u3))) * u4^(-1 - alpha) -
(1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) * (1 +
4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * (1 + 9 * alpha) *
u1^(-1 - alpha) * u10^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * (u4^(-1 - alpha) * log(u4))) *
u5^(-1 - alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 +
3 * alpha) * (1 + 4 * alpha) * (1 + 5 * alpha) * (1 +
6 * alpha) * (1 + 7 * alpha) * (1 + 8 * alpha) * (1 +
9 * alpha) * u1^(-1 - alpha) * u10^(-1 - alpha) * u2^(-1 -
alpha) * u3^(-1 - alpha) * u4^(-1 - alpha) * (u5^(-1 -
alpha) * log(u5))) * u6^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
(1 + 8 * alpha) * (1 + 9 * alpha) * u1^(-1 - alpha) *
u10^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) *
u4^(-1 - alpha) * u5^(-1 - alpha) * (u6^(-1 - alpha) *
log(u6))) * u7^(-1 - alpha) - (1 + alpha) * (1 + 2 *
alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) * (1 + 5 *
alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) * (1 + 8 *
alpha) * (1 + 9 * alpha) * u1^(-1 - alpha) * u10^(-1 -
alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) * u4^(-1 -
alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) * (u7^(-1 -
alpha) * log(u7))) * u8^(-1 - alpha) - (1 + alpha) *
(1 + 2 * alpha) * (1 + 3 * alpha) * (1 + 4 * alpha) *
(1 + 5 * alpha) * (1 + 6 * alpha) * (1 + 7 * alpha) *
(1 + 8 * alpha) * (1 + 9 * alpha) * u1^(-1 - alpha) *
u10^(-1 - alpha) * u2^(-1 - alpha) * u3^(-1 - alpha) *
u4^(-1 - alpha) * u5^(-1 - alpha) * u6^(-1 - alpha) *
u7^(-1 - alpha) * (u8^(-1 - alpha) * log(u8))) * u9^(-1 -
alpha) - (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * (1 + 9 * alpha) *
u1^(-1 - alpha) * u10^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) *
u6^(-1 - alpha) * u7^(-1 - alpha) * u8^(-1 - alpha) *
(u9^(-1 - alpha) * log(u9))) * (-9 + u1^(-alpha) + u10^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-10 -
1/alpha) + (1 + alpha) * (1 + 2 * alpha) * (1 + 3 * alpha) *
(1 + 4 * alpha) * (1 + 5 * alpha) * (1 + 6 * alpha) *
(1 + 7 * alpha) * (1 + 8 * alpha) * (1 + 9 * alpha) *
u1^(-1 - alpha) * u10^(-1 - alpha) * u2^(-1 - alpha) *
u3^(-1 - alpha) * u4^(-1 - alpha) * u5^(-1 - alpha) *
u6^(-1 - alpha) * u7^(-1 - alpha) * u8^(-1 - alpha) *
u9^(-1 - alpha) * ((-9 + u1^(-alpha) + u10^(-alpha) +
u2^(-alpha) + u3^(-alpha) + u4^(-alpha) + u5^(-alpha) +
u6^(-alpha) + u7^(-alpha) + u8^(-alpha) + u9^(-alpha))^(-10 -
1/alpha) * (log((-9 + u1^(-alpha) + u10^(-alpha) + u2^(-alpha) +
u3^(-alpha) + u4^(-alpha) + u5^(-alpha) + u6^(-alpha) +
u7^(-alpha) + u8^(-alpha) + u9^(-alpha))) * (1/alpha^2)) -
(-9 + u1^(-alpha) + u10^(-alpha) + u2^(-alpha) + u3^(-alpha) +
u4^(-alpha) + u5^(-alpha) + u6^(-alpha) + u7^(-alpha) +
u8^(-alpha) + u9^(-alpha))^((-10 - 1/alpha) - 1) *
((-10 - 1/alpha) * (u9^(-alpha) * log(u9) + (u8^(-alpha) *
log(u8) + (u7^(-alpha) * log(u7) + (u6^(-alpha) *
log(u6) + (u5^(-alpha) * log(u5) + (u4^(-alpha) *
log(u4) + (u3^(-alpha) * log(u3) + (u2^(-alpha) *
log(u2) + (u10^(-alpha) * log(u10) + u1^(-alpha) *
log(u1)))))))))))))
`claytonCopula.pdfDerWrtPar.algr` <-
expression({
.expr2 <- u1^-alpha
.expr3 <- 1/alpha
.expr4 <- .expr2^.expr3
.value <- u1 * (.expr4 * (log(.expr2) * (1/alpha^2)) + .expr2^(.expr3 -
1) * (.expr3 * (.expr2 * log(u1))))/(u1 * .expr4)^2
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr1 <- -1
.expr2 <- .expr1 - alpha
.expr3 <- u1^.expr2
.expr4 <- 1 + alpha
.expr5 <- log(u1)
.expr9 <- u2^.expr2
.expr11 <- .expr4 * .expr3
.expr12 <- log(u2)
.expr16 <- -alpha
.expr17 <- u1^.expr16
.expr19 <- u2^.expr16
.expr20 <- .expr1 + .expr17 + .expr19
.expr23 <- -2 - 1/alpha
.expr24 <- .expr20^.expr23
.value <- ((.expr3 - .expr4 * (.expr3 * .expr5)) * .expr9 -
.expr11 * (.expr9 * .expr12)) * .expr24 + .expr11 * .expr9 *
(.expr24 * (log(.expr20) * (1/alpha^2)) - .expr20^(.expr23 -
1) * (.expr23 * (.expr19 * .expr12 + .expr17 * .expr5)))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- -1 - alpha
.expr8 <- u1^.expr7
.expr10 <- .expr3 * .expr2
.expr11 <- log(u1)
.expr15 <- u2^.expr7
.expr17 <- .expr10 * .expr8
.expr18 <- log(u2)
.expr22 <- u3^.expr7
.expr24 <- .expr17 * .expr15
.expr25 <- log(u3)
.expr30 <- -alpha
.expr31 <- u1^.expr30
.expr33 <- u2^.expr30
.expr35 <- u3^.expr30
.expr36 <- -2 + .expr31 + .expr33 + .expr35
.expr39 <- -3 - 1/alpha
.expr40 <- .expr36^.expr39
.value <- ((((.expr2 + .expr3 * 2) * .expr8 - .expr10 * (.expr8 *
.expr11)) * .expr15 - .expr17 * (.expr15 * .expr18)) *
.expr22 - .expr24 * (.expr22 * .expr25)) * .expr40 +
.expr24 * .expr22 * (.expr40 * (log(.expr36) * (1/alpha^2)) -
.expr36^(.expr39 - 1) * (.expr39 * (.expr35 * .expr25 +
(.expr33 * .expr18 + .expr31 * .expr11))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- 1 + 3 * alpha
.expr9 <- .expr3 * .expr2
.expr13 <- -1 - alpha
.expr14 <- u1^.expr13
.expr16 <- .expr9 * .expr7
.expr17 <- log(u1)
.expr21 <- u2^.expr13
.expr23 <- .expr16 * .expr14
.expr24 <- log(u2)
.expr28 <- u3^.expr13
.expr30 <- .expr23 * .expr21
.expr31 <- log(u3)
.expr35 <- u4^.expr13
.expr37 <- .expr30 * .expr28
.expr38 <- log(u4)
.expr43 <- -alpha
.expr44 <- u1^.expr43
.expr46 <- u2^.expr43
.expr48 <- u3^.expr43
.expr50 <- u4^.expr43
.expr51 <- -3 + .expr44 + .expr46 + .expr48 + .expr50
.expr54 <- -4 - 1/alpha
.expr55 <- .expr51^.expr54
.value <- ((((((.expr2 + .expr3 * 2) * .expr7 + .expr9 *
3) * .expr14 - .expr16 * (.expr14 * .expr17)) * .expr21 -
.expr23 * (.expr21 * .expr24)) * .expr28 - .expr30 *
(.expr28 * .expr31)) * .expr35 - .expr37 * (.expr35 *
.expr38)) * .expr55 + .expr37 * .expr35 * (.expr55 *
(log(.expr51) * (1/alpha^2)) - .expr51^(.expr54 - 1) *
(.expr54 * (.expr50 * .expr38 + (.expr48 * .expr31 +
(.expr46 * .expr24 + .expr44 * .expr17)))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- 1 + 3 * alpha
.expr9 <- .expr3 * .expr2
.expr13 <- 1 + 4 * alpha
.expr15 <- .expr9 * .expr7
.expr19 <- -1 - alpha
.expr20 <- u1^.expr19
.expr22 <- .expr15 * .expr13
.expr23 <- log(u1)
.expr27 <- u2^.expr19
.expr29 <- .expr22 * .expr20
.expr30 <- log(u2)
.expr34 <- u3^.expr19
.expr36 <- .expr29 * .expr27
.expr37 <- log(u3)
.expr41 <- u4^.expr19
.expr43 <- .expr36 * .expr34
.expr44 <- log(u4)
.expr48 <- u5^.expr19
.expr50 <- .expr43 * .expr41
.expr51 <- log(u5)
.expr56 <- -alpha
.expr57 <- u1^.expr56
.expr59 <- u2^.expr56
.expr61 <- u3^.expr56
.expr63 <- u4^.expr56
.expr65 <- u5^.expr56
.expr66 <- -4 + .expr57 + .expr59 + .expr61 + .expr63 + .expr65
.expr69 <- -5 - 1/alpha
.expr70 <- .expr66^.expr69
.value <- ((((((((.expr2 + .expr3 * 2) * .expr7 + .expr9 *
3) * .expr13 + .expr15 * 4) * .expr20 - .expr22 * (.expr20 *
.expr23)) * .expr27 - .expr29 * (.expr27 * .expr30)) *
.expr34 - .expr36 * (.expr34 * .expr37)) * .expr41 -
.expr43 * (.expr41 * .expr44)) * .expr48 - .expr50 *
(.expr48 * .expr51)) * .expr70 + .expr50 * .expr48 *
(.expr70 * (log(.expr66) * (1/alpha^2)) - .expr66^(.expr69 -
1) * (.expr69 * (.expr65 * .expr51 + (.expr63 * .expr44 +
(.expr61 * .expr37 + (.expr59 * .expr30 + .expr57 *
.expr23))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- 1 + 3 * alpha
.expr9 <- .expr3 * .expr2
.expr13 <- 1 + 4 * alpha
.expr15 <- .expr9 * .expr7
.expr19 <- 1 + 5 * alpha
.expr21 <- .expr15 * .expr13
.expr25 <- -1 - alpha
.expr26 <- u1^.expr25
.expr28 <- .expr21 * .expr19
.expr29 <- log(u1)
.expr33 <- u2^.expr25
.expr35 <- .expr28 * .expr26
.expr36 <- log(u2)
.expr40 <- u3^.expr25
.expr42 <- .expr35 * .expr33
.expr43 <- log(u3)
.expr47 <- u4^.expr25
.expr49 <- .expr42 * .expr40
.expr50 <- log(u4)
.expr54 <- u5^.expr25
.expr56 <- .expr49 * .expr47
.expr57 <- log(u5)
.expr61 <- u6^.expr25
.expr63 <- .expr56 * .expr54
.expr64 <- log(u6)
.expr69 <- -alpha
.expr70 <- u1^.expr69
.expr72 <- u2^.expr69
.expr74 <- u3^.expr69
.expr76 <- u4^.expr69
.expr78 <- u5^.expr69
.expr80 <- u6^.expr69
.expr81 <- -5 + .expr70 + .expr72 + .expr74 + .expr76 + .expr78 +
.expr80
.expr84 <- -6 - 1/alpha
.expr85 <- .expr81^.expr84
.value <- ((((((((((.expr2 + .expr3 * 2) * .expr7 + .expr9 *
3) * .expr13 + .expr15 * 4) * .expr19 + .expr21 * 5) *
.expr26 - .expr28 * (.expr26 * .expr29)) * .expr33 -
.expr35 * (.expr33 * .expr36)) * .expr40 - .expr42 *
(.expr40 * .expr43)) * .expr47 - .expr49 * (.expr47 *
.expr50)) * .expr54 - .expr56 * (.expr54 * .expr57)) *
.expr61 - .expr63 * (.expr61 * .expr64)) * .expr85 +
.expr63 * .expr61 * (.expr85 * (log(.expr81) * (1/alpha^2)) -
.expr81^(.expr84 - 1) * (.expr84 * (.expr80 * .expr64 +
(.expr78 * .expr57 + (.expr76 * .expr50 + (.expr74 *
.expr43 + (.expr72 * .expr36 + .expr70 * .expr29)))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- 1 + 3 * alpha
.expr9 <- .expr3 * .expr2
.expr13 <- 1 + 4 * alpha
.expr15 <- .expr9 * .expr7
.expr19 <- 1 + 5 * alpha
.expr21 <- .expr15 * .expr13
.expr25 <- 1 + 6 * alpha
.expr27 <- .expr21 * .expr19
.expr31 <- -1 - alpha
.expr32 <- u1^.expr31
.expr34 <- .expr27 * .expr25
.expr35 <- log(u1)
.expr39 <- u2^.expr31
.expr41 <- .expr34 * .expr32
.expr42 <- log(u2)
.expr46 <- u3^.expr31
.expr48 <- .expr41 * .expr39
.expr49 <- log(u3)
.expr53 <- u4^.expr31
.expr55 <- .expr48 * .expr46
.expr56 <- log(u4)
.expr60 <- u5^.expr31
.expr62 <- .expr55 * .expr53
.expr63 <- log(u5)
.expr67 <- u6^.expr31
.expr69 <- .expr62 * .expr60
.expr70 <- log(u6)
.expr74 <- u7^.expr31
.expr76 <- .expr69 * .expr67
.expr77 <- log(u7)
.expr82 <- -alpha
.expr83 <- u1^.expr82
.expr85 <- u2^.expr82
.expr87 <- u3^.expr82
.expr89 <- u4^.expr82
.expr91 <- u5^.expr82
.expr93 <- u6^.expr82
.expr95 <- u7^.expr82
.expr96 <- -6 + .expr83 + .expr85 + .expr87 + .expr89 + .expr91 +
.expr93 + .expr95
.expr99 <- -7 - 1/alpha
.expr100 <- .expr96^.expr99
.value <- ((((((((((((.expr2 + .expr3 * 2) * .expr7 + .expr9 *
3) * .expr13 + .expr15 * 4) * .expr19 + .expr21 * 5) *
.expr25 + .expr27 * 6) * .expr32 - .expr34 * (.expr32 *
.expr35)) * .expr39 - .expr41 * (.expr39 * .expr42)) *
.expr46 - .expr48 * (.expr46 * .expr49)) * .expr53 -
.expr55 * (.expr53 * .expr56)) * .expr60 - .expr62 *
(.expr60 * .expr63)) * .expr67 - .expr69 * (.expr67 *
.expr70)) * .expr74 - .expr76 * (.expr74 * .expr77)) *
.expr100 + .expr76 * .expr74 * (.expr100 * (log(.expr96) *
(1/alpha^2)) - .expr96^(.expr99 - 1) * (.expr99 * (.expr95 *
.expr77 + (.expr93 * .expr70 + (.expr91 * .expr63 + (.expr89 *
.expr56 + (.expr87 * .expr49 + (.expr85 * .expr42 + .expr83 *
.expr35))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- 1 + 3 * alpha
.expr9 <- .expr3 * .expr2
.expr13 <- 1 + 4 * alpha
.expr15 <- .expr9 * .expr7
.expr19 <- 1 + 5 * alpha
.expr21 <- .expr15 * .expr13
.expr25 <- 1 + 6 * alpha
.expr27 <- .expr21 * .expr19
.expr31 <- 1 + 7 * alpha
.expr33 <- .expr27 * .expr25
.expr37 <- -1 - alpha
.expr38 <- u1^.expr37
.expr40 <- .expr33 * .expr31
.expr41 <- log(u1)
.expr45 <- u2^.expr37
.expr47 <- .expr40 * .expr38
.expr48 <- log(u2)
.expr52 <- u3^.expr37
.expr54 <- .expr47 * .expr45
.expr55 <- log(u3)
.expr59 <- u4^.expr37
.expr61 <- .expr54 * .expr52
.expr62 <- log(u4)
.expr66 <- u5^.expr37
.expr68 <- .expr61 * .expr59
.expr69 <- log(u5)
.expr73 <- u6^.expr37
.expr75 <- .expr68 * .expr66
.expr76 <- log(u6)
.expr80 <- u7^.expr37
.expr82 <- .expr75 * .expr73
.expr83 <- log(u7)
.expr87 <- u8^.expr37
.expr89 <- .expr82 * .expr80
.expr90 <- log(u8)
.expr95 <- -alpha
.expr96 <- u1^.expr95
.expr98 <- u2^.expr95
.expr100 <- u3^.expr95
.expr102 <- u4^.expr95
.expr104 <- u5^.expr95
.expr106 <- u6^.expr95
.expr108 <- u7^.expr95
.expr110 <- u8^.expr95
.expr111 <- -7 + .expr96 + .expr98 + .expr100 + .expr102 +
.expr104 + .expr106 + .expr108 + .expr110
.expr114 <- -8 - 1/alpha
.expr115 <- .expr111^.expr114
.value <- ((((((((((((((.expr2 + .expr3 * 2) * .expr7 + .expr9 *
3) * .expr13 + .expr15 * 4) * .expr19 + .expr21 * 5) *
.expr25 + .expr27 * 6) * .expr31 + .expr33 * 7) * .expr38 -
.expr40 * (.expr38 * .expr41)) * .expr45 - .expr47 *
(.expr45 * .expr48)) * .expr52 - .expr54 * (.expr52 *
.expr55)) * .expr59 - .expr61 * (.expr59 * .expr62)) *
.expr66 - .expr68 * (.expr66 * .expr69)) * .expr73 -
.expr75 * (.expr73 * .expr76)) * .expr80 - .expr82 *
(.expr80 * .expr83)) * .expr87 - .expr89 * (.expr87 *
.expr90)) * .expr115 + .expr89 * .expr87 * (.expr115 *
(log(.expr111) * (1/alpha^2)) - .expr111^(.expr114 -
1) * (.expr114 * (.expr110 * .expr90 + (.expr108 * .expr83 +
(.expr106 * .expr76 + (.expr104 * .expr69 + (.expr102 *
.expr62 + (.expr100 * .expr55 + (.expr98 * .expr48 +
.expr96 * .expr41)))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- 1 + 3 * alpha
.expr9 <- .expr3 * .expr2
.expr13 <- 1 + 4 * alpha
.expr15 <- .expr9 * .expr7
.expr19 <- 1 + 5 * alpha
.expr21 <- .expr15 * .expr13
.expr25 <- 1 + 6 * alpha
.expr27 <- .expr21 * .expr19
.expr31 <- 1 + 7 * alpha
.expr33 <- .expr27 * .expr25
.expr37 <- 1 + 8 * alpha
.expr39 <- .expr33 * .expr31
.expr43 <- -1 - alpha
.expr44 <- u1^.expr43
.expr46 <- .expr39 * .expr37
.expr47 <- log(u1)
.expr51 <- u2^.expr43
.expr53 <- .expr46 * .expr44
.expr54 <- log(u2)
.expr58 <- u3^.expr43
.expr60 <- .expr53 * .expr51
.expr61 <- log(u3)
.expr65 <- u4^.expr43
.expr67 <- .expr60 * .expr58
.expr68 <- log(u4)
.expr72 <- u5^.expr43
.expr74 <- .expr67 * .expr65
.expr75 <- log(u5)
.expr79 <- u6^.expr43
.expr81 <- .expr74 * .expr72
.expr82 <- log(u6)
.expr86 <- u7^.expr43
.expr88 <- .expr81 * .expr79
.expr89 <- log(u7)
.expr93 <- u8^.expr43
.expr95 <- .expr88 * .expr86
.expr96 <- log(u8)
.expr100 <- u9^.expr43
.expr102 <- .expr95 * .expr93
.expr103 <- log(u9)
.expr108 <- -alpha
.expr109 <- u1^.expr108
.expr111 <- u2^.expr108
.expr113 <- u3^.expr108
.expr115 <- u4^.expr108
.expr117 <- u5^.expr108
.expr119 <- u6^.expr108
.expr121 <- u7^.expr108
.expr123 <- u8^.expr108
.expr125 <- u9^.expr108
.expr126 <- -8 + .expr109 + .expr111 + .expr113 + .expr115 +
.expr117 + .expr119 + .expr121 + .expr123 + .expr125
.expr129 <- -9 - 1/alpha
.expr130 <- .expr126^.expr129
.value <- ((((((((((((((((.expr2 + .expr3 * 2) * .expr7 +
.expr9 * 3) * .expr13 + .expr15 * 4) * .expr19 + .expr21 *
5) * .expr25 + .expr27 * 6) * .expr31 + .expr33 * 7) *
.expr37 + .expr39 * 8) * .expr44 - .expr46 * (.expr44 *
.expr47)) * .expr51 - .expr53 * (.expr51 * .expr54)) *
.expr58 - .expr60 * (.expr58 * .expr61)) * .expr65 -
.expr67 * (.expr65 * .expr68)) * .expr72 - .expr74 *
(.expr72 * .expr75)) * .expr79 - .expr81 * (.expr79 *
.expr82)) * .expr86 - .expr88 * (.expr86 * .expr89)) *
.expr93 - .expr95 * (.expr93 * .expr96)) * .expr100 -
.expr102 * (.expr100 * .expr103)) * .expr130 + .expr102 *
.expr100 * (.expr130 * (log(.expr126) * (1/alpha^2)) -
.expr126^(.expr129 - 1) * (.expr129 * (.expr125 * .expr103 +
(.expr123 * .expr96 + (.expr121 * .expr89 + (.expr119 *
.expr82 + (.expr117 * .expr75 + (.expr115 * .expr68 +
(.expr113 * .expr61 + (.expr111 * .expr54 + .expr109 *
.expr47))))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.expr2 <- 1 + 2 * alpha
.expr3 <- 1 + alpha
.expr7 <- 1 + 3 * alpha
.expr9 <- .expr3 * .expr2
.expr13 <- 1 + 4 * alpha
.expr15 <- .expr9 * .expr7
.expr19 <- 1 + 5 * alpha
.expr21 <- .expr15 * .expr13
.expr25 <- 1 + 6 * alpha
.expr27 <- .expr21 * .expr19
.expr31 <- 1 + 7 * alpha
.expr33 <- .expr27 * .expr25
.expr37 <- 1 + 8 * alpha
.expr39 <- .expr33 * .expr31
.expr43 <- 1 + 9 * alpha
.expr45 <- .expr39 * .expr37
.expr49 <- -1 - alpha
.expr50 <- u1^.expr49
.expr52 <- .expr45 * .expr43
.expr53 <- log(u1)
.expr57 <- u10^.expr49
.expr59 <- .expr52 * .expr50
.expr60 <- log(u10)
.expr64 <- u2^.expr49
.expr66 <- .expr59 * .expr57
.expr67 <- log(u2)
.expr71 <- u3^.expr49
.expr73 <- .expr66 * .expr64
.expr74 <- log(u3)
.expr78 <- u4^.expr49
.expr80 <- .expr73 * .expr71
.expr81 <- log(u4)
.expr85 <- u5^.expr49
.expr87 <- .expr80 * .expr78
.expr88 <- log(u5)
.expr92 <- u6^.expr49
.expr94 <- .expr87 * .expr85
.expr95 <- log(u6)
.expr99 <- u7^.expr49
.expr101 <- .expr94 * .expr92
.expr102 <- log(u7)
.expr106 <- u8^.expr49
.expr108 <- .expr101 * .expr99
.expr109 <- log(u8)
.expr113 <- u9^.expr49
.expr115 <- .expr108 * .expr106
.expr116 <- log(u9)
.expr121 <- -alpha
.expr122 <- u1^.expr121
.expr124 <- u10^.expr121
.expr126 <- u2^.expr121
.expr128 <- u3^.expr121
.expr130 <- u4^.expr121
.expr132 <- u5^.expr121
.expr134 <- u6^.expr121
.expr136 <- u7^.expr121
.expr138 <- u8^.expr121
.expr140 <- u9^.expr121
.expr141 <- -9 + .expr122 + .expr124 + .expr126 + .expr128 +
.expr130 + .expr132 + .expr134 + .expr136 + .expr138 +
.expr140
.expr144 <- -10 - 1/alpha
.expr145 <- .expr141^.expr144
.value <- ((((((((((((((((((.expr2 + .expr3 * 2) * .expr7 +
.expr9 * 3) * .expr13 + .expr15 * 4) * .expr19 + .expr21 *
5) * .expr25 + .expr27 * 6) * .expr31 + .expr33 * 7) *
.expr37 + .expr39 * 8) * .expr43 + .expr45 * 9) * .expr50 -
.expr52 * (.expr50 * .expr53)) * .expr57 - .expr59 *
(.expr57 * .expr60)) * .expr64 - .expr66 * (.expr64 *
.expr67)) * .expr71 - .expr73 * (.expr71 * .expr74)) *
.expr78 - .expr80 * (.expr78 * .expr81)) * .expr85 -
.expr87 * (.expr85 * .expr88)) * .expr92 - .expr94 *
(.expr92 * .expr95)) * .expr99 - .expr101 * (.expr99 *
.expr102)) * .expr106 - .expr108 * (.expr106 * .expr109)) *
.expr113 - .expr115 * (.expr113 * .expr116)) * .expr145 +
.expr115 * .expr113 * (.expr145 * (log(.expr141) * (1/alpha^2)) -
.expr141^(.expr144 - 1) * (.expr144 * (.expr140 *
.expr116 + (.expr138 * .expr109 + (.expr136 *
.expr102 + (.expr134 * .expr95 + (.expr132 *
.expr88 + (.expr130 * .expr81 + (.expr128 * .expr74 +
(.expr126 * .expr67 + (.expr124 * .expr60 + .expr122 *
.expr53)))))))))))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
})
`claytonCopula.genfunDer.expr` <-
expression(-(alpha * u^(-1 - alpha)), alpha * (1 + alpha) * u^(-2 -
alpha))
`claytonCopula.genfunDer.algr` <-
expression({
.value <- -(alpha * u^(-1 - alpha))
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
}, {
.value <- alpha * (1 + alpha) * u^(-2 - alpha)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("s")))
.grad[, "s"] <- 0
attr(.value, "gradient") <- .grad
.value
})
Try the copula package in your browser
Any scripts or data that you put into this service are public.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.