mvrcop

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

# S4 method for 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)
#> Warning: margins correct? Currently, have no function(s) named: ptriangle
#> Warning: margins correct? Currently, have no function(s) named: dtriangle
#> Error in qtriangle(x, a = 50, b = 100, c = 58.5): could not find function "qtriangle"
splom(data.frame(t(params(vbObj)@.Data)), pch=".")