## print.summary.nscosinor.R
#' printing a summary of an nscosinor
#' @name print.summary.nscosinor
#' @method print summary.nscosinor
#' @param \dots further arguments passed to or from other methods.
#' @param x a \code{summary.nsCosinor} object produced by
#' \code{summary.nsCosinor}
#' @export
print.summary.nscosinor <- function(x, ...){
if (class(x)!="summary.nsCosinor"){
stop("Object must be of class 'summary.nsCosinor'")
}
k<-length(x$cycles);
cat("Statistics for non-stationary cosinor based on MCMC chains\n")
cat("Number of MCMC samples = ",x$niters-x$burnin+1,"\n",sep="", ...)
cat("\nStandard deviations\n")
cat("Residual, mean=",x$stats$errorstats[1],", 95% CI [",x$stats$errorstats[2],", ",x$stats$errorstats[3],"]\n",sep="", ...)
if(k==1){
cat("Cycle=",x$cycles,"\n",sep="", ...)
cat("Season, mean=",x$stats$wstats[1],", 95% CI [",x$stats$wstats[2],", ",x$stats$wstats[3],"]\n",sep="", ...)
cat("\nPhase and amplitude\n")
cat("Cycle=",x$cycles,"\n",sep="", ...)
cat("Amplitude, mean=",x$stats$ampstats[1],", 95% CI [",x$stats$ampstats[2],", ",x$stats$ampstats[3],"]\n",sep="", ...)
cat("Phase (radians), mean=",x$stats$phasestats[1],", 95% CI [",x$stats$phasestats[2],", ",x$stats$phasestats[3],"]\n",sep="", ...)
} # end of k==1
if(k>=2){
for (index in 1:k){
cat("Cycle=",x$cycles[index],"\n",sep="", ...)
cat("Season, mean=",x$stats$wstats[index,1],", 95% CI [",x$stats$wstats[index,2],", ",x$stats$wstats[index,3],"]\n",sep="", ...)
}
cat("\nPhase and amplitude\n")
for (index in 1:k){
cat("Cycle=",x$cycles[index],"\n",sep="", ...)
cat("Amplitude, mean=",x$stats$ampstats[index,1],", 95% CI [",x$stats$ampstats[index,2],", ",x$stats$ampstats[index,3],"]\n",sep="", ...)
cat("Phase (radians), mean=",x$stats$phasestats[index,1],", 95% CI [",x$stats$phasestats[index,2],", ",x$stats$phasestats[index,3],"]\n",sep="", ...)
}
} # end of k>=2
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.