R/estSigmaR.R

Defines functions estSigmaR

Documented in estSigmaR

estSigmaR <- function(model, digits=2)
{
  ## 1  Parse args
  x <- if(class(model)=="scape") model$Dev else model  # allow data frame
  if(is.null(x))
    stop("element 'Dev' not found")

  ## 2  Calculate sigmaR
  ss <- sapply(x[c("Initial","Annual")], function(x) sum(x^2))
  n <- sapply(x[c("Initial","Annual")], length)
  output <- sqrt(ss/(n-1))
  if(!is.null(digits))
    output <- round(output, digits=digits)

  output
}

Try the scape package in your browser

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

scape documentation built on Nov. 23, 2020, 5:08 p.m.