Displaying 1 result from an estimated 1 matches for "numknots".
Did you mean:
nknots
2013 Jan 18
0
problem that arises after using the new version of "BRugs"
...}
samplesSet(parametersToSave)
iterIncSize <- round(nIter/(nThin * 100))
for (i in 1:100) {
modelUpdate(iterIncSize, thin = nThin)
}
sims <- samplesStats("*")
return(list(Stats = sims))
}
require(BRugs)
numKnots <- 25
range.x <- c(0,202)
ox <- order(x)
x <- x[ox]
delta <- delta[ox]
n <- length(x)
xTil <- as.numeric(names(table(x)))
cnts <- as.numeric(table(x))
nTil <- length(xTil)
deltaTil <- rep(NA, nTil)
for (iTil in 1:nTil)...