getpars: Extract parameters from Kfunc lme models

In robertbagchi/ReplicatedPointPatterns: Analysis of replicated point patterns with covariates

Description

Used internally to extract the parameters from a kfunclme object.

Usage

getpars(mod, lin.comb.Ct)

Arguments

mod	Model to get parameters from
lin.comb.Ct	Linear combination to be used for predictions

Author(s)

Robert Bagchi Maintainer: Robert Bagchi <[e-mail]>

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (mod, lin.comb.Ct) 
{
    if (!is.null(mod)) {
        beta.r <- fixef(mod)
        vcov.r <- vcov(mod)
        est.Kmean.r <- as.vector(lin.comb.Ct %*% fixef(mod))
        est.Kse.r <- sqrt(diag(lin.comb.Ct %*% vcov(mod) %*% 
            t(lin.comb.Ct)))
    }
    else {
        beta.r <- rep(NA, ncol(lin.comb.Ct))
        vcov.r <- matrix(NA, ncol = length(beta.r), nrow = length(beta.r))
        est.Kmean.r <- rep(NA, nrow(lin.comb.Ct))
        est.Kse.r <- rep(NA, nrow(lin.comb.Ct))
    }
    return(list(pred.r = est.Kmean.r, se.pred.r = est.Kse.r, 
        beta.r = beta.r, vcov.r = vcov.r))
  }

robertbagchi/ReplicatedPointPatterns documentation built on May 27, 2019, 10:32 a.m.