Skip to content

Commit e7d9403

Browse files
committed
Day 2
1 parent 30edf84 commit e7d9403

10 files changed

+3035
-10
lines changed

README.md

+2-1
Original file line numberDiff line numberDiff line change
@@ -33,7 +33,8 @@ Each day will consist of 3 sessions, roughly one hour each, with short breaks in
3333
| **Day 2** | **Behavioral complexities** | |
3434
| | 1. Statistical assumptions and nuisance variables | [Slides](https://peter.solymos.org/qpad-workshop/day2-1-intro.pdf) |
3535
| | 2. Behavioral complexities | [Notes](https://peter.solymos.org/qpad-workshop/day2-2-behavior.html) |
36-
| | 3. Removal modeling techniques | [Notes](https://peter.solymos.org/qpad-workshop/day2-3-removel.html) |
36+
| | 3. Removal modeling techniques 1. | [Notes](https://peter.solymos.org/qpad-workshop/day2-3-removal.html) |
37+
| | 3. Removal modeling techniques 2. | [Notes](https://peter.solymos.org/qpad-workshop/day2-3-mixtures.html) |
3738
| **Day 3** | **The detection process** | |
3839
| | 1. Distance sampling | |
3940
| | 2. Calibrating population density | |

day2-1-intro.Rmd

+14-4
Original file line numberDiff line numberDiff line change
@@ -283,15 +283,15 @@ And a lot more that we'll cover later
283283

284284
# Why do we need simulation
285285

286-
Probabilistic simulation is useful to test how well a stat method is working _if_ the assumptions are met.
286+
Statistical techniques are expected to provide unbiased estimates when the data generation process follows the assumptions of the model. __Testing this is standard practice for any new technique.__
287287

288-
It is not quite as useful for assessing how robust the method is when the assumptions are _not_ met.
288+
But testing if the __model assumptions are realistic__ is rarely evaluated with the same rigor.
289289

290-
You need to link to biological mechanisms to do sensitivity analysis.
290+
For that we need to link to biological mechanisms to do sensitivity analysis. This is a lot more complex than a coin toss.
291291

292292
***
293293

294-
# It is a good pass time
294+
# Capturing complexity
295295

296296
\centering
297297
_"I've yet to see any problem, however complicated, which when you looked at it the right way didn’t become still more complicated."_ --- Poul Anderson, Call Me Joe
@@ -309,6 +309,16 @@ _"Love the simulation we’re dreaming in"_ --- Dua Lipa, Physical
309309

310310
***
311311

312+
# Some history
313+
314+
Movement trajectories of bird flocks to mitigate mortality near airports and wind farms.
315+
316+
Beth Rigby ([2016](https://conservancy.umn.edu/handle/11299/193439)): _Simulating Effects of Imperfect Detectability in Bird Surveys_
317+
318+
But not a lot more...
319+
320+
***
321+
312322
# bSims goals
313323

314324
- Allow easy __testing of statistical assumptions__ and explore effects of violating these assumptions

day2-3-removal.Rmd

+1-1
Original file line numberDiff line numberDiff line change
@@ -211,7 +211,7 @@ Singing rates of birds vary with time of day, time of year, breeding status, and
211211

212212
Let's fit a couple of time-varying models using `DAY` and `TSSR` as covariates:
213213

214-
```{r beh-Me,cache=TRUE}
214+
```{r beh-Me}
215215
Me1 <- cmulti(Y | D ~ DAY, X, type="rem")
216216
Me2 <- cmulti(Y | D ~ TSSR, X, type="rem")
217217
```

day2-4-mixtures.Rmd

+4-3
Original file line numberDiff line numberDiff line change
@@ -10,7 +10,7 @@ load("data/josm-data.rda") # JOSM data
1010
set.seed(1)
1111
```
1212

13-
```{r eval=FALSE}
13+
```{r eval=TRUE}
1414
## data analysis from day2-3 file
1515
yall <- Xtab(~ SiteID + Dur + SpeciesID,
1616
josm$counts[josm$counts$DetectType1 != "V",])
@@ -385,8 +385,9 @@ xall <- bsims_all(
385385
density=c(1, 1, 0),
386386
vocal_rate=phi,
387387
move_rate=1, movement=0.2,
388-
tau=tau,
389-
tint=tbr, rint=rbr)
388+
tau=0.8,
389+
tint=c(3,5,10),
390+
rint=c(0.5, 1, 1.5))
390391
xall
391392
```
392393

0 commit comments

Comments
 (0)