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