Displaying 3 results from an estimated 3 matches for "stn_pos".
Did you mean:
src_pos
2010 Jun 09
0
fitting t copula
...help given.
myCop.t <- ellipCopula(family = "t", dim = 2, dispstr = "toep", param = 0.5, df = 8)
myCop.t
myMvd <- mvdc(copula = myCop.t, margins = c("gamma", "gamma"), paramMargins = list(list(mean = 0, sd = 2), list(mean = 0,
sd = 1)))
myMvd
dat <- stn_pos ## observed data
loglikMvdc(c(0, 2, 0, 1, 0.5,8), dat, myMvd) ## loglikelihood
mm <- apply(dat, 2, mean)
vv <- apply(dat, 2, var)
rho <- rcorr(stn_pos,type="spearman")[[1]]; round(rho,2)
rbind(mm,vv)
b1.0 <- c(mm[1]^2/vv[1], vv[1]/mm[1])
b2.0 <- c(mm[2]^2/vv[2], vv[2]/mm[...
2010 Jun 10
0
error message fitting tcopula
...Jun Yan in “Enjoy the joy of copula” but I’m not sure how to correct the error based on the error message. I tried so many ways but still could not get it working.
loglik.marg <- function(b, x) sum(dgamma(x, shape = b[1], scale = b[2], log = TRUE))
ctrl <- list(fnscale = -1)
#dat <- stn_pos[,1:2] ## observed data
myCop.t <- ellipCopula(family = "t", param = 0.5,dim = 2, dispstr = "un", df = 8)
myCop.t
myMvd <- mvdc(copula = myCop.t, margins = c("gamma", "gamma"), paramMargins = list(list(shape = 1.5, scale = 38),
list(shape = 1.7, scale...
2010 Jun 10
0
error message in fitting tcopula
...oy the joy of copula" but I'm not sure
how to correct the error based on the error message. I tried so many
ways but still could not get it working.
loglik.marg <- function(b, x) sum(dgamma(x, shape = b[1], scale = b[2],
log = TRUE))
ctrl <- list(fnscale = -1)
#dat <- stn_pos[,1:2] ## observed data
myCop.t <- ellipCopula(family = "t", param = 0.5,dim = 2, dispstr =
"un", df = 8)
myCop.t
myMvd <- mvdc(copula = myCop.t, margins = c("gamma", "gamma"),
paramMargins = list(list(shape = 1.5, scale = 38),
list(shape = 1.7,...