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

wuhailongbbcc

4328

收藏 2013-07-28

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

扫码加我拉你入群

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

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

全部回复

DM小菜鸟

2015-2-14 15:27:55

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改掉，自己赋值

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群