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
|
### R code from vignette source 'DiscreteBayes.Rnw'
###################################################
### code chunk number 1: DiscreteBayes.Rnw:21-23
###################################################
p <- seq(0, 1, by = 0.01)
prior <- 1 / 101 + 0 * p
###################################################
### code chunk number 2: DiscreteBayes.Rnw:25-28
###################################################
plot(p, prior,
type="h",
main="Prior Distribution")
###################################################
### code chunk number 3: DiscreteBayes.Rnw:35-37
###################################################
library(LearnBayes)
post <- pdisc(p, prior, c(20, 12))
###################################################
### code chunk number 4: DiscreteBayes.Rnw:39-42
###################################################
plot(p, post,
type="h",
main="Posterior Distribution")
###################################################
### code chunk number 5: DiscreteBayes.Rnw:47-48
###################################################
discint(cbind(p, post), 0.90)
###################################################
### code chunk number 6: DiscreteBayes.Rnw:57-60
###################################################
n <- 20
s <- 0:20
pred.probs <- pdiscp(p, post, n, s)
###################################################
### code chunk number 7: DiscreteBayes.Rnw:63-66
###################################################
plot(s, pred.probs,
type="h",
main="Predictive Distribution")
###################################################
### code chunk number 8: DiscreteBayes.Rnw:72-74
###################################################
prior <- rep(1/11, 11)
names(prior) <- 20:30
###################################################
### code chunk number 9: DiscreteBayes.Rnw:78-79
###################################################
y <- c(24, 25, 31, 31, 22, 21, 26, 20, 16, 22)
###################################################
### code chunk number 10: DiscreteBayes.Rnw:83-84
###################################################
post <- discrete.bayes(dpois, prior, y)
###################################################
### code chunk number 11: DiscreteBayes.Rnw:89-90
###################################################
print(post)
###################################################
### code chunk number 12: DiscreteBayes.Rnw:93-94
###################################################
plot(post)
###################################################
### code chunk number 13: DiscreteBayes.Rnw:97-98
###################################################
summary(post)
|
