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
82
83
84
85
86
87
88
89
90
91
92
93
94
|
"chanames" <-
function (x, allow.null = TRUE)
{
if (is.mcmc.list(x)) {
if (is.null(names(x)))
if (allow.null)
NULL
else paste("chain", 1:length(x), sep = "")
else names(x)
}
else NULL
}
"chanames<-" <-
function (x, value)
{
if (is.mcmc.list(x))
names(x) <- value
else stop("Not an mcmc.list object")
x
}
"varnames" <-
function (x, allow.null = TRUE)
{
if (!is.mcmc(x) && !is.mcmc.list(x))
return(NULL)
y <- if (is.mcmc(x))
dimnames(x)[[2]]
else if (is.mcmc.list(x))
dimnames(x[[1]])[[2]]
if (is.null(y) && !allow.null)
y <- paste("var", 1:nvar(x), sep = "")
return(y)
}
"varnames<-" <-
function (x, value)
{
if (is.mcmc(x)) {
dimnames(x)[[2]] <- value
}
else if (is.mcmc.list(x))
for (i in 1:nchain(x)) varnames(x[[i]]) <- value
else stop("Not an mcmc or mcmc.list object")
x
}
"nchain" <-
function (x)
{
if (is.mcmc(x))
1
else if (is.mcmc.list(x))
length(x)
else NULL
}
"nvar" <-
function (x)
{
if (is.mcmc(x)) {
if (is.matrix(x)) ncol(x) else 1
}
else if (is.mcmc.list(x)) {
if (is.matrix(x[[1]])) ncol(x[[1]]) else 1
}
else NULL
}
"niter" <-
function (x)
{
if (is.mcmc(x)) {
if (is.matrix(x)) nrow(x) else length(x)
}
else if (is.mcmc.list(x)) {
if (is.matrix(x[[1]])) nrow(x[[1]]) else length(x[[1]])
}
else NULL
}
|
