getGibbsSummary {phenoCDM} | R Documentation |
Summarize Output of the CDM Model
Description
This function return a summary of the output from the Gibbs-Sampling of the CDM model.
Usage
getGibbsSummary(ssOut, burnin = NULL, colNames = NULL,
sigmaPerSeason = TRUE)
Arguments
ssOut |
CDM output list. |
burnin |
Number of burnin itterations . |
colNames |
vector of charachters includes names of each variable in the output. |
sigmaPerSeason |
logical value indicating whether each site/season has a separate process error |
Examples
#Summarize CDM Model Ouput
ssSim <- phenoSim(nSites = 2, #number of sites
nTSet = 30, #number of Time steps
beta = c(1, 2), #beta coefficients
sig = .01, #process error
tau = .1, #observation error
plotFlag = TRUE, #whether plot the data or not
miss = 0.05, #fraction of missing data
ymax = c(6, 3) #maximum of saturation trajectory
)
ssOut <- fitCDM(x = ssSim$x, #predictors
nGibbs = 200,
nBurnin = 100,
z = ssSim$z,#response
connect = ssSim$connect, #connectivity of time data
quiet=TRUE)
summ <- getGibbsSummary(ssOut, burnin = 100, sigmaPerSeason = FALSE)
colMeans(summ$ymax)
colMeans(summ$betas)
colMeans(summ$tau)
colMeans(summ$sigma)
[Package phenoCDM version 0.1.3 Index]