Nothing
"samplesStats" <-
function(node, beg = samplesGetBeg(), end = samplesGetEnd(),
firstChain = samplesGetFirstChain(), lastChain = samplesGetLastChain(),
thin = samplesGetThin())
# Calculate statistics for monitored node
{
oldBeg <- samplesGetBeg()
oldEnd <- samplesGetEnd()
oldFirstChain <- samplesGetFirstChain()
oldLastChain <- samplesGetLastChain()
oldThin <- samplesGetThin()
on.exit({
samplesSetBeg(oldBeg)
samplesSetEnd(oldEnd)
samplesSetFirstChain(oldFirstChain)
samplesSetLastChain(oldLastChain)
samplesSetThin(oldThin)
})
samplesSetBeg(beg)
samplesSetEnd(end)
samplesSetFirstChain(firstChain)
samplesSetLastChain(lastChain)
thin <- max(c(thin, 1))
samplesSetThin(thin)
if (is.R()){
result <- data.frame(mean=NULL, sd=NULL, MC_error = NULL, val2.5pc=NULL,
median=NULL, val97.5pc=NULL, start = NULL, sample=NULL)
} else {
result <- data.frame(mean=numeric(), sd=numeric(), MC.error = numeric(),
val2.5pc=numeric(), median=numeric(), val97.5pc=numeric(),
start = numeric(), sample=numeric())
}
for(i in seq(along=node)){
command <- paste(.SamplesGlobalsCmd(node[i]), "SamplesEmbed.StatsGuard;SamplesEmbed.Stats")
.CmdInterpreter(command)
buffer <- file.path(tempdir(), "buffer.txt")
rlb <- readLines(buffer)
len <- length(rlb)
if (len > 1)
result <- rbind(result, read.table(buffer))
else{
if(length(grep("val97.5pc", rlb)))
message("Variable ", node[i], " has probably not been updated")
else
message("Variable ", node[i], ": ", rlb)
}
}
return(result)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.