residuals-methods: residuals method

Description Usage Arguments Value Author(s) Examples

Description

Provides the residuals of the fitted FMRs model from an fmrsfit-class

Usage

1
2
3
4
residuals(object, ...)

## S4 method for signature 'fmrsfit'
residuals(object, ...)

Arguments

object

An fmrsfit-class

...

Other possible arguments

Value

A numeric array of dimension-nObs-nComp

Author(s)

Farhad Shokoohi <shokoohi@icloud.com>

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
set.seed(1980)
nComp = 2
nCov = 10
nObs = 500
dispersion = c(1, 1)
mixProp = c(0.4, 0.6)
rho = 0.5
coeff1 = c( 2,  2, -1, -2, 1, 2, 0, 0,  0, 0,  0)
coeff2 = c(-1, -1,  1,  2, 0, 0, 0, 0, -1, 2, -2)
umax = 40

dat <- fmrs.gendata(nObs = nObs, nComp = nComp, nCov = nCov,
 coeff = c(coeff1, coeff2), dispersion = dispersion,
 mixProp = mixProp, rho = rho, umax = umax,
 disFamily = 'lnorm')

res.mle <- fmrs.mle(y = dat$y, x = dat$x, delta = dat$delta,
 nComp = nComp, disFamily = 'lnorm',
 initCoeff = rnorm(nComp*nCov+nComp),
 initDispersion = rep(1, nComp),
 initmixProp = rep(1/nComp, nComp))
head(residuals(res.mle))

fmrs documentation built on Nov. 8, 2020, 5:50 p.m.