R/variances.R

Defines functions variances offdiags

Documented in offdiags variances

### Return position of variance elements in the parameter vector (without mean parameters)
### Optimization constraints are needed on these parameters
##' @export
variances <- function(x,mean=FALSE) {
##  if (is.null(x$parpos))
##    x$parpos <- parpos(x)
  x$parpos <- parpos(Model(x),mean=TRUE) 
  res <- diag(x$parpos$P)[diag(index(x)$P0)==1]
  if (!mean) {
    return(res - index(x)$npar.mean)    
  }
  return(res)
}
## And the off-diagonal (covariance) parameters
##' @export
offdiags <- function(x,mean=FALSE) {
  parpos <- parpos(x,mean=mean)
  pp <- parpos$P
  pp[lower.tri(pp)][(index(x)$P0)[lower.tri(pp)]==1]
}

Try the lava package in your browser

Any scripts or data that you put into this service are public.

lava documentation built on May 2, 2019, 4:49 p.m.