请教 一下tCopula 中自由度df的下界

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

Copula Package 中的tCopula函数要求其参数df的下界大于1，但是作图时，我想看一下df趋于零（小于1)的情形，请问如何实现，谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Copula opula TCO 自由度 package 自由度

tCopula <- function(param = NA_real_, dim = 2L, dispstr = "ex",
                    df = 4, df.fixed = FALSE)
{
    dim <- as.integer(dim)
    stopifnot((pdim <- length(param)) >= 1, is.numeric(param))
    if(pdim == 1 && is.na(param)) ## extend it (rho only!)
        pdim <- length(param <- rep(param, length.out = npar.ellip(dim, dispstr)))
    parameters <- param
    param.names <- paste("rho", seq_len(pdim), sep=".")
    param.lowbnd <- rep.int(-1, pdim)
    param.upbnd         <- rep.int( 1, pdim)
    if (!df.fixed) { ## df is another parameter __at end__
        parameters <- c(parameters, df)
        param.names <- c(param.names, "df")
        param.lowbnd <- c(param.lowbnd, 1e-6)
        param.upbnd <- c(param.upbnd, Inf)
    }

    new("tCopula",
        dispstr = dispstr,
        dimension = dim,
        parameters = parameters,
        df = df,
        df.fixed = df.fixed,
        param.names = param.names,
        param.lowbnd = param.lowbnd,
        param.upbnd = param.upbnd,
        fullname = paste("t copula family", if(df.fixed) paste("df fixed at", df)),
        getRho = function(obj) {
            par <- obj@parameters
            if (obj@df.fixed) par else par[-length(par)]
        }
        )
}
把上面那句df改掉，自己赋值