Work at SourceForge, help us to make it a better place! We have an immediate need for a Support Technician in our San Francisco or Denver office.

Close

[87be2a]: R / jags.object.R Maximize Restore History

Download this file

jags.object.R    82 lines (71 with data), 2.2 kB

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
update.jags <- function(object, n.iter = 1, by, progress.bar, ...)
{
if (!is.numeric(n.iter) || n.iter < 1) {
stop("Invalid n.iter")
}
adapting <- .Call("is_adapting", object$ptr(), PACKAGE="rjags")
on.exit(object$sync())
if (missing(progress.bar)) {
progress.bar <- getOption("jags.pb")
}
if (!is.null(progress.bar)) {
match.arg(progress.bar, c("text","gui","none"))
if (progress.bar=="none")
progress.bar <- NULL
}
do.pb <- interactive() && !is.null(progress.bar) && n.iter >= 100
if (do.pb) {
start.iter <- object$iter()
end.iter <- start.iter + n.iter
pb <- switch(progress.bar,
"text" = txtProgressBar(start.iter, end.iter,
initial = start.iter, style=3, width=50,
char=ifelse(adapting,"+","*")),
"gui" = updatePB(start.iter, end.iter, adapting))
}
## Set refresh frequency for progress bar
if (missing(by) || by <= 0) {
by <- min(ceiling(n.iter/50), 100)
}
else {
by <- ceiling(by)
}
## Do updates
n <- n.iter
while (n > 0) {
.Call("update", object$ptr(), min(n,by), PACKAGE="rjags")
if (do.pb) {
switch(progress.bar,
"text" = setTxtProgressBar(pb, object$iter()),
"gui" = setPB(pb, object$iter()))
}
n <- n - by
}
if (do.pb) {
close(pb)
}
invisible(NULL)
}
adapt <- function(object, n.iter, end.adaptation = FALSE, ...)
{
if(.Call("is_adapting", object$ptr(), PACKAGE="rjags")) {
update(object, n.iter, ...)
ok <- .Call("check_adaptation", object$ptr(), PACKAGE="rjags")
if (end.adaptation) {
.Call("adapt_off", object$ptr(), PACKAGE="rjags")
}
return(ok)
}
else {
return(TRUE)
}
}
coef.jags <- function(object, chain = 1, ...) {
if (!is.numeric(chain) || chain < 1 || chain > object$nchain()) {
stop("Invalid chain")
}
object$state(internal=FALSE)[[chain]]
}
variable.names.jags <- function(object, ...) {
.Call("get_variable_names", object$ptr(), PACKAGE="rjags")
}