1### R code from vignette source 'examples.rnw'
2
3###################################################
4### code chunk number 1: preliminaries
5###################################################
6library(MSwM)
7set.seed(2341)
8
9
10###################################################
11### code chunk number 2: examples.rnw:50-51
12###################################################
13data(example)
14
15
16###################################################
17### code chunk number 3: examples.rnw:55-56
18###################################################
19plot(ts(example))
20
21
22###################################################
23### code chunk number 4: examples.rnw:63-65
24###################################################
25mod=lm(y~x,example)
26summary(mod)
27
28
29###################################################
30### code chunk number 5: fig1_3
31###################################################
32qqnorm(resid(mod))
33qqline(resid(mod))
34
35
36###################################################
37### code chunk number 6: fig1_4
38###################################################
39acf(resid(mod))
40
41
42###################################################
43### code chunk number 7: examples.rnw:88-90
44###################################################
45mod.mswm=msmFit(mod,k=2,p=1,sw=c(TRUE,TRUE,TRUE,TRUE),control=list(parallel=FALSE))
46summary(mod.mswm)
47
48
49###################################################
50### code chunk number 8: examples.rnw:97-98
51###################################################
52plot(mod.mswm)
53
54
55###################################################
56### code chunk number 9: fig1_6
57###################################################
58plotDiag(mod.mswm,which=2)
59
60
61###################################################
62### code chunk number 10: fig1_7
63###################################################
64plotDiag(mod.mswm,which=3)
65
66
67###################################################
68### code chunk number 11: pp1m3
69###################################################
70   plotProb(mod.mswm,which=1)
71
72
73###################################################
74### code chunk number 12: pp2m3
75###################################################
76   plotProb(mod.mswm,which=2)
77
78
79###################################################
80### code chunk number 13: pp3m3
81###################################################
82   plotReg(mod.mswm,expl="x")
83
84
85###################################################
86### code chunk number 14: examples.rnw:156-158
87###################################################
88data(traffic)
89plot(ts(traffic[2:4]),main="Traffic")
90
91
92###################################################
93### code chunk number 15: examples.rnw:165-167
94###################################################
95model=glm(NDead~Temp+Prec,traffic,family="poisson")
96summary(model)
97
98
99###################################################
100### code chunk number 16: examples.rnw:172-174
101###################################################
102m1=msmFit(model,k=2,sw=c(TRUE,TRUE,TRUE),family="poisson",control=list(parallel=FALSE))
103summary(m1)
104
105
106###################################################
107### code chunk number 17: examples.rnw:177-178
108###################################################
109intervals(m1)
110
111
112###################################################
113### code chunk number 18: examples.rnw:182-183
114###################################################
115plot(m1)
116
117
118###################################################
119### code chunk number 19: fig2_3
120###################################################
121plotDiag(m1,which=2)
122
123
124###################################################
125### code chunk number 20: fig2_4
126###################################################
127plotDiag(m1,which=3)
128
129
130###################################################
131### code chunk number 21: pp1m1
132###################################################
133   plotProb(m1,which=2)
134
135
136