SimInf practice
Alfredo Acosta
2023-06-22
1 SIR practice model
# Parameters
beta <- 4.0706 # Calculated nls montly
gama <- 2.7941 # Calculated nls montly
#Number of animals in the model
n <- 1000
# Initial state of the population compartments
u0 <- data.frame(S=rep(999,n), I=rep(1,n), R=rep(0,n))
# time
tspan <- seq(1,180,1)1.1 Run the model
set.seed(555)
model <- SIR(u0=u0, tspan = tspan, beta = beta, gamma = gama)
result <- run(model=model, threads = 1)
result## Model: SIR
## Number of nodes: 1000
## Number of transitions: 2
## Number of scheduled events: 0
##
## Local data
## ----------
## Parameter Value
## beta 4.0706
## gamma 2.7941
##
## Compartments
## ------------
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## S 302.00 491.00 997.00 825.53 999.00 999.00
## I 0.00 0.00 0.00 0.35 0.00 119.00
## R 0.00 1.00 3.00 174.12 503.00 698.001.2 Plot results
plot(result)
plot(result, range = FALSE, cex=0.1)
plot(prevalence(result, I~S+I+R), cex=0.1, col="red")
## Running on differente time scale (days)
# Parameters
beta <- 4.0706/30 # Calculated nls montly
gama <- 2.7941/30 # Calculated nls montly
#Number of animals in the model
n <- 1000
# Initial state of the population compartments
u0 <- data.frame(S=rep(999,n), I=rep(1,n), R=rep(0,n))
# time
tspan <- seq(1,180*30,30)1.3 Running on days
set.seed(555)
model <- SIR(u0=u0, tspan = tspan, beta = beta, gamma = gama)
result <- run(model=model, threads = 1)
result## Model: SIR
## Number of nodes: 1000
## Number of transitions: 2
## Number of scheduled events: 0
##
## Local data
## ----------
## Parameter Value
## beta 0.13568667
## gamma 0.09313667
##
## Compartments
## ------------
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## S 317.000 497.000 998.000 830.207 999.000 999.000
## I 0.000 0.000 0.000 0.347 0.000 104.000
## R 0.000 1.000 2.000 169.445 503.000 683.0001.4 Ploting on days
plot(result)
plot(result, range = FALSE, cex=0.1)
plot(prevalence(result, I~S+I+R), cex=0.1, col="red")
2 Second model SIR
# Example 2 Giving te equations (transitions)
model2 <- mparse(transitions = c("S -> beta*S*I/(S+I+R) -> I", "I -> gamma*I -> R"),
compartments = c("S", "I", "R"),
gdata = c(beta = 0.16, gamma = 0.077),
u0 = data.frame(S = 100, I = 1, R = 0),
tspan = 1:1000)2.1 Run and plot
set.seed(123)
result2 <- run(model=model2, threads = 1)
result2## Model: SimInf_model
## Number of nodes: 1
## Number of transitions: 2
## Number of scheduled events: 0
##
## Global data
## -----------
## Parameter Value
## beta 0.160
## gamma 0.077
##
## Compartments
## ------------
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## S 13.00 13.00 13.00 18.06 13.00 100.00
## I 0.00 0.00 0.00 1.24 0.00 21.00
## R 0.00 88.00 88.00 81.70 88.00 88.00plot(result2)
plot(result2, node = 1:10, range = FALSE)## Warning in plot.window(...): "node" is not a graphical parameter## Warning in plot.xy(xy, type, ...): "node" is not a graphical parameter## Warning in axis(side = side, at = at, labels = labels, ...): "node" is not a
## graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "node" is not a
## graphical parameter## Warning in title(...): "node" is not a graphical parameter## Warning in plot.xy(xy.coords(x, y), type = type, ...): "node" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "node" is not a
## graphical parameter
plot(prevalence(result2, I~.), cex=0.3, col="red")
3 Third model SEIR
# Example 2 Giving te equations (transitions)
model3 <- mparse(transitions = c("S -> alpha*S -> E",
"E -> tau*E -> I",
"I -> beta*S*I/(S+I+R) -> R",
"R -> gamma*R -> S"),
compartments = c("S", "E", "I", "R"),
gdata = c(alpha = 0.05, tau= 0.002, beta = 0.16, gamma = 0.077),
u0 = data.frame(S = 100, E= 20, I = 5, R = 0),
tspan = 1:1000)
result3 <- run(model=model3, threads = 1)
result3## Model: SimInf_model
## Number of nodes: 1
## Number of transitions: 4
## Number of scheduled events: 0
##
## Global data
## -----------
## Parameter Value
## alpha 0.050
## tau 0.002
## beta 0.160
## gamma 0.077
##
## Compartments
## ------------
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## S 0.00 0.00 2.00 4.83 5.00 96.00
## E 24.00 95.00 109.00 100.99 115.00 121.00
## I 0.00 3.00 7.00 17.33 24.00 71.00
## R 0.00 0.00 2.00 1.85 3.00 8.00plot(result3)
plot(prevalence(result3, I~.), cex=0.3, col="red")
t <- data.frame(trajectory(result3))
t$N <- t$S + t$E + t$I + t$R
t## node time S E I R N
## 1 1 1 96 24 5 0 125
## 2 1 2 93 27 3 2 125
## 3 1 3 91 29 3 2 125
## 4 1 4 89 31 3 2 125
## 5 1 5 85 35 3 2 125
## 6 1 6 85 36 3 1 125
## 7 1 7 79 41 4 1 125
## 8 1 8 70 50 4 1 125
## 9 1 9 67 53 3 2 125
## 10 1 10 65 56 2 2 125
## 11 1 11 60 60 2 3 125
## 12 1 12 56 63 3 3 125
## 13 1 13 55 64 3 3 125
## 14 1 14 51 68 2 4 125
## 15 1 15 49 71 2 3 125
## 16 1 16 46 74 2 3 125
## 17 1 17 43 78 2 2 125
## 18 1 18 42 79 2 2 125
## 19 1 19 41 80 2 2 125
## 20 1 20 42 80 2 1 125
## 21 1 21 41 81 2 1 125
## 22 1 22 41 81 2 1 125
## 23 1 23 39 84 1 1 125
## 24 1 24 38 85 1 1 125
## 25 1 25 35 88 1 1 125
## 26 1 26 34 88 1 2 125
## 27 1 27 32 90 1 2 125
## 28 1 28 31 91 1 2 125
## 29 1 29 30 92 1 2 125
## 30 1 30 28 94 1 2 125
## 31 1 31 25 97 1 2 125
## 32 1 32 23 98 2 2 125
## 33 1 33 22 99 2 2 125
## 34 1 34 23 98 3 1 125
## 35 1 35 24 98 3 0 125
## 36 1 36 24 98 3 0 125
## 37 1 37 24 98 3 0 125
## 38 1 38 22 100 1 2 125
## 39 1 39 22 100 1 2 125
## 40 1 40 20 102 1 2 125
## 41 1 41 18 103 0 4 125
## 42 1 42 16 105 0 4 125
## 43 1 43 16 106 0 3 125
## 44 1 44 15 107 0 3 125
## 45 1 45 15 107 0 3 125
## 46 1 46 14 108 0 3 125
## 47 1 47 14 108 0 3 125
## 48 1 48 13 109 0 3 125
## 49 1 49 14 109 0 2 125
## 50 1 50 13 109 1 2 125
## 51 1 51 11 110 2 2 125
## 52 1 52 11 110 2 2 125
## 53 1 53 11 111 2 1 125
## 54 1 54 10 112 2 1 125
## 55 1 55 7 114 3 1 125
## 56 1 56 5 116 3 1 125
## 57 1 57 4 118 3 0 125
## 58 1 58 4 118 3 0 125
## 59 1 59 4 118 3 0 125
## 60 1 60 3 119 3 0 125
## 61 1 61 3 119 3 0 125
## 62 1 62 3 118 4 0 125
## 63 1 63 3 118 3 1 125
## 64 1 64 3 118 3 1 125
## 65 1 65 3 118 3 1 125
## 66 1 66 3 118 3 1 125
## 67 1 67 3 117 4 1 125
## 68 1 68 3 117 3 2 125
## 69 1 69 3 117 3 2 125
## 70 1 70 3 117 3 2 125
## 71 1 71 3 117 3 2 125
## 72 1 72 3 118 3 1 125
## 73 1 73 3 118 3 1 125
## 74 1 74 2 119 2 2 125
## 75 1 75 2 119 2 2 125
## 76 1 76 1 120 2 2 125
## 77 1 77 2 120 2 1 125
## 78 1 78 2 119 3 1 125
## 79 1 79 2 119 3 1 125
## 80 1 80 2 119 3 1 125
## 81 1 81 2 119 3 1 125
## 82 1 82 2 119 3 1 125
## 83 1 83 2 119 3 1 125
## 84 1 84 2 119 2 2 125
## 85 1 85 2 119 2 2 125
## 86 1 86 2 119 2 2 125
## 87 1 87 2 117 4 2 125
## 88 1 88 2 117 4 2 125
## 89 1 89 2 117 4 2 125
## 90 1 90 2 117 4 2 125
## 91 1 91 2 117 4 2 125
## 92 1 92 2 117 4 2 125
## 93 1 93 2 116 5 2 125
## 94 1 94 2 116 5 2 125
## 95 1 95 2 116 5 2 125
## 96 1 96 2 116 5 2 125
## 97 1 97 2 116 5 2 125
## 98 1 98 2 116 5 2 125
## 99 1 99 2 115 6 2 125
## 100 1 100 2 115 6 2 125
## 101 1 101 1 116 5 3 125
## 102 1 102 1 116 4 4 125
## 103 1 103 3 116 2 4 125
## 104 1 104 3 117 2 3 125
## 105 1 105 3 117 0 5 125
## 106 1 106 3 117 0 5 125
## 107 1 107 3 117 0 5 125
## 108 1 108 3 117 0 5 125
## 109 1 109 3 117 0 5 125
## 110 1 110 3 117 0 5 125
## 111 1 111 3 116 1 5 125
## 112 1 112 3 114 3 5 125
## 113 1 113 2 115 3 5 125
## 114 1 114 2 115 3 5 125
## 115 1 115 2 114 4 5 125
## 116 1 116 2 114 4 5 125
## 117 1 117 2 114 4 5 125
## 118 1 118 2 114 4 5 125
## 119 1 119 2 114 4 5 125
## 120 1 120 2 114 4 5 125
## 121 1 121 2 114 4 5 125
## 122 1 122 2 114 4 5 125
## 123 1 123 2 114 4 5 125
## 124 1 124 3 114 4 4 125
## 125 1 125 3 114 4 4 125
## 126 1 126 2 114 5 4 125
## 127 1 127 2 113 6 4 125
## 128 1 128 2 113 6 4 125
## 129 1 129 3 113 5 4 125
## 130 1 130 4 113 3 5 125
## 131 1 131 4 113 2 6 125
## 132 1 132 6 112 3 4 125
## 133 1 133 5 113 3 4 125
## 134 1 134 3 115 3 4 125
## 135 1 135 5 115 3 2 125
## 136 1 136 5 115 3 2 125
## 137 1 137 5 115 3 2 125
## 138 1 138 6 115 3 1 125
## 139 1 139 6 114 3 2 125
## 140 1 140 6 114 3 2 125
## 141 1 141 6 114 2 3 125
## 142 1 142 5 116 2 2 125
## 143 1 143 5 116 2 2 125
## 144 1 144 5 115 3 2 125
## 145 1 145 5 115 3 2 125
## 146 1 146 5 115 3 2 125
## 147 1 147 4 115 4 2 125
## 148 1 148 4 115 4 2 125
## 149 1 149 3 116 4 2 125
## 150 1 150 3 116 3 3 125
## 151 1 151 3 116 3 3 125
## 152 1 152 2 117 3 3 125
## 153 1 153 4 117 3 1 125
## 154 1 154 4 116 4 1 125
## 155 1 155 4 116 4 1 125
## 156 1 156 4 115 5 1 125
## 157 1 157 5 114 6 0 125
## 158 1 158 5 114 5 1 125
## 159 1 159 5 113 6 1 125
## 160 1 160 5 113 6 1 125
## 161 1 161 5 113 5 2 125
## 162 1 162 5 113 4 3 125
## 163 1 163 5 110 7 3 125
## 164 1 164 4 110 8 3 125
## 165 1 165 3 111 8 3 125
## 166 1 166 3 111 8 3 125
## 167 1 167 3 110 9 3 125
## 168 1 168 3 110 9 3 125
## 169 1 169 3 110 7 5 125
## 170 1 170 3 110 6 6 125
## 171 1 171 4 110 5 6 125
## 172 1 172 4 110 5 6 125
## 173 1 173 4 110 5 6 125
## 174 1 174 5 110 5 5 125
## 175 1 175 7 110 5 3 125
## 176 1 176 8 110 5 2 125
## 177 1 177 8 110 3 4 125
## 178 1 178 7 111 3 4 125
## 179 1 179 8 111 1 5 125
## 180 1 180 8 111 1 5 125
## 181 1 181 7 112 1 5 125
## 182 1 182 6 113 1 5 125
## 183 1 183 6 113 1 5 125
## 184 1 184 5 113 2 5 125
## 185 1 185 6 113 2 4 125
## 186 1 186 7 112 3 3 125
## 187 1 187 7 112 3 3 125
## 188 1 188 6 112 4 3 125
## 189 1 189 6 111 5 3 125
## 190 1 190 5 112 4 4 125
## 191 1 191 5 111 5 4 125
## 192 1 192 5 111 5 4 125
## 193 1 193 4 111 6 4 125
## 194 1 194 5 111 6 3 125
## 195 1 195 5 111 6 3 125
## 196 1 196 6 111 6 2 125
## 197 1 197 6 111 6 2 125
## 198 1 198 6 110 6 3 125
## 199 1 199 5 111 6 3 125
## 200 1 200 6 110 6 3 125
## 201 1 201 6 110 6 3 125
## 202 1 202 5 111 5 4 125
## 203 1 203 5 110 6 4 125
## 204 1 204 5 110 6 4 125
## 205 1 205 5 110 6 4 125
## 206 1 206 5 110 4 6 125
## 207 1 207 6 109 5 5 125
## 208 1 208 6 109 4 6 125
## 209 1 209 7 109 4 5 125
## 210 1 210 7 109 4 5 125
## 211 1 211 7 109 3 6 125
## 212 1 212 6 110 2 7 125
## 213 1 213 7 109 3 6 125
## 214 1 214 6 110 3 6 125
## 215 1 215 5 111 3 6 125
## 216 1 216 4 112 3 6 125
## 217 1 217 3 113 3 6 125
## 218 1 218 3 113 3 6 125
## 219 1 219 4 113 3 5 125
## 220 1 220 5 112 4 4 125
## 221 1 221 4 111 6 4 125
## 222 1 222 4 111 6 4 125
## 223 1 223 4 111 6 4 125
## 224 1 224 4 111 6 4 125
## 225 1 225 4 111 6 4 125
## 226 1 226 5 110 7 3 125
## 227 1 227 5 110 7 3 125
## 228 1 228 5 111 6 3 125
## 229 1 229 5 111 5 4 125
## 230 1 230 6 111 5 3 125
## 231 1 231 6 111 5 3 125
## 232 1 232 7 111 5 2 125
## 233 1 233 7 111 4 3 125
## 234 1 234 7 112 4 2 125
## 235 1 235 7 112 3 3 125
## 236 1 236 7 112 2 4 125
## 237 1 237 8 112 2 3 125
## 238 1 238 9 112 3 1 125
## 239 1 239 8 113 3 1 125
## 240 1 240 7 114 3 1 125
## 241 1 241 7 114 3 1 125
## 242 1 242 6 115 1 3 125
## 243 1 243 6 115 1 3 125
## 244 1 244 6 115 1 3 125
## 245 1 245 6 115 1 3 125
## 246 1 246 7 115 1 2 125
## 247 1 247 7 115 1 2 125
## 248 1 248 7 115 1 2 125
## 249 1 249 7 115 1 2 125
## 250 1 250 6 117 1 1 125
## 251 1 251 5 118 1 1 125
## 252 1 252 5 118 1 1 125
## 253 1 253 5 118 1 1 125
## 254 1 254 5 118 1 1 125
## 255 1 255 5 118 1 1 125
## 256 1 256 5 118 1 1 125
## 257 1 257 5 119 0 1 125
## 258 1 258 5 118 1 1 125
## 259 1 259 5 118 1 1 125
## 260 1 260 5 118 1 1 125
## 261 1 261 5 117 1 2 125
## 262 1 262 5 117 0 3 125
## 263 1 263 5 117 0 3 125
## 264 1 264 5 117 0 3 125
## 265 1 265 4 118 0 3 125
## 266 1 266 4 118 0 3 125
## 267 1 267 3 119 0 3 125
## 268 1 268 3 119 0 3 125
## 269 1 269 2 120 0 3 125
## 270 1 270 1 121 0 3 125
## 271 1 271 1 121 0 3 125
## 272 1 272 1 121 0 3 125
## 273 1 273 1 121 0 3 125
## 274 1 274 2 121 0 2 125
## 275 1 275 2 121 0 2 125
## 276 1 276 2 120 1 2 125
## 277 1 277 2 120 1 2 125
## 278 1 278 2 120 1 2 125
## 279 1 279 2 120 1 2 125
## 280 1 280 2 120 1 2 125
## 281 1 281 2 120 1 2 125
## 282 1 282 3 120 1 1 125
## 283 1 283 3 120 1 1 125
## 284 1 284 2 120 2 1 125
## 285 1 285 2 120 2 1 125
## 286 1 286 2 119 3 1 125
## 287 1 287 2 119 3 1 125
## 288 1 288 1 120 3 1 125
## 289 1 289 1 120 3 1 125
## 290 1 290 1 120 3 1 125
## 291 1 291 1 120 3 1 125
## 292 1 292 1 120 3 1 125
## 293 1 293 1 120 3 1 125
## 294 1 294 1 120 3 1 125
## 295 1 295 2 120 3 0 125
## 296 1 296 2 120 3 0 125
## 297 1 297 2 120 2 1 125
## 298 1 298 1 121 2 1 125
## 299 1 299 1 121 2 1 125
## 300 1 300 2 121 2 0 125
## 301 1 301 2 121 1 1 125
## 302 1 302 1 121 2 1 125
## 303 1 303 1 121 2 1 125
## 304 1 304 1 121 2 1 125
## 305 1 305 1 120 3 1 125
## 306 1 306 1 120 3 1 125
## 307 1 307 1 120 3 1 125
## 308 1 308 1 120 3 1 125
## 309 1 309 1 120 3 1 125
## 310 1 310 1 119 4 1 125
## 311 1 311 1 119 4 1 125
## 312 1 312 1 119 4 1 125
## 313 1 313 1 118 5 1 125
## 314 1 314 1 118 5 1 125
## 315 1 315 1 118 5 1 125
## 316 1 316 1 118 5 1 125
## 317 1 317 1 117 6 1 125
## 318 1 318 1 117 6 1 125
## 319 1 319 1 116 7 1 125
## 320 1 320 1 115 8 1 125
## 321 1 321 1 115 8 1 125
## 322 1 322 1 115 8 1 125
## 323 1 323 1 115 8 1 125
## 324 1 324 1 114 9 1 125
## 325 1 325 1 112 11 1 125
## 326 1 326 1 113 11 0 125
## 327 1 327 1 112 12 0 125
## 328 1 328 1 111 13 0 125
## 329 1 329 1 111 13 0 125
## 330 1 330 1 111 13 0 125
## 331 1 331 1 111 13 0 125
## 332 1 332 1 111 13 0 125
## 333 1 333 1 111 13 0 125
## 334 1 334 1 111 13 0 125
## 335 1 335 1 111 13 0 125
## 336 1 336 1 111 12 1 125
## 337 1 337 2 111 12 0 125
## 338 1 338 2 111 12 0 125
## 339 1 339 2 111 12 0 125
## 340 1 340 2 110 13 0 125
## 341 1 341 2 109 13 1 125
## 342 1 342 2 109 13 1 125
## 343 1 343 2 109 13 1 125
## 344 1 344 2 109 13 1 125
## 345 1 345 2 109 13 1 125
## 346 1 346 2 109 13 1 125
## 347 1 347 2 109 13 1 125
## 348 1 348 3 108 14 0 125
## 349 1 349 3 108 13 1 125
## 350 1 350 3 108 12 2 125
## 351 1 351 3 108 12 2 125
## 352 1 352 2 108 13 2 125
## 353 1 353 2 108 13 2 125
## 354 1 354 2 108 13 2 125
## 355 1 355 2 107 14 2 125
## 356 1 356 2 107 14 2 125
## 357 1 357 2 107 14 2 125
## 358 1 358 2 107 14 2 125
## 359 1 359 3 107 13 2 125
## 360 1 360 3 107 13 2 125
## 361 1 361 3 107 13 2 125
## 362 1 362 3 107 13 2 125
## 363 1 363 3 107 12 3 125
## 364 1 364 3 107 12 3 125
## 365 1 365 3 107 12 3 125
## 366 1 366 3 107 12 3 125
## 367 1 367 3 107 12 3 125
## 368 1 368 4 107 12 2 125
## 369 1 369 4 106 12 3 125
## 370 1 370 4 106 12 3 125
## 371 1 371 4 107 10 4 125
## 372 1 372 4 106 11 4 125
## 373 1 373 4 105 11 5 125
## 374 1 374 4 105 10 6 125
## 375 1 375 4 105 10 6 125
## 376 1 376 4 105 10 6 125
## 377 1 377 4 104 11 6 125
## 378 1 378 4 104 9 8 125
## 379 1 379 6 103 10 6 125
## 380 1 380 6 103 10 6 125
## 381 1 381 6 103 9 7 125
## 382 1 382 6 103 8 8 125
## 383 1 383 7 102 9 7 125
## 384 1 384 8 101 9 7 125
## 385 1 385 8 101 9 7 125
## 386 1 386 10 101 8 6 125
## 387 1 387 9 102 8 6 125
## 388 1 388 9 102 7 7 125
## 389 1 389 8 103 7 7 125
## 390 1 390 8 103 7 7 125
## 391 1 391 8 103 7 7 125
## 392 1 392 8 104 7 6 125
## 393 1 393 8 104 7 6 125
## 394 1 394 7 104 8 6 125
## 395 1 395 8 104 9 4 125
## 396 1 396 7 105 9 4 125
## 397 1 397 8 104 10 3 125
## 398 1 398 8 104 9 4 125
## 399 1 399 7 105 9 4 125
## 400 1 400 7 105 9 4 125
## 401 1 401 7 105 9 4 125
## 402 1 402 6 105 10 4 125
## 403 1 403 6 105 10 4 125
## 404 1 404 7 104 8 6 125
## 405 1 405 7 104 8 6 125
## 406 1 406 6 105 8 6 125
## 407 1 407 6 106 8 5 125
## 408 1 408 7 106 7 5 125
## 409 1 409 7 105 8 5 125
## 410 1 410 7 106 8 4 125
## 411 1 411 8 106 7 4 125
## 412 1 412 8 106 7 4 125
## 413 1 413 8 106 6 5 125
## 414 1 414 8 106 6 5 125
## 415 1 415 9 106 5 5 125
## 416 1 416 8 107 5 5 125
## 417 1 417 9 107 4 5 125
## 418 1 418 10 105 5 5 125
## 419 1 419 10 105 5 5 125
## 420 1 420 10 105 4 6 125
## 421 1 421 10 106 4 5 125
## 422 1 422 9 107 4 5 125
## 423 1 423 8 107 5 5 125
## 424 1 424 9 107 4 5 125
## 425 1 425 8 107 5 5 125
## 426 1 426 9 108 4 4 125
## 427 1 427 8 109 2 6 125
## 428 1 428 11 109 2 3 125
## 429 1 429 11 109 2 3 125
## 430 1 430 11 109 2 3 125
## 431 1 431 11 109 2 3 125
## 432 1 432 11 109 1 4 125
## 433 1 433 12 109 1 3 125
## 434 1 434 11 110 1 3 125
## 435 1 435 11 110 1 3 125
## 436 1 436 13 109 2 1 125
## 437 1 437 13 109 2 1 125
## 438 1 438 14 109 2 0 125
## 439 1 439 13 110 2 0 125
## 440 1 440 12 111 1 1 125
## 441 1 441 11 112 1 1 125
## 442 1 442 12 112 1 0 125
## 443 1 443 11 113 1 0 125
## 444 1 444 10 113 2 0 125
## 445 1 445 10 113 1 1 125
## 446 1 446 10 113 0 2 125
## 447 1 447 11 113 0 1 125
## 448 1 448 9 115 0 1 125
## 449 1 449 9 115 0 1 125
## 450 1 450 8 115 1 1 125
## 451 1 451 8 115 1 1 125
## 452 1 452 7 116 0 2 125
## 453 1 453 7 116 0 2 125
## 454 1 454 7 114 2 2 125
## 455 1 455 7 114 2 2 125
## 456 1 456 7 115 1 2 125
## 457 1 457 7 115 1 2 125
## 458 1 458 7 115 1 2 125
## 459 1 459 7 115 1 2 125
## 460 1 460 7 115 1 2 125
## 461 1 461 8 115 1 1 125
## 462 1 462 8 115 1 1 125
## 463 1 463 7 116 0 2 125
## 464 1 464 7 116 0 2 125
## 465 1 465 7 116 0 2 125
## 466 1 466 6 116 1 2 125
## 467 1 467 6 116 1 2 125
## 468 1 468 5 117 1 2 125
## 469 1 469 4 118 1 2 125
## 470 1 470 4 117 2 2 125
## 471 1 471 4 117 1 3 125
## 472 1 472 4 117 1 3 125
## 473 1 473 4 117 1 3 125
## 474 1 474 4 117 1 3 125
## 475 1 475 5 116 2 2 125
## 476 1 476 6 116 2 1 125
## 477 1 477 6 116 2 1 125
## 478 1 478 5 116 3 1 125
## 479 1 479 6 116 2 1 125
## 480 1 480 5 117 2 1 125
## 481 1 481 5 116 3 1 125
## 482 1 482 5 115 4 1 125
## 483 1 483 6 115 4 0 125
## 484 1 484 6 114 3 2 125
## 485 1 485 6 114 3 2 125
## 486 1 486 6 114 3 2 125
## 487 1 487 6 114 3 2 125
## 488 1 488 6 115 3 1 125
## 489 1 489 6 115 3 1 125
## 490 1 490 6 114 4 1 125
## 491 1 491 6 114 4 1 125
## 492 1 492 6 114 4 1 125
## 493 1 493 5 115 4 1 125
## 494 1 494 5 115 4 1 125
## 495 1 495 4 116 4 1 125
## 496 1 496 4 116 4 1 125
## 497 1 497 5 116 3 1 125
## 498 1 498 6 116 3 0 125
## 499 1 499 6 115 4 0 125
## 500 1 500 5 116 4 0 125
## 501 1 501 5 116 4 0 125
## 502 1 502 5 116 4 0 125
## 503 1 503 5 116 4 0 125
## 504 1 504 5 116 4 0 125
## 505 1 505 4 117 4 0 125
## 506 1 506 4 116 5 0 125
## 507 1 507 4 116 4 1 125
## 508 1 508 3 117 3 2 125
## 509 1 509 3 116 4 2 125
## 510 1 510 4 116 3 2 125
## 511 1 511 4 115 4 2 125
## 512 1 512 4 115 4 2 125
## 513 1 513 4 114 5 2 125
## 514 1 514 4 113 6 2 125
## 515 1 515 4 113 6 2 125
## 516 1 516 4 113 5 3 125
## 517 1 517 3 114 5 3 125
## 518 1 518 3 113 5 4 125
## 519 1 519 3 113 5 4 125
## 520 1 520 3 113 5 4 125
## 521 1 521 4 113 5 3 125
## 522 1 522 4 113 5 3 125
## 523 1 523 3 113 6 3 125
## 524 1 524 2 114 6 3 125
## 525 1 525 2 114 6 3 125
## 526 1 526 2 114 6 3 125
## 527 1 527 2 114 6 3 125
## 528 1 528 2 114 6 3 125
## 529 1 529 2 114 6 3 125
## 530 1 530 2 114 6 3 125
## 531 1 531 2 114 5 4 125
## 532 1 532 2 113 6 4 125
## 533 1 533 3 112 7 3 125
## 534 1 534 3 112 7 3 125
## 535 1 535 4 112 7 2 125
## 536 1 536 4 111 8 2 125
## 537 1 537 4 111 8 2 125
## 538 1 538 4 111 7 3 125
## 539 1 539 3 111 8 3 125
## 540 1 540 3 111 7 4 125
## 541 1 541 3 111 7 4 125
## 542 1 542 3 111 7 4 125
## 543 1 543 3 111 6 5 125
## 544 1 544 3 111 6 5 125
## 545 1 545 5 111 6 3 125
## 546 1 546 5 110 6 4 125
## 547 1 547 4 111 6 4 125
## 548 1 548 5 111 6 3 125
## 549 1 549 5 111 6 3 125
## 550 1 550 5 111 6 3 125
## 551 1 551 4 112 6 3 125
## 552 1 552 4 112 6 3 125
## 553 1 553 4 111 7 3 125
## 554 1 554 4 111 7 3 125
## 555 1 555 4 111 7 3 125
## 556 1 556 3 112 7 3 125
## 557 1 557 3 111 8 3 125
## 558 1 558 2 112 8 3 125
## 559 1 559 2 111 9 3 125
## 560 1 560 2 110 10 3 125
## 561 1 561 3 110 10 2 125
## 562 1 562 2 110 10 3 125
## 563 1 563 2 109 11 3 125
## 564 1 564 2 108 12 3 125
## 565 1 565 2 108 11 4 125
## 566 1 566 3 108 11 3 125
## 567 1 567 3 108 11 3 125
## 568 1 568 3 107 11 4 125
## 569 1 569 3 107 10 5 125
## 570 1 570 3 107 10 5 125
## 571 1 571 4 107 9 5 125
## 572 1 572 4 108 7 6 125
## 573 1 573 6 108 7 4 125
## 574 1 574 6 108 7 4 125
## 575 1 575 5 109 7 4 125
## 576 1 576 4 109 7 5 125
## 577 1 577 4 109 7 5 125
## 578 1 578 5 109 7 4 125
## 579 1 579 4 110 7 4 125
## 580 1 580 4 110 7 4 125
## 581 1 581 4 110 7 4 125
## 582 1 582 6 110 6 3 125
## 583 1 583 5 111 7 2 125
## 584 1 584 5 111 6 3 125
## 585 1 585 5 111 6 3 125
## 586 1 586 5 110 7 3 125
## 587 1 587 6 109 8 2 125
## 588 1 588 6 109 7 3 125
## 589 1 589 6 109 5 5 125
## 590 1 590 5 110 5 5 125
## 591 1 591 5 109 6 5 125
## 592 1 592 6 108 7 4 125
## 593 1 593 6 107 8 4 125
## 594 1 594 7 106 8 4 125
## 595 1 595 8 106 7 4 125
## 596 1 596 6 107 8 4 125
## 597 1 597 7 107 8 3 125
## 598 1 598 6 108 8 3 125
## 599 1 599 5 109 8 3 125
## 600 1 600 3 111 8 3 125
## 601 1 601 3 111 8 3 125
## 602 1 602 3 111 8 3 125
## 603 1 603 4 111 8 2 125
## 604 1 604 5 111 8 1 125
## 605 1 605 5 111 7 2 125
## 606 1 606 5 111 7 2 125
## 607 1 607 6 109 9 1 125
## 608 1 608 6 109 7 3 125
## 609 1 609 7 109 7 2 125
## 610 1 610 5 111 7 2 125
## 611 1 611 6 111 7 1 125
## 612 1 612 6 111 6 2 125
## 613 1 613 6 111 5 3 125
## 614 1 614 7 111 5 2 125
## 615 1 615 7 111 5 2 125
## 616 1 616 7 109 7 2 125
## 617 1 617 6 110 6 3 125
## 618 1 618 6 109 5 5 125
## 619 1 619 6 109 4 6 125
## 620 1 620 6 109 4 6 125
## 621 1 621 7 109 4 5 125
## 622 1 622 7 108 5 5 125
## 623 1 623 7 107 5 6 125
## 624 1 624 6 108 5 6 125
## 625 1 625 6 109 4 6 125
## 626 1 626 6 109 4 6 125
## 627 1 627 6 109 4 6 125
## 628 1 628 7 109 4 5 125
## 629 1 629 4 112 4 5 125
## 630 1 630 5 112 4 4 125
## 631 1 631 5 111 5 4 125
## 632 1 632 4 112 5 4 125
## 633 1 633 4 112 4 5 125
## 634 1 634 3 113 4 5 125
## 635 1 635 3 112 5 5 125
## 636 1 636 1 114 5 5 125
## 637 1 637 4 114 5 2 125
## 638 1 638 4 114 5 2 125
## 639 1 639 4 114 5 2 125
## 640 1 640 4 114 4 3 125
## 641 1 641 3 115 4 3 125
## 642 1 642 2 116 4 3 125
## 643 1 643 2 116 4 3 125
## 644 1 644 1 117 4 3 125
## 645 1 645 1 117 4 3 125
## 646 1 646 1 116 5 3 125
## 647 1 647 1 116 5 3 125
## 648 1 648 1 116 5 3 125
## 649 1 649 1 116 5 3 125
## 650 1 650 1 117 5 2 125
## 651 1 651 1 117 5 2 125
## 652 1 652 2 117 5 1 125
## 653 1 653 2 116 6 1 125
## 654 1 654 2 116 6 1 125
## 655 1 655 3 116 6 0 125
## 656 1 656 2 117 6 0 125
## 657 1 657 2 117 6 0 125
## 658 1 658 2 116 7 0 125
## 659 1 659 2 116 7 0 125
## 660 1 660 2 116 7 0 125
## 661 1 661 2 116 7 0 125
## 662 1 662 2 116 7 0 125
## 663 1 663 2 115 7 1 125
## 664 1 664 2 114 8 1 125
## 665 1 665 2 114 8 1 125
## 666 1 666 1 115 8 1 125
## 667 1 667 1 115 7 2 125
## 668 1 668 1 115 7 2 125
## 669 1 669 1 115 7 2 125
## 670 1 670 1 115 7 2 125
## 671 1 671 1 115 7 2 125
## 672 1 672 1 115 7 2 125
## 673 1 673 1 115 7 2 125
## 674 1 674 1 114 8 2 125
## 675 1 675 0 115 8 2 125
## 676 1 676 0 115 8 2 125
## 677 1 677 0 115 8 2 125
## 678 1 678 1 115 8 1 125
## 679 1 679 1 115 8 1 125
## 680 1 680 1 115 8 1 125
## 681 1 681 1 114 9 1 125
## 682 1 682 1 114 8 2 125
## 683 1 683 1 114 8 2 125
## 684 1 684 1 114 8 2 125
## 685 1 685 1 114 8 2 125
## 686 1 686 1 113 9 2 125
## 687 1 687 1 112 10 2 125
## 688 1 688 1 112 10 2 125
## 689 1 689 1 112 10 2 125
## 690 1 690 1 111 11 2 125
## 691 1 691 1 111 11 2 125
## 692 1 692 1 111 11 2 125
## 693 1 693 1 111 11 2 125
## 694 1 694 1 109 13 2 125
## 695 1 695 1 109 13 2 125
## 696 1 696 1 108 14 2 125
## 697 1 697 1 108 14 2 125
## 698 1 698 1 108 14 2 125
## 699 1 699 1 109 14 1 125
## 700 1 700 1 109 14 1 125
## 701 1 701 1 109 14 1 125
## 702 1 702 1 108 15 1 125
## 703 1 703 2 107 16 0 125
## 704 1 704 2 107 16 0 125
## 705 1 705 2 107 16 0 125
## 706 1 706 2 106 17 0 125
## 707 1 707 2 105 18 0 125
## 708 1 708 2 105 18 0 125
## 709 1 709 1 106 18 0 125
## 710 1 710 1 106 18 0 125
## 711 1 711 1 105 19 0 125
## 712 1 712 1 105 19 0 125
## 713 1 713 1 105 19 0 125
## 714 1 714 0 105 20 0 125
## 715 1 715 0 104 21 0 125
## 716 1 716 0 104 21 0 125
## 717 1 717 0 103 22 0 125
## 718 1 718 0 103 22 0 125
## 719 1 719 0 103 22 0 125
## 720 1 720 0 103 22 0 125
## 721 1 721 0 103 22 0 125
## 722 1 722 0 103 22 0 125
## 723 1 723 0 103 22 0 125
## 724 1 724 0 103 22 0 125
## 725 1 725 0 102 23 0 125
## 726 1 726 0 102 23 0 125
## 727 1 727 0 102 23 0 125
## 728 1 728 0 102 23 0 125
## 729 1 729 0 102 23 0 125
## 730 1 730 0 102 23 0 125
## 731 1 731 0 102 23 0 125
## 732 1 732 0 102 23 0 125
## 733 1 733 0 102 23 0 125
## 734 1 734 0 102 23 0 125
## 735 1 735 0 102 23 0 125
## 736 1 736 0 102 23 0 125
## 737 1 737 0 102 23 0 125
## 738 1 738 0 101 24 0 125
## 739 1 739 0 101 24 0 125
## 740 1 740 0 101 24 0 125
## 741 1 741 0 101 24 0 125
## 742 1 742 0 101 24 0 125
## 743 1 743 0 101 24 0 125
## 744 1 744 0 101 24 0 125
## 745 1 745 0 101 24 0 125
## 746 1 746 0 101 24 0 125
## 747 1 747 0 101 24 0 125
## 748 1 748 0 101 24 0 125
## 749 1 749 0 101 24 0 125
## 750 1 750 0 101 24 0 125
## 751 1 751 0 101 24 0 125
## 752 1 752 0 100 25 0 125
## 753 1 753 0 100 25 0 125
## 754 1 754 0 100 25 0 125
## 755 1 755 0 100 25 0 125
## 756 1 756 0 100 25 0 125
## 757 1 757 0 100 25 0 125
## 758 1 758 0 100 25 0 125
## 759 1 759 0 99 26 0 125
## 760 1 760 0 99 26 0 125
## 761 1 761 0 99 26 0 125
## 762 1 762 0 99 26 0 125
## 763 1 763 0 99 26 0 125
## 764 1 764 0 99 26 0 125
## 765 1 765 0 99 26 0 125
## 766 1 766 0 99 26 0 125
## 767 1 767 0 99 26 0 125
## 768 1 768 0 98 27 0 125
## 769 1 769 0 98 27 0 125
## 770 1 770 0 97 28 0 125
## 771 1 771 0 97 28 0 125
## 772 1 772 0 96 29 0 125
## 773 1 773 0 96 29 0 125
## 774 1 774 0 96 29 0 125
## 775 1 775 0 95 30 0 125
## 776 1 776 0 95 30 0 125
## 777 1 777 0 95 30 0 125
## 778 1 778 0 95 30 0 125
## 779 1 779 0 95 30 0 125
## 780 1 780 0 95 30 0 125
## 781 1 781 0 95 30 0 125
## 782 1 782 0 94 31 0 125
## 783 1 783 0 94 31 0 125
## 784 1 784 0 94 31 0 125
## 785 1 785 0 93 32 0 125
## 786 1 786 0 93 32 0 125
## 787 1 787 0 93 32 0 125
## 788 1 788 0 93 32 0 125
## 789 1 789 0 93 32 0 125
## 790 1 790 0 93 32 0 125
## 791 1 791 0 93 32 0 125
## 792 1 792 0 93 32 0 125
## 793 1 793 0 93 32 0 125
## 794 1 794 0 93 32 0 125
## 795 1 795 0 93 32 0 125
## 796 1 796 0 93 32 0 125
## 797 1 797 0 93 32 0 125
## 798 1 798 0 91 34 0 125
## 799 1 799 0 91 34 0 125
## 800 1 800 0 89 36 0 125
## 801 1 801 0 89 36 0 125
## 802 1 802 0 88 37 0 125
## 803 1 803 0 88 37 0 125
## 804 1 804 0 88 37 0 125
## 805 1 805 0 88 37 0 125
## 806 1 806 0 88 37 0 125
## 807 1 807 0 87 38 0 125
## 808 1 808 0 87 38 0 125
## 809 1 809 0 86 39 0 125
## 810 1 810 0 85 40 0 125
## 811 1 811 0 85 40 0 125
## 812 1 812 0 85 40 0 125
## 813 1 813 0 84 41 0 125
## 814 1 814 0 84 41 0 125
## 815 1 815 0 84 41 0 125
## 816 1 816 0 84 41 0 125
## 817 1 817 0 83 42 0 125
## 818 1 818 0 82 43 0 125
## 819 1 819 0 82 43 0 125
## 820 1 820 0 82 43 0 125
## 821 1 821 0 82 43 0 125
## 822 1 822 0 82 43 0 125
## 823 1 823 0 82 43 0 125
## 824 1 824 0 82 43 0 125
## 825 1 825 0 82 43 0 125
## 826 1 826 0 82 43 0 125
## 827 1 827 0 82 43 0 125
## 828 1 828 0 82 43 0 125
## 829 1 829 0 82 43 0 125
## 830 1 830 0 82 43 0 125
## 831 1 831 0 82 43 0 125
## 832 1 832 0 82 43 0 125
## 833 1 833 0 82 43 0 125
## 834 1 834 0 81 44 0 125
## 835 1 835 0 80 45 0 125
## 836 1 836 0 80 45 0 125
## 837 1 837 0 80 45 0 125
## 838 1 838 0 80 45 0 125
## 839 1 839 0 80 45 0 125
## 840 1 840 0 80 45 0 125
## 841 1 841 0 80 45 0 125
## 842 1 842 0 79 46 0 125
## 843 1 843 0 79 46 0 125
## 844 1 844 0 79 46 0 125
## 845 1 845 0 79 46 0 125
## 846 1 846 0 79 46 0 125
## 847 1 847 0 78 47 0 125
## 848 1 848 0 78 47 0 125
## 849 1 849 0 78 47 0 125
## 850 1 850 0 78 47 0 125
## 851 1 851 0 78 47 0 125
## 852 1 852 0 77 48 0 125
## 853 1 853 0 76 49 0 125
## 854 1 854 0 76 49 0 125
## 855 1 855 0 76 49 0 125
## 856 1 856 0 76 49 0 125
## 857 1 857 0 76 49 0 125
## 858 1 858 0 75 50 0 125
## 859 1 859 0 75 50 0 125
## 860 1 860 0 75 50 0 125
## 861 1 861 0 75 50 0 125
## 862 1 862 0 75 50 0 125
## 863 1 863 0 75 50 0 125
## 864 1 864 0 75 50 0 125
## 865 1 865 0 75 50 0 125
## 866 1 866 0 75 50 0 125
## 867 1 867 0 75 50 0 125
## 868 1 868 0 75 50 0 125
## 869 1 869 0 74 51 0 125
## 870 1 870 0 73 52 0 125
## 871 1 871 0 72 53 0 125
## 872 1 872 0 72 53 0 125
## 873 1 873 0 71 54 0 125
## 874 1 874 0 70 55 0 125
## 875 1 875 0 70 55 0 125
## 876 1 876 0 70 55 0 125
## 877 1 877 0 69 56 0 125
## 878 1 878 0 69 56 0 125
## 879 1 879 0 69 56 0 125
## 880 1 880 0 69 56 0 125
## 881 1 881 0 69 56 0 125
## 882 1 882 0 69 56 0 125
## 883 1 883 0 69 56 0 125
## 884 1 884 0 69 56 0 125
## 885 1 885 0 69 56 0 125
## 886 1 886 0 68 57 0 125
## 887 1 887 0 68 57 0 125
## 888 1 888 0 68 57 0 125
## 889 1 889 0 68 57 0 125
## 890 1 890 0 68 57 0 125
## 891 1 891 0 68 57 0 125
## 892 1 892 0 68 57 0 125
## 893 1 893 0 68 57 0 125
## 894 1 894 0 66 59 0 125
## 895 1 895 0 66 59 0 125
## 896 1 896 0 66 59 0 125
## 897 1 897 0 66 59 0 125
## 898 1 898 0 66 59 0 125
## 899 1 899 0 66 59 0 125
## 900 1 900 0 66 59 0 125
## 901 1 901 0 66 59 0 125
## 902 1 902 0 66 59 0 125
## 903 1 903 0 66 59 0 125
## 904 1 904 0 66 59 0 125
## 905 1 905 0 66 59 0 125
## 906 1 906 0 66 59 0 125
## 907 1 907 0 66 59 0 125
## 908 1 908 0 66 59 0 125
## 909 1 909 0 66 59 0 125
## 910 1 910 0 66 59 0 125
## 911 1 911 0 65 60 0 125
## 912 1 912 0 65 60 0 125
## 913 1 913 0 65 60 0 125
## 914 1 914 0 64 61 0 125
## 915 1 915 0 64 61 0 125
## 916 1 916 0 64 61 0 125
## 917 1 917 0 64 61 0 125
## 918 1 918 0 64 61 0 125
## 919 1 919 0 64 61 0 125
## 920 1 920 0 64 61 0 125
## 921 1 921 0 64 61 0 125
## 922 1 922 0 64 61 0 125
## 923 1 923 0 63 62 0 125
## 924 1 924 0 62 63 0 125
## 925 1 925 0 62 63 0 125
## 926 1 926 0 62 63 0 125
## 927 1 927 0 62 63 0 125
## 928 1 928 0 62 63 0 125
## 929 1 929 0 61 64 0 125
## 930 1 930 0 61 64 0 125
## 931 1 931 0 61 64 0 125
## 932 1 932 0 61 64 0 125
## 933 1 933 0 61 64 0 125
## 934 1 934 0 61 64 0 125
## 935 1 935 0 61 64 0 125
## 936 1 936 0 60 65 0 125
## 937 1 937 0 60 65 0 125
## 938 1 938 0 60 65 0 125
## 939 1 939 0 60 65 0 125
## 940 1 940 0 60 65 0 125
## 941 1 941 0 60 65 0 125
## 942 1 942 0 60 65 0 125
## 943 1 943 0 60 65 0 125
## 944 1 944 0 60 65 0 125
## 945 1 945 0 60 65 0 125
## 946 1 946 0 60 65 0 125
## 947 1 947 0 60 65 0 125
## 948 1 948 0 60 65 0 125
## 949 1 949 0 60 65 0 125
## 950 1 950 0 60 65 0 125
## 951 1 951 0 60 65 0 125
## 952 1 952 0 60 65 0 125
## 953 1 953 0 60 65 0 125
## 954 1 954 0 60 65 0 125
## 955 1 955 0 60 65 0 125
## 956 1 956 0 60 65 0 125
## 957 1 957 0 60 65 0 125
## 958 1 958 0 60 65 0 125
## 959 1 959 0 60 65 0 125
## 960 1 960 0 60 65 0 125
## 961 1 961 0 60 65 0 125
## 962 1 962 0 60 65 0 125
## 963 1 963 0 60 65 0 125
## 964 1 964 0 60 65 0 125
## 965 1 965 0 60 65 0 125
## 966 1 966 0 60 65 0 125
## 967 1 967 0 60 65 0 125
## 968 1 968 0 60 65 0 125
## 969 1 969 0 60 65 0 125
## 970 1 970 0 59 66 0 125
## 971 1 971 0 59 66 0 125
## 972 1 972 0 59 66 0 125
## 973 1 973 0 59 66 0 125
## 974 1 974 0 59 66 0 125
## 975 1 975 0 59 66 0 125
## 976 1 976 0 59 66 0 125
## 977 1 977 0 59 66 0 125
## 978 1 978 0 58 67 0 125
## 979 1 979 0 58 67 0 125
## 980 1 980 0 58 67 0 125
## 981 1 981 0 57 68 0 125
## 982 1 982 0 57 68 0 125
## 983 1 983 0 57 68 0 125
## 984 1 984 0 57 68 0 125
## 985 1 985 0 57 68 0 125
## 986 1 986 0 57 68 0 125
## 987 1 987 0 57 68 0 125
## 988 1 988 0 57 68 0 125
## 989 1 989 0 57 68 0 125
## 990 1 990 0 57 68 0 125
## 991 1 991 0 56 69 0 125
## 992 1 992 0 56 69 0 125
## 993 1 993 0 56 69 0 125
## 994 1 994 0 55 70 0 125
## 995 1 995 0 54 71 0 125
## 996 1 996 0 54 71 0 125
## 997 1 997 0 54 71 0 125
## 998 1 998 0 54 71 0 125
## 999 1 999 0 54 71 0 125
## 1000 1 1000 0 54 71 0 1253.1 Trajectory
3.2 ? What the index is refering to
trajectory(result3, compartments = "I", index = 1)## node time I
## 1 1 1 5
## 2 1 2 3
## 3 1 3 3
## 4 1 4 3
## 5 1 5 3
## 6 1 6 3
## 7 1 7 4
## 8 1 8 4
## 9 1 9 3
## 10 1 10 2
## 11 1 11 2
## 12 1 12 3
## 13 1 13 3
## 14 1 14 2
## 15 1 15 2
## 16 1 16 2
## 17 1 17 2
## 18 1 18 2
## 19 1 19 2
## 20 1 20 2
## 21 1 21 2
## 22 1 22 2
## 23 1 23 1
## 24 1 24 1
## 25 1 25 1
## 26 1 26 1
## 27 1 27 1
## 28 1 28 1
## 29 1 29 1
## 30 1 30 1
## 31 1 31 1
## 32 1 32 2
## 33 1 33 2
## 34 1 34 3
## 35 1 35 3
## 36 1 36 3
## 37 1 37 3
## 38 1 38 1
## 39 1 39 1
## 40 1 40 1
## 41 1 41 0
## 42 1 42 0
## 43 1 43 0
## 44 1 44 0
## 45 1 45 0
## 46 1 46 0
## 47 1 47 0
## 48 1 48 0
## 49 1 49 0
## 50 1 50 1
## 51 1 51 2
## 52 1 52 2
## 53 1 53 2
## 54 1 54 2
## 55 1 55 3
## 56 1 56 3
## 57 1 57 3
## 58 1 58 3
## 59 1 59 3
## 60 1 60 3
## 61 1 61 3
## 62 1 62 4
## 63 1 63 3
## 64 1 64 3
## 65 1 65 3
## 66 1 66 3
## 67 1 67 4
## 68 1 68 3
## 69 1 69 3
## 70 1 70 3
## 71 1 71 3
## 72 1 72 3
## 73 1 73 3
## 74 1 74 2
## 75 1 75 2
## 76 1 76 2
## 77 1 77 2
## 78 1 78 3
## 79 1 79 3
## 80 1 80 3
## 81 1 81 3
## 82 1 82 3
## 83 1 83 3
## 84 1 84 2
## 85 1 85 2
## 86 1 86 2
## 87 1 87 4
## 88 1 88 4
## 89 1 89 4
## 90 1 90 4
## 91 1 91 4
## 92 1 92 4
## 93 1 93 5
## 94 1 94 5
## 95 1 95 5
## 96 1 96 5
## 97 1 97 5
## 98 1 98 5
## 99 1 99 6
## 100 1 100 6
## 101 1 101 5
## 102 1 102 4
## 103 1 103 2
## 104 1 104 2
## 105 1 105 0
## 106 1 106 0
## 107 1 107 0
## 108 1 108 0
## 109 1 109 0
## 110 1 110 0
## 111 1 111 1
## 112 1 112 3
## 113 1 113 3
## 114 1 114 3
## 115 1 115 4
## 116 1 116 4
## 117 1 117 4
## 118 1 118 4
## 119 1 119 4
## 120 1 120 4
## 121 1 121 4
## 122 1 122 4
## 123 1 123 4
## 124 1 124 4
## 125 1 125 4
## 126 1 126 5
## 127 1 127 6
## 128 1 128 6
## 129 1 129 5
## 130 1 130 3
## 131 1 131 2
## 132 1 132 3
## 133 1 133 3
## 134 1 134 3
## 135 1 135 3
## 136 1 136 3
## 137 1 137 3
## 138 1 138 3
## 139 1 139 3
## 140 1 140 3
## 141 1 141 2
## 142 1 142 2
## 143 1 143 2
## 144 1 144 3
## 145 1 145 3
## 146 1 146 3
## 147 1 147 4
## 148 1 148 4
## 149 1 149 4
## 150 1 150 3
## 151 1 151 3
## 152 1 152 3
## 153 1 153 3
## 154 1 154 4
## 155 1 155 4
## 156 1 156 5
## 157 1 157 6
## 158 1 158 5
## 159 1 159 6
## 160 1 160 6
## 161 1 161 5
## 162 1 162 4
## 163 1 163 7
## 164 1 164 8
## 165 1 165 8
## 166 1 166 8
## 167 1 167 9
## 168 1 168 9
## 169 1 169 7
## 170 1 170 6
## 171 1 171 5
## 172 1 172 5
## 173 1 173 5
## 174 1 174 5
## 175 1 175 5
## 176 1 176 5
## 177 1 177 3
## 178 1 178 3
## 179 1 179 1
## 180 1 180 1
## 181 1 181 1
## 182 1 182 1
## 183 1 183 1
## 184 1 184 2
## 185 1 185 2
## 186 1 186 3
## 187 1 187 3
## 188 1 188 4
## 189 1 189 5
## 190 1 190 4
## 191 1 191 5
## 192 1 192 5
## 193 1 193 6
## 194 1 194 6
## 195 1 195 6
## 196 1 196 6
## 197 1 197 6
## 198 1 198 6
## 199 1 199 6
## 200 1 200 6
## 201 1 201 6
## 202 1 202 5
## 203 1 203 6
## 204 1 204 6
## 205 1 205 6
## 206 1 206 4
## 207 1 207 5
## 208 1 208 4
## 209 1 209 4
## 210 1 210 4
## 211 1 211 3
## 212 1 212 2
## 213 1 213 3
## 214 1 214 3
## 215 1 215 3
## 216 1 216 3
## 217 1 217 3
## 218 1 218 3
## 219 1 219 3
## 220 1 220 4
## 221 1 221 6
## 222 1 222 6
## 223 1 223 6
## 224 1 224 6
## 225 1 225 6
## 226 1 226 7
## 227 1 227 7
## 228 1 228 6
## 229 1 229 5
## 230 1 230 5
## 231 1 231 5
## 232 1 232 5
## 233 1 233 4
## 234 1 234 4
## 235 1 235 3
## 236 1 236 2
## 237 1 237 2
## 238 1 238 3
## 239 1 239 3
## 240 1 240 3
## 241 1 241 3
## 242 1 242 1
## 243 1 243 1
## 244 1 244 1
## 245 1 245 1
## 246 1 246 1
## 247 1 247 1
## 248 1 248 1
## 249 1 249 1
## 250 1 250 1
## 251 1 251 1
## 252 1 252 1
## 253 1 253 1
## 254 1 254 1
## 255 1 255 1
## 256 1 256 1
## 257 1 257 0
## 258 1 258 1
## 259 1 259 1
## 260 1 260 1
## 261 1 261 1
## 262 1 262 0
## 263 1 263 0
## 264 1 264 0
## 265 1 265 0
## 266 1 266 0
## 267 1 267 0
## 268 1 268 0
## 269 1 269 0
## 270 1 270 0
## 271 1 271 0
## 272 1 272 0
## 273 1 273 0
## 274 1 274 0
## 275 1 275 0
## 276 1 276 1
## 277 1 277 1
## 278 1 278 1
## 279 1 279 1
## 280 1 280 1
## 281 1 281 1
## 282 1 282 1
## 283 1 283 1
## 284 1 284 2
## 285 1 285 2
## 286 1 286 3
## 287 1 287 3
## 288 1 288 3
## 289 1 289 3
## 290 1 290 3
## 291 1 291 3
## 292 1 292 3
## 293 1 293 3
## 294 1 294 3
## 295 1 295 3
## 296 1 296 3
## 297 1 297 2
## 298 1 298 2
## 299 1 299 2
## 300 1 300 2
## 301 1 301 1
## 302 1 302 2
## 303 1 303 2
## 304 1 304 2
## 305 1 305 3
## 306 1 306 3
## 307 1 307 3
## 308 1 308 3
## 309 1 309 3
## 310 1 310 4
## 311 1 311 4
## 312 1 312 4
## 313 1 313 5
## 314 1 314 5
## 315 1 315 5
## 316 1 316 5
## 317 1 317 6
## 318 1 318 6
## 319 1 319 7
## 320 1 320 8
## 321 1 321 8
## 322 1 322 8
## 323 1 323 8
## 324 1 324 9
## 325 1 325 11
## 326 1 326 11
## 327 1 327 12
## 328 1 328 13
## 329 1 329 13
## 330 1 330 13
## 331 1 331 13
## 332 1 332 13
## 333 1 333 13
## 334 1 334 13
## 335 1 335 13
## 336 1 336 12
## 337 1 337 12
## 338 1 338 12
## 339 1 339 12
## 340 1 340 13
## 341 1 341 13
## 342 1 342 13
## 343 1 343 13
## 344 1 344 13
## 345 1 345 13
## 346 1 346 13
## 347 1 347 13
## 348 1 348 14
## 349 1 349 13
## 350 1 350 12
## 351 1 351 12
## 352 1 352 13
## 353 1 353 13
## 354 1 354 13
## 355 1 355 14
## 356 1 356 14
## 357 1 357 14
## 358 1 358 14
## 359 1 359 13
## 360 1 360 13
## 361 1 361 13
## 362 1 362 13
## 363 1 363 12
## 364 1 364 12
## 365 1 365 12
## 366 1 366 12
## 367 1 367 12
## 368 1 368 12
## 369 1 369 12
## 370 1 370 12
## 371 1 371 10
## 372 1 372 11
## 373 1 373 11
## 374 1 374 10
## 375 1 375 10
## 376 1 376 10
## 377 1 377 11
## 378 1 378 9
## 379 1 379 10
## 380 1 380 10
## 381 1 381 9
## 382 1 382 8
## 383 1 383 9
## 384 1 384 9
## 385 1 385 9
## 386 1 386 8
## 387 1 387 8
## 388 1 388 7
## 389 1 389 7
## 390 1 390 7
## 391 1 391 7
## 392 1 392 7
## 393 1 393 7
## 394 1 394 8
## 395 1 395 9
## 396 1 396 9
## 397 1 397 10
## 398 1 398 9
## 399 1 399 9
## 400 1 400 9
## 401 1 401 9
## 402 1 402 10
## 403 1 403 10
## 404 1 404 8
## 405 1 405 8
## 406 1 406 8
## 407 1 407 8
## 408 1 408 7
## 409 1 409 8
## 410 1 410 8
## 411 1 411 7
## 412 1 412 7
## 413 1 413 6
## 414 1 414 6
## 415 1 415 5
## 416 1 416 5
## 417 1 417 4
## 418 1 418 5
## 419 1 419 5
## 420 1 420 4
## 421 1 421 4
## 422 1 422 4
## 423 1 423 5
## 424 1 424 4
## 425 1 425 5
## 426 1 426 4
## 427 1 427 2
## 428 1 428 2
## 429 1 429 2
## 430 1 430 2
## 431 1 431 2
## 432 1 432 1
## 433 1 433 1
## 434 1 434 1
## 435 1 435 1
## 436 1 436 2
## 437 1 437 2
## 438 1 438 2
## 439 1 439 2
## 440 1 440 1
## 441 1 441 1
## 442 1 442 1
## 443 1 443 1
## 444 1 444 2
## 445 1 445 1
## 446 1 446 0
## 447 1 447 0
## 448 1 448 0
## 449 1 449 0
## 450 1 450 1
## 451 1 451 1
## 452 1 452 0
## 453 1 453 0
## 454 1 454 2
## 455 1 455 2
## 456 1 456 1
## 457 1 457 1
## 458 1 458 1
## 459 1 459 1
## 460 1 460 1
## 461 1 461 1
## 462 1 462 1
## 463 1 463 0
## 464 1 464 0
## 465 1 465 0
## 466 1 466 1
## 467 1 467 1
## 468 1 468 1
## 469 1 469 1
## 470 1 470 2
## 471 1 471 1
## 472 1 472 1
## 473 1 473 1
## 474 1 474 1
## 475 1 475 2
## 476 1 476 2
## 477 1 477 2
## 478 1 478 3
## 479 1 479 2
## 480 1 480 2
## 481 1 481 3
## 482 1 482 4
## 483 1 483 4
## 484 1 484 3
## 485 1 485 3
## 486 1 486 3
## 487 1 487 3
## 488 1 488 3
## 489 1 489 3
## 490 1 490 4
## 491 1 491 4
## 492 1 492 4
## 493 1 493 4
## 494 1 494 4
## 495 1 495 4
## 496 1 496 4
## 497 1 497 3
## 498 1 498 3
## 499 1 499 4
## 500 1 500 4
## 501 1 501 4
## 502 1 502 4
## 503 1 503 4
## 504 1 504 4
## 505 1 505 4
## 506 1 506 5
## 507 1 507 4
## 508 1 508 3
## 509 1 509 4
## 510 1 510 3
## 511 1 511 4
## 512 1 512 4
## 513 1 513 5
## 514 1 514 6
## 515 1 515 6
## 516 1 516 5
## 517 1 517 5
## 518 1 518 5
## 519 1 519 5
## 520 1 520 5
## 521 1 521 5
## 522 1 522 5
## 523 1 523 6
## 524 1 524 6
## 525 1 525 6
## 526 1 526 6
## 527 1 527 6
## 528 1 528 6
## 529 1 529 6
## 530 1 530 6
## 531 1 531 5
## 532 1 532 6
## 533 1 533 7
## 534 1 534 7
## 535 1 535 7
## 536 1 536 8
## 537 1 537 8
## 538 1 538 7
## 539 1 539 8
## 540 1 540 7
## 541 1 541 7
## 542 1 542 7
## 543 1 543 6
## 544 1 544 6
## 545 1 545 6
## 546 1 546 6
## 547 1 547 6
## 548 1 548 6
## 549 1 549 6
## 550 1 550 6
## 551 1 551 6
## 552 1 552 6
## 553 1 553 7
## 554 1 554 7
## 555 1 555 7
## 556 1 556 7
## 557 1 557 8
## 558 1 558 8
## 559 1 559 9
## 560 1 560 10
## 561 1 561 10
## 562 1 562 10
## 563 1 563 11
## 564 1 564 12
## 565 1 565 11
## 566 1 566 11
## 567 1 567 11
## 568 1 568 11
## 569 1 569 10
## 570 1 570 10
## 571 1 571 9
## 572 1 572 7
## 573 1 573 7
## 574 1 574 7
## 575 1 575 7
## 576 1 576 7
## 577 1 577 7
## 578 1 578 7
## 579 1 579 7
## 580 1 580 7
## 581 1 581 7
## 582 1 582 6
## 583 1 583 7
## 584 1 584 6
## 585 1 585 6
## 586 1 586 7
## 587 1 587 8
## 588 1 588 7
## 589 1 589 5
## 590 1 590 5
## 591 1 591 6
## 592 1 592 7
## 593 1 593 8
## 594 1 594 8
## 595 1 595 7
## 596 1 596 8
## 597 1 597 8
## 598 1 598 8
## 599 1 599 8
## 600 1 600 8
## 601 1 601 8
## 602 1 602 8
## 603 1 603 8
## 604 1 604 8
## 605 1 605 7
## 606 1 606 7
## 607 1 607 9
## 608 1 608 7
## 609 1 609 7
## 610 1 610 7
## 611 1 611 7
## 612 1 612 6
## 613 1 613 5
## 614 1 614 5
## 615 1 615 5
## 616 1 616 7
## 617 1 617 6
## 618 1 618 5
## 619 1 619 4
## 620 1 620 4
## 621 1 621 4
## 622 1 622 5
## 623 1 623 5
## 624 1 624 5
## 625 1 625 4
## 626 1 626 4
## 627 1 627 4
## 628 1 628 4
## 629 1 629 4
## 630 1 630 4
## 631 1 631 5
## 632 1 632 5
## 633 1 633 4
## 634 1 634 4
## 635 1 635 5
## 636 1 636 5
## 637 1 637 5
## 638 1 638 5
## 639 1 639 5
## 640 1 640 4
## 641 1 641 4
## 642 1 642 4
## 643 1 643 4
## 644 1 644 4
## 645 1 645 4
## 646 1 646 5
## 647 1 647 5
## 648 1 648 5
## 649 1 649 5
## 650 1 650 5
## 651 1 651 5
## 652 1 652 5
## 653 1 653 6
## 654 1 654 6
## 655 1 655 6
## 656 1 656 6
## 657 1 657 6
## 658 1 658 7
## 659 1 659 7
## 660 1 660 7
## 661 1 661 7
## 662 1 662 7
## 663 1 663 7
## 664 1 664 8
## 665 1 665 8
## 666 1 666 8
## 667 1 667 7
## 668 1 668 7
## 669 1 669 7
## 670 1 670 7
## 671 1 671 7
## 672 1 672 7
## 673 1 673 7
## 674 1 674 8
## 675 1 675 8
## 676 1 676 8
## 677 1 677 8
## 678 1 678 8
## 679 1 679 8
## 680 1 680 8
## 681 1 681 9
## 682 1 682 8
## 683 1 683 8
## 684 1 684 8
## 685 1 685 8
## 686 1 686 9
## 687 1 687 10
## 688 1 688 10
## 689 1 689 10
## 690 1 690 11
## 691 1 691 11
## 692 1 692 11
## 693 1 693 11
## 694 1 694 13
## 695 1 695 13
## 696 1 696 14
## 697 1 697 14
## 698 1 698 14
## 699 1 699 14
## 700 1 700 14
## 701 1 701 14
## 702 1 702 15
## 703 1 703 16
## 704 1 704 16
## 705 1 705 16
## 706 1 706 17
## 707 1 707 18
## 708 1 708 18
## 709 1 709 18
## 710 1 710 18
## 711 1 711 19
## 712 1 712 19
## 713 1 713 19
## 714 1 714 20
## 715 1 715 21
## 716 1 716 21
## 717 1 717 22
## 718 1 718 22
## 719 1 719 22
## 720 1 720 22
## 721 1 721 22
## 722 1 722 22
## 723 1 723 22
## 724 1 724 22
## 725 1 725 23
## 726 1 726 23
## 727 1 727 23
## 728 1 728 23
## 729 1 729 23
## 730 1 730 23
## 731 1 731 23
## 732 1 732 23
## 733 1 733 23
## 734 1 734 23
## 735 1 735 23
## 736 1 736 23
## 737 1 737 23
## 738 1 738 24
## 739 1 739 24
## 740 1 740 24
## 741 1 741 24
## 742 1 742 24
## 743 1 743 24
## 744 1 744 24
## 745 1 745 24
## 746 1 746 24
## 747 1 747 24
## 748 1 748 24
## 749 1 749 24
## 750 1 750 24
## 751 1 751 24
## 752 1 752 25
## 753 1 753 25
## 754 1 754 25
## 755 1 755 25
## 756 1 756 25
## 757 1 757 25
## 758 1 758 25
## 759 1 759 26
## 760 1 760 26
## 761 1 761 26
## 762 1 762 26
## 763 1 763 26
## 764 1 764 26
## 765 1 765 26
## 766 1 766 26
## 767 1 767 26
## 768 1 768 27
## 769 1 769 27
## 770 1 770 28
## 771 1 771 28
## 772 1 772 29
## 773 1 773 29
## 774 1 774 29
## 775 1 775 30
## 776 1 776 30
## 777 1 777 30
## 778 1 778 30
## 779 1 779 30
## 780 1 780 30
## 781 1 781 30
## 782 1 782 31
## 783 1 783 31
## 784 1 784 31
## 785 1 785 32
## 786 1 786 32
## 787 1 787 32
## 788 1 788 32
## 789 1 789 32
## 790 1 790 32
## 791 1 791 32
## 792 1 792 32
## 793 1 793 32
## 794 1 794 32
## 795 1 795 32
## 796 1 796 32
## 797 1 797 32
## 798 1 798 34
## 799 1 799 34
## 800 1 800 36
## 801 1 801 36
## 802 1 802 37
## 803 1 803 37
## 804 1 804 37
## 805 1 805 37
## 806 1 806 37
## 807 1 807 38
## 808 1 808 38
## 809 1 809 39
## 810 1 810 40
## 811 1 811 40
## 812 1 812 40
## 813 1 813 41
## 814 1 814 41
## 815 1 815 41
## 816 1 816 41
## 817 1 817 42
## 818 1 818 43
## 819 1 819 43
## 820 1 820 43
## 821 1 821 43
## 822 1 822 43
## 823 1 823 43
## 824 1 824 43
## 825 1 825 43
## 826 1 826 43
## 827 1 827 43
## 828 1 828 43
## 829 1 829 43
## 830 1 830 43
## 831 1 831 43
## 832 1 832 43
## 833 1 833 43
## 834 1 834 44
## 835 1 835 45
## 836 1 836 45
## 837 1 837 45
## 838 1 838 45
## 839 1 839 45
## 840 1 840 45
## 841 1 841 45
## 842 1 842 46
## 843 1 843 46
## 844 1 844 46
## 845 1 845 46
## 846 1 846 46
## 847 1 847 47
## 848 1 848 47
## 849 1 849 47
## 850 1 850 47
## 851 1 851 47
## 852 1 852 48
## 853 1 853 49
## 854 1 854 49
## 855 1 855 49
## 856 1 856 49
## 857 1 857 49
## 858 1 858 50
## 859 1 859 50
## 860 1 860 50
## 861 1 861 50
## 862 1 862 50
## 863 1 863 50
## 864 1 864 50
## 865 1 865 50
## 866 1 866 50
## 867 1 867 50
## 868 1 868 50
## 869 1 869 51
## 870 1 870 52
## 871 1 871 53
## 872 1 872 53
## 873 1 873 54
## 874 1 874 55
## 875 1 875 55
## 876 1 876 55
## 877 1 877 56
## 878 1 878 56
## 879 1 879 56
## 880 1 880 56
## 881 1 881 56
## 882 1 882 56
## 883 1 883 56
## 884 1 884 56
## 885 1 885 56
## 886 1 886 57
## 887 1 887 57
## 888 1 888 57
## 889 1 889 57
## 890 1 890 57
## 891 1 891 57
## 892 1 892 57
## 893 1 893 57
## 894 1 894 59
## 895 1 895 59
## 896 1 896 59
## 897 1 897 59
## 898 1 898 59
## 899 1 899 59
## 900 1 900 59
## 901 1 901 59
## 902 1 902 59
## 903 1 903 59
## 904 1 904 59
## 905 1 905 59
## 906 1 906 59
## 907 1 907 59
## 908 1 908 59
## 909 1 909 59
## 910 1 910 59
## 911 1 911 60
## 912 1 912 60
## 913 1 913 60
## 914 1 914 61
## 915 1 915 61
## 916 1 916 61
## 917 1 917 61
## 918 1 918 61
## 919 1 919 61
## 920 1 920 61
## 921 1 921 61
## 922 1 922 61
## 923 1 923 62
## 924 1 924 63
## 925 1 925 63
## 926 1 926 63
## 927 1 927 63
## 928 1 928 63
## 929 1 929 64
## 930 1 930 64
## 931 1 931 64
## 932 1 932 64
## 933 1 933 64
## 934 1 934 64
## 935 1 935 64
## 936 1 936 65
## 937 1 937 65
## 938 1 938 65
## 939 1 939 65
## 940 1 940 65
## 941 1 941 65
## 942 1 942 65
## 943 1 943 65
## 944 1 944 65
## 945 1 945 65
## 946 1 946 65
## 947 1 947 65
## 948 1 948 65
## 949 1 949 65
## 950 1 950 65
## 951 1 951 65
## 952 1 952 65
## 953 1 953 65
## 954 1 954 65
## 955 1 955 65
## 956 1 956 65
## 957 1 957 65
## 958 1 958 65
## 959 1 959 65
## 960 1 960 65
## 961 1 961 65
## 962 1 962 65
## 963 1 963 65
## 964 1 964 65
## 965 1 965 65
## 966 1 966 65
## 967 1 967 65
## 968 1 968 65
## 969 1 969 65
## 970 1 970 66
## 971 1 971 66
## 972 1 972 66
## 973 1 973 66
## 974 1 974 66
## 975 1 975 66
## 976 1 976 66
## 977 1 977 66
## 978 1 978 67
## 979 1 979 67
## 980 1 980 67
## 981 1 981 68
## 982 1 982 68
## 983 1 983 68
## 984 1 984 68
## 985 1 985 68
## 986 1 986 68
## 987 1 987 68
## 988 1 988 68
## 989 1 989 68
## 990 1 990 68
## 991 1 991 69
## 992 1 992 69
## 993 1 993 69
## 994 1 994 70
## 995 1 995 71
## 996 1 996 71
## 997 1 997 71
## 998 1 998 71
## 999 1 999 71
## 1000 1 1000 713.3 Prevalence
prevalence(result3, I~.)## time prevalence
## 1 1 0.040
## 2 2 0.024
## 3 3 0.024
## 4 4 0.024
## 5 5 0.024
## 6 6 0.024
## 7 7 0.032
## 8 8 0.032
## 9 9 0.024
## 10 10 0.016
## 11 11 0.016
## 12 12 0.024
## 13 13 0.024
## 14 14 0.016
## 15 15 0.016
## 16 16 0.016
## 17 17 0.016
## 18 18 0.016
## 19 19 0.016
## 20 20 0.016
## 21 21 0.016
## 22 22 0.016
## 23 23 0.008
## 24 24 0.008
## 25 25 0.008
## 26 26 0.008
## 27 27 0.008
## 28 28 0.008
## 29 29 0.008
## 30 30 0.008
## 31 31 0.008
## 32 32 0.016
## 33 33 0.016
## 34 34 0.024
## 35 35 0.024
## 36 36 0.024
## 37 37 0.024
## 38 38 0.008
## 39 39 0.008
## 40 40 0.008
## 41 41 0.000
## 42 42 0.000
## 43 43 0.000
## 44 44 0.000
## 45 45 0.000
## 46 46 0.000
## 47 47 0.000
## 48 48 0.000
## 49 49 0.000
## 50 50 0.008
## 51 51 0.016
## 52 52 0.016
## 53 53 0.016
## 54 54 0.016
## 55 55 0.024
## 56 56 0.024
## 57 57 0.024
## 58 58 0.024
## 59 59 0.024
## 60 60 0.024
## 61 61 0.024
## 62 62 0.032
## 63 63 0.024
## 64 64 0.024
## 65 65 0.024
## 66 66 0.024
## 67 67 0.032
## 68 68 0.024
## 69 69 0.024
## 70 70 0.024
## 71 71 0.024
## 72 72 0.024
## 73 73 0.024
## 74 74 0.016
## 75 75 0.016
## 76 76 0.016
## 77 77 0.016
## 78 78 0.024
## 79 79 0.024
## 80 80 0.024
## 81 81 0.024
## 82 82 0.024
## 83 83 0.024
## 84 84 0.016
## 85 85 0.016
## 86 86 0.016
## 87 87 0.032
## 88 88 0.032
## 89 89 0.032
## 90 90 0.032
## 91 91 0.032
## 92 92 0.032
## 93 93 0.040
## 94 94 0.040
## 95 95 0.040
## 96 96 0.040
## 97 97 0.040
## 98 98 0.040
## 99 99 0.048
## 100 100 0.048
## 101 101 0.040
## 102 102 0.032
## 103 103 0.016
## 104 104 0.016
## 105 105 0.000
## 106 106 0.000
## 107 107 0.000
## 108 108 0.000
## 109 109 0.000
## 110 110 0.000
## 111 111 0.008
## 112 112 0.024
## 113 113 0.024
## 114 114 0.024
## 115 115 0.032
## 116 116 0.032
## 117 117 0.032
## 118 118 0.032
## 119 119 0.032
## 120 120 0.032
## 121 121 0.032
## 122 122 0.032
## 123 123 0.032
## 124 124 0.032
## 125 125 0.032
## 126 126 0.040
## 127 127 0.048
## 128 128 0.048
## 129 129 0.040
## 130 130 0.024
## 131 131 0.016
## 132 132 0.024
## 133 133 0.024
## 134 134 0.024
## 135 135 0.024
## 136 136 0.024
## 137 137 0.024
## 138 138 0.024
## 139 139 0.024
## 140 140 0.024
## 141 141 0.016
## 142 142 0.016
## 143 143 0.016
## 144 144 0.024
## 145 145 0.024
## 146 146 0.024
## 147 147 0.032
## 148 148 0.032
## 149 149 0.032
## 150 150 0.024
## 151 151 0.024
## 152 152 0.024
## 153 153 0.024
## 154 154 0.032
## 155 155 0.032
## 156 156 0.040
## 157 157 0.048
## 158 158 0.040
## 159 159 0.048
## 160 160 0.048
## 161 161 0.040
## 162 162 0.032
## 163 163 0.056
## 164 164 0.064
## 165 165 0.064
## 166 166 0.064
## 167 167 0.072
## 168 168 0.072
## 169 169 0.056
## 170 170 0.048
## 171 171 0.040
## 172 172 0.040
## 173 173 0.040
## 174 174 0.040
## 175 175 0.040
## 176 176 0.040
## 177 177 0.024
## 178 178 0.024
## 179 179 0.008
## 180 180 0.008
## 181 181 0.008
## 182 182 0.008
## 183 183 0.008
## 184 184 0.016
## 185 185 0.016
## 186 186 0.024
## 187 187 0.024
## 188 188 0.032
## 189 189 0.040
## 190 190 0.032
## 191 191 0.040
## 192 192 0.040
## 193 193 0.048
## 194 194 0.048
## 195 195 0.048
## 196 196 0.048
## 197 197 0.048
## 198 198 0.048
## 199 199 0.048
## 200 200 0.048
## 201 201 0.048
## 202 202 0.040
## 203 203 0.048
## 204 204 0.048
## 205 205 0.048
## 206 206 0.032
## 207 207 0.040
## 208 208 0.032
## 209 209 0.032
## 210 210 0.032
## 211 211 0.024
## 212 212 0.016
## 213 213 0.024
## 214 214 0.024
## 215 215 0.024
## 216 216 0.024
## 217 217 0.024
## 218 218 0.024
## 219 219 0.024
## 220 220 0.032
## 221 221 0.048
## 222 222 0.048
## 223 223 0.048
## 224 224 0.048
## 225 225 0.048
## 226 226 0.056
## 227 227 0.056
## 228 228 0.048
## 229 229 0.040
## 230 230 0.040
## 231 231 0.040
## 232 232 0.040
## 233 233 0.032
## 234 234 0.032
## 235 235 0.024
## 236 236 0.016
## 237 237 0.016
## 238 238 0.024
## 239 239 0.024
## 240 240 0.024
## 241 241 0.024
## 242 242 0.008
## 243 243 0.008
## 244 244 0.008
## 245 245 0.008
## 246 246 0.008
## 247 247 0.008
## 248 248 0.008
## 249 249 0.008
## 250 250 0.008
## 251 251 0.008
## 252 252 0.008
## 253 253 0.008
## 254 254 0.008
## 255 255 0.008
## 256 256 0.008
## 257 257 0.000
## 258 258 0.008
## 259 259 0.008
## 260 260 0.008
## 261 261 0.008
## 262 262 0.000
## 263 263 0.000
## 264 264 0.000
## 265 265 0.000
## 266 266 0.000
## 267 267 0.000
## 268 268 0.000
## 269 269 0.000
## 270 270 0.000
## 271 271 0.000
## 272 272 0.000
## 273 273 0.000
## 274 274 0.000
## 275 275 0.000
## 276 276 0.008
## 277 277 0.008
## 278 278 0.008
## 279 279 0.008
## 280 280 0.008
## 281 281 0.008
## 282 282 0.008
## 283 283 0.008
## 284 284 0.016
## 285 285 0.016
## 286 286 0.024
## 287 287 0.024
## 288 288 0.024
## 289 289 0.024
## 290 290 0.024
## 291 291 0.024
## 292 292 0.024
## 293 293 0.024
## 294 294 0.024
## 295 295 0.024
## 296 296 0.024
## 297 297 0.016
## 298 298 0.016
## 299 299 0.016
## 300 300 0.016
## 301 301 0.008
## 302 302 0.016
## 303 303 0.016
## 304 304 0.016
## 305 305 0.024
## 306 306 0.024
## 307 307 0.024
## 308 308 0.024
## 309 309 0.024
## 310 310 0.032
## 311 311 0.032
## 312 312 0.032
## 313 313 0.040
## 314 314 0.040
## 315 315 0.040
## 316 316 0.040
## 317 317 0.048
## 318 318 0.048
## 319 319 0.056
## 320 320 0.064
## 321 321 0.064
## 322 322 0.064
## 323 323 0.064
## 324 324 0.072
## 325 325 0.088
## 326 326 0.088
## 327 327 0.096
## 328 328 0.104
## 329 329 0.104
## 330 330 0.104
## 331 331 0.104
## 332 332 0.104
## 333 333 0.104
## 334 334 0.104
## 335 335 0.104
## 336 336 0.096
## 337 337 0.096
## 338 338 0.096
## 339 339 0.096
## 340 340 0.104
## 341 341 0.104
## 342 342 0.104
## 343 343 0.104
## 344 344 0.104
## 345 345 0.104
## 346 346 0.104
## 347 347 0.104
## 348 348 0.112
## 349 349 0.104
## 350 350 0.096
## 351 351 0.096
## 352 352 0.104
## 353 353 0.104
## 354 354 0.104
## 355 355 0.112
## 356 356 0.112
## 357 357 0.112
## 358 358 0.112
## 359 359 0.104
## 360 360 0.104
## 361 361 0.104
## 362 362 0.104
## 363 363 0.096
## 364 364 0.096
## 365 365 0.096
## 366 366 0.096
## 367 367 0.096
## 368 368 0.096
## 369 369 0.096
## 370 370 0.096
## 371 371 0.080
## 372 372 0.088
## 373 373 0.088
## 374 374 0.080
## 375 375 0.080
## 376 376 0.080
## 377 377 0.088
## 378 378 0.072
## 379 379 0.080
## 380 380 0.080
## 381 381 0.072
## 382 382 0.064
## 383 383 0.072
## 384 384 0.072
## 385 385 0.072
## 386 386 0.064
## 387 387 0.064
## 388 388 0.056
## 389 389 0.056
## 390 390 0.056
## 391 391 0.056
## 392 392 0.056
## 393 393 0.056
## 394 394 0.064
## 395 395 0.072
## 396 396 0.072
## 397 397 0.080
## 398 398 0.072
## 399 399 0.072
## 400 400 0.072
## 401 401 0.072
## 402 402 0.080
## 403 403 0.080
## 404 404 0.064
## 405 405 0.064
## 406 406 0.064
## 407 407 0.064
## 408 408 0.056
## 409 409 0.064
## 410 410 0.064
## 411 411 0.056
## 412 412 0.056
## 413 413 0.048
## 414 414 0.048
## 415 415 0.040
## 416 416 0.040
## 417 417 0.032
## 418 418 0.040
## 419 419 0.040
## 420 420 0.032
## 421 421 0.032
## 422 422 0.032
## 423 423 0.040
## 424 424 0.032
## 425 425 0.040
## 426 426 0.032
## 427 427 0.016
## 428 428 0.016
## 429 429 0.016
## 430 430 0.016
## 431 431 0.016
## 432 432 0.008
## 433 433 0.008
## 434 434 0.008
## 435 435 0.008
## 436 436 0.016
## 437 437 0.016
## 438 438 0.016
## 439 439 0.016
## 440 440 0.008
## 441 441 0.008
## 442 442 0.008
## 443 443 0.008
## 444 444 0.016
## 445 445 0.008
## 446 446 0.000
## 447 447 0.000
## 448 448 0.000
## 449 449 0.000
## 450 450 0.008
## 451 451 0.008
## 452 452 0.000
## 453 453 0.000
## 454 454 0.016
## 455 455 0.016
## 456 456 0.008
## 457 457 0.008
## 458 458 0.008
## 459 459 0.008
## 460 460 0.008
## 461 461 0.008
## 462 462 0.008
## 463 463 0.000
## 464 464 0.000
## 465 465 0.000
## 466 466 0.008
## 467 467 0.008
## 468 468 0.008
## 469 469 0.008
## 470 470 0.016
## 471 471 0.008
## 472 472 0.008
## 473 473 0.008
## 474 474 0.008
## 475 475 0.016
## 476 476 0.016
## 477 477 0.016
## 478 478 0.024
## 479 479 0.016
## 480 480 0.016
## 481 481 0.024
## 482 482 0.032
## 483 483 0.032
## 484 484 0.024
## 485 485 0.024
## 486 486 0.024
## 487 487 0.024
## 488 488 0.024
## 489 489 0.024
## 490 490 0.032
## 491 491 0.032
## 492 492 0.032
## 493 493 0.032
## 494 494 0.032
## 495 495 0.032
## 496 496 0.032
## 497 497 0.024
## 498 498 0.024
## 499 499 0.032
## 500 500 0.032
## 501 501 0.032
## 502 502 0.032
## 503 503 0.032
## 504 504 0.032
## 505 505 0.032
## 506 506 0.040
## 507 507 0.032
## 508 508 0.024
## 509 509 0.032
## 510 510 0.024
## 511 511 0.032
## 512 512 0.032
## 513 513 0.040
## 514 514 0.048
## 515 515 0.048
## 516 516 0.040
## 517 517 0.040
## 518 518 0.040
## 519 519 0.040
## 520 520 0.040
## 521 521 0.040
## 522 522 0.040
## 523 523 0.048
## 524 524 0.048
## 525 525 0.048
## 526 526 0.048
## 527 527 0.048
## 528 528 0.048
## 529 529 0.048
## 530 530 0.048
## 531 531 0.040
## 532 532 0.048
## 533 533 0.056
## 534 534 0.056
## 535 535 0.056
## 536 536 0.064
## 537 537 0.064
## 538 538 0.056
## 539 539 0.064
## 540 540 0.056
## 541 541 0.056
## 542 542 0.056
## 543 543 0.048
## 544 544 0.048
## 545 545 0.048
## 546 546 0.048
## 547 547 0.048
## 548 548 0.048
## 549 549 0.048
## 550 550 0.048
## 551 551 0.048
## 552 552 0.048
## 553 553 0.056
## 554 554 0.056
## 555 555 0.056
## 556 556 0.056
## 557 557 0.064
## 558 558 0.064
## 559 559 0.072
## 560 560 0.080
## 561 561 0.080
## 562 562 0.080
## 563 563 0.088
## 564 564 0.096
## 565 565 0.088
## 566 566 0.088
## 567 567 0.088
## 568 568 0.088
## 569 569 0.080
## 570 570 0.080
## 571 571 0.072
## 572 572 0.056
## 573 573 0.056
## 574 574 0.056
## 575 575 0.056
## 576 576 0.056
## 577 577 0.056
## 578 578 0.056
## 579 579 0.056
## 580 580 0.056
## 581 581 0.056
## 582 582 0.048
## 583 583 0.056
## 584 584 0.048
## 585 585 0.048
## 586 586 0.056
## 587 587 0.064
## 588 588 0.056
## 589 589 0.040
## 590 590 0.040
## 591 591 0.048
## 592 592 0.056
## 593 593 0.064
## 594 594 0.064
## 595 595 0.056
## 596 596 0.064
## 597 597 0.064
## 598 598 0.064
## 599 599 0.064
## 600 600 0.064
## 601 601 0.064
## 602 602 0.064
## 603 603 0.064
## 604 604 0.064
## 605 605 0.056
## 606 606 0.056
## 607 607 0.072
## 608 608 0.056
## 609 609 0.056
## 610 610 0.056
## 611 611 0.056
## 612 612 0.048
## 613 613 0.040
## 614 614 0.040
## 615 615 0.040
## 616 616 0.056
## 617 617 0.048
## 618 618 0.040
## 619 619 0.032
## 620 620 0.032
## 621 621 0.032
## 622 622 0.040
## 623 623 0.040
## 624 624 0.040
## 625 625 0.032
## 626 626 0.032
## 627 627 0.032
## 628 628 0.032
## 629 629 0.032
## 630 630 0.032
## 631 631 0.040
## 632 632 0.040
## 633 633 0.032
## 634 634 0.032
## 635 635 0.040
## 636 636 0.040
## 637 637 0.040
## 638 638 0.040
## 639 639 0.040
## 640 640 0.032
## 641 641 0.032
## 642 642 0.032
## 643 643 0.032
## 644 644 0.032
## 645 645 0.032
## 646 646 0.040
## 647 647 0.040
## 648 648 0.040
## 649 649 0.040
## 650 650 0.040
## 651 651 0.040
## 652 652 0.040
## 653 653 0.048
## 654 654 0.048
## 655 655 0.048
## 656 656 0.048
## 657 657 0.048
## 658 658 0.056
## 659 659 0.056
## 660 660 0.056
## 661 661 0.056
## 662 662 0.056
## 663 663 0.056
## 664 664 0.064
## 665 665 0.064
## 666 666 0.064
## 667 667 0.056
## 668 668 0.056
## 669 669 0.056
## 670 670 0.056
## 671 671 0.056
## 672 672 0.056
## 673 673 0.056
## 674 674 0.064
## 675 675 0.064
## 676 676 0.064
## 677 677 0.064
## 678 678 0.064
## 679 679 0.064
## 680 680 0.064
## 681 681 0.072
## 682 682 0.064
## 683 683 0.064
## 684 684 0.064
## 685 685 0.064
## 686 686 0.072
## 687 687 0.080
## 688 688 0.080
## 689 689 0.080
## 690 690 0.088
## 691 691 0.088
## 692 692 0.088
## 693 693 0.088
## 694 694 0.104
## 695 695 0.104
## 696 696 0.112
## 697 697 0.112
## 698 698 0.112
## 699 699 0.112
## 700 700 0.112
## 701 701 0.112
## 702 702 0.120
## 703 703 0.128
## 704 704 0.128
## 705 705 0.128
## 706 706 0.136
## 707 707 0.144
## 708 708 0.144
## 709 709 0.144
## 710 710 0.144
## 711 711 0.152
## 712 712 0.152
## 713 713 0.152
## 714 714 0.160
## 715 715 0.168
## 716 716 0.168
## 717 717 0.176
## 718 718 0.176
## 719 719 0.176
## 720 720 0.176
## 721 721 0.176
## 722 722 0.176
## 723 723 0.176
## 724 724 0.176
## 725 725 0.184
## 726 726 0.184
## 727 727 0.184
## 728 728 0.184
## 729 729 0.184
## 730 730 0.184
## 731 731 0.184
## 732 732 0.184
## 733 733 0.184
## 734 734 0.184
## 735 735 0.184
## 736 736 0.184
## 737 737 0.184
## 738 738 0.192
## 739 739 0.192
## 740 740 0.192
## 741 741 0.192
## 742 742 0.192
## 743 743 0.192
## 744 744 0.192
## 745 745 0.192
## 746 746 0.192
## 747 747 0.192
## 748 748 0.192
## 749 749 0.192
## 750 750 0.192
## 751 751 0.192
## 752 752 0.200
## 753 753 0.200
## 754 754 0.200
## 755 755 0.200
## 756 756 0.200
## 757 757 0.200
## 758 758 0.200
## 759 759 0.208
## 760 760 0.208
## 761 761 0.208
## 762 762 0.208
## 763 763 0.208
## 764 764 0.208
## 765 765 0.208
## 766 766 0.208
## 767 767 0.208
## 768 768 0.216
## 769 769 0.216
## 770 770 0.224
## 771 771 0.224
## 772 772 0.232
## 773 773 0.232
## 774 774 0.232
## 775 775 0.240
## 776 776 0.240
## 777 777 0.240
## 778 778 0.240
## 779 779 0.240
## 780 780 0.240
## 781 781 0.240
## 782 782 0.248
## 783 783 0.248
## 784 784 0.248
## 785 785 0.256
## 786 786 0.256
## 787 787 0.256
## 788 788 0.256
## 789 789 0.256
## 790 790 0.256
## 791 791 0.256
## 792 792 0.256
## 793 793 0.256
## 794 794 0.256
## 795 795 0.256
## 796 796 0.256
## 797 797 0.256
## 798 798 0.272
## 799 799 0.272
## 800 800 0.288
## 801 801 0.288
## 802 802 0.296
## 803 803 0.296
## 804 804 0.296
## 805 805 0.296
## 806 806 0.296
## 807 807 0.304
## 808 808 0.304
## 809 809 0.312
## 810 810 0.320
## 811 811 0.320
## 812 812 0.320
## 813 813 0.328
## 814 814 0.328
## 815 815 0.328
## 816 816 0.328
## 817 817 0.336
## 818 818 0.344
## 819 819 0.344
## 820 820 0.344
## 821 821 0.344
## 822 822 0.344
## 823 823 0.344
## 824 824 0.344
## 825 825 0.344
## 826 826 0.344
## 827 827 0.344
## 828 828 0.344
## 829 829 0.344
## 830 830 0.344
## 831 831 0.344
## 832 832 0.344
## 833 833 0.344
## 834 834 0.352
## 835 835 0.360
## 836 836 0.360
## 837 837 0.360
## 838 838 0.360
## 839 839 0.360
## 840 840 0.360
## 841 841 0.360
## 842 842 0.368
## 843 843 0.368
## 844 844 0.368
## 845 845 0.368
## 846 846 0.368
## 847 847 0.376
## 848 848 0.376
## 849 849 0.376
## 850 850 0.376
## 851 851 0.376
## 852 852 0.384
## 853 853 0.392
## 854 854 0.392
## 855 855 0.392
## 856 856 0.392
## 857 857 0.392
## 858 858 0.400
## 859 859 0.400
## 860 860 0.400
## 861 861 0.400
## 862 862 0.400
## 863 863 0.400
## 864 864 0.400
## 865 865 0.400
## 866 866 0.400
## 867 867 0.400
## 868 868 0.400
## 869 869 0.408
## 870 870 0.416
## 871 871 0.424
## 872 872 0.424
## 873 873 0.432
## 874 874 0.440
## 875 875 0.440
## 876 876 0.440
## 877 877 0.448
## 878 878 0.448
## 879 879 0.448
## 880 880 0.448
## 881 881 0.448
## 882 882 0.448
## 883 883 0.448
## 884 884 0.448
## 885 885 0.448
## 886 886 0.456
## 887 887 0.456
## 888 888 0.456
## 889 889 0.456
## 890 890 0.456
## 891 891 0.456
## 892 892 0.456
## 893 893 0.456
## 894 894 0.472
## 895 895 0.472
## 896 896 0.472
## 897 897 0.472
## 898 898 0.472
## 899 899 0.472
## 900 900 0.472
## 901 901 0.472
## 902 902 0.472
## 903 903 0.472
## 904 904 0.472
## 905 905 0.472
## 906 906 0.472
## 907 907 0.472
## 908 908 0.472
## 909 909 0.472
## 910 910 0.472
## 911 911 0.480
## 912 912 0.480
## 913 913 0.480
## 914 914 0.488
## 915 915 0.488
## 916 916 0.488
## 917 917 0.488
## 918 918 0.488
## 919 919 0.488
## 920 920 0.488
## 921 921 0.488
## 922 922 0.488
## 923 923 0.496
## 924 924 0.504
## 925 925 0.504
## 926 926 0.504
## 927 927 0.504
## 928 928 0.504
## 929 929 0.512
## 930 930 0.512
## 931 931 0.512
## 932 932 0.512
## 933 933 0.512
## 934 934 0.512
## 935 935 0.512
## 936 936 0.520
## 937 937 0.520
## 938 938 0.520
## 939 939 0.520
## 940 940 0.520
## 941 941 0.520
## 942 942 0.520
## 943 943 0.520
## 944 944 0.520
## 945 945 0.520
## 946 946 0.520
## 947 947 0.520
## 948 948 0.520
## 949 949 0.520
## 950 950 0.520
## 951 951 0.520
## 952 952 0.520
## 953 953 0.520
## 954 954 0.520
## 955 955 0.520
## 956 956 0.520
## 957 957 0.520
## 958 958 0.520
## 959 959 0.520
## 960 960 0.520
## 961 961 0.520
## 962 962 0.520
## 963 963 0.520
## 964 964 0.520
## 965 965 0.520
## 966 966 0.520
## 967 967 0.520
## 968 968 0.520
## 969 969 0.520
## 970 970 0.528
## 971 971 0.528
## 972 972 0.528
## 973 973 0.528
## 974 974 0.528
## 975 975 0.528
## 976 976 0.528
## 977 977 0.528
## 978 978 0.536
## 979 979 0.536
## 980 980 0.536
## 981 981 0.544
## 982 982 0.544
## 983 983 0.544
## 984 984 0.544
## 985 985 0.544
## 986 986 0.544
## 987 987 0.544
## 988 988 0.544
## 989 989 0.544
## 990 990 0.544
## 991 991 0.552
## 992 992 0.552
## 993 993 0.552
## 994 994 0.560
## 995 995 0.568
## 996 996 0.568
## 997 997 0.568
## 998 998 0.568
## 999 999 0.568
## 1000 1000 0.568# Each node
prevalence(result3, I~., level=3)## node time prevalence
## 1 1 1 0.040
## 2 1 2 0.024
## 3 1 3 0.024
## 4 1 4 0.024
## 5 1 5 0.024
## 6 1 6 0.024
## 7 1 7 0.032
## 8 1 8 0.032
## 9 1 9 0.024
## 10 1 10 0.016
## 11 1 11 0.016
## 12 1 12 0.024
## 13 1 13 0.024
## 14 1 14 0.016
## 15 1 15 0.016
## 16 1 16 0.016
## 17 1 17 0.016
## 18 1 18 0.016
## 19 1 19 0.016
## 20 1 20 0.016
## 21 1 21 0.016
## 22 1 22 0.016
## 23 1 23 0.008
## 24 1 24 0.008
## 25 1 25 0.008
## 26 1 26 0.008
## 27 1 27 0.008
## 28 1 28 0.008
## 29 1 29 0.008
## 30 1 30 0.008
## 31 1 31 0.008
## 32 1 32 0.016
## 33 1 33 0.016
## 34 1 34 0.024
## 35 1 35 0.024
## 36 1 36 0.024
## 37 1 37 0.024
## 38 1 38 0.008
## 39 1 39 0.008
## 40 1 40 0.008
## 41 1 41 0.000
## 42 1 42 0.000
## 43 1 43 0.000
## 44 1 44 0.000
## 45 1 45 0.000
## 46 1 46 0.000
## 47 1 47 0.000
## 48 1 48 0.000
## 49 1 49 0.000
## 50 1 50 0.008
## 51 1 51 0.016
## 52 1 52 0.016
## 53 1 53 0.016
## 54 1 54 0.016
## 55 1 55 0.024
## 56 1 56 0.024
## 57 1 57 0.024
## 58 1 58 0.024
## 59 1 59 0.024
## 60 1 60 0.024
## 61 1 61 0.024
## 62 1 62 0.032
## 63 1 63 0.024
## 64 1 64 0.024
## 65 1 65 0.024
## 66 1 66 0.024
## 67 1 67 0.032
## 68 1 68 0.024
## 69 1 69 0.024
## 70 1 70 0.024
## 71 1 71 0.024
## 72 1 72 0.024
## 73 1 73 0.024
## 74 1 74 0.016
## 75 1 75 0.016
## 76 1 76 0.016
## 77 1 77 0.016
## 78 1 78 0.024
## 79 1 79 0.024
## 80 1 80 0.024
## 81 1 81 0.024
## 82 1 82 0.024
## 83 1 83 0.024
## 84 1 84 0.016
## 85 1 85 0.016
## 86 1 86 0.016
## 87 1 87 0.032
## 88 1 88 0.032
## 89 1 89 0.032
## 90 1 90 0.032
## 91 1 91 0.032
## 92 1 92 0.032
## 93 1 93 0.040
## 94 1 94 0.040
## 95 1 95 0.040
## 96 1 96 0.040
## 97 1 97 0.040
## 98 1 98 0.040
## 99 1 99 0.048
## 100 1 100 0.048
## 101 1 101 0.040
## 102 1 102 0.032
## 103 1 103 0.016
## 104 1 104 0.016
## 105 1 105 0.000
## 106 1 106 0.000
## 107 1 107 0.000
## 108 1 108 0.000
## 109 1 109 0.000
## 110 1 110 0.000
## 111 1 111 0.008
## 112 1 112 0.024
## 113 1 113 0.024
## 114 1 114 0.024
## 115 1 115 0.032
## 116 1 116 0.032
## 117 1 117 0.032
## 118 1 118 0.032
## 119 1 119 0.032
## 120 1 120 0.032
## 121 1 121 0.032
## 122 1 122 0.032
## 123 1 123 0.032
## 124 1 124 0.032
## 125 1 125 0.032
## 126 1 126 0.040
## 127 1 127 0.048
## 128 1 128 0.048
## 129 1 129 0.040
## 130 1 130 0.024
## 131 1 131 0.016
## 132 1 132 0.024
## 133 1 133 0.024
## 134 1 134 0.024
## 135 1 135 0.024
## 136 1 136 0.024
## 137 1 137 0.024
## 138 1 138 0.024
## 139 1 139 0.024
## 140 1 140 0.024
## 141 1 141 0.016
## 142 1 142 0.016
## 143 1 143 0.016
## 144 1 144 0.024
## 145 1 145 0.024
## 146 1 146 0.024
## 147 1 147 0.032
## 148 1 148 0.032
## 149 1 149 0.032
## 150 1 150 0.024
## 151 1 151 0.024
## 152 1 152 0.024
## 153 1 153 0.024
## 154 1 154 0.032
## 155 1 155 0.032
## 156 1 156 0.040
## 157 1 157 0.048
## 158 1 158 0.040
## 159 1 159 0.048
## 160 1 160 0.048
## 161 1 161 0.040
## 162 1 162 0.032
## 163 1 163 0.056
## 164 1 164 0.064
## 165 1 165 0.064
## 166 1 166 0.064
## 167 1 167 0.072
## 168 1 168 0.072
## 169 1 169 0.056
## 170 1 170 0.048
## 171 1 171 0.040
## 172 1 172 0.040
## 173 1 173 0.040
## 174 1 174 0.040
## 175 1 175 0.040
## 176 1 176 0.040
## 177 1 177 0.024
## 178 1 178 0.024
## 179 1 179 0.008
## 180 1 180 0.008
## 181 1 181 0.008
## 182 1 182 0.008
## 183 1 183 0.008
## 184 1 184 0.016
## 185 1 185 0.016
## 186 1 186 0.024
## 187 1 187 0.024
## 188 1 188 0.032
## 189 1 189 0.040
## 190 1 190 0.032
## 191 1 191 0.040
## 192 1 192 0.040
## 193 1 193 0.048
## 194 1 194 0.048
## 195 1 195 0.048
## 196 1 196 0.048
## 197 1 197 0.048
## 198 1 198 0.048
## 199 1 199 0.048
## 200 1 200 0.048
## 201 1 201 0.048
## 202 1 202 0.040
## 203 1 203 0.048
## 204 1 204 0.048
## 205 1 205 0.048
## 206 1 206 0.032
## 207 1 207 0.040
## 208 1 208 0.032
## 209 1 209 0.032
## 210 1 210 0.032
## 211 1 211 0.024
## 212 1 212 0.016
## 213 1 213 0.024
## 214 1 214 0.024
## 215 1 215 0.024
## 216 1 216 0.024
## 217 1 217 0.024
## 218 1 218 0.024
## 219 1 219 0.024
## 220 1 220 0.032
## 221 1 221 0.048
## 222 1 222 0.048
## 223 1 223 0.048
## 224 1 224 0.048
## 225 1 225 0.048
## 226 1 226 0.056
## 227 1 227 0.056
## 228 1 228 0.048
## 229 1 229 0.040
## 230 1 230 0.040
## 231 1 231 0.040
## 232 1 232 0.040
## 233 1 233 0.032
## 234 1 234 0.032
## 235 1 235 0.024
## 236 1 236 0.016
## 237 1 237 0.016
## 238 1 238 0.024
## 239 1 239 0.024
## 240 1 240 0.024
## 241 1 241 0.024
## 242 1 242 0.008
## 243 1 243 0.008
## 244 1 244 0.008
## 245 1 245 0.008
## 246 1 246 0.008
## 247 1 247 0.008
## 248 1 248 0.008
## 249 1 249 0.008
## 250 1 250 0.008
## 251 1 251 0.008
## 252 1 252 0.008
## 253 1 253 0.008
## 254 1 254 0.008
## 255 1 255 0.008
## 256 1 256 0.008
## 257 1 257 0.000
## 258 1 258 0.008
## 259 1 259 0.008
## 260 1 260 0.008
## 261 1 261 0.008
## 262 1 262 0.000
## 263 1 263 0.000
## 264 1 264 0.000
## 265 1 265 0.000
## 266 1 266 0.000
## 267 1 267 0.000
## 268 1 268 0.000
## 269 1 269 0.000
## 270 1 270 0.000
## 271 1 271 0.000
## 272 1 272 0.000
## 273 1 273 0.000
## 274 1 274 0.000
## 275 1 275 0.000
## 276 1 276 0.008
## 277 1 277 0.008
## 278 1 278 0.008
## 279 1 279 0.008
## 280 1 280 0.008
## 281 1 281 0.008
## 282 1 282 0.008
## 283 1 283 0.008
## 284 1 284 0.016
## 285 1 285 0.016
## 286 1 286 0.024
## 287 1 287 0.024
## 288 1 288 0.024
## 289 1 289 0.024
## 290 1 290 0.024
## 291 1 291 0.024
## 292 1 292 0.024
## 293 1 293 0.024
## 294 1 294 0.024
## 295 1 295 0.024
## 296 1 296 0.024
## 297 1 297 0.016
## 298 1 298 0.016
## 299 1 299 0.016
## 300 1 300 0.016
## 301 1 301 0.008
## 302 1 302 0.016
## 303 1 303 0.016
## 304 1 304 0.016
## 305 1 305 0.024
## 306 1 306 0.024
## 307 1 307 0.024
## 308 1 308 0.024
## 309 1 309 0.024
## 310 1 310 0.032
## 311 1 311 0.032
## 312 1 312 0.032
## 313 1 313 0.040
## 314 1 314 0.040
## 315 1 315 0.040
## 316 1 316 0.040
## 317 1 317 0.048
## 318 1 318 0.048
## 319 1 319 0.056
## 320 1 320 0.064
## 321 1 321 0.064
## 322 1 322 0.064
## 323 1 323 0.064
## 324 1 324 0.072
## 325 1 325 0.088
## 326 1 326 0.088
## 327 1 327 0.096
## 328 1 328 0.104
## 329 1 329 0.104
## 330 1 330 0.104
## 331 1 331 0.104
## 332 1 332 0.104
## 333 1 333 0.104
## 334 1 334 0.104
## 335 1 335 0.104
## 336 1 336 0.096
## 337 1 337 0.096
## 338 1 338 0.096
## 339 1 339 0.096
## 340 1 340 0.104
## 341 1 341 0.104
## 342 1 342 0.104
## 343 1 343 0.104
## 344 1 344 0.104
## 345 1 345 0.104
## 346 1 346 0.104
## 347 1 347 0.104
## 348 1 348 0.112
## 349 1 349 0.104
## 350 1 350 0.096
## 351 1 351 0.096
## 352 1 352 0.104
## 353 1 353 0.104
## 354 1 354 0.104
## 355 1 355 0.112
## 356 1 356 0.112
## 357 1 357 0.112
## 358 1 358 0.112
## 359 1 359 0.104
## 360 1 360 0.104
## 361 1 361 0.104
## 362 1 362 0.104
## 363 1 363 0.096
## 364 1 364 0.096
## 365 1 365 0.096
## 366 1 366 0.096
## 367 1 367 0.096
## 368 1 368 0.096
## 369 1 369 0.096
## 370 1 370 0.096
## 371 1 371 0.080
## 372 1 372 0.088
## 373 1 373 0.088
## 374 1 374 0.080
## 375 1 375 0.080
## 376 1 376 0.080
## 377 1 377 0.088
## 378 1 378 0.072
## 379 1 379 0.080
## 380 1 380 0.080
## 381 1 381 0.072
## 382 1 382 0.064
## 383 1 383 0.072
## 384 1 384 0.072
## 385 1 385 0.072
## 386 1 386 0.064
## 387 1 387 0.064
## 388 1 388 0.056
## 389 1 389 0.056
## 390 1 390 0.056
## 391 1 391 0.056
## 392 1 392 0.056
## 393 1 393 0.056
## 394 1 394 0.064
## 395 1 395 0.072
## 396 1 396 0.072
## 397 1 397 0.080
## 398 1 398 0.072
## 399 1 399 0.072
## 400 1 400 0.072
## 401 1 401 0.072
## 402 1 402 0.080
## 403 1 403 0.080
## 404 1 404 0.064
## 405 1 405 0.064
## 406 1 406 0.064
## 407 1 407 0.064
## 408 1 408 0.056
## 409 1 409 0.064
## 410 1 410 0.064
## 411 1 411 0.056
## 412 1 412 0.056
## 413 1 413 0.048
## 414 1 414 0.048
## 415 1 415 0.040
## 416 1 416 0.040
## 417 1 417 0.032
## 418 1 418 0.040
## 419 1 419 0.040
## 420 1 420 0.032
## 421 1 421 0.032
## 422 1 422 0.032
## 423 1 423 0.040
## 424 1 424 0.032
## 425 1 425 0.040
## 426 1 426 0.032
## 427 1 427 0.016
## 428 1 428 0.016
## 429 1 429 0.016
## 430 1 430 0.016
## 431 1 431 0.016
## 432 1 432 0.008
## 433 1 433 0.008
## 434 1 434 0.008
## 435 1 435 0.008
## 436 1 436 0.016
## 437 1 437 0.016
## 438 1 438 0.016
## 439 1 439 0.016
## 440 1 440 0.008
## 441 1 441 0.008
## 442 1 442 0.008
## 443 1 443 0.008
## 444 1 444 0.016
## 445 1 445 0.008
## 446 1 446 0.000
## 447 1 447 0.000
## 448 1 448 0.000
## 449 1 449 0.000
## 450 1 450 0.008
## 451 1 451 0.008
## 452 1 452 0.000
## 453 1 453 0.000
## 454 1 454 0.016
## 455 1 455 0.016
## 456 1 456 0.008
## 457 1 457 0.008
## 458 1 458 0.008
## 459 1 459 0.008
## 460 1 460 0.008
## 461 1 461 0.008
## 462 1 462 0.008
## 463 1 463 0.000
## 464 1 464 0.000
## 465 1 465 0.000
## 466 1 466 0.008
## 467 1 467 0.008
## 468 1 468 0.008
## 469 1 469 0.008
## 470 1 470 0.016
## 471 1 471 0.008
## 472 1 472 0.008
## 473 1 473 0.008
## 474 1 474 0.008
## 475 1 475 0.016
## 476 1 476 0.016
## 477 1 477 0.016
## 478 1 478 0.024
## 479 1 479 0.016
## 480 1 480 0.016
## 481 1 481 0.024
## 482 1 482 0.032
## 483 1 483 0.032
## 484 1 484 0.024
## 485 1 485 0.024
## 486 1 486 0.024
## 487 1 487 0.024
## 488 1 488 0.024
## 489 1 489 0.024
## 490 1 490 0.032
## 491 1 491 0.032
## 492 1 492 0.032
## 493 1 493 0.032
## 494 1 494 0.032
## 495 1 495 0.032
## 496 1 496 0.032
## 497 1 497 0.024
## 498 1 498 0.024
## 499 1 499 0.032
## 500 1 500 0.032
## 501 1 501 0.032
## 502 1 502 0.032
## 503 1 503 0.032
## 504 1 504 0.032
## 505 1 505 0.032
## 506 1 506 0.040
## 507 1 507 0.032
## 508 1 508 0.024
## 509 1 509 0.032
## 510 1 510 0.024
## 511 1 511 0.032
## 512 1 512 0.032
## 513 1 513 0.040
## 514 1 514 0.048
## 515 1 515 0.048
## 516 1 516 0.040
## 517 1 517 0.040
## 518 1 518 0.040
## 519 1 519 0.040
## 520 1 520 0.040
## 521 1 521 0.040
## 522 1 522 0.040
## 523 1 523 0.048
## 524 1 524 0.048
## 525 1 525 0.048
## 526 1 526 0.048
## 527 1 527 0.048
## 528 1 528 0.048
## 529 1 529 0.048
## 530 1 530 0.048
## 531 1 531 0.040
## 532 1 532 0.048
## 533 1 533 0.056
## 534 1 534 0.056
## 535 1 535 0.056
## 536 1 536 0.064
## 537 1 537 0.064
## 538 1 538 0.056
## 539 1 539 0.064
## 540 1 540 0.056
## 541 1 541 0.056
## 542 1 542 0.056
## 543 1 543 0.048
## 544 1 544 0.048
## 545 1 545 0.048
## 546 1 546 0.048
## 547 1 547 0.048
## 548 1 548 0.048
## 549 1 549 0.048
## 550 1 550 0.048
## 551 1 551 0.048
## 552 1 552 0.048
## 553 1 553 0.056
## 554 1 554 0.056
## 555 1 555 0.056
## 556 1 556 0.056
## 557 1 557 0.064
## 558 1 558 0.064
## 559 1 559 0.072
## 560 1 560 0.080
## 561 1 561 0.080
## 562 1 562 0.080
## 563 1 563 0.088
## 564 1 564 0.096
## 565 1 565 0.088
## 566 1 566 0.088
## 567 1 567 0.088
## 568 1 568 0.088
## 569 1 569 0.080
## 570 1 570 0.080
## 571 1 571 0.072
## 572 1 572 0.056
## 573 1 573 0.056
## 574 1 574 0.056
## 575 1 575 0.056
## 576 1 576 0.056
## 577 1 577 0.056
## 578 1 578 0.056
## 579 1 579 0.056
## 580 1 580 0.056
## 581 1 581 0.056
## 582 1 582 0.048
## 583 1 583 0.056
## 584 1 584 0.048
## 585 1 585 0.048
## 586 1 586 0.056
## 587 1 587 0.064
## 588 1 588 0.056
## 589 1 589 0.040
## 590 1 590 0.040
## 591 1 591 0.048
## 592 1 592 0.056
## 593 1 593 0.064
## 594 1 594 0.064
## 595 1 595 0.056
## 596 1 596 0.064
## 597 1 597 0.064
## 598 1 598 0.064
## 599 1 599 0.064
## 600 1 600 0.064
## 601 1 601 0.064
## 602 1 602 0.064
## 603 1 603 0.064
## 604 1 604 0.064
## 605 1 605 0.056
## 606 1 606 0.056
## 607 1 607 0.072
## 608 1 608 0.056
## 609 1 609 0.056
## 610 1 610 0.056
## 611 1 611 0.056
## 612 1 612 0.048
## 613 1 613 0.040
## 614 1 614 0.040
## 615 1 615 0.040
## 616 1 616 0.056
## 617 1 617 0.048
## 618 1 618 0.040
## 619 1 619 0.032
## 620 1 620 0.032
## 621 1 621 0.032
## 622 1 622 0.040
## 623 1 623 0.040
## 624 1 624 0.040
## 625 1 625 0.032
## 626 1 626 0.032
## 627 1 627 0.032
## 628 1 628 0.032
## 629 1 629 0.032
## 630 1 630 0.032
## 631 1 631 0.040
## 632 1 632 0.040
## 633 1 633 0.032
## 634 1 634 0.032
## 635 1 635 0.040
## 636 1 636 0.040
## 637 1 637 0.040
## 638 1 638 0.040
## 639 1 639 0.040
## 640 1 640 0.032
## 641 1 641 0.032
## 642 1 642 0.032
## 643 1 643 0.032
## 644 1 644 0.032
## 645 1 645 0.032
## 646 1 646 0.040
## 647 1 647 0.040
## 648 1 648 0.040
## 649 1 649 0.040
## 650 1 650 0.040
## 651 1 651 0.040
## 652 1 652 0.040
## 653 1 653 0.048
## 654 1 654 0.048
## 655 1 655 0.048
## 656 1 656 0.048
## 657 1 657 0.048
## 658 1 658 0.056
## 659 1 659 0.056
## 660 1 660 0.056
## 661 1 661 0.056
## 662 1 662 0.056
## 663 1 663 0.056
## 664 1 664 0.064
## 665 1 665 0.064
## 666 1 666 0.064
## 667 1 667 0.056
## 668 1 668 0.056
## 669 1 669 0.056
## 670 1 670 0.056
## 671 1 671 0.056
## 672 1 672 0.056
## 673 1 673 0.056
## 674 1 674 0.064
## 675 1 675 0.064
## 676 1 676 0.064
## 677 1 677 0.064
## 678 1 678 0.064
## 679 1 679 0.064
## 680 1 680 0.064
## 681 1 681 0.072
## 682 1 682 0.064
## 683 1 683 0.064
## 684 1 684 0.064
## 685 1 685 0.064
## 686 1 686 0.072
## 687 1 687 0.080
## 688 1 688 0.080
## 689 1 689 0.080
## 690 1 690 0.088
## 691 1 691 0.088
## 692 1 692 0.088
## 693 1 693 0.088
## 694 1 694 0.104
## 695 1 695 0.104
## 696 1 696 0.112
## 697 1 697 0.112
## 698 1 698 0.112
## 699 1 699 0.112
## 700 1 700 0.112
## 701 1 701 0.112
## 702 1 702 0.120
## 703 1 703 0.128
## 704 1 704 0.128
## 705 1 705 0.128
## 706 1 706 0.136
## 707 1 707 0.144
## 708 1 708 0.144
## 709 1 709 0.144
## 710 1 710 0.144
## 711 1 711 0.152
## 712 1 712 0.152
## 713 1 713 0.152
## 714 1 714 0.160
## 715 1 715 0.168
## 716 1 716 0.168
## 717 1 717 0.176
## 718 1 718 0.176
## 719 1 719 0.176
## 720 1 720 0.176
## 721 1 721 0.176
## 722 1 722 0.176
## 723 1 723 0.176
## 724 1 724 0.176
## 725 1 725 0.184
## 726 1 726 0.184
## 727 1 727 0.184
## 728 1 728 0.184
## 729 1 729 0.184
## 730 1 730 0.184
## 731 1 731 0.184
## 732 1 732 0.184
## 733 1 733 0.184
## 734 1 734 0.184
## 735 1 735 0.184
## 736 1 736 0.184
## 737 1 737 0.184
## 738 1 738 0.192
## 739 1 739 0.192
## 740 1 740 0.192
## 741 1 741 0.192
## 742 1 742 0.192
## 743 1 743 0.192
## 744 1 744 0.192
## 745 1 745 0.192
## 746 1 746 0.192
## 747 1 747 0.192
## 748 1 748 0.192
## 749 1 749 0.192
## 750 1 750 0.192
## 751 1 751 0.192
## 752 1 752 0.200
## 753 1 753 0.200
## 754 1 754 0.200
## 755 1 755 0.200
## 756 1 756 0.200
## 757 1 757 0.200
## 758 1 758 0.200
## 759 1 759 0.208
## 760 1 760 0.208
## 761 1 761 0.208
## 762 1 762 0.208
## 763 1 763 0.208
## 764 1 764 0.208
## 765 1 765 0.208
## 766 1 766 0.208
## 767 1 767 0.208
## 768 1 768 0.216
## 769 1 769 0.216
## 770 1 770 0.224
## 771 1 771 0.224
## 772 1 772 0.232
## 773 1 773 0.232
## 774 1 774 0.232
## 775 1 775 0.240
## 776 1 776 0.240
## 777 1 777 0.240
## 778 1 778 0.240
## 779 1 779 0.240
## 780 1 780 0.240
## 781 1 781 0.240
## 782 1 782 0.248
## 783 1 783 0.248
## 784 1 784 0.248
## 785 1 785 0.256
## 786 1 786 0.256
## 787 1 787 0.256
## 788 1 788 0.256
## 789 1 789 0.256
## 790 1 790 0.256
## 791 1 791 0.256
## 792 1 792 0.256
## 793 1 793 0.256
## 794 1 794 0.256
## 795 1 795 0.256
## 796 1 796 0.256
## 797 1 797 0.256
## 798 1 798 0.272
## 799 1 799 0.272
## 800 1 800 0.288
## 801 1 801 0.288
## 802 1 802 0.296
## 803 1 803 0.296
## 804 1 804 0.296
## 805 1 805 0.296
## 806 1 806 0.296
## 807 1 807 0.304
## 808 1 808 0.304
## 809 1 809 0.312
## 810 1 810 0.320
## 811 1 811 0.320
## 812 1 812 0.320
## 813 1 813 0.328
## 814 1 814 0.328
## 815 1 815 0.328
## 816 1 816 0.328
## 817 1 817 0.336
## 818 1 818 0.344
## 819 1 819 0.344
## 820 1 820 0.344
## 821 1 821 0.344
## 822 1 822 0.344
## 823 1 823 0.344
## 824 1 824 0.344
## 825 1 825 0.344
## 826 1 826 0.344
## 827 1 827 0.344
## 828 1 828 0.344
## 829 1 829 0.344
## 830 1 830 0.344
## 831 1 831 0.344
## 832 1 832 0.344
## 833 1 833 0.344
## 834 1 834 0.352
## 835 1 835 0.360
## 836 1 836 0.360
## 837 1 837 0.360
## 838 1 838 0.360
## 839 1 839 0.360
## 840 1 840 0.360
## 841 1 841 0.360
## 842 1 842 0.368
## 843 1 843 0.368
## 844 1 844 0.368
## 845 1 845 0.368
## 846 1 846 0.368
## 847 1 847 0.376
## 848 1 848 0.376
## 849 1 849 0.376
## 850 1 850 0.376
## 851 1 851 0.376
## 852 1 852 0.384
## 853 1 853 0.392
## 854 1 854 0.392
## 855 1 855 0.392
## 856 1 856 0.392
## 857 1 857 0.392
## 858 1 858 0.400
## 859 1 859 0.400
## 860 1 860 0.400
## 861 1 861 0.400
## 862 1 862 0.400
## 863 1 863 0.400
## 864 1 864 0.400
## 865 1 865 0.400
## 866 1 866 0.400
## 867 1 867 0.400
## 868 1 868 0.400
## 869 1 869 0.408
## 870 1 870 0.416
## 871 1 871 0.424
## 872 1 872 0.424
## 873 1 873 0.432
## 874 1 874 0.440
## 875 1 875 0.440
## 876 1 876 0.440
## 877 1 877 0.448
## 878 1 878 0.448
## 879 1 879 0.448
## 880 1 880 0.448
## 881 1 881 0.448
## 882 1 882 0.448
## 883 1 883 0.448
## 884 1 884 0.448
## 885 1 885 0.448
## 886 1 886 0.456
## 887 1 887 0.456
## 888 1 888 0.456
## 889 1 889 0.456
## 890 1 890 0.456
## 891 1 891 0.456
## 892 1 892 0.456
## 893 1 893 0.456
## 894 1 894 0.472
## 895 1 895 0.472
## 896 1 896 0.472
## 897 1 897 0.472
## 898 1 898 0.472
## 899 1 899 0.472
## 900 1 900 0.472
## 901 1 901 0.472
## 902 1 902 0.472
## 903 1 903 0.472
## 904 1 904 0.472
## 905 1 905 0.472
## 906 1 906 0.472
## 907 1 907 0.472
## 908 1 908 0.472
## 909 1 909 0.472
## 910 1 910 0.472
## 911 1 911 0.480
## 912 1 912 0.480
## 913 1 913 0.480
## 914 1 914 0.488
## 915 1 915 0.488
## 916 1 916 0.488
## 917 1 917 0.488
## 918 1 918 0.488
## 919 1 919 0.488
## 920 1 920 0.488
## 921 1 921 0.488
## 922 1 922 0.488
## 923 1 923 0.496
## 924 1 924 0.504
## 925 1 925 0.504
## 926 1 926 0.504
## 927 1 927 0.504
## 928 1 928 0.504
## 929 1 929 0.512
## 930 1 930 0.512
## 931 1 931 0.512
## 932 1 932 0.512
## 933 1 933 0.512
## 934 1 934 0.512
## 935 1 935 0.512
## 936 1 936 0.520
## 937 1 937 0.520
## 938 1 938 0.520
## 939 1 939 0.520
## 940 1 940 0.520
## 941 1 941 0.520
## 942 1 942 0.520
## 943 1 943 0.520
## 944 1 944 0.520
## 945 1 945 0.520
## 946 1 946 0.520
## 947 1 947 0.520
## 948 1 948 0.520
## 949 1 949 0.520
## 950 1 950 0.520
## 951 1 951 0.520
## 952 1 952 0.520
## 953 1 953 0.520
## 954 1 954 0.520
## 955 1 955 0.520
## 956 1 956 0.520
## 957 1 957 0.520
## 958 1 958 0.520
## 959 1 959 0.520
## 960 1 960 0.520
## 961 1 961 0.520
## 962 1 962 0.520
## 963 1 963 0.520
## 964 1 964 0.520
## 965 1 965 0.520
## 966 1 966 0.520
## 967 1 967 0.520
## 968 1 968 0.520
## 969 1 969 0.520
## 970 1 970 0.528
## 971 1 971 0.528
## 972 1 972 0.528
## 973 1 973 0.528
## 974 1 974 0.528
## 975 1 975 0.528
## 976 1 976 0.528
## 977 1 977 0.528
## 978 1 978 0.536
## 979 1 979 0.536
## 980 1 980 0.536
## 981 1 981 0.544
## 982 1 982 0.544
## 983 1 983 0.544
## 984 1 984 0.544
## 985 1 985 0.544
## 986 1 986 0.544
## 987 1 987 0.544
## 988 1 988 0.544
## 989 1 989 0.544
## 990 1 990 0.544
## 991 1 991 0.552
## 992 1 992 0.552
## 993 1 993 0.552
## 994 1 994 0.560
## 995 1 995 0.568
## 996 1 996 0.568
## 997 1 997 0.568
## 998 1 998 0.568
## 999 1 999 0.568
## 1000 1 1000 0.5683.4 Plot compartments
plot(result3, range=0.95, type = "p", cex=0.051)
3.5 Plot prevalence
# plot(result3, index= 1:5, range=0.95)
plot(result3, I~.)
4 SIV model translation from dif eq
4.1 Parameters and transitions
# Rate of change
# dS <- (1 - fc(t)*e) * mu * (S + I + V) - (beta * S * I) / (S + I + V) + w * mu * V - mu * S
# dV <- fc(t)*e * mu * (S + I + V) - w * mu * V - mu * V
# dI <- (beta * S * I) / (S + I + V) - gama * I - mu * I
transitions <- c(
"@ -> (1-c)*e*mu*(S+I+V)-> S",
"S -> (beta*S*I)/(S+I+V)-> I",
"V -> w*mu*V -> S",
"S -> mu*S -> @",
"@ -> c*e*mu*(S+I+V)-> V",
"V -> w*mu*V-> S",
"V -> mu*V -> @",
"S -> (beta*S*I)/(S+I+V)-> I",
"I -> gama*I -> @",
"I -> mu*I -> @")
compartments <- c("S", "V", "I")
gdata <- c(beta = 4.0706,
gama= 2.7941,
e = 0.95,
mu = (1/(350/30)),
w = 1/(180/30),
c = 0.84)
u0 <- data.frame(S = 1418861, I = 1604, V = 579535)
n <- 2000000
# Initial state of the population compartments
# u0 <- data.frame(S=rep(999,n), I=rep(1,n), R=rep(0,n))
# time
tspan <- seq(1,180,1)4.2 Define the model
SIV <- mparse(transitions = transitions,
compartments = compartments,
u0 = u0,
gdata = gdata,
tspan = tspan)4.3 Run the SIV
set.seed(555)
result4 <- run(model=SIV, threads = 1)4.4 Plots of SIV
co <-c("#E69F00", "#009E73", "red")
result4## Model: SimInf_model
## Number of nodes: 1
## Number of transitions: 10
## Number of scheduled events: 0
##
## Global data
## -----------
## Parameter Value
## beta 4.07060000
## gama 2.79410000
## e 0.95000000
## mu 0.08571429
## w 0.16666667
## c 0.84000000
##
## Compartments
## ------------
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## S 129175 170384 203582 213492 240030 1299321
## V 234653 286188 345384 367673 420747 713061
## I 0 0 0 2309 0 171052plot(result4, col=co)
plot(result4, range= FALSE)
plot(prevalence(result4, I~.), cex=0.3, col="red")
trajectory(result4)## node time S V I
## 1 1 1 1299321 645088 23649
## 2 1 2 860523 696390 171052
## 3 1 3 261719 713061 169490
## 4 1 4 134315 698894 40956
## 5 1 5 129175 677605 8077
## 6 1 6 141692 657189 1717
## 7 1 7 156478 637785 452
## 8 1 8 170408 619490 125
## 9 1 9 182925 603121 42
## 10 1 10 193836 588479 30
## 11 1 11 203438 575377 25
## 12 1 12 211797 563531 11
## 13 1 13 219303 552919 15
## 14 1 14 225806 543459 10
## 15 1 15 231618 534616 1
## 16 1 16 236819 526670 0
## 17 1 17 241059 519103 0
## 18 1 18 244940 512099 0
## 19 1 19 248281 505250 0
## 20 1 20 251243 499075 0
## 21 1 21 253513 493362 0
## 22 1 22 255229 488530 0
## 23 1 23 257169 483028 0
## 24 1 24 259064 478262 0
## 25 1 25 260211 474354 0
## 26 1 26 261095 470819 0
## 27 1 27 261889 466865 0
## 28 1 28 262068 463568 0
## 29 1 29 262295 460512 0
## 30 1 30 262285 457481 0
## 31 1 31 262083 454454 0
## 32 1 32 261844 452092 0
## 33 1 33 261453 449394 0
## 34 1 34 260689 446579 0
## 35 1 35 260207 443930 0
## 36 1 36 259691 441419 0
## 37 1 37 258865 439210 0
## 38 1 38 257986 437031 0
## 39 1 39 257677 434928 0
## 40 1 40 256960 432808 0
## 41 1 41 256257 429865 0
## 42 1 42 255705 427883 0
## 43 1 43 254828 426001 0
## 44 1 44 253637 424497 0
## 45 1 45 253108 422665 0
## 46 1 46 252239 420108 0
## 47 1 47 251269 417828 0
## 48 1 48 250225 415306 0
## 49 1 49 249524 413136 0
## 50 1 50 248566 411412 0
## 51 1 51 247404 409196 0
## 52 1 52 246152 407355 0
## 53 1 53 245245 405595 0
## 54 1 54 244186 403654 0
## 55 1 55 243375 401909 0
## 56 1 56 242536 399761 0
## 57 1 57 241732 397467 0
## 58 1 58 240532 396325 0
## 59 1 59 239863 394422 0
## 60 1 60 238698 393203 0
## 61 1 61 237642 392085 0
## 62 1 62 236596 390192 0
## 63 1 63 236034 388005 0
## 64 1 64 234850 386020 0
## 65 1 65 233857 384538 0
## 66 1 66 232992 383014 0
## 67 1 67 231834 381593 0
## 68 1 68 230822 380371 0
## 69 1 69 229856 379106 0
## 70 1 70 229066 377304 0
## 71 1 71 228082 375563 0
## 72 1 72 227038 374331 0
## 73 1 73 226134 372745 0
## 74 1 74 225451 371146 0
## 75 1 75 224422 369397 0
## 76 1 76 223279 367962 0
## 77 1 77 222569 366498 0
## 78 1 78 221814 365544 0
## 79 1 79 221107 363617 0
## 80 1 80 220430 362028 0
## 81 1 81 219421 361082 0
## 82 1 82 218347 359240 0
## 83 1 83 216774 357777 0
## 84 1 84 215842 356105 0
## 85 1 85 214840 354904 0
## 86 1 86 214001 352969 0
## 87 1 87 213206 351266 0
## 88 1 88 212783 349694 0
## 89 1 89 211608 347842 0
## 90 1 90 210504 346161 0
## 91 1 91 209635 344606 0
## 92 1 92 208535 343247 0
## 93 1 93 207830 342108 0
## 94 1 94 207069 340923 0
## 95 1 95 206275 339254 0
## 96 1 96 205506 337757 0
## 97 1 97 204634 336439 0
## 98 1 98 203727 335024 0
## 99 1 99 202711 333653 0
## 100 1 100 201926 332438 0
## 101 1 101 200475 331049 0
## 102 1 102 199771 329692 0
## 103 1 103 199045 328321 0
## 104 1 104 198356 326753 0
## 105 1 105 197753 325424 0
## 106 1 106 196900 324251 0
## 107 1 107 196035 323158 0
## 108 1 108 195600 321547 0
## 109 1 109 194767 320085 0
## 110 1 110 194002 318737 0
## 111 1 111 193406 316854 0
## 112 1 112 192684 315373 0
## 113 1 113 191541 314113 0
## 114 1 114 190801 313169 0
## 115 1 115 189859 312303 0
## 116 1 116 189060 310881 0
## 117 1 117 188301 309180 0
## 118 1 118 187320 307726 0
## 119 1 119 186255 306891 0
## 120 1 120 185342 305554 0
## 121 1 121 184718 304001 0
## 122 1 122 183554 303124 0
## 123 1 123 182896 301853 0
## 124 1 124 182161 300049 0
## 125 1 125 181532 298944 0
## 126 1 126 180445 297469 0
## 127 1 127 179750 296264 0
## 128 1 128 178955 295181 0
## 129 1 129 178242 293936 0
## 130 1 130 177283 292823 0
## 131 1 131 176589 291738 0
## 132 1 132 175892 290607 0
## 133 1 133 175284 289204 0
## 134 1 134 174691 287693 0
## 135 1 135 173948 286490 0
## 136 1 136 173096 285283 0
## 137 1 137 172473 283860 0
## 138 1 138 171649 282437 0
## 139 1 139 170778 281624 0
## 140 1 140 170312 280150 0
## 141 1 141 169547 278638 0
## 142 1 142 168960 277728 0
## 143 1 143 168024 276540 0
## 144 1 144 167436 275617 0
## 145 1 145 166813 274282 0
## 146 1 146 166220 272954 0
## 147 1 147 165553 271842 0
## 148 1 148 164704 270150 0
## 149 1 149 164111 269173 0
## 150 1 150 163372 268368 0
## 151 1 151 162655 267053 0
## 152 1 152 161633 266001 0
## 153 1 153 160701 264742 0
## 154 1 154 160104 263830 0
## 155 1 155 159388 262589 0
## 156 1 156 158827 261238 0
## 157 1 157 158295 260104 0
## 158 1 158 157952 259010 0
## 159 1 159 157412 257830 0
## 160 1 160 156720 257009 0
## 161 1 161 155710 256429 0
## 162 1 162 155095 255574 0
## 163 1 163 154339 254461 0
## 164 1 164 153498 253242 0
## 165 1 165 152805 252122 0
## 166 1 166 151886 250900 0
## 167 1 167 151435 249652 0
## 168 1 168 150664 248537 0
## 169 1 169 150133 247325 0
## 170 1 170 149724 246100 0
## 171 1 171 148932 245384 0
## 172 1 172 148388 244000 0
## 173 1 173 147905 242890 0
## 174 1 174 147457 241843 0
## 175 1 175 146622 240619 0
## 176 1 176 146044 239529 0
## 177 1 177 145337 237919 0
## 178 1 178 144607 236678 0
## 179 1 179 143683 235772 0
## 180 1 180 143173 234653 0