R/CVer.R

In stats4sd/LSDer: Fishers LSD & %CV from Mixed Models

#' Calculating % coefficient of variation from mixed model
#' @param model An object of class \code{merMod}
#' @keywords CV lmer
#' @export
#' @examples

CVer<-function(model){
  if(class(model)=="merMod"|class(model)=="merModLmerTest"|class(model)=="lmerMod"|class(model)=="lmerModLmerTest"){
   
    GM<-mean(model@frame[,1])
    MS<-sigma(model)**2
return(paste(round(100*sqrt(MS)/GM,2),"%",sep=""))
  }
  else(stop("model not of class merMod or merModLmerTest"))
}