simOcc {AHMbook} | R Documentation |

Function to simulate static occupancy measurements replicated at M sites during J occasions. Population closure is assumed for each site. Expected occurrence may be affected by elevation (elev), forest cover (forest) and their interaction. Expected detection probability may be affected by elevation, wind speed (wind) and their interaction.

simOcc(M = 267, J = 3, mean.occupancy = 0.6, beta1 = -2, beta2 = 2, beta3 = 1, mean.detection = 0.3, time.effects = c(-1, 1), alpha1 = -1, alpha2 = -3, alpha3 = 0, sd.lp = 0.5, b = 2, show.plot = TRUE)

`M` |
Number of spatial replicates (sites) |

`J` |
Number of temporal replicates (occasions) |

`mean.occupancy` |
Mean occurrence at value 0 of occurrence covariates |

`beta1` |
Main effect of elevation on occurrence |

`beta2` |
Main effect of forest cover on occurrence |

`beta3` |
Interaction effect on occurrence of elevation and forest cover |

`mean.detection` |
Mean detection prob. at value 0 of detection covariates |

`time.effects` |
bounds (on logit scale) for uniform distribution from which time effects gamma will be drawn |

`alpha1` |
Main effect of elevation on detection probability |

`alpha2` |
Main effect of wind speed on detection probability |

`alpha3` |
Interaction effect on detection of elevation and wind speed |

`sd.lp` |
standard deviation of random site effects (on logit scale) |

`b` |
constant value of 'behavioral response' leading to 'trap-happiness' (if b > 0) or 'trap shyness' (if b < 0) |

`show.plot` |
if TRUE, plots of the data will be displayed; set to FALSE if you are running simulations |

A list with the values of the input arguments and the following additional elements:

`gamma ` |
The time effects, a vector of length J |

`eps ` |
Individual random effects, a vector of length M |

`elev ` |
Elevation, a vector of length M |

`forest ` |
Forest cover, a vector of length M |

`wind ` |
wind speed, a M x J matrix |

`psi ` |
Probability of occurrence, a vector of length M |

`z ` |
Realized occurrence (0/1), a vector of length M |

`p ` |
probability of capture, possibly with a behavioral effect, a M x J matrix |

`p0 ` |
probability of capture when not captured on previous occasion, a M x J matrix |

`p1 ` |
probability of capture when captured on previous occasion, a M x J matrix |

`y ` |
simulated capture history, a M x J matrix |

`sumZ ` |
True number of occupied sites |

`sumZ.obs ` |
Number of sites observed to be occupied |

`psi.fs.true ` |
True proportion of occupied sites in sample (sumZ/N) |

`psi.fs.obs ` |
Proportion of sites observed to be occupied (sumZ.obs/N) |

Marc Kéry & Andy Royle

Kéry, M. & Royle, J.A. (2016) *Applied Hierarchical Modeling in Ecology* AHM1 - 10.8.

# Generate data with the default arguments and look at the structure tmp <- simOcc() str(tmp) # Simplest possible occupancy model, with constant occupancy and detection str(simOcc(mean.occ=0.6, beta1=0, beta2=0, beta3=0, mean.det=0.3, time.effects=c(0, 0), alpha1=0, alpha2=0, alpha3=0, sd.lp=0, b=0)) # psi = 1 (i.e., species occurs at every site) str(simOcc(mean.occ=1)) # p = 1 (i.e., species is always detected when it occurs) str(simOcc(mean.det=1)) # Other potentially interesting settings include these: str(simOcc(J = 2)) # Only 2 surveys str(simOcc(M = 1, J = 100)) # No spatial replicates, but 100 measurements str(simOcc(beta3 = 1)) # Including interaction elev-wind on p str(simOcc(mean.occ = 0.96)) # A really common species str(simOcc(mean.occ = 0.05)) # A really rare species str(simOcc(mean.det = 0.96)) # A really easy species str(simOcc(mean.det = 0.05)) # A really hard species str(simOcc(mean.det = 0)) # The dreaded invisible species str(simOcc(alpha1=-2, beta1=2)) # Opposing effects of elev on psi and p str(simOcc(J = 10, time.effects = c(-5, 5))) # Huge time effects on p str(simOcc(sd.lp = 10)) # Huge (random) site effects on p str(simOcc(J = 10, b = 0)) # No behavioural response in p str(simOcc(J = 10, b = 2)) # Trap happiness str(simOcc(J = 10, b = -2)) # Trap shyness

[Package *AHMbook* version 0.2.3 Index]