calculateCC <- function(models, CC = 50, raw = FALSE ){
sdCC <- lapply(models, function(replicates){
suppressMessages(suppressWarnings(sd(unlist(nplr::getEstimates(replicates,1-(CC/100))[c(2,3,4)]))))
})
CC <- lapply(models, function(replicates){
suppressMessages(suppressWarnings(nplr::getEstimates(replicates,1-(CC/100))[3]))
})
seCC <- lapply(models, function(replicates){
replicates@stdErr[1]
})
aCC <- lapply(models,function(replicates){
replicates@goodness$gof
})
if (length(models)==0){
return (c(CC="NR", SD="NR", SE="NR",mADJ=paste("<")))
}else{
if (raw == TRUE){
cc <- lapply(models, function(replicates){
suppressMessages(suppressWarnings(nplr::getEstimates(replicates,1-(CC/100))))
})
} else{
cc <- return(c(CC = ifelse(test = is.nan(round(mean(unlist(CC),na.rm = TRUE),6)),yes = ">",no = ifelse(round(mean(unlist(CC),na.rm = TRUE),6)<0,yes = "<",no = round(mean(unlist(CC),na.rm = TRUE),6))), SD = round(mean(unlist(sdCC),na.rm = TRUE),6), SE = round(mean(unlist(seCC),na.rm = TRUE),6), mADJ = round(mean(unlist(aCC),na.rm = TRUE),6)))
}}
return(cc)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.