mvrcop-a4aGr: mvrcop

In flr/FLa4a: A Simple and Robust Statistical Catch at Age Model

mvrcop

Description

Method to generate multivariate parameters with user-defined copulas and marginals for a4aGr objects.

Usage

## S4 method for signature 'numeric,a4aGr'
mvrcop(n = 1, mvdc, ...)

Arguments

n	the number of iterations
mvdc	the a4aGr object
...	arguments to be passed to the rMvdc and copula methods

Value

an FLModelSim object with n groups of parameters

Examples

mm <- matrix(NA, ncol=3, nrow=3)
diag(mm) <- c(50, 0.001,0.001)
mm[upper.tri(mm)] <- mm[lower.tri(mm)] <- c(0.1,0.01,0.00004)
md <- ~linf*(1-exp(-k*(t-t0)))
imd <- ~t0-1/k*log(1-len/linf)
prs <- FLPar(linf=58.5, k=0.086, t0=0.001, units=c("cm","yr^-1","yr"))
vbObj <- a4aGr(grMod=md, grInvMod=imd, params=prs, vcov=mm, distr="norm")
pars <- list(list(a=50, b=100, c=58.5), list(a=0.06, b=0.2, c=0.086), list(a=0, b=0.005, c=0.001))
#In the following, the third, fourth and fifth arguments refer to the copula,
#  while the final two arguments refer to the marginal distributions:
vbObj <- mvrcop(10000, vbObj, copula="archmCopula", family="clayton", param=2, 
   margins="triangle", paramMargins=pars)
splom(data.frame(t(params(vbObj)@.Data)), pch=".")

flr/FLa4a documentation built on Dec. 21, 2024, 9:28 p.m.