# Analysis of mallard data illustrating basic RMark features
# How to fit a nest level categorical variable (habitat)
# 2019-05-01 CJS Initial code
# This is the mallard data that ships with RMark
library(ggplot2)
## Registered S3 methods overwritten by 'ggplot2':
## method from
## [.quosures rlang
## c.quosures rlang
## print.quosures rlang
library(readxl)
library(RMark)
## This is RMark 2.2.6
## Documentation available at http://www.phidot.org/software/mark/rmark/RMarkDocumentation.zip
# The dataframe must contain the following fields with the following names
#
# FirstFound: day the nest was first found
# LastPresent: last day that a chick was present in the nest
# LastChecked: last day the nest was checked
# Fate: fate of the nest; 0=hatch an
# Freq: number of nests with this data
#
# Also contains the following fields
# Robel - Measurement of Robel pole of visibility of nest
# PpnGrass - proportion of grassland cover on the 10.4 km2 study site t
# AgeFound - Age of nest when found
# AgeDay1 - Age of nest on day 1 of study (can be negative)
# Habitat - N=Native; P=Planted; W=Wetland; R=roadside right of way
#
# In this example, the multiple visits to a nest have been collapsed
# to a single record for each nest.
# In more complex examples, you may have multple records per nest
# as shown in the mallard example.
#
malldata <- readxl::read_excel(file.path("..","mallard.xlsx"),
sheet="mallard")
head(malldata)
## # A tibble: 6 x 10
## FirstFound LastPresent LastChecked Fate Freq Robel PpnGrass AgeFound
## <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 73 89 89 0 1 6 0.800 13
## 2 63 90 90 0 1 3.75 0.668 5
## 3 70 70 76 1 1 3.75 0.800 13
## 4 63 81 85 1 1 3.12 0.668 6
## 5 61 61 66 1 1 4.5 0.668 4
## 6 57 57 61 1 1 4.25 0.668 1
## # ... with 2 more variables: AgeDay1 <dbl>, Habitat <chr>
malldata <- as.data.frame(malldata)
malldata$Habitat <- factor(malldata$Habitat) # RMark wants categorical variable to be factors
# what are the parameters of the model
# There is only one parameter, the daily survival probality (S)
setup.parameters("Nest", check=TRUE)
## [1] "S"
# 1. Process the data.
# The nocc variable the number of sampling occasions in the data
mall.proc <- process.data(malldata, model="Nest",
group="Habitat",
nocc=max(malldata$LastChecked))
mall.proc
## $data
## FirstFound LastPresent LastChecked Fate freq Robel PpnGrass AgeFound
## 1 73 89 89 0 1 6.000 0.8002 13
## 2 63 90 90 0 1 3.750 0.6684 5
## 3 70 70 76 1 1 3.750 0.8002 13
## 4 63 81 85 1 1 3.125 0.6684 6
## 5 61 61 66 1 1 4.500 0.6684 4
## 6 57 57 61 1 1 4.250 0.6684 1
## 7 67 88 88 0 1 8.000 0.6684 12
## 8 58 87 87 0 1 5.875 0.4651 3
## 9 74 86 86 0 1 4.500 0.8002 20
## 10 65 65 70 1 1 3.625 0.8002 11
## 11 63 88 88 0 1 5.250 0.8666 10
## 12 61 66 70 1 1 3.625 0.6684 8
## 13 63 84 84 0 1 3.750 0.8002 11
## 14 74 74 78 1 1 5.625 0.8002 22
## 15 56 56 63 1 1 4.875 0.6684 4
## 16 57 87 87 0 1 3.500 0.6684 5
## 17 74 86 86 0 1 3.375 0.8002 22
## 18 68 87 87 0 1 4.875 0.8002 16
## 19 56 56 63 1 1 3.875 0.6684 4
## 20 66 66 73 1 1 3.375 0.8002 15
## 21 81 84 84 0 1 5.750 0.8002 30
## 22 63 63 68 1 1 5.250 0.8666 13
## 23 56 56 63 1 1 5.875 0.6684 7
## 24 56 63 69 1 1 5.750 0.6684 7
## 25 74 84 84 0 1 4.125 0.8002 25
## 26 55 64 69 1 1 5.875 0.4651 7
## 27 56 80 80 0 1 5.000 0.6684 8
## 28 56 56 63 1 1 2.250 0.6684 8
## 29 50 60 64 1 1 3.750 0.7281 2
## 30 58 58 63 1 1 3.375 0.6684 10
## 31 58 73 77 1 1 1.875 0.8666 10
## 32 57 80 80 0 1 2.250 0.8666 9
## 33 63 63 68 1 1 4.625 0.8002 16
## 34 60 79 79 0 1 6.125 0.6684 13
## 35 61 61 66 1 1 4.750 0.6684 14
## 36 56 69 74 1 1 4.125 0.6684 9
## 37 65 81 81 0 1 5.000 0.8002 18
## 38 57 68 73 1 1 3.000 0.4651 11
## 39 73 76 76 0 1 4.875 0.8002 27
## 40 54 64 68 1 1 5.000 0.4651 8
## 41 75 79 79 0 1 3.750 0.8002 29
## 42 55 55 59 1 1 3.000 0.6684 9
## 43 50 56 60 1 1 5.500 0.7281 4
## 44 58 58 64 1 1 5.625 0.8666 12
## 45 62 80 80 0 1 4.875 0.6684 16
## 46 54 74 74 0 1 4.625 0.4651 9
## 47 61 66 70 1 1 3.250 0.8666 16
## 48 73 76 76 0 1 6.000 0.8002 28
## 49 52 52 57 1 1 4.250 0.7281 7
## 50 53 75 75 0 1 4.750 0.7281 9
## 51 52 61 66 1 1 3.375 0.7281 8
## 52 63 78 78 0 1 4.750 0.8666 19
## 53 63 77 77 0 1 3.625 0.8002 19
## 54 70 70 76 1 1 4.375 0.8002 26
## 55 73 79 79 0 1 3.750 0.8002 29
## 56 59 77 77 0 1 4.125 0.8666 15
## 57 50 50 56 1 1 3.625 0.7281 7
## 58 63 73 73 0 1 2.750 0.8666 20
## 59 63 73 73 0 1 3.500 0.8002 20
## 60 49 49 55 1 1 3.750 0.4651 6
## 61 56 63 69 1 1 6.250 0.6684 14
## 62 47 53 59 1 1 5.250 0.4651 5
## 63 44 44 48 1 1 4.125 0.4651 2
## 64 55 74 74 0 1 5.125 0.4651 14
## 65 45 56 61 1 1 4.750 0.4651 4
## 66 47 53 59 1 1 6.875 0.4651 7
## 67 61 61 66 1 1 4.500 0.6684 21
## 68 61 61 66 1 1 4.625 0.6684 21
## 69 67 70 70 0 1 5.625 0.6684 27
## 70 55 68 68 0 1 4.375 0.4651 15
## 71 43 46 50 1 1 4.750 0.8666 3
## 72 47 67 67 0 1 3.875 0.8666 7
## 73 64 73 73 0 1 3.500 0.8002 24
## 74 50 60 64 1 1 1.875 0.7281 11
## 75 63 63 68 1 1 2.500 0.8002 24
## 76 54 68 74 1 1 4.625 0.4651 15
## 77 64 64 69 1 1 3.500 0.8002 25
## 78 43 70 70 0 1 4.625 0.5079 4
## 79 42 46 46 0 1 5.750 0.4651 4
## 80 60 60 65 1 1 3.750 0.8666 22
## 81 43 46 46 0 1 3.875 0.8666 5
## 82 61 69 69 0 1 4.750 0.6684 23
## 83 61 66 70 1 1 4.125 0.6684 23
## 84 68 72 72 0 1 4.875 0.8002 30
## 85 52 65 69 1 1 4.125 0.7281 14
## 86 47 71 71 0 1 3.625 0.8666 9
## 87 65 70 72 1 1 4.375 0.8002 28
## 88 43 43 49 1 1 4.000 0.4651 6
## 89 41 46 46 0 1 2.375 0.3344 4
## 90 43 46 50 1 1 5.125 0.8666 6
## 91 56 72 72 0 1 5.500 0.6684 19
## 92 42 68 72 1 1 5.875 0.4651 5
## 93 63 67 67 0 1 4.125 0.8002 26
## 94 62 70 70 0 1 4.500 0.8666 25
## 95 55 66 66 0 1 3.750 0.4651 18
## 96 44 72 72 0 1 3.750 0.8666 7
## 97 67 68 68 0 1 9.250 0.6684 31
## 98 61 65 65 0 1 4.625 0.6684 25
## 99 62 67 67 0 1 3.500 0.6684 26
## 100 42 57 60 1 1 5.375 0.3344 6
## 101 52 52 57 1 1 2.875 0.7281 16
## 102 50 64 69 1 1 4.000 0.7281 14
## 103 48 55 60 1 1 5.875 0.7281 12
## 104 63 67 67 0 1 4.625 0.8666 27
## 105 62 67 67 0 1 5.000 0.6684 26
## 106 63 67 67 0 1 4.375 0.6684 28
## 107 43 55 60 1 1 3.375 0.5079 8
## 108 40 50 56 1 1 2.375 0.3344 5
## 109 60 60 65 1 1 5.500 0.6684 25
## 110 42 64 68 1 1 1.750 0.5079 7
## 111 59 70 70 0 1 4.500 0.6684 24
## 112 42 66 66 0 1 2.875 0.3344 7
## 113 40 40 44 1 1 4.625 0.7265 6
## 114 46 46 52 1 1 5.250 0.4651 12
## 115 58 63 63 0 1 4.875 0.8666 24
## 116 40 60 64 1 1 3.375 0.3344 6
## 117 40 50 56 1 1 2.875 0.3344 6
## 118 42 46 46 0 1 5.375 0.4651 8
## 119 40 44 50 1 1 2.875 0.3344 6
## 120 42 42 48 1 1 2.625 0.3344 9
## 121 57 66 68 1 1 4.750 0.6684 24
## 122 63 67 67 0 1 2.750 0.8666 30
## 123 40 44 50 1 1 4.625 0.7265 7
## 124 47 67 67 0 1 4.375 0.7281 14
## 125 57 57 61 1 1 4.250 0.6684 24
## 126 40 44 50 1 1 4.000 0.3344 7
## 127 38 55 65 1 1 4.875 0.3344 5
## 128 56 56 63 1 1 4.500 0.6684 24
## 129 58 63 63 0 1 5.875 0.8666 26
## 130 60 66 66 0 1 5.250 0.6684 28
## 131 40 66 66 0 1 2.875 0.3344 8
## 132 40 50 56 1 1 3.000 0.3344 8
## 133 40 63 63 0 1 3.250 0.3344 8
## 134 47 63 63 0 1 4.250 0.4651 15
## 135 37 65 65 0 1 4.750 0.4651 5
## 136 60 66 66 0 1 5.125 0.8666 28
## 137 38 38 44 1 1 5.500 0.3344 7
## 138 40 40 44 1 1 3.000 0.3344 9
## 139 41 63 63 0 1 5.000 0.5079 10
## 140 44 65 65 0 1 3.500 0.8666 13
## 141 42 46 46 0 1 5.125 0.4651 11
## 142 55 63 63 0 1 5.125 0.4651 24
## 143 47 66 66 0 1 2.750 0.4651 16
## 144 41 57 62 1 1 3.875 0.7265 10
## 145 48 66 66 0 1 6.250 0.7281 17
## 146 40 40 44 1 1 3.375 0.3344 9
## 147 55 63 63 0 1 4.875 0.4651 25
## 148 40 40 44 1 1 3.500 0.3344 10
## 149 41 62 65 1 1 4.000 0.5079 11
## 150 38 38 43 1 1 2.000 0.3344 8
## 151 58 58 64 1 1 5.125 0.8666 28
## 152 41 45 51 1 1 3.375 0.3344 11
## 153 36 36 41 1 1 2.750 0.3344 6
## 154 40 40 44 1 1 2.875 0.3344 10
## 155 62 63 63 0 1 2.500 0.6684 32
## 156 38 38 44 1 1 4.250 0.3344 8
## 157 35 40 45 1 1 4.125 0.3344 5
## 158 36 64 64 0 1 4.250 0.3344 6
## 159 35 35 40 1 1 2.875 0.3344 5
## 160 36 41 46 1 1 2.750 0.3344 7
## 161 40 61 61 0 1 3.000 0.3344 11
## 162 47 63 63 0 1 3.000 0.4651 18
## 163 31 31 36 1 1 4.125 0.6587 2
## 164 34 37 42 1 1 1.750 0.1876 5
## 165 55 63 63 0 1 5.250 0.4651 26
## 166 38 43 49 1 1 3.875 0.3344 9
## 167 35 35 40 1 1 3.750 0.3344 6
## 168 40 64 64 0 1 2.500 0.4651 11
## 169 36 36 42 1 1 2.000 0.3344 7
## 170 31 60 60 0 1 1.250 0.1876 3
## 171 38 55 60 1 1 3.625 0.3344 10
## 172 35 57 57 0 1 3.250 0.3344 7
## 173 43 63 63 0 1 3.500 0.4651 15
## 174 43 47 55 1 1 3.625 0.5079 15
## 175 29 29 34 1 1 3.375 0.1876 2
## 176 43 59 59 0 1 3.875 0.5079 16
## 177 40 55 55 0 1 2.000 0.7265 13
## 178 41 45 51 1 1 2.500 0.3344 14
## 179 54 62 62 0 1 4.500 0.4651 27
## 180 36 51 51 0 1 4.625 0.3344 9
## 181 37 42 48 1 1 4.000 0.3344 10
## 182 43 59 59 0 1 3.500 0.5079 16
## 183 41 62 62 0 1 3.500 0.5079 14
## 184 40 55 61 1 1 6.750 0.7265 13
## 185 30 42 42 0 1 4.875 0.6587 3
## 186 36 59 59 0 1 4.125 0.3344 10
## 187 31 31 36 1 1 3.375 0.6587 5
## 188 54 59 61 1 1 4.625 0.4651 28
## 189 35 52 56 1 1 3.000 0.3344 9
## 190 42 46 46 0 1 5.125 0.4651 16
## 191 38 55 60 1 1 2.750 0.3344 12
## 192 54 61 61 0 1 2.875 0.4651 28
## 193 54 61 61 0 1 3.875 0.4651 28
## 194 27 31 31 0 1 3.400 0.6587 1
## 195 34 58 58 0 1 5.125 0.9671 8
## 196 40 50 56 1 1 3.500 0.3344 14
## 197 48 55 60 1 1 3.250 0.7281 23
## 198 31 31 36 1 1 2.875 0.6587 6
## 199 31 31 36 1 1 3.750 0.6587 6
## 200 31 42 46 1 1 4.000 0.6587 6
## 201 31 46 54 1 1 3.625 0.6587 6
## 202 31 58 58 0 1 4.875 0.6587 6
## 203 34 34 38 1 1 4.500 0.9671 9
## 204 31 58 58 0 1 3.625 0.9671 6
## 205 31 42 46 1 1 2.500 0.9671 6
## 206 40 44 50 1 1 3.250 0.3344 15
## 207 43 59 59 0 1 4.875 0.5079 18
## 208 41 50 50 0 1 3.000 0.3344 17
## 209 29 29 34 1 1 3.125 0.1876 5
## 210 36 36 40 1 1 2.875 0.3344 12
## 211 41 45 51 1 1 3.000 0.3344 17
## 212 25 25 29 1 1 4.250 0.9671 1
## 213 36 52 58 1 1 4.500 0.3344 12
## 214 31 42 46 1 1 3.375 0.9671 7
## 215 31 58 58 0 1 4.375 0.6587 7
## 216 31 31 42 1 1 4.750 0.1876 7
## 217 31 56 56 0 1 2.000 0.1876 7
## 218 31 31 36 1 1 1.750 0.1876 7
## 219 31 46 54 1 1 4.000 0.6587 7
## 220 30 56 56 0 1 1.125 0.6587 6
## 221 25 30 35 1 1 3.375 0.6587 2
## 222 31 31 36 1 1 3.750 0.6587 8
## 223 32 34 34 0 1 4.250 0.9671 9
## 224 47 47 51 1 1 4.250 0.7281 24
## 225 29 29 34 1 1 3.125 0.1876 6
## 226 30 30 41 1 1 3.250 0.6587 7
## 227 26 35 35 0 1 2.250 0.9671 3
## 228 35 41 45 1 1 4.625 0.3344 12
## 229 52 58 58 0 1 4.250 0.4651 29
## 230 26 30 30 0 1 1.750 0.9671 3
## 231 26 48 48 0 1 2.125 0.6587 4
## 232 31 55 55 0 1 4.375 0.6587 9
## 233 28 35 40 1 1 5.125 0.1876 6
## 234 34 34 37 1 1 3.000 0.1876 12
## 235 31 31 36 1 1 3.500 0.6587 9
## 236 31 42 46 1 1 2.625 0.9671 9
## 237 29 40 45 1 1 3.625 0.1876 8
## 238 41 41 45 1 1 2.875 0.3344 20
## 239 35 56 56 0 1 3.375 0.3344 14
## 240 27 30 34 1 1 0.750 0.9671 6
## 241 29 29 34 1 1 3.125 0.1876 8
## 242 25 25 32 1 1 3.125 0.6587 4
## 243 31 31 36 1 1 4.000 0.6587 10
## 244 31 53 53 0 1 3.500 0.9671 10
## 245 22 45 53 1 1 1.000 0.6587 1
## 246 40 40 44 1 1 2.625 0.3344 20
## 247 26 40 45 1 1 3.750 0.6587 6
## 248 26 26 30 1 1 1.375 0.6587 6
## 249 29 34 37 1 1 3.500 0.1876 9
## 250 41 55 55 0 1 3.875 0.7265 21
## 251 40 54 54 0 1 5.625 0.7265 20
## 252 40 40 44 1 1 3.000 0.3344 20
## 253 29 51 51 0 1 2.250 0.1876 10
## 254 31 46 54 1 1 3.625 0.6587 12
## 255 27 34 40 1 1 2.375 0.9671 8
## 256 26 30 35 1 1 2.875 0.9671 7
## 257 27 30 35 1 1 4.000 0.9671 8
## 258 25 49 54 1 1 4.375 0.9671 6
## 259 23 51 51 0 1 3.125 0.9378 4
## 260 23 23 27 1 1 2.250 0.6587 4
## 261 32 53 53 0 1 4.167 0.9671 13
## 262 45 50 54 1 1 3.875 0.4651 26
## 263 27 37 42 1 1 4.000 0.6587 8
## 264 28 54 54 0 1 2.875 0.1876 9
## 265 29 40 45 1 1 3.000 0.1876 10
## 266 26 35 41 1 1 2.500 0.9671 8
## 267 26 40 40 0 1 1.500 0.9671 8
## 268 29 29 34 1 1 2.375 0.1876 11
## 269 36 53 53 0 1 4.375 0.3344 18
## 270 35 53 53 0 1 2.625 0.3344 17
## 271 47 52 52 0 1 4.375 0.8666 29
## 272 31 46 53 1 1 3.375 0.9671 13
## 273 25 25 29 1 1 2.750 0.9671 7
## 274 31 52 52 0 1 1.625 0.9671 13
## 275 46 51 51 0 1 5.500 0.4651 28
## 276 36 40 45 1 1 3.750 0.3344 18
## 277 30 30 41 1 1 2.250 0.6587 12
## 278 26 30 30 0 1 2.250 0.9671 9
## 279 24 52 52 0 1 3.625 0.6587 7
## 280 41 50 50 0 1 2.000 0.7265 24
## 281 40 49 49 0 1 2.375 0.7265 23
## 282 29 52 52 0 1 4.375 0.6587 12
## 283 36 46 46 0 1 1.750 0.3344 20
## 284 31 51 51 0 1 3.500 0.9671 15
## 285 43 50 51 1 1 4.500 0.8666 27
## 286 40 51 51 0 1 5.625 0.3344 24
## 287 40 49 49 0 1 5.125 0.7265 24
## 288 32 44 44 0 1 5.750 0.1876 16
## 289 20 31 31 0 1 2.625 0.9378 5
## 290 19 19 21 1 1 3.250 0.9378 4
## 291 23 32 37 1 1 3.125 0.6587 8
## 292 47 48 48 0 1 3.000 0.8666 32
## 293 44 48 48 0 1 3.000 0.8666 29
## 294 21 30 30 0 1 3.625 0.9248 6
## 295 22 22 27 1 1 1.375 0.9378 7
## 296 25 46 46 0 1 5.375 0.6587 10
## 297 29 44 50 1 1 4.500 0.6587 14
## 298 29 29 34 1 1 2.625 0.1876 14
## 299 41 50 50 0 1 3.750 0.3344 26
## 300 38 38 43 1 1 3.750 0.3344 23
## 301 40 46 46 0 1 2.625 0.3344 26
## 302 21 30 41 1 1 1.000 0.9248 7
## 303 21 25 30 1 1 3.125 0.9248 7
## 304 40 49 49 0 1 3.375 0.3344 26
## 305 36 49 49 0 1 6.250 0.3344 22
## 306 40 49 49 0 1 6.750 0.7265 26
## 307 43 44 44 0 1 3.500 0.5079 29
## 308 23 46 46 0 1 5.000 0.9378 9
## 309 22 37 42 1 1 2.750 0.9378 8
## 310 31 49 49 0 1 3.250 0.6587 17
## 311 23 46 46 0 1 3.625 0.9378 9
## 312 44 44 48 1 1 5.125 0.4651 30
## 313 40 48 48 0 1 3.250 0.3344 26
## 314 17 17 23 1 1 2.625 0.9248 4
## 315 47 48 48 0 1 3.500 0.8666 34
## 316 35 46 46 0 1 3.500 0.1876 22
## 317 22 47 47 0 1 2.250 0.6587 9
## 318 16 22 22 0 1 3.125 0.9248 3
## 319 20 28 31 1 1 1.000 0.9378 7
## 320 17 23 29 1 1 5.875 0.9248 4
## 321 29 29 34 1 1 3.500 0.1876 16
## 322 32 43 47 1 1 4.500 0.9671 19
## 323 29 40 45 1 1 3.125 0.1876 16
## 324 23 47 47 0 1 3.250 0.6587 11
## 325 19 46 46 0 1 2.625 0.6607 7
## 326 20 36 36 0 1 1.000 0.9378 8
## 327 25 29 32 1 1 3.375 0.9671 13
## 328 27 31 37 1 1 2.375 0.6587 15
## 329 21 25 29 1 1 2.125 0.9378 9
## 330 21 29 34 1 1 1.750 0.9378 9
## 331 35 47 47 0 1 3.625 0.3344 23
## 332 21 31 36 1 1 2.625 0.9378 9
## 333 42 47 47 0 1 4.250 0.5079 30
## 334 22 37 42 1 1 1.375 0.9378 10
## 335 21 44 44 0 1 1.625 0.9378 10
## 336 29 34 37 1 1 2.000 0.1876 18
## 337 20 42 46 1 1 3.375 0.9248 9
## 338 20 20 24 1 1 2.375 0.9378 9
## 339 12 21 25 1 1 0.875 0.8765 1
## 340 30 46 46 0 1 2.375 0.6587 19
## 341 14 24 29 1 1 5.000 0.6607 3
## 342 14 21 25 1 1 4.375 0.9248 3
## 343 17 17 22 1 1 2.375 0.6607 6
## 344 13 13 17 1 1 3.375 0.8765 2
## 345 13 13 17 1 1 3.250 0.8765 2
## 346 13 21 21 0 1 3.375 0.8765 2
## 347 17 17 22 1 1 2.625 0.6607 6
## 348 23 42 42 0 1 2.875 0.9378 12
## 349 21 46 46 0 1 2.000 0.9378 10
## 350 21 46 46 0 1 2.250 0.9378 10
## 351 32 44 44 0 1 4.750 0.9671 22
## 352 23 27 32 1 1 3.250 0.6587 13
## 353 16 22 22 0 1 4.875 0.9248 6
## 354 15 24 27 1 1 3.625 0.9248 5
## 355 21 29 34 1 1 1.875 0.9378 11
## 356 36 41 41 0 1 3.125 0.3344 26
## 357 12 45 45 0 1 2.500 0.2290 2
## 358 20 20 24 1 1 1.875 0.9378 11
## 359 35 41 41 0 1 2.500 0.1876 26
## 360 25 43 43 0 1 5.000 0.9671 16
## 361 26 44 44 0 1 2.750 0.9671 17
## 362 31 44 44 0 1 3.500 0.9671 22
## 363 16 35 41 1 1 3.500 0.6607 7
## 364 22 25 31 1 1 6.000 0.6587 13
## 365 27 31 37 1 1 2.125 0.6587 18
## 366 31 44 44 0 1 3.500 0.6587 22
## 367 9 9 13 1 1 5.125 0.9197 1
## 368 16 34 34 0 1 4.250 0.6607 8
## 369 13 43 43 0 1 2.250 0.6607 5
## 370 11 15 21 1 1 3.375 0.2290 3
## 371 21 36 36 0 1 3.500 0.9378 13
## 372 16 16 21 1 1 2.875 0.6607 8
## 373 13 13 17 1 1 2.250 0.8765 5
## 374 15 15 20 1 1 3.125 0.9248 7
## 375 31 43 43 0 1 3.125 0.6587 23
## 376 13 17 17 0 1 2.750 0.8765 6
## 377 8 8 12 1 1 3.625 0.9378 1
## 378 15 20 25 1 1 6.375 0.6607 8
## 379 22 25 31 1 1 2.625 0.6587 15
## 380 24 37 37 0 1 2.250 0.9671 17
## 381 16 25 29 1 1 2.875 0.6607 9
## 382 24 39 39 0 1 3.625 0.9671 17
## 383 25 41 41 0 1 1.625 0.9671 18
## 384 10 42 42 0 1 4.250 0.2290 3
## 385 23 27 32 1 1 2.750 0.6587 16
## 386 16 35 35 0 1 3.500 0.6607 9
## 387 17 41 41 0 1 1.625 0.9248 10
## 388 26 39 39 0 1 5.500 0.6587 20
## 389 7 12 17 1 1 4.375 0.2290 1
## 390 17 34 40 1 1 4.375 0.9248 11
## 391 16 40 40 0 1 4.125 0.6607 10
## 392 25 32 33 1 1 3.375 0.9671 19
## 393 11 15 21 1 1 3.125 0.8765 5
## 394 14 20 20 0 1 4.125 0.6607 8
## 395 15 25 29 1 1 0.875 0.6607 9
## 396 16 39 39 0 1 3.750 0.6607 10
## 397 20 24 28 1 1 2.250 0.9378 15
## 398 9 40 40 0 1 1.375 0.9197 4
## 399 11 34 38 1 1 2.750 0.8765 6
## 400 23 36 36 0 1 2.250 0.9378 18
## 401 14 37 37 0 1 0.625 0.9248 9
## 402 22 22 27 1 1 1.000 0.9378 17
## 403 22 40 40 0 1 4.875 0.9378 17
## 404 21 40 40 0 1 2.125 0.9248 16
## 405 10 10 15 1 1 3.250 0.2290 5
## 406 23 23 27 1 1 2.750 0.6587 18
## 407 10 10 14 1 1 3.000 0.9197 5
## 408 28 32 38 1 1 3.625 0.6587 23
## 409 16 34 34 0 1 3.375 0.6607 11
## 410 34 38 38 0 1 4.250 0.9671 29
## 411 13 32 32 0 1 1.875 0.6607 8
## 412 13 28 32 1 1 3.250 0.6607 8
## 413 13 28 32 1 1 1.500 0.2290 8
## 414 12 16 22 1 1 3.250 0.2290 7
## 415 11 35 35 0 1 1.250 0.2290 6
## 416 10 24 27 1 1 3.625 0.9197 5
## 417 6 15 20 1 1 4.000 0.2290 1
## 418 13 21 21 0 1 3.125 0.8765 9
## 419 10 10 15 1 1 3.500 0.2290 6
## 420 10 15 21 1 1 3.250 0.2290 6
## 421 10 25 29 1 1 3.375 0.2290 6
## 422 10 15 21 1 1 5.875 0.2290 6
## 423 10 10 15 1 1 3.250 0.2290 6
## 424 10 10 15 1 1 5.000 0.2290 6
## 425 11 15 21 1 1 3.375 0.2290 7
## 426 13 23 28 1 1 2.500 0.2290 9
## 427 13 32 37 1 1 2.625 0.2290 9
## 428 13 20 23 1 1 3.125 0.6607 9
## 429 15 25 29 1 1 2.125 0.6607 11
## 430 8 8 12 1 1 1.750 0.9197 4
## 431 9 31 31 0 1 3.500 0.9197 5
## 432 9 25 31 1 1 3.125 0.9197 5
## 433 10 20 24 1 1 2.125 0.9197 6
## 434 10 27 27 0 1 3.375 0.9197 6
## 435 11 21 24 1 1 4.125 0.8765 7
## 436 21 21 25 1 1 4.250 0.9248 17
## 437 13 30 30 0 1 2.750 0.8765 9
## 438 13 39 39 0 1 2.000 0.8765 9
## 439 13 13 17 1 1 2.875 0.8765 9
## 440 15 37 39 1 1 2.500 0.9248 11
## 441 10 21 25 1 1 4.375 0.2290 6
## 442 25 25 30 1 1 2.500 0.6587 21
## 443 30 38 38 0 1 2.375 0.6587 26
## 444 8 22 27 1 1 3.875 0.2290 5
## 445 8 17 23 1 1 5.625 0.9378 5
## 446 9 38 38 0 1 3.375 0.2290 6
## 447 23 38 38 0 1 2.875 0.9378 20
## 448 12 22 27 1 1 3.000 0.2290 9
## 449 13 36 36 0 1 2.500 0.2290 10
## 450 6 29 32 1 1 4.000 0.2290 3
## 451 10 10 15 1 1 2.750 0.2290 7
## 452 10 15 21 1 1 4.250 0.2290 7
## 453 30 38 38 0 1 2.375 0.6587 27
## 454 20 38 38 0 1 5.875 0.9248 17
## 455 10 37 37 0 1 2.875 0.9197 7
## 456 13 38 38 0 1 2.750 0.6607 10
## 457 10 21 25 1 1 4.750 0.2290 7
## 458 13 36 36 0 1 1.250 0.6607 10
## 459 10 14 20 1 1 3.000 0.9197 7
## 460 30 38 38 0 1 2.125 0.6587 27
## 461 16 38 38 0 1 3.125 0.6607 13
## 462 17 22 27 1 1 4.000 0.6607 14
## 463 20 35 35 0 1 1.750 0.9378 17
## 464 10 10 14 1 1 3.750 0.9197 8
## 465 9 9 13 1 1 5.750 0.9197 7
## 466 8 12 17 1 1 4.250 0.9378 6
## 467 7 30 35 1 1 2.750 0.2290 5
## 468 7 17 22 1 1 4.500 0.2290 5
## 469 10 15 21 1 1 3.000 0.2290 8
## 470 11 26 31 1 1 5.750 0.2290 9
## 471 11 37 37 0 1 3.375 0.2290 9
## 472 30 30 37 1 1 4.125 0.6587 28
## 473 14 37 37 0 1 4.000 0.9248 12
## 474 14 30 37 1 1 3.125 0.9248 12
## 475 11 31 36 1 1 3.250 0.2290 9
## 476 11 26 31 1 1 2.500 0.2290 9
## 477 10 10 14 1 1 3.750 0.9197 8
## 478 10 37 37 0 1 3.375 0.9197 8
## 479 5 20 24 1 1 3.250 0.2290 3
## 480 10 36 36 0 1 2.750 0.9197 8
## 481 11 15 21 1 1 4.000 0.2290 9
## 482 10 14 20 1 1 2.125 0.9197 8
## 483 6 6 11 1 1 2.625 0.2290 4
## 484 9 36 36 0 1 1.750 0.9197 8
## 485 24 24 29 1 1 3.625 0.9671 23
## 486 11 21 24 1 1 3.750 0.8765 10
## 487 28 34 34 0 1 3.625 0.9671 27
## 488 9 36 36 0 1 4.725 0.9197 8
## 489 10 14 20 1 1 4.000 0.9197 9
## 490 9 25 25 0 1 5.750 0.9197 8
## 491 11 35 35 0 1 3.625 0.2290 10
## 492 10 10 15 1 1 4.875 0.2290 9
## 493 10 33 33 0 1 4.750 0.2290 9
## 494 10 35 35 0 1 5.875 0.2290 9
## 495 10 33 33 0 1 4.500 0.2290 9
## 496 10 36 36 0 1 5.375 0.2290 9
## 497 7 16 21 1 1 4.250 0.2290 6
## 498 6 35 35 0 1 3.750 0.2290 5
## 499 6 36 36 0 1 4.250 0.2290 5
## 500 23 32 35 1 1 2.375 0.6587 23
## 501 10 10 15 1 1 3.375 0.2290 10
## 502 23 27 32 1 1 2.000 0.6587 23
## 503 1 35 35 0 1 4.500 0.9616 1
## 504 4 9 14 1 1 3.875 0.9378 4
## 505 24 35 35 0 1 2.375 0.9671 24
## 506 9 9 13 1 1 6.000 0.9197 9
## 507 8 23 28 1 1 3.125 0.2290 8
## 508 23 23 28 1 1 3.125 0.9378 23
## 509 7 34 34 0 1 3.375 0.2290 7
## 510 8 22 27 1 1 3.750 0.2290 9
## 511 9 13 20 1 1 2.375 0.9197 10
## 512 16 34 34 0 1 3.250 0.6607 17
## 513 10 14 21 1 1 3.625 0.9197 11
## 514 6 20 25 1 1 3.875 0.2290 7
## 515 26 34 34 0 1 2.500 0.9671 27
## 516 23 23 27 1 1 3.500 0.6587 24
## 517 8 8 13 1 1 3.375 0.2290 9
## 518 7 34 34 0 1 3.625 0.2290 8
## 519 30 34 34 0 1 1.750 0.6587 31
## 520 8 12 17 1 1 5.625 0.9197 9
## 521 3 26 30 1 1 2.750 0.9378 4
## 522 10 27 27 0 1 3.375 0.9197 12
## 523 3 13 17 1 1 3.750 0.9378 5
## 524 8 13 17 1 1 4.375 0.2290 10
## 525 2 7 12 1 1 1.625 0.9378 4
## 526 4 9 14 1 1 3.875 0.9378 6
## 527 8 23 29 1 1 5.000 0.9197 10
## 528 7 12 17 1 1 5.375 0.2290 9
## 529 1 31 33 1 1 3.375 0.9616 3
## 530 7 17 22 1 1 0.875 0.2290 9
## 531 14 33 33 0 1 3.125 0.8765 16
## 532 6 11 15 1 1 3.125 0.2290 8
## 533 1 15 21 1 1 0.875 0.9616 3
## 534 1 11 15 1 1 2.750 0.9616 3
## 535 2 2 7 1 1 1.875 0.9616 4
## 536 7 12 12 0 1 3.000 0.2290 10
## 537 17 32 32 0 1 5.875 0.9248 20
## 538 7 7 12 1 1 2.750 0.2290 10
## 539 30 31 31 0 1 1.875 0.6587 33
## 540 2 7 12 1 1 2.750 0.9616 5
## 541 14 32 32 0 1 3.750 0.9248 17
## 542 23 27 32 1 1 3.375 0.6587 26
## 543 8 22 27 1 1 4.000 0.2290 11
## 544 10 30 30 0 1 4.125 0.9197 14
## 545 3 31 31 0 1 4.375 0.9616 7
## 546 23 29 29 0 1 2.625 0.9378 27
## 547 4 30 31 1 1 2.375 0.9378 8
## 548 7 12 17 1 1 8.000 0.2290 11
## 549 10 15 21 1 1 4.500 0.2290 14
## 550 12 27 30 1 1 2.000 0.2290 17
## 551 15 20 24 1 1 3.000 0.9248 20
## 552 21 30 30 0 1 3.125 0.9248 26
## 553 13 30 30 0 1 3.375 0.8765 18
## 554 4 19 23 1 1 4.125 0.2290 9
## 555 14 30 30 0 1 5.250 0.9248 19
## 556 21 30 30 0 1 4.250 0.9248 26
## 557 4 30 30 0 1 2.000 0.9378 9
## 558 10 14 20 1 1 2.250 0.9197 16
## 559 4 9 14 1 1 1.625 0.9378 10
## 560 26 29 29 0 1 2.500 0.9671 32
## 561 10 15 21 1 1 4.250 0.2290 16
## 562 15 20 25 1 1 4.625 0.6607 21
## 563 13 13 17 1 1 2.250 0.8765 20
## 564 16 26 26 0 1 3.250 0.6607 25
## 565 16 16 21 1 1 4.000 0.6607 25
## AgeDay1 Habitat group
## 1 -59 R 3
## 2 -57 P 2
## 3 -56 P 2
## 4 -56 N 1
## 5 -56 P 2
## 6 -55 P 2
## 7 -54 P 2
## 8 -54 N 1
## 9 -53 P 2
## 10 -53 N 1
## 11 -52 P 2
## 12 -52 P 2
## 13 -51 P 2
## 14 -51 N 1
## 15 -51 N 1
## 16 -51 P 2
## 17 -51 P 2
## 18 -51 W 4
## 19 -51 N 1
## 20 -50 P 2
## 21 -50 N 1
## 22 -49 P 2
## 23 -48 N 1
## 24 -48 N 1
## 25 -48 R 3
## 26 -47 P 2
## 27 -47 N 1
## 28 -47 R 3
## 29 -47 R 3
## 30 -47 P 2
## 31 -47 N 1
## 32 -47 N 1
## 33 -46 P 2
## 34 -46 P 2
## 35 -46 P 2
## 36 -46 N 1
## 37 -46 P 2
## 38 -45 N 1
## 39 -45 P 2
## 40 -45 N 1
## 41 -45 P 2
## 42 -45 N 1
## 43 -45 R 3
## 44 -45 P 2
## 45 -45 N 1
## 46 -44 N 1
## 47 -44 N 1
## 48 -44 W 4
## 49 -44 P 2
## 50 -43 R 3
## 51 -43 N 1
## 52 -43 P 2
## 53 -43 P 2
## 54 -43 P 2
## 55 -43 P 2
## 56 -43 P 2
## 57 -42 P 2
## 58 -42 P 2
## 59 -42 P 2
## 60 -42 W 4
## 61 -41 N 1
## 62 -41 P 2
## 63 -41 N 1
## 64 -40 P 2
## 65 -40 N 1
## 66 -39 N 1
## 67 -39 P 2
## 68 -39 P 2
## 69 -39 P 2
## 70 -39 P 2
## 71 -39 N 1
## 72 -39 P 2
## 73 -39 W 4
## 74 -38 P 2
## 75 -38 P 2
## 76 -38 P 2
## 77 -38 P 2
## 78 -38 P 2
## 79 -37 P 2
## 80 -37 N 1
## 81 -37 R 3
## 82 -37 P 2
## 83 -37 P 2
## 84 -37 P 2
## 85 -37 R 3
## 86 -37 P 2
## 87 -36 P 2
## 88 -36 N 1
## 89 -36 P 2
## 90 -36 P 2
## 91 -36 R 3
## 92 -36 P 2
## 93 -36 P 2
## 94 -36 R 3
## 95 -36 P 2
## 96 -36 P 2
## 97 -35 P 2
## 98 -35 P 2
## 99 -35 W 4
## 100 -35 R 3
## 101 -35 W 4
## 102 -35 N 1
## 103 -35 R 3
## 104 -35 P 2
## 105 -35 N 1
## 106 -34 P 2
## 107 -34 P 2
## 108 -34 P 2
## 109 -34 N 1
## 110 -34 N 1
## 111 -34 P 2
## 112 -34 W 4
## 113 -33 N 1
## 114 -33 P 2
## 115 -33 P 2
## 116 -33 P 2
## 117 -33 P 2
## 118 -33 P 2
## 119 -33 P 2
## 120 -32 W 4
## 121 -32 P 2
## 122 -32 P 2
## 123 -32 N 1
## 124 -32 P 2
## 125 -32 P 2
## 126 -32 P 2
## 127 -32 N 1
## 128 -31 R 3
## 129 -31 N 1
## 130 -31 N 1
## 131 -31 N 1
## 132 -31 P 2
## 133 -31 P 2
## 134 -31 P 2
## 135 -31 P 2
## 136 -31 R 3
## 137 -30 N 1
## 138 -30 P 2
## 139 -30 P 2
## 140 -30 P 2
## 141 -30 P 2
## 142 -30 P 2
## 143 -30 P 2
## 144 -30 P 2
## 145 -30 R 3
## 146 -30 P 2
## 147 -29 P 2
## 148 -29 P 2
## 149 -29 P 2
## 150 -29 N 1
## 151 -29 P 2
## 152 -29 P 2
## 153 -29 R 3
## 154 -29 P 2
## 155 -29 N 1
## 156 -29 R 3
## 157 -29 P 2
## 158 -29 P 2
## 159 -29 N 1
## 160 -28 P 2
## 161 -28 P 2
## 162 -28 P 2
## 163 -28 P 2
## 164 -28 W 4
## 165 -28 P 2
## 166 -28 N 1
## 167 -28 N 1
## 168 -28 P 2
## 169 -28 N 1
## 170 -27 N 1
## 171 -27 N 1
## 172 -27 P 2
## 173 -27 P 2
## 174 -27 P 2
## 175 -26 P 2
## 176 -26 P 2
## 177 -26 W 4
## 178 -26 P 2
## 179 -26 P 2
## 180 -26 P 2
## 181 -26 N 1
## 182 -26 P 2
## 183 -26 P 2
## 184 -26 N 1
## 185 -26 N 1
## 186 -25 P 2
## 187 -25 P 2
## 188 -25 N 1
## 189 -25 P 2
## 190 -25 P 2
## 191 -25 N 1
## 192 -25 P 2
## 193 -25 P 2
## 194 -25 N 1
## 195 -25 N 1
## 196 -25 P 2
## 197 -24 N 1
## 198 -24 P 2
## 199 -24 P 2
## 200 -24 P 2
## 201 -24 P 2
## 202 -24 P 2
## 203 -24 N 1
## 204 -24 P 2
## 205 -24 P 2
## 206 -24 P 2
## 207 -24 P 2
## 208 -23 P 2
## 209 -23 P 2
## 210 -23 P 2
## 211 -23 P 2
## 212 -23 N 1
## 213 -23 P 2
## 214 -23 R 3
## 215 -23 P 2
## 216 -23 N 1
## 217 -23 N 1
## 218 -23 N 1
## 219 -23 P 2
## 220 -23 P 2
## 221 -22 N 1
## 222 -22 P 2
## 223 -22 N 1
## 224 -22 P 2
## 225 -22 P 2
## 226 -22 P 2
## 227 -22 P 2
## 228 -22 P 2
## 229 -22 N 1
## 230 -22 N 1
## 231 -21 P 2
## 232 -21 P 2
## 233 -21 P 2
## 234 -21 P 2
## 235 -21 P 2
## 236 -21 P 2
## 237 -20 P 2
## 238 -20 P 2
## 239 -20 P 2
## 240 -20 N 1
## 241 -20 P 2
## 242 -20 N 1
## 243 -20 P 2
## 244 -20 P 2
## 245 -20 R 3
## 246 -19 N 1
## 247 -19 N 1
## 248 -19 N 1
## 249 -19 P 2
## 250 -19 P 2
## 251 -19 N 1
## 252 -19 P 2
## 253 -18 P 2
## 254 -18 P 2
## 255 -18 P 2
## 256 -18 P 2
## 257 -18 P 2
## 258 -18 N 1
## 259 -18 P 2
## 260 -18 P 2
## 261 -18 W 4
## 262 -18 N 1
## 263 -18 N 1
## 264 -18 R 3
## 265 -18 P 2
## 266 -17 P 2
## 267 -17 P 2
## 268 -17 P 2
## 269 -17 P 2
## 270 -17 P 2
## 271 -17 P 2
## 272 -17 P 2
## 273 -17 R 3
## 274 -17 N 1
## 275 -17 P 2
## 276 -17 P 2
## 277 -17 P 2
## 278 -16 N 1
## 279 -16 N 1
## 280 -16 P 2
## 281 -16 W 4
## 282 -16 N 1
## 283 -15 N 1
## 284 -15 P 2
## 285 -15 P 2
## 286 -15 N 1
## 287 -15 N 1
## 288 -15 N 1
## 289 -14 N 1
## 290 -14 N 1
## 291 -14 P 2
## 292 -14 P 2
## 293 -14 P 2
## 294 -14 N 1
## 295 -14 W 4
## 296 -14 N 1
## 297 -14 N 1
## 298 -14 P 2
## 299 -14 P 2
## 300 -14 N 1
## 301 -13 P 2
## 302 -13 N 1
## 303 -13 N 1
## 304 -13 P 2
## 305 -13 P 2
## 306 -13 N 1
## 307 -13 P 2
## 308 -13 P 2
## 309 -13 R 3
## 310 -13 P 2
## 311 -13 P 2
## 312 -13 N 1
## 313 -13 P 2
## 314 -12 N 1
## 315 -12 P 2
## 316 -12 P 2
## 317 -12 P 2
## 318 -12 P 2
## 319 -12 N 1
## 320 -12 N 1
## 321 -12 P 2
## 322 -12 N 1
## 323 -12 P 2
## 324 -11 P 2
## 325 -11 W 4
## 326 -11 N 1
## 327 -11 N 1
## 328 -11 N 1
## 329 -11 P 2
## 330 -11 R 3
## 331 -11 P 2
## 332 -11 P 2
## 333 -11 N 1
## 334 -11 W 4
## 335 -10 R 3
## 336 -10 P 2
## 337 -10 N 1
## 338 -10 W 4
## 339 -10 R 3
## 340 -10 P 2
## 341 -10 N 1
## 342 -10 N 1
## 343 -10 P 2
## 344 -10 P 2
## 345 -10 P 2
## 346 -10 P 2
## 347 -10 P 2
## 348 -10 P 2
## 349 -10 R 3
## 350 -10 P 2
## 351 -9 N 1
## 352 -9 P 2
## 353 -9 P 2
## 354 -9 P 2
## 355 -9 R 3
## 356 -9 P 2
## 357 -9 W 4
## 358 -8 R 3
## 359 -8 P 2
## 360 -8 N 1
## 361 -8 P 2
## 362 -8 N 1
## 363 -8 P 2
## 364 -8 N 1
## 365 -8 N 1
## 366 -8 P 2
## 367 -7 P 2
## 368 -7 P 2
## 369 -7 P 2
## 370 -7 P 2
## 371 -7 R 3
## 372 -7 P 2
## 373 -7 W 4
## 374 -7 P 2
## 375 -7 P 2
## 376 -6 P 2
## 377 -6 N 1
## 378 -6 N 1
## 379 -6 N 1
## 380 -6 N 1
## 381 -6 P 2
## 382 -6 W 4
## 383 -6 R 3
## 384 -6 P 2
## 385 -6 P 2
## 386 -6 P 2
## 387 -6 N 1
## 388 -5 N 1
## 389 -5 W 4
## 390 -5 N 1
## 391 -5 P 2
## 392 -5 N 1
## 393 -5 N 1
## 394 -5 N 1
## 395 -5 N 1
## 396 -5 P 2
## 397 -4 R 3
## 398 -4 N 1
## 399 -4 N 1
## 400 -4 P 2
## 401 -4 N 1
## 402 -4 N 1
## 403 -4 N 1
## 404 -4 N 1
## 405 -4 P 2
## 406 -4 P 2
## 407 -4 P 2
## 408 -4 N 1
## 409 -4 P 2
## 410 -4 N 1
## 411 -4 P 2
## 412 -4 R 3
## 413 -4 R 3
## 414 -4 R 3
## 415 -4 N 1
## 416 -4 P 2
## 417 -4 P 2
## 418 -3 P 2
## 419 -3 P 2
## 420 -3 P 2
## 421 -3 P 2
## 422 -3 W 4
## 423 -3 P 2
## 424 -3 P 2
## 425 -3 W 4
## 426 -3 N 1
## 427 -3 W 4
## 428 -3 P 2
## 429 -3 N 1
## 430 -3 N 1
## 431 -3 W 4
## 432 -3 N 1
## 433 -3 P 2
## 434 -3 P 2
## 435 -3 N 1
## 436 -3 N 1
## 437 -3 P 2
## 438 -3 P 2
## 439 -3 P 2
## 440 -3 P 2
## 441 -3 P 2
## 442 -3 N 1
## 443 -3 P 2
## 444 -2 P 2
## 445 -2 N 1
## 446 -2 N 1
## 447 -2 P 2
## 448 -2 W 4
## 449 -2 R 3
## 450 -2 N 1
## 451 -2 P 2
## 452 -2 W 4
## 453 -2 R 3
## 454 -2 N 1
## 455 -2 P 2
## 456 -2 P 2
## 457 -2 P 2
## 458 -2 P 2
## 459 -2 P 2
## 460 -2 P 2
## 461 -2 P 2
## 462 -2 N 1
## 463 -2 N 1
## 464 -1 P 2
## 465 -1 N 1
## 466 -1 N 1
## 467 -1 P 2
## 468 -1 P 2
## 469 -1 W 4
## 470 -1 W 4
## 471 -1 R 3
## 472 -1 P 2
## 473 -1 N 1
## 474 -1 N 1
## 475 -1 N 1
## 476 -1 N 1
## 477 -1 P 2
## 478 -1 P 2
## 479 -1 N 1
## 480 -1 P 2
## 481 -1 W 4
## 482 -1 P 2
## 483 -1 W 4
## 484 0 N 1
## 485 0 N 1
## 486 0 P 2
## 487 0 P 2
## 488 0 N 1
## 489 0 P 2
## 490 0 P 2
## 491 0 P 2
## 492 0 P 2
## 493 0 W 4
## 494 0 P 2
## 495 0 P 2
## 496 0 P 2
## 497 0 P 2
## 498 0 P 2
## 499 0 N 1
## 500 1 P 2
## 501 1 P 2
## 502 1 P 2
## 503 1 P 2
## 504 1 N 1
## 505 1 W 4
## 506 1 P 2
## 507 1 N 1
## 508 1 P 2
## 509 1 P 2
## 510 2 P 2
## 511 2 N 1
## 512 2 P 2
## 513 2 P 2
## 514 2 W 4
## 515 2 P 2
## 516 2 P 2
## 517 2 P 2
## 518 2 P 2
## 519 2 P 2
## 520 2 N 1
## 521 2 P 2
## 522 3 R 3
## 523 3 P 2
## 524 3 P 2
## 525 3 N 1
## 526 3 N 1
## 527 3 N 1
## 528 3 P 2
## 529 3 P 2
## 530 3 P 2
## 531 3 N 1
## 532 3 N 1
## 533 3 P 2
## 534 3 P 2
## 535 3 R 3
## 536 4 P 2
## 537 4 N 1
## 538 4 W 4
## 539 4 P 2
## 540 4 N 1
## 541 4 N 1
## 542 4 P 2
## 543 4 P 2
## 544 5 P 2
## 545 5 P 2
## 546 5 P 2
## 547 5 N 1
## 548 5 P 2
## 549 5 P 2
## 550 6 R 3
## 551 6 P 2
## 552 6 N 1
## 553 6 P 2
## 554 6 P 2
## 555 6 W 4
## 556 6 N 1
## 557 6 N 1
## 558 7 P 2
## 559 7 N 1
## 560 7 P 2
## 561 7 W 4
## 562 7 N 1
## 563 8 P 2
## 564 10 P 2
## 565 10 P 2
##
## $model
## [1] "Nest"
##
## $mixtures
## [1] 1
##
## $freq
## HabitatN HabitatP HabitatR HabitatW
## 1 0 0 1 0
## 2 0 1 0 0
## 3 0 1 0 0
## 4 1 0 0 0
## 5 0 1 0 0
## 6 0 1 0 0
## 7 0 1 0 0
## 8 1 0 0 0
## 9 0 1 0 0
## 10 1 0 0 0
## 11 0 1 0 0
## 12 0 1 0 0
## 13 0 1 0 0
## 14 1 0 0 0
## 15 1 0 0 0
## 16 0 1 0 0
## 17 0 1 0 0
## 18 0 0 0 1
## 19 1 0 0 0
## 20 0 1 0 0
## 21 1 0 0 0
## 22 0 1 0 0
## 23 1 0 0 0
## 24 1 0 0 0
## 25 0 0 1 0
## 26 0 1 0 0
## 27 1 0 0 0
## 28 0 0 1 0
## 29 0 0 1 0
## 30 0 1 0 0
## 31 1 0 0 0
## 32 1 0 0 0
## 33 0 1 0 0
## 34 0 1 0 0
## 35 0 1 0 0
## 36 1 0 0 0
## 37 0 1 0 0
## 38 1 0 0 0
## 39 0 1 0 0
## 40 1 0 0 0
## 41 0 1 0 0
## 42 1 0 0 0
## 43 0 0 1 0
## 44 0 1 0 0
## 45 1 0 0 0
## 46 1 0 0 0
## 47 1 0 0 0
## 48 0 0 0 1
## 49 0 1 0 0
## 50 0 0 1 0
## 51 1 0 0 0
## 52 0 1 0 0
## 53 0 1 0 0
## 54 0 1 0 0
## 55 0 1 0 0
## 56 0 1 0 0
## 57 0 1 0 0
## 58 0 1 0 0
## 59 0 1 0 0
## 60 0 0 0 1
## 61 1 0 0 0
## 62 0 1 0 0
## 63 1 0 0 0
## 64 0 1 0 0
## 65 1 0 0 0
## 66 1 0 0 0
## 67 0 1 0 0
## 68 0 1 0 0
## 69 0 1 0 0
## 70 0 1 0 0
## 71 1 0 0 0
## 72 0 1 0 0
## 73 0 0 0 1
## 74 0 1 0 0
## 75 0 1 0 0
## 76 0 1 0 0
## 77 0 1 0 0
## 78 0 1 0 0
## 79 0 1 0 0
## 80 1 0 0 0
## 81 0 0 1 0
## 82 0 1 0 0
## 83 0 1 0 0
## 84 0 1 0 0
## 85 0 0 1 0
## 86 0 1 0 0
## 87 0 1 0 0
## 88 1 0 0 0
## 89 0 1 0 0
## 90 0 1 0 0
## 91 0 0 1 0
## 92 0 1 0 0
## 93 0 1 0 0
## 94 0 0 1 0
## 95 0 1 0 0
## 96 0 1 0 0
## 97 0 1 0 0
## 98 0 1 0 0
## 99 0 0 0 1
## 100 0 0 1 0
## 101 0 0 0 1
## 102 1 0 0 0
## 103 0 0 1 0
## 104 0 1 0 0
## 105 1 0 0 0
## 106 0 1 0 0
## 107 0 1 0 0
## 108 0 1 0 0
## 109 1 0 0 0
## 110 1 0 0 0
## 111 0 1 0 0
## 112 0 0 0 1
## 113 1 0 0 0
## 114 0 1 0 0
## 115 0 1 0 0
## 116 0 1 0 0
## 117 0 1 0 0
## 118 0 1 0 0
## 119 0 1 0 0
## 120 0 0 0 1
## 121 0 1 0 0
## 122 0 1 0 0
## 123 1 0 0 0
## 124 0 1 0 0
## 125 0 1 0 0
## 126 0 1 0 0
## 127 1 0 0 0
## 128 0 0 1 0
## 129 1 0 0 0
## 130 1 0 0 0
## 131 1 0 0 0
## 132 0 1 0 0
## 133 0 1 0 0
## 134 0 1 0 0
## 135 0 1 0 0
## 136 0 0 1 0
## 137 1 0 0 0
## 138 0 1 0 0
## 139 0 1 0 0
## 140 0 1 0 0
## 141 0 1 0 0
## 142 0 1 0 0
## 143 0 1 0 0
## 144 0 1 0 0
## 145 0 0 1 0
## 146 0 1 0 0
## 147 0 1 0 0
## 148 0 1 0 0
## 149 0 1 0 0
## 150 1 0 0 0
## 151 0 1 0 0
## 152 0 1 0 0
## 153 0 0 1 0
## 154 0 1 0 0
## 155 1 0 0 0
## 156 0 0 1 0
## 157 0 1 0 0
## 158 0 1 0 0
## 159 1 0 0 0
## 160 0 1 0 0
## 161 0 1 0 0
## 162 0 1 0 0
## 163 0 1 0 0
## 164 0 0 0 1
## 165 0 1 0 0
## 166 1 0 0 0
## 167 1 0 0 0
## 168 0 1 0 0
## 169 1 0 0 0
## 170 1 0 0 0
## 171 1 0 0 0
## 172 0 1 0 0
## 173 0 1 0 0
## 174 0 1 0 0
## 175 0 1 0 0
## 176 0 1 0 0
## 177 0 0 0 1
## 178 0 1 0 0
## 179 0 1 0 0
## 180 0 1 0 0
## 181 1 0 0 0
## 182 0 1 0 0
## 183 0 1 0 0
## 184 1 0 0 0
## 185 1 0 0 0
## 186 0 1 0 0
## 187 0 1 0 0
## 188 1 0 0 0
## 189 0 1 0 0
## 190 0 1 0 0
## 191 1 0 0 0
## 192 0 1 0 0
## 193 0 1 0 0
## 194 1 0 0 0
## 195 1 0 0 0
## 196 0 1 0 0
## 197 1 0 0 0
## 198 0 1 0 0
## 199 0 1 0 0
## 200 0 1 0 0
## 201 0 1 0 0
## 202 0 1 0 0
## 203 1 0 0 0
## 204 0 1 0 0
## 205 0 1 0 0
## 206 0 1 0 0
## 207 0 1 0 0
## 208 0 1 0 0
## 209 0 1 0 0
## 210 0 1 0 0
## 211 0 1 0 0
## 212 1 0 0 0
## 213 0 1 0 0
## 214 0 0 1 0
## 215 0 1 0 0
## 216 1 0 0 0
## 217 1 0 0 0
## 218 1 0 0 0
## 219 0 1 0 0
## 220 0 1 0 0
## 221 1 0 0 0
## 222 0 1 0 0
## 223 1 0 0 0
## 224 0 1 0 0
## 225 0 1 0 0
## 226 0 1 0 0
## 227 0 1 0 0
## 228 0 1 0 0
## 229 1 0 0 0
## 230 1 0 0 0
## 231 0 1 0 0
## 232 0 1 0 0
## 233 0 1 0 0
## 234 0 1 0 0
## 235 0 1 0 0
## 236 0 1 0 0
## 237 0 1 0 0
## 238 0 1 0 0
## 239 0 1 0 0
## 240 1 0 0 0
## 241 0 1 0 0
## 242 1 0 0 0
## 243 0 1 0 0
## 244 0 1 0 0
## 245 0 0 1 0
## 246 1 0 0 0
## 247 1 0 0 0
## 248 1 0 0 0
## 249 0 1 0 0
## 250 0 1 0 0
## 251 1 0 0 0
## 252 0 1 0 0
## 253 0 1 0 0
## 254 0 1 0 0
## 255 0 1 0 0
## 256 0 1 0 0
## 257 0 1 0 0
## 258 1 0 0 0
## 259 0 1 0 0
## 260 0 1 0 0
## 261 0 0 0 1
## 262 1 0 0 0
## 263 1 0 0 0
## 264 0 0 1 0
## 265 0 1 0 0
## 266 0 1 0 0
## 267 0 1 0 0
## 268 0 1 0 0
## 269 0 1 0 0
## 270 0 1 0 0
## 271 0 1 0 0
## 272 0 1 0 0
## 273 0 0 1 0
## 274 1 0 0 0
## 275 0 1 0 0
## 276 0 1 0 0
## 277 0 1 0 0
## 278 1 0 0 0
## 279 1 0 0 0
## 280 0 1 0 0
## 281 0 0 0 1
## 282 1 0 0 0
## 283 1 0 0 0
## 284 0 1 0 0
## 285 0 1 0 0
## 286 1 0 0 0
## 287 1 0 0 0
## 288 1 0 0 0
## 289 1 0 0 0
## 290 1 0 0 0
## 291 0 1 0 0
## 292 0 1 0 0
## 293 0 1 0 0
## 294 1 0 0 0
## 295 0 0 0 1
## 296 1 0 0 0
## 297 1 0 0 0
## 298 0 1 0 0
## 299 0 1 0 0
## 300 1 0 0 0
## 301 0 1 0 0
## 302 1 0 0 0
## 303 1 0 0 0
## 304 0 1 0 0
## 305 0 1 0 0
## 306 1 0 0 0
## 307 0 1 0 0
## 308 0 1 0 0
## 309 0 0 1 0
## 310 0 1 0 0
## 311 0 1 0 0
## 312 1 0 0 0
## 313 0 1 0 0
## 314 1 0 0 0
## 315 0 1 0 0
## 316 0 1 0 0
## 317 0 1 0 0
## 318 0 1 0 0
## 319 1 0 0 0
## 320 1 0 0 0
## 321 0 1 0 0
## 322 1 0 0 0
## 323 0 1 0 0
## 324 0 1 0 0
## 325 0 0 0 1
## 326 1 0 0 0
## 327 1 0 0 0
## 328 1 0 0 0
## 329 0 1 0 0
## 330 0 0 1 0
## 331 0 1 0 0
## 332 0 1 0 0
## 333 1 0 0 0
## 334 0 0 0 1
## 335 0 0 1 0
## 336 0 1 0 0
## 337 1 0 0 0
## 338 0 0 0 1
## 339 0 0 1 0
## 340 0 1 0 0
## 341 1 0 0 0
## 342 1 0 0 0
## 343 0 1 0 0
## 344 0 1 0 0
## 345 0 1 0 0
## 346 0 1 0 0
## 347 0 1 0 0
## 348 0 1 0 0
## 349 0 0 1 0
## 350 0 1 0 0
## 351 1 0 0 0
## 352 0 1 0 0
## 353 0 1 0 0
## 354 0 1 0 0
## 355 0 0 1 0
## 356 0 1 0 0
## 357 0 0 0 1
## 358 0 0 1 0
## 359 0 1 0 0
## 360 1 0 0 0
## 361 0 1 0 0
## 362 1 0 0 0
## 363 0 1 0 0
## 364 1 0 0 0
## 365 1 0 0 0
## 366 0 1 0 0
## 367 0 1 0 0
## 368 0 1 0 0
## 369 0 1 0 0
## 370 0 1 0 0
## 371 0 0 1 0
## 372 0 1 0 0
## 373 0 0 0 1
## 374 0 1 0 0
## 375 0 1 0 0
## 376 0 1 0 0
## 377 1 0 0 0
## 378 1 0 0 0
## 379 1 0 0 0
## 380 1 0 0 0
## 381 0 1 0 0
## 382 0 0 0 1
## 383 0 0 1 0
## 384 0 1 0 0
## 385 0 1 0 0
## 386 0 1 0 0
## 387 1 0 0 0
## 388 1 0 0 0
## 389 0 0 0 1
## 390 1 0 0 0
## 391 0 1 0 0
## 392 1 0 0 0
## 393 1 0 0 0
## 394 1 0 0 0
## 395 1 0 0 0
## 396 0 1 0 0
## 397 0 0 1 0
## 398 1 0 0 0
## 399 1 0 0 0
## 400 0 1 0 0
## 401 1 0 0 0
## 402 1 0 0 0
## 403 1 0 0 0
## 404 1 0 0 0
## 405 0 1 0 0
## 406 0 1 0 0
## 407 0 1 0 0
## 408 1 0 0 0
## 409 0 1 0 0
## 410 1 0 0 0
## 411 0 1 0 0
## 412 0 0 1 0
## 413 0 0 1 0
## 414 0 0 1 0
## 415 1 0 0 0
## 416 0 1 0 0
## 417 0 1 0 0
## 418 0 1 0 0
## 419 0 1 0 0
## 420 0 1 0 0
## 421 0 1 0 0
## 422 0 0 0 1
## 423 0 1 0 0
## 424 0 1 0 0
## 425 0 0 0 1
## 426 1 0 0 0
## 427 0 0 0 1
## 428 0 1 0 0
## 429 1 0 0 0
## 430 1 0 0 0
## 431 0 0 0 1
## 432 1 0 0 0
## 433 0 1 0 0
## 434 0 1 0 0
## 435 1 0 0 0
## 436 1 0 0 0
## 437 0 1 0 0
## 438 0 1 0 0
## 439 0 1 0 0
## 440 0 1 0 0
## 441 0 1 0 0
## 442 1 0 0 0
## 443 0 1 0 0
## 444 0 1 0 0
## 445 1 0 0 0
## 446 1 0 0 0
## 447 0 1 0 0
## 448 0 0 0 1
## 449 0 0 1 0
## 450 1 0 0 0
## 451 0 1 0 0
## 452 0 0 0 1
## 453 0 0 1 0
## 454 1 0 0 0
## 455 0 1 0 0
## 456 0 1 0 0
## 457 0 1 0 0
## 458 0 1 0 0
## 459 0 1 0 0
## 460 0 1 0 0
## 461 0 1 0 0
## 462 1 0 0 0
## 463 1 0 0 0
## 464 0 1 0 0
## 465 1 0 0 0
## 466 1 0 0 0
## 467 0 1 0 0
## 468 0 1 0 0
## 469 0 0 0 1
## 470 0 0 0 1
## 471 0 0 1 0
## 472 0 1 0 0
## 473 1 0 0 0
## 474 1 0 0 0
## 475 1 0 0 0
## 476 1 0 0 0
## 477 0 1 0 0
## 478 0 1 0 0
## 479 1 0 0 0
## 480 0 1 0 0
## 481 0 0 0 1
## 482 0 1 0 0
## 483 0 0 0 1
## 484 1 0 0 0
## 485 1 0 0 0
## 486 0 1 0 0
## 487 0 1 0 0
## 488 1 0 0 0
## 489 0 1 0 0
## 490 0 1 0 0
## 491 0 1 0 0
## 492 0 1 0 0
## 493 0 0 0 1
## 494 0 1 0 0
## 495 0 1 0 0
## 496 0 1 0 0
## 497 0 1 0 0
## 498 0 1 0 0
## 499 1 0 0 0
## 500 0 1 0 0
## 501 0 1 0 0
## 502 0 1 0 0
## 503 0 1 0 0
## 504 1 0 0 0
## 505 0 0 0 1
## 506 0 1 0 0
## 507 1 0 0 0
## 508 0 1 0 0
## 509 0 1 0 0
## 510 0 1 0 0
## 511 1 0 0 0
## 512 0 1 0 0
## 513 0 1 0 0
## 514 0 0 0 1
## 515 0 1 0 0
## 516 0 1 0 0
## 517 0 1 0 0
## 518 0 1 0 0
## 519 0 1 0 0
## 520 1 0 0 0
## 521 0 1 0 0
## 522 0 0 1 0
## 523 0 1 0 0
## 524 0 1 0 0
## 525 1 0 0 0
## 526 1 0 0 0
## 527 1 0 0 0
## 528 0 1 0 0
## 529 0 1 0 0
## 530 0 1 0 0
## 531 1 0 0 0
## 532 1 0 0 0
## 533 0 1 0 0
## 534 0 1 0 0
## 535 0 0 1 0
## 536 0 1 0 0
## 537 1 0 0 0
## 538 0 0 0 1
## 539 0 1 0 0
## 540 1 0 0 0
## 541 1 0 0 0
## 542 0 1 0 0
## 543 0 1 0 0
## 544 0 1 0 0
## 545 0 1 0 0
## 546 0 1 0 0
## 547 1 0 0 0
## 548 0 1 0 0
## 549 0 1 0 0
## 550 0 0 1 0
## 551 0 1 0 0
## 552 1 0 0 0
## 553 0 1 0 0
## 554 0 1 0 0
## 555 0 0 0 1
## 556 1 0 0 0
## 557 1 0 0 0
## 558 0 1 0 0
## 559 1 0 0 0
## 560 0 1 0 0
## 561 0 0 0 1
## 562 1 0 0 0
## 563 0 1 0 0
## 564 0 1 0 0
## 565 0 1 0 0
##
## $nocc
## [1] 90
##
## $nocc.secondary
## NULL
##
## $time.intervals
## [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
## [36] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
## [71] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##
## $begin.time
## [1] 1
##
## $age.unit
## [1] 1
##
## $initial.ages
## [1] 0 0 0 0
##
## $group.covariates
## Habitat
## 1 N
## 2 P
## 3 R
## 4 W
##
## $nstrata
## [1] 1
##
## $strata.labels
## [1] ""
##
## $counts
## NULL
##
## $reverse
## [1] FALSE
##
## $areas
## NULL
##
## $events
## NULL
# 2. Examine and/or modify the ddl. (Not done here)
mall.ddl <- make.design.data(mall.proc)
str(mall.ddl)
## List of 2
## $ S :'data.frame': 356 obs. of 8 variables:
## ..$ par.index : int [1:356] 1 2 3 4 5 6 7 8 9 10 ...
## ..$ model.index: num [1:356] 1 2 3 4 5 6 7 8 9 10 ...
## ..$ group : Factor w/ 4 levels "N","P","R","W": 1 1 1 1 1 1 1 1 1 1 ...
## ..$ age : Factor w/ 89 levels "0","1","2","3",..: 1 2 3 4 5 6 7 8 9 10 ...
## ..$ time : Factor w/ 89 levels "1","2","3","4",..: 1 2 3 4 5 6 7 8 9 10 ...
## ..$ Age : num [1:356] 0 1 2 3 4 5 6 7 8 9 ...
## ..$ Time : num [1:356] 0 1 2 3 4 5 6 7 8 9 ...
## ..$ Habitat : Factor w/ 4 levels "N","P","R","W": 1 1 1 1 1 1 1 1 1 1 ...
## $ pimtypes:List of 1
## ..$ S:List of 1
## .. ..$ pim.type: chr "all"
mall.ddl
## $S
## par.index model.index group age time Age Time Habitat
## 1 1 1 N 0 1 0 0 N
## 2 2 2 N 1 2 1 1 N
## 3 3 3 N 2 3 2 2 N
## 4 4 4 N 3 4 3 3 N
## 5 5 5 N 4 5 4 4 N
## 6 6 6 N 5 6 5 5 N
## 7 7 7 N 6 7 6 6 N
## 8 8 8 N 7 8 7 7 N
## 9 9 9 N 8 9 8 8 N
## 10 10 10 N 9 10 9 9 N
## 11 11 11 N 10 11 10 10 N
## 12 12 12 N 11 12 11 11 N
## 13 13 13 N 12 13 12 12 N
## 14 14 14 N 13 14 13 13 N
## 15 15 15 N 14 15 14 14 N
## 16 16 16 N 15 16 15 15 N
## 17 17 17 N 16 17 16 16 N
## 18 18 18 N 17 18 17 17 N
## 19 19 19 N 18 19 18 18 N
## 20 20 20 N 19 20 19 19 N
## 21 21 21 N 20 21 20 20 N
## 22 22 22 N 21 22 21 21 N
## 23 23 23 N 22 23 22 22 N
## 24 24 24 N 23 24 23 23 N
## 25 25 25 N 24 25 24 24 N
## 26 26 26 N 25 26 25 25 N
## 27 27 27 N 26 27 26 26 N
## 28 28 28 N 27 28 27 27 N
## 29 29 29 N 28 29 28 28 N
## 30 30 30 N 29 30 29 29 N
## 31 31 31 N 30 31 30 30 N
## 32 32 32 N 31 32 31 31 N
## 33 33 33 N 32 33 32 32 N
## 34 34 34 N 33 34 33 33 N
## 35 35 35 N 34 35 34 34 N
## 36 36 36 N 35 36 35 35 N
## 37 37 37 N 36 37 36 36 N
## 38 38 38 N 37 38 37 37 N
## 39 39 39 N 38 39 38 38 N
## 40 40 40 N 39 40 39 39 N
## 41 41 41 N 40 41 40 40 N
## 42 42 42 N 41 42 41 41 N
## 43 43 43 N 42 43 42 42 N
## 44 44 44 N 43 44 43 43 N
## 45 45 45 N 44 45 44 44 N
## 46 46 46 N 45 46 45 45 N
## 47 47 47 N 46 47 46 46 N
## 48 48 48 N 47 48 47 47 N
## 49 49 49 N 48 49 48 48 N
## 50 50 50 N 49 50 49 49 N
## 51 51 51 N 50 51 50 50 N
## 52 52 52 N 51 52 51 51 N
## 53 53 53 N 52 53 52 52 N
## 54 54 54 N 53 54 53 53 N
## 55 55 55 N 54 55 54 54 N
## 56 56 56 N 55 56 55 55 N
## 57 57 57 N 56 57 56 56 N
## 58 58 58 N 57 58 57 57 N
## 59 59 59 N 58 59 58 58 N
## 60 60 60 N 59 60 59 59 N
## 61 61 61 N 60 61 60 60 N
## 62 62 62 N 61 62 61 61 N
## 63 63 63 N 62 63 62 62 N
## 64 64 64 N 63 64 63 63 N
## 65 65 65 N 64 65 64 64 N
## 66 66 66 N 65 66 65 65 N
## 67 67 67 N 66 67 66 66 N
## 68 68 68 N 67 68 67 67 N
## 69 69 69 N 68 69 68 68 N
## 70 70 70 N 69 70 69 69 N
## 71 71 71 N 70 71 70 70 N
## 72 72 72 N 71 72 71 71 N
## 73 73 73 N 72 73 72 72 N
## 74 74 74 N 73 74 73 73 N
## 75 75 75 N 74 75 74 74 N
## 76 76 76 N 75 76 75 75 N
## 77 77 77 N 76 77 76 76 N
## 78 78 78 N 77 78 77 77 N
## 79 79 79 N 78 79 78 78 N
## 80 80 80 N 79 80 79 79 N
## 81 81 81 N 80 81 80 80 N
## 82 82 82 N 81 82 81 81 N
## 83 83 83 N 82 83 82 82 N
## 84 84 84 N 83 84 83 83 N
## 85 85 85 N 84 85 84 84 N
## 86 86 86 N 85 86 85 85 N
## 87 87 87 N 86 87 86 86 N
## 88 88 88 N 87 88 87 87 N
## 89 89 89 N 88 89 88 88 N
## 90 90 90 P 0 1 0 0 P
## 91 91 91 P 1 2 1 1 P
## 92 92 92 P 2 3 2 2 P
## 93 93 93 P 3 4 3 3 P
## 94 94 94 P 4 5 4 4 P
## 95 95 95 P 5 6 5 5 P
## 96 96 96 P 6 7 6 6 P
## 97 97 97 P 7 8 7 7 P
## 98 98 98 P 8 9 8 8 P
## 99 99 99 P 9 10 9 9 P
## 100 100 100 P 10 11 10 10 P
## 101 101 101 P 11 12 11 11 P
## 102 102 102 P 12 13 12 12 P
## 103 103 103 P 13 14 13 13 P
## 104 104 104 P 14 15 14 14 P
## 105 105 105 P 15 16 15 15 P
## 106 106 106 P 16 17 16 16 P
## 107 107 107 P 17 18 17 17 P
## 108 108 108 P 18 19 18 18 P
## 109 109 109 P 19 20 19 19 P
## 110 110 110 P 20 21 20 20 P
## 111 111 111 P 21 22 21 21 P
## 112 112 112 P 22 23 22 22 P
## 113 113 113 P 23 24 23 23 P
## 114 114 114 P 24 25 24 24 P
## 115 115 115 P 25 26 25 25 P
## 116 116 116 P 26 27 26 26 P
## 117 117 117 P 27 28 27 27 P
## 118 118 118 P 28 29 28 28 P
## 119 119 119 P 29 30 29 29 P
## 120 120 120 P 30 31 30 30 P
## 121 121 121 P 31 32 31 31 P
## 122 122 122 P 32 33 32 32 P
## 123 123 123 P 33 34 33 33 P
## 124 124 124 P 34 35 34 34 P
## 125 125 125 P 35 36 35 35 P
## 126 126 126 P 36 37 36 36 P
## 127 127 127 P 37 38 37 37 P
## 128 128 128 P 38 39 38 38 P
## 129 129 129 P 39 40 39 39 P
## 130 130 130 P 40 41 40 40 P
## 131 131 131 P 41 42 41 41 P
## 132 132 132 P 42 43 42 42 P
## 133 133 133 P 43 44 43 43 P
## 134 134 134 P 44 45 44 44 P
## 135 135 135 P 45 46 45 45 P
## 136 136 136 P 46 47 46 46 P
## 137 137 137 P 47 48 47 47 P
## 138 138 138 P 48 49 48 48 P
## 139 139 139 P 49 50 49 49 P
## 140 140 140 P 50 51 50 50 P
## 141 141 141 P 51 52 51 51 P
## 142 142 142 P 52 53 52 52 P
## 143 143 143 P 53 54 53 53 P
## 144 144 144 P 54 55 54 54 P
## 145 145 145 P 55 56 55 55 P
## 146 146 146 P 56 57 56 56 P
## 147 147 147 P 57 58 57 57 P
## 148 148 148 P 58 59 58 58 P
## 149 149 149 P 59 60 59 59 P
## 150 150 150 P 60 61 60 60 P
## 151 151 151 P 61 62 61 61 P
## 152 152 152 P 62 63 62 62 P
## 153 153 153 P 63 64 63 63 P
## 154 154 154 P 64 65 64 64 P
## 155 155 155 P 65 66 65 65 P
## 156 156 156 P 66 67 66 66 P
## 157 157 157 P 67 68 67 67 P
## 158 158 158 P 68 69 68 68 P
## 159 159 159 P 69 70 69 69 P
## 160 160 160 P 70 71 70 70 P
## 161 161 161 P 71 72 71 71 P
## 162 162 162 P 72 73 72 72 P
## 163 163 163 P 73 74 73 73 P
## 164 164 164 P 74 75 74 74 P
## 165 165 165 P 75 76 75 75 P
## 166 166 166 P 76 77 76 76 P
## 167 167 167 P 77 78 77 77 P
## 168 168 168 P 78 79 78 78 P
## 169 169 169 P 79 80 79 79 P
## 170 170 170 P 80 81 80 80 P
## 171 171 171 P 81 82 81 81 P
## 172 172 172 P 82 83 82 82 P
## 173 173 173 P 83 84 83 83 P
## 174 174 174 P 84 85 84 84 P
## 175 175 175 P 85 86 85 85 P
## 176 176 176 P 86 87 86 86 P
## 177 177 177 P 87 88 87 87 P
## 178 178 178 P 88 89 88 88 P
## 179 179 179 R 0 1 0 0 R
## 180 180 180 R 1 2 1 1 R
## 181 181 181 R 2 3 2 2 R
## 182 182 182 R 3 4 3 3 R
## 183 183 183 R 4 5 4 4 R
## 184 184 184 R 5 6 5 5 R
## 185 185 185 R 6 7 6 6 R
## 186 186 186 R 7 8 7 7 R
## 187 187 187 R 8 9 8 8 R
## 188 188 188 R 9 10 9 9 R
## 189 189 189 R 10 11 10 10 R
## 190 190 190 R 11 12 11 11 R
## 191 191 191 R 12 13 12 12 R
## 192 192 192 R 13 14 13 13 R
## 193 193 193 R 14 15 14 14 R
## 194 194 194 R 15 16 15 15 R
## 195 195 195 R 16 17 16 16 R
## 196 196 196 R 17 18 17 17 R
## 197 197 197 R 18 19 18 18 R
## 198 198 198 R 19 20 19 19 R
## 199 199 199 R 20 21 20 20 R
## 200 200 200 R 21 22 21 21 R
## 201 201 201 R 22 23 22 22 R
## 202 202 202 R 23 24 23 23 R
## 203 203 203 R 24 25 24 24 R
## 204 204 204 R 25 26 25 25 R
## 205 205 205 R 26 27 26 26 R
## 206 206 206 R 27 28 27 27 R
## 207 207 207 R 28 29 28 28 R
## 208 208 208 R 29 30 29 29 R
## 209 209 209 R 30 31 30 30 R
## 210 210 210 R 31 32 31 31 R
## 211 211 211 R 32 33 32 32 R
## 212 212 212 R 33 34 33 33 R
## 213 213 213 R 34 35 34 34 R
## 214 214 214 R 35 36 35 35 R
## 215 215 215 R 36 37 36 36 R
## 216 216 216 R 37 38 37 37 R
## 217 217 217 R 38 39 38 38 R
## 218 218 218 R 39 40 39 39 R
## 219 219 219 R 40 41 40 40 R
## 220 220 220 R 41 42 41 41 R
## 221 221 221 R 42 43 42 42 R
## 222 222 222 R 43 44 43 43 R
## 223 223 223 R 44 45 44 44 R
## 224 224 224 R 45 46 45 45 R
## 225 225 225 R 46 47 46 46 R
## 226 226 226 R 47 48 47 47 R
## 227 227 227 R 48 49 48 48 R
## 228 228 228 R 49 50 49 49 R
## 229 229 229 R 50 51 50 50 R
## 230 230 230 R 51 52 51 51 R
## 231 231 231 R 52 53 52 52 R
## 232 232 232 R 53 54 53 53 R
## 233 233 233 R 54 55 54 54 R
## 234 234 234 R 55 56 55 55 R
## 235 235 235 R 56 57 56 56 R
## 236 236 236 R 57 58 57 57 R
## 237 237 237 R 58 59 58 58 R
## 238 238 238 R 59 60 59 59 R
## 239 239 239 R 60 61 60 60 R
## 240 240 240 R 61 62 61 61 R
## 241 241 241 R 62 63 62 62 R
## 242 242 242 R 63 64 63 63 R
## 243 243 243 R 64 65 64 64 R
## 244 244 244 R 65 66 65 65 R
## 245 245 245 R 66 67 66 66 R
## 246 246 246 R 67 68 67 67 R
## 247 247 247 R 68 69 68 68 R
## 248 248 248 R 69 70 69 69 R
## 249 249 249 R 70 71 70 70 R
## 250 250 250 R 71 72 71 71 R
## 251 251 251 R 72 73 72 72 R
## 252 252 252 R 73 74 73 73 R
## 253 253 253 R 74 75 74 74 R
## 254 254 254 R 75 76 75 75 R
## 255 255 255 R 76 77 76 76 R
## 256 256 256 R 77 78 77 77 R
## 257 257 257 R 78 79 78 78 R
## 258 258 258 R 79 80 79 79 R
## 259 259 259 R 80 81 80 80 R
## 260 260 260 R 81 82 81 81 R
## 261 261 261 R 82 83 82 82 R
## 262 262 262 R 83 84 83 83 R
## 263 263 263 R 84 85 84 84 R
## 264 264 264 R 85 86 85 85 R
## 265 265 265 R 86 87 86 86 R
## 266 266 266 R 87 88 87 87 R
## 267 267 267 R 88 89 88 88 R
## 268 268 268 W 0 1 0 0 W
## 269 269 269 W 1 2 1 1 W
## 270 270 270 W 2 3 2 2 W
## 271 271 271 W 3 4 3 3 W
## 272 272 272 W 4 5 4 4 W
## 273 273 273 W 5 6 5 5 W
## 274 274 274 W 6 7 6 6 W
## 275 275 275 W 7 8 7 7 W
## 276 276 276 W 8 9 8 8 W
## 277 277 277 W 9 10 9 9 W
## 278 278 278 W 10 11 10 10 W
## 279 279 279 W 11 12 11 11 W
## 280 280 280 W 12 13 12 12 W
## 281 281 281 W 13 14 13 13 W
## 282 282 282 W 14 15 14 14 W
## 283 283 283 W 15 16 15 15 W
## 284 284 284 W 16 17 16 16 W
## 285 285 285 W 17 18 17 17 W
## 286 286 286 W 18 19 18 18 W
## 287 287 287 W 19 20 19 19 W
## 288 288 288 W 20 21 20 20 W
## 289 289 289 W 21 22 21 21 W
## 290 290 290 W 22 23 22 22 W
## 291 291 291 W 23 24 23 23 W
## 292 292 292 W 24 25 24 24 W
## 293 293 293 W 25 26 25 25 W
## 294 294 294 W 26 27 26 26 W
## 295 295 295 W 27 28 27 27 W
## 296 296 296 W 28 29 28 28 W
## 297 297 297 W 29 30 29 29 W
## 298 298 298 W 30 31 30 30 W
## 299 299 299 W 31 32 31 31 W
## 300 300 300 W 32 33 32 32 W
## 301 301 301 W 33 34 33 33 W
## 302 302 302 W 34 35 34 34 W
## 303 303 303 W 35 36 35 35 W
## 304 304 304 W 36 37 36 36 W
## 305 305 305 W 37 38 37 37 W
## 306 306 306 W 38 39 38 38 W
## 307 307 307 W 39 40 39 39 W
## 308 308 308 W 40 41 40 40 W
## 309 309 309 W 41 42 41 41 W
## 310 310 310 W 42 43 42 42 W
## 311 311 311 W 43 44 43 43 W
## 312 312 312 W 44 45 44 44 W
## 313 313 313 W 45 46 45 45 W
## 314 314 314 W 46 47 46 46 W
## 315 315 315 W 47 48 47 47 W
## 316 316 316 W 48 49 48 48 W
## 317 317 317 W 49 50 49 49 W
## 318 318 318 W 50 51 50 50 W
## 319 319 319 W 51 52 51 51 W
## 320 320 320 W 52 53 52 52 W
## 321 321 321 W 53 54 53 53 W
## 322 322 322 W 54 55 54 54 W
## 323 323 323 W 55 56 55 55 W
## 324 324 324 W 56 57 56 56 W
## 325 325 325 W 57 58 57 57 W
## 326 326 326 W 58 59 58 58 W
## 327 327 327 W 59 60 59 59 W
## 328 328 328 W 60 61 60 60 W
## 329 329 329 W 61 62 61 61 W
## 330 330 330 W 62 63 62 62 W
## 331 331 331 W 63 64 63 63 W
## 332 332 332 W 64 65 64 64 W
## 333 333 333 W 65 66 65 65 W
## 334 334 334 W 66 67 66 66 W
## 335 335 335 W 67 68 67 67 W
## 336 336 336 W 68 69 68 68 W
## 337 337 337 W 69 70 69 69 W
## 338 338 338 W 70 71 70 70 W
## 339 339 339 W 71 72 71 71 W
## 340 340 340 W 72 73 72 72 W
## 341 341 341 W 73 74 73 73 W
## 342 342 342 W 74 75 74 74 W
## 343 343 343 W 75 76 75 75 W
## 344 344 344 W 76 77 76 76 W
## 345 345 345 W 77 78 77 77 W
## 346 346 346 W 78 79 78 78 W
## 347 347 347 W 79 80 79 79 W
## 348 348 348 W 80 81 80 80 W
## 349 349 349 W 81 82 81 81 W
## 350 350 350 W 82 83 82 82 W
## 351 351 351 W 83 84 83 83 W
## 352 352 352 W 84 85 84 84 W
## 353 353 353 W 85 86 85 85 W
## 354 354 354 W 86 87 86 86 W
## 355 355 355 W 87 88 87 87 W
## 356 356 356 W 88 89 88 88 W
##
## $pimtypes
## $pimtypes$S
## $pimtypes$S$pim.type
## [1] "all"
# 3. Fit a particular model
# This is a model with S varying by habitat type
mod.hab <- RMark::mark(mall.proc, ddl=mall.ddl,
model="Nest",
model.parameters=list(
S =list(formula=~Habitat)
)
)
##
## Output summary for Nest model
## Name : S(~Habitat)
##
## Npar : 4
## -2lnL: 1563.951
## AICc : 1571.957
##
## Beta
## estimate se lcl ucl
## S:(Intercept) 2.8629313 0.0992682 2.6683656 3.0574970
## S:HabitatP 0.2226790 0.1273492 -0.0269255 0.4722835
## S:HabitatR 0.2111142 0.2356356 -0.2507317 0.6729600
## S:HabitatW 0.0929137 0.2454492 -0.3881667 0.5739941
##
##
## Real Parameter S
## 1 2 3 4 5 6
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 7 8 9 10 11 12
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 13 14 15 16 17 18
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 19 20 21 22 23 24
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 25 26 27 28 29 30
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 31 32 33 34 35 36
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 37 38 39 40 41 42
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 43 44 45 46 47 48
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 49 50 51 52 53 54
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 55 56 57 58 59 60
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 61 62 63 64 65 66
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 67 68 69 70 71 72
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 73 74 75 76 77 78
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 79 80 81 82 83 84
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 85 86 87 88 89
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
summary(mod.hab)
## Output summary for Nest model
## Name : S(~Habitat)
##
## Npar : 4
## -2lnL: 1563.951
## AICc : 1571.957
##
## Beta
## estimate se lcl ucl
## S:(Intercept) 2.8629313 0.0992682 2.6683656 3.0574970
## S:HabitatP 0.2226790 0.1273492 -0.0269255 0.4722835
## S:HabitatR 0.2111142 0.2356356 -0.2507317 0.6729600
## S:HabitatW 0.0929137 0.2454492 -0.3881667 0.5739941
##
##
## Real Parameter S
## 1 2 3 4 5 6
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 7 8 9 10 11 12
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 13 14 15 16 17 18
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 19 20 21 22 23 24
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 25 26 27 28 29 30
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 31 32 33 34 35 36
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 37 38 39 40 41 42
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 43 44 45 46 47 48
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 49 50 51 52 53 54
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 55 56 57 58 59 60
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 61 62 63 64 65 66
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 67 68 69 70 71 72
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 73 74 75 76 77 78
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 79 80 81 82 83 84
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
## 85 86 87 88 89
## Group:HabitatN 0.9459833 0.9459833 0.9459833 0.9459833 0.9459833
## Group:HabitatP 0.9562953 0.9562953 0.9562953 0.9562953 0.9562953
## Group:HabitatR 0.9558094 0.9558094 0.9558094 0.9558094 0.9558094
## Group:HabitatW 0.9505390 0.9505390 0.9505390 0.9505390 0.9505390
# Look the objects returned in more details
names(mod.hab)
## [1] "data" "model" "title"
## [4] "model.name" "links" "mixtures"
## [7] "call" "parameters" "time.intervals"
## [10] "number.of.groups" "group.labels" "nocc"
## [13] "begin.time" "covariates" "fixed"
## [16] "design.matrix" "pims" "design.data"
## [19] "strata.labels" "mlogit.list" "profile.int"
## [22] "simplify" "model.parameters" "results"
## [25] "output"
names(mod.hab$results)
## [1] "lnl" "deviance" "deviance.df"
## [4] "npar" "n" "AICc"
## [7] "beta" "real" "beta.vcv"
## [10] "derived" "derived.vcv" "covariate.values"
## [13] "singular" "real.vcv"
# look at estimates on beta and original scale
mod.hab$results$beta # on the logit scale
## estimate se lcl ucl
## S:(Intercept) 2.8629313 0.0992682 2.6683656 3.0574970
## S:HabitatP 0.2226790 0.1273492 -0.0269255 0.4722835
## S:HabitatR 0.2111142 0.2356356 -0.2507317 0.6729600
## S:HabitatW 0.0929137 0.2454492 -0.3881667 0.5739941
mod.hab$results$real# on the regular 0-1 scale for each habitat
## estimate se lcl ucl fixed note
## S gN a0 t1 0.9459833 0.0050725 0.9351340 0.9551051
## S gP a0 t1 0.9562953 0.0033341 0.9492739 0.9623833
## S gR a0 t1 0.9558094 0.0090265 0.9343297 0.9704853
## S gW a0 t1 0.9505390 0.0105538 0.9252465 0.9675738
# derived variabldes is the nest survival probability over the (nocc) days
names(mod.hab$results$derived)
## [1] "S Overall Survival"
mod.hab$results$derived$"S Overall Survival"
## estimate se lcl ucl
## 1 0.007138846 0.003406925 0.002794926 0.01811154
## 2 0.018736500 0.005813810 0.010169356 0.03427112
## 3 0.017907864 0.015051503 0.003395653 0.08890875
## 4 0.010947760 0.010818212 0.001559456 0.07273834
# alternatively
get.real(mod.hab, "S", se=TRUE)
## all.diff.index par.index estimate se lcl
## S gN a0 t1 1 1 0.9459833 0.0050725 0.9351340
## S gN a1 t2 2 1 0.9459833 0.0050725 0.9351340
## S gN a2 t3 3 1 0.9459833 0.0050725 0.9351340
## S gN a3 t4 4 1 0.9459833 0.0050725 0.9351340
## S gN a4 t5 5 1 0.9459833 0.0050725 0.9351340
## S gN a5 t6 6 1 0.9459833 0.0050725 0.9351340
## S gN a6 t7 7 1 0.9459833 0.0050725 0.9351340
## S gN a7 t8 8 1 0.9459833 0.0050725 0.9351340
## S gN a8 t9 9 1 0.9459833 0.0050725 0.9351340
## S gN a9 t10 10 1 0.9459833 0.0050725 0.9351340
## S gN a10 t11 11 1 0.9459833 0.0050725 0.9351340
## S gN a11 t12 12 1 0.9459833 0.0050725 0.9351340
## S gN a12 t13 13 1 0.9459833 0.0050725 0.9351340
## S gN a13 t14 14 1 0.9459833 0.0050725 0.9351340
## S gN a14 t15 15 1 0.9459833 0.0050725 0.9351340
## S gN a15 t16 16 1 0.9459833 0.0050725 0.9351340
## S gN a16 t17 17 1 0.9459833 0.0050725 0.9351340
## S gN a17 t18 18 1 0.9459833 0.0050725 0.9351340
## S gN a18 t19 19 1 0.9459833 0.0050725 0.9351340
## S gN a19 t20 20 1 0.9459833 0.0050725 0.9351340
## S gN a20 t21 21 1 0.9459833 0.0050725 0.9351340
## S gN a21 t22 22 1 0.9459833 0.0050725 0.9351340
## S gN a22 t23 23 1 0.9459833 0.0050725 0.9351340
## S gN a23 t24 24 1 0.9459833 0.0050725 0.9351340
## S gN a24 t25 25 1 0.9459833 0.0050725 0.9351340
## S gN a25 t26 26 1 0.9459833 0.0050725 0.9351340
## S gN a26 t27 27 1 0.9459833 0.0050725 0.9351340
## S gN a27 t28 28 1 0.9459833 0.0050725 0.9351340
## S gN a28 t29 29 1 0.9459833 0.0050725 0.9351340
## S gN a29 t30 30 1 0.9459833 0.0050725 0.9351340
## S gN a30 t31 31 1 0.9459833 0.0050725 0.9351340
## S gN a31 t32 32 1 0.9459833 0.0050725 0.9351340
## S gN a32 t33 33 1 0.9459833 0.0050725 0.9351340
## S gN a33 t34 34 1 0.9459833 0.0050725 0.9351340
## S gN a34 t35 35 1 0.9459833 0.0050725 0.9351340
## S gN a35 t36 36 1 0.9459833 0.0050725 0.9351340
## S gN a36 t37 37 1 0.9459833 0.0050725 0.9351340
## S gN a37 t38 38 1 0.9459833 0.0050725 0.9351340
## S gN a38 t39 39 1 0.9459833 0.0050725 0.9351340
## S gN a39 t40 40 1 0.9459833 0.0050725 0.9351340
## S gN a40 t41 41 1 0.9459833 0.0050725 0.9351340
## S gN a41 t42 42 1 0.9459833 0.0050725 0.9351340
## S gN a42 t43 43 1 0.9459833 0.0050725 0.9351340
## S gN a43 t44 44 1 0.9459833 0.0050725 0.9351340
## S gN a44 t45 45 1 0.9459833 0.0050725 0.9351340
## S gN a45 t46 46 1 0.9459833 0.0050725 0.9351340
## S gN a46 t47 47 1 0.9459833 0.0050725 0.9351340
## S gN a47 t48 48 1 0.9459833 0.0050725 0.9351340
## S gN a48 t49 49 1 0.9459833 0.0050725 0.9351340
## S gN a49 t50 50 1 0.9459833 0.0050725 0.9351340
## S gN a50 t51 51 1 0.9459833 0.0050725 0.9351340
## S gN a51 t52 52 1 0.9459833 0.0050725 0.9351340
## S gN a52 t53 53 1 0.9459833 0.0050725 0.9351340
## S gN a53 t54 54 1 0.9459833 0.0050725 0.9351340
## S gN a54 t55 55 1 0.9459833 0.0050725 0.9351340
## S gN a55 t56 56 1 0.9459833 0.0050725 0.9351340
## S gN a56 t57 57 1 0.9459833 0.0050725 0.9351340
## S gN a57 t58 58 1 0.9459833 0.0050725 0.9351340
## S gN a58 t59 59 1 0.9459833 0.0050725 0.9351340
## S gN a59 t60 60 1 0.9459833 0.0050725 0.9351340
## S gN a60 t61 61 1 0.9459833 0.0050725 0.9351340
## S gN a61 t62 62 1 0.9459833 0.0050725 0.9351340
## S gN a62 t63 63 1 0.9459833 0.0050725 0.9351340
## S gN a63 t64 64 1 0.9459833 0.0050725 0.9351340
## S gN a64 t65 65 1 0.9459833 0.0050725 0.9351340
## S gN a65 t66 66 1 0.9459833 0.0050725 0.9351340
## S gN a66 t67 67 1 0.9459833 0.0050725 0.9351340
## S gN a67 t68 68 1 0.9459833 0.0050725 0.9351340
## S gN a68 t69 69 1 0.9459833 0.0050725 0.9351340
## S gN a69 t70 70 1 0.9459833 0.0050725 0.9351340
## S gN a70 t71 71 1 0.9459833 0.0050725 0.9351340
## S gN a71 t72 72 1 0.9459833 0.0050725 0.9351340
## S gN a72 t73 73 1 0.9459833 0.0050725 0.9351340
## S gN a73 t74 74 1 0.9459833 0.0050725 0.9351340
## S gN a74 t75 75 1 0.9459833 0.0050725 0.9351340
## S gN a75 t76 76 1 0.9459833 0.0050725 0.9351340
## S gN a76 t77 77 1 0.9459833 0.0050725 0.9351340
## S gN a77 t78 78 1 0.9459833 0.0050725 0.9351340
## S gN a78 t79 79 1 0.9459833 0.0050725 0.9351340
## S gN a79 t80 80 1 0.9459833 0.0050725 0.9351340
## S gN a80 t81 81 1 0.9459833 0.0050725 0.9351340
## S gN a81 t82 82 1 0.9459833 0.0050725 0.9351340
## S gN a82 t83 83 1 0.9459833 0.0050725 0.9351340
## S gN a83 t84 84 1 0.9459833 0.0050725 0.9351340
## S gN a84 t85 85 1 0.9459833 0.0050725 0.9351340
## S gN a85 t86 86 1 0.9459833 0.0050725 0.9351340
## S gN a86 t87 87 1 0.9459833 0.0050725 0.9351340
## S gN a87 t88 88 1 0.9459833 0.0050725 0.9351340
## S gN a88 t89 89 1 0.9459833 0.0050725 0.9351340
## S gP a0 t1 90 2 0.9562953 0.0033341 0.9492739
## S gP a1 t2 91 2 0.9562953 0.0033341 0.9492739
## S gP a2 t3 92 2 0.9562953 0.0033341 0.9492739
## S gP a3 t4 93 2 0.9562953 0.0033341 0.9492739
## S gP a4 t5 94 2 0.9562953 0.0033341 0.9492739
## S gP a5 t6 95 2 0.9562953 0.0033341 0.9492739
## S gP a6 t7 96 2 0.9562953 0.0033341 0.9492739
## S gP a7 t8 97 2 0.9562953 0.0033341 0.9492739
## S gP a8 t9 98 2 0.9562953 0.0033341 0.9492739
## S gP a9 t10 99 2 0.9562953 0.0033341 0.9492739
## S gP a10 t11 100 2 0.9562953 0.0033341 0.9492739
## S gP a11 t12 101 2 0.9562953 0.0033341 0.9492739
## S gP a12 t13 102 2 0.9562953 0.0033341 0.9492739
## S gP a13 t14 103 2 0.9562953 0.0033341 0.9492739
## S gP a14 t15 104 2 0.9562953 0.0033341 0.9492739
## S gP a15 t16 105 2 0.9562953 0.0033341 0.9492739
## S gP a16 t17 106 2 0.9562953 0.0033341 0.9492739
## S gP a17 t18 107 2 0.9562953 0.0033341 0.9492739
## S gP a18 t19 108 2 0.9562953 0.0033341 0.9492739
## S gP a19 t20 109 2 0.9562953 0.0033341 0.9492739
## S gP a20 t21 110 2 0.9562953 0.0033341 0.9492739
## S gP a21 t22 111 2 0.9562953 0.0033341 0.9492739
## S gP a22 t23 112 2 0.9562953 0.0033341 0.9492739
## S gP a23 t24 113 2 0.9562953 0.0033341 0.9492739
## S gP a24 t25 114 2 0.9562953 0.0033341 0.9492739
## S gP a25 t26 115 2 0.9562953 0.0033341 0.9492739
## S gP a26 t27 116 2 0.9562953 0.0033341 0.9492739
## S gP a27 t28 117 2 0.9562953 0.0033341 0.9492739
## S gP a28 t29 118 2 0.9562953 0.0033341 0.9492739
## S gP a29 t30 119 2 0.9562953 0.0033341 0.9492739
## S gP a30 t31 120 2 0.9562953 0.0033341 0.9492739
## S gP a31 t32 121 2 0.9562953 0.0033341 0.9492739
## S gP a32 t33 122 2 0.9562953 0.0033341 0.9492739
## S gP a33 t34 123 2 0.9562953 0.0033341 0.9492739
## S gP a34 t35 124 2 0.9562953 0.0033341 0.9492739
## S gP a35 t36 125 2 0.9562953 0.0033341 0.9492739
## S gP a36 t37 126 2 0.9562953 0.0033341 0.9492739
## S gP a37 t38 127 2 0.9562953 0.0033341 0.9492739
## S gP a38 t39 128 2 0.9562953 0.0033341 0.9492739
## S gP a39 t40 129 2 0.9562953 0.0033341 0.9492739
## S gP a40 t41 130 2 0.9562953 0.0033341 0.9492739
## S gP a41 t42 131 2 0.9562953 0.0033341 0.9492739
## S gP a42 t43 132 2 0.9562953 0.0033341 0.9492739
## S gP a43 t44 133 2 0.9562953 0.0033341 0.9492739
## S gP a44 t45 134 2 0.9562953 0.0033341 0.9492739
## S gP a45 t46 135 2 0.9562953 0.0033341 0.9492739
## S gP a46 t47 136 2 0.9562953 0.0033341 0.9492739
## S gP a47 t48 137 2 0.9562953 0.0033341 0.9492739
## S gP a48 t49 138 2 0.9562953 0.0033341 0.9492739
## S gP a49 t50 139 2 0.9562953 0.0033341 0.9492739
## S gP a50 t51 140 2 0.9562953 0.0033341 0.9492739
## S gP a51 t52 141 2 0.9562953 0.0033341 0.9492739
## S gP a52 t53 142 2 0.9562953 0.0033341 0.9492739
## S gP a53 t54 143 2 0.9562953 0.0033341 0.9492739
## S gP a54 t55 144 2 0.9562953 0.0033341 0.9492739
## S gP a55 t56 145 2 0.9562953 0.0033341 0.9492739
## S gP a56 t57 146 2 0.9562953 0.0033341 0.9492739
## S gP a57 t58 147 2 0.9562953 0.0033341 0.9492739
## S gP a58 t59 148 2 0.9562953 0.0033341 0.9492739
## S gP a59 t60 149 2 0.9562953 0.0033341 0.9492739
## S gP a60 t61 150 2 0.9562953 0.0033341 0.9492739
## S gP a61 t62 151 2 0.9562953 0.0033341 0.9492739
## S gP a62 t63 152 2 0.9562953 0.0033341 0.9492739
## S gP a63 t64 153 2 0.9562953 0.0033341 0.9492739
## S gP a64 t65 154 2 0.9562953 0.0033341 0.9492739
## S gP a65 t66 155 2 0.9562953 0.0033341 0.9492739
## S gP a66 t67 156 2 0.9562953 0.0033341 0.9492739
## S gP a67 t68 157 2 0.9562953 0.0033341 0.9492739
## S gP a68 t69 158 2 0.9562953 0.0033341 0.9492739
## S gP a69 t70 159 2 0.9562953 0.0033341 0.9492739
## S gP a70 t71 160 2 0.9562953 0.0033341 0.9492739
## S gP a71 t72 161 2 0.9562953 0.0033341 0.9492739
## S gP a72 t73 162 2 0.9562953 0.0033341 0.9492739
## S gP a73 t74 163 2 0.9562953 0.0033341 0.9492739
## S gP a74 t75 164 2 0.9562953 0.0033341 0.9492739
## S gP a75 t76 165 2 0.9562953 0.0033341 0.9492739
## S gP a76 t77 166 2 0.9562953 0.0033341 0.9492739
## S gP a77 t78 167 2 0.9562953 0.0033341 0.9492739
## S gP a78 t79 168 2 0.9562953 0.0033341 0.9492739
## S gP a79 t80 169 2 0.9562953 0.0033341 0.9492739
## S gP a80 t81 170 2 0.9562953 0.0033341 0.9492739
## S gP a81 t82 171 2 0.9562953 0.0033341 0.9492739
## S gP a82 t83 172 2 0.9562953 0.0033341 0.9492739
## S gP a83 t84 173 2 0.9562953 0.0033341 0.9492739
## S gP a84 t85 174 2 0.9562953 0.0033341 0.9492739
## S gP a85 t86 175 2 0.9562953 0.0033341 0.9492739
## S gP a86 t87 176 2 0.9562953 0.0033341 0.9492739
## S gP a87 t88 177 2 0.9562953 0.0033341 0.9492739
## S gP a88 t89 178 2 0.9562953 0.0033341 0.9492739
## S gR a0 t1 179 3 0.9558094 0.0090265 0.9343297
## S gR a1 t2 180 3 0.9558094 0.0090265 0.9343297
## S gR a2 t3 181 3 0.9558094 0.0090265 0.9343297
## S gR a3 t4 182 3 0.9558094 0.0090265 0.9343297
## S gR a4 t5 183 3 0.9558094 0.0090265 0.9343297
## S gR a5 t6 184 3 0.9558094 0.0090265 0.9343297
## S gR a6 t7 185 3 0.9558094 0.0090265 0.9343297
## S gR a7 t8 186 3 0.9558094 0.0090265 0.9343297
## S gR a8 t9 187 3 0.9558094 0.0090265 0.9343297
## S gR a9 t10 188 3 0.9558094 0.0090265 0.9343297
## S gR a10 t11 189 3 0.9558094 0.0090265 0.9343297
## S gR a11 t12 190 3 0.9558094 0.0090265 0.9343297
## S gR a12 t13 191 3 0.9558094 0.0090265 0.9343297
## S gR a13 t14 192 3 0.9558094 0.0090265 0.9343297
## S gR a14 t15 193 3 0.9558094 0.0090265 0.9343297
## S gR a15 t16 194 3 0.9558094 0.0090265 0.9343297
## S gR a16 t17 195 3 0.9558094 0.0090265 0.9343297
## S gR a17 t18 196 3 0.9558094 0.0090265 0.9343297
## S gR a18 t19 197 3 0.9558094 0.0090265 0.9343297
## S gR a19 t20 198 3 0.9558094 0.0090265 0.9343297
## S gR a20 t21 199 3 0.9558094 0.0090265 0.9343297
## S gR a21 t22 200 3 0.9558094 0.0090265 0.9343297
## S gR a22 t23 201 3 0.9558094 0.0090265 0.9343297
## S gR a23 t24 202 3 0.9558094 0.0090265 0.9343297
## S gR a24 t25 203 3 0.9558094 0.0090265 0.9343297
## S gR a25 t26 204 3 0.9558094 0.0090265 0.9343297
## S gR a26 t27 205 3 0.9558094 0.0090265 0.9343297
## S gR a27 t28 206 3 0.9558094 0.0090265 0.9343297
## S gR a28 t29 207 3 0.9558094 0.0090265 0.9343297
## S gR a29 t30 208 3 0.9558094 0.0090265 0.9343297
## S gR a30 t31 209 3 0.9558094 0.0090265 0.9343297
## S gR a31 t32 210 3 0.9558094 0.0090265 0.9343297
## S gR a32 t33 211 3 0.9558094 0.0090265 0.9343297
## S gR a33 t34 212 3 0.9558094 0.0090265 0.9343297
## S gR a34 t35 213 3 0.9558094 0.0090265 0.9343297
## S gR a35 t36 214 3 0.9558094 0.0090265 0.9343297
## S gR a36 t37 215 3 0.9558094 0.0090265 0.9343297
## S gR a37 t38 216 3 0.9558094 0.0090265 0.9343297
## S gR a38 t39 217 3 0.9558094 0.0090265 0.9343297
## S gR a39 t40 218 3 0.9558094 0.0090265 0.9343297
## S gR a40 t41 219 3 0.9558094 0.0090265 0.9343297
## S gR a41 t42 220 3 0.9558094 0.0090265 0.9343297
## S gR a42 t43 221 3 0.9558094 0.0090265 0.9343297
## S gR a43 t44 222 3 0.9558094 0.0090265 0.9343297
## S gR a44 t45 223 3 0.9558094 0.0090265 0.9343297
## S gR a45 t46 224 3 0.9558094 0.0090265 0.9343297
## S gR a46 t47 225 3 0.9558094 0.0090265 0.9343297
## S gR a47 t48 226 3 0.9558094 0.0090265 0.9343297
## S gR a48 t49 227 3 0.9558094 0.0090265 0.9343297
## S gR a49 t50 228 3 0.9558094 0.0090265 0.9343297
## S gR a50 t51 229 3 0.9558094 0.0090265 0.9343297
## S gR a51 t52 230 3 0.9558094 0.0090265 0.9343297
## S gR a52 t53 231 3 0.9558094 0.0090265 0.9343297
## S gR a53 t54 232 3 0.9558094 0.0090265 0.9343297
## S gR a54 t55 233 3 0.9558094 0.0090265 0.9343297
## S gR a55 t56 234 3 0.9558094 0.0090265 0.9343297
## S gR a56 t57 235 3 0.9558094 0.0090265 0.9343297
## S gR a57 t58 236 3 0.9558094 0.0090265 0.9343297
## S gR a58 t59 237 3 0.9558094 0.0090265 0.9343297
## S gR a59 t60 238 3 0.9558094 0.0090265 0.9343297
## S gR a60 t61 239 3 0.9558094 0.0090265 0.9343297
## S gR a61 t62 240 3 0.9558094 0.0090265 0.9343297
## S gR a62 t63 241 3 0.9558094 0.0090265 0.9343297
## S gR a63 t64 242 3 0.9558094 0.0090265 0.9343297
## S gR a64 t65 243 3 0.9558094 0.0090265 0.9343297
## S gR a65 t66 244 3 0.9558094 0.0090265 0.9343297
## S gR a66 t67 245 3 0.9558094 0.0090265 0.9343297
## S gR a67 t68 246 3 0.9558094 0.0090265 0.9343297
## S gR a68 t69 247 3 0.9558094 0.0090265 0.9343297
## S gR a69 t70 248 3 0.9558094 0.0090265 0.9343297
## S gR a70 t71 249 3 0.9558094 0.0090265 0.9343297
## S gR a71 t72 250 3 0.9558094 0.0090265 0.9343297
## S gR a72 t73 251 3 0.9558094 0.0090265 0.9343297
## S gR a73 t74 252 3 0.9558094 0.0090265 0.9343297
## S gR a74 t75 253 3 0.9558094 0.0090265 0.9343297
## S gR a75 t76 254 3 0.9558094 0.0090265 0.9343297
## S gR a76 t77 255 3 0.9558094 0.0090265 0.9343297
## S gR a77 t78 256 3 0.9558094 0.0090265 0.9343297
## S gR a78 t79 257 3 0.9558094 0.0090265 0.9343297
## S gR a79 t80 258 3 0.9558094 0.0090265 0.9343297
## S gR a80 t81 259 3 0.9558094 0.0090265 0.9343297
## S gR a81 t82 260 3 0.9558094 0.0090265 0.9343297
## S gR a82 t83 261 3 0.9558094 0.0090265 0.9343297
## S gR a83 t84 262 3 0.9558094 0.0090265 0.9343297
## S gR a84 t85 263 3 0.9558094 0.0090265 0.9343297
## S gR a85 t86 264 3 0.9558094 0.0090265 0.9343297
## S gR a86 t87 265 3 0.9558094 0.0090265 0.9343297
## S gR a87 t88 266 3 0.9558094 0.0090265 0.9343297
## S gR a88 t89 267 3 0.9558094 0.0090265 0.9343297
## S gW a0 t1 268 4 0.9505390 0.0105538 0.9252465
## S gW a1 t2 269 4 0.9505390 0.0105538 0.9252465
## S gW a2 t3 270 4 0.9505390 0.0105538 0.9252465
## S gW a3 t4 271 4 0.9505390 0.0105538 0.9252465
## S gW a4 t5 272 4 0.9505390 0.0105538 0.9252465
## S gW a5 t6 273 4 0.9505390 0.0105538 0.9252465
## S gW a6 t7 274 4 0.9505390 0.0105538 0.9252465
## S gW a7 t8 275 4 0.9505390 0.0105538 0.9252465
## S gW a8 t9 276 4 0.9505390 0.0105538 0.9252465
## S gW a9 t10 277 4 0.9505390 0.0105538 0.9252465
## S gW a10 t11 278 4 0.9505390 0.0105538 0.9252465
## S gW a11 t12 279 4 0.9505390 0.0105538 0.9252465
## S gW a12 t13 280 4 0.9505390 0.0105538 0.9252465
## S gW a13 t14 281 4 0.9505390 0.0105538 0.9252465
## S gW a14 t15 282 4 0.9505390 0.0105538 0.9252465
## S gW a15 t16 283 4 0.9505390 0.0105538 0.9252465
## S gW a16 t17 284 4 0.9505390 0.0105538 0.9252465
## S gW a17 t18 285 4 0.9505390 0.0105538 0.9252465
## S gW a18 t19 286 4 0.9505390 0.0105538 0.9252465
## S gW a19 t20 287 4 0.9505390 0.0105538 0.9252465
## S gW a20 t21 288 4 0.9505390 0.0105538 0.9252465
## S gW a21 t22 289 4 0.9505390 0.0105538 0.9252465
## S gW a22 t23 290 4 0.9505390 0.0105538 0.9252465
## S gW a23 t24 291 4 0.9505390 0.0105538 0.9252465
## S gW a24 t25 292 4 0.9505390 0.0105538 0.9252465
## S gW a25 t26 293 4 0.9505390 0.0105538 0.9252465
## S gW a26 t27 294 4 0.9505390 0.0105538 0.9252465
## S gW a27 t28 295 4 0.9505390 0.0105538 0.9252465
## S gW a28 t29 296 4 0.9505390 0.0105538 0.9252465
## S gW a29 t30 297 4 0.9505390 0.0105538 0.9252465
## S gW a30 t31 298 4 0.9505390 0.0105538 0.9252465
## S gW a31 t32 299 4 0.9505390 0.0105538 0.9252465
## S gW a32 t33 300 4 0.9505390 0.0105538 0.9252465
## S gW a33 t34 301 4 0.9505390 0.0105538 0.9252465
## S gW a34 t35 302 4 0.9505390 0.0105538 0.9252465
## S gW a35 t36 303 4 0.9505390 0.0105538 0.9252465
## S gW a36 t37 304 4 0.9505390 0.0105538 0.9252465
## S gW a37 t38 305 4 0.9505390 0.0105538 0.9252465
## S gW a38 t39 306 4 0.9505390 0.0105538 0.9252465
## S gW a39 t40 307 4 0.9505390 0.0105538 0.9252465
## S gW a40 t41 308 4 0.9505390 0.0105538 0.9252465
## S gW a41 t42 309 4 0.9505390 0.0105538 0.9252465
## S gW a42 t43 310 4 0.9505390 0.0105538 0.9252465
## S gW a43 t44 311 4 0.9505390 0.0105538 0.9252465
## S gW a44 t45 312 4 0.9505390 0.0105538 0.9252465
## S gW a45 t46 313 4 0.9505390 0.0105538 0.9252465
## S gW a46 t47 314 4 0.9505390 0.0105538 0.9252465
## S gW a47 t48 315 4 0.9505390 0.0105538 0.9252465
## S gW a48 t49 316 4 0.9505390 0.0105538 0.9252465
## S gW a49 t50 317 4 0.9505390 0.0105538 0.9252465
## S gW a50 t51 318 4 0.9505390 0.0105538 0.9252465
## S gW a51 t52 319 4 0.9505390 0.0105538 0.9252465
## S gW a52 t53 320 4 0.9505390 0.0105538 0.9252465
## S gW a53 t54 321 4 0.9505390 0.0105538 0.9252465
## S gW a54 t55 322 4 0.9505390 0.0105538 0.9252465
## S gW a55 t56 323 4 0.9505390 0.0105538 0.9252465
## S gW a56 t57 324 4 0.9505390 0.0105538 0.9252465
## S gW a57 t58 325 4 0.9505390 0.0105538 0.9252465
## S gW a58 t59 326 4 0.9505390 0.0105538 0.9252465
## S gW a59 t60 327 4 0.9505390 0.0105538 0.9252465
## S gW a60 t61 328 4 0.9505390 0.0105538 0.9252465
## S gW a61 t62 329 4 0.9505390 0.0105538 0.9252465
## S gW a62 t63 330 4 0.9505390 0.0105538 0.9252465
## S gW a63 t64 331 4 0.9505390 0.0105538 0.9252465
## S gW a64 t65 332 4 0.9505390 0.0105538 0.9252465
## S gW a65 t66 333 4 0.9505390 0.0105538 0.9252465
## S gW a66 t67 334 4 0.9505390 0.0105538 0.9252465
## S gW a67 t68 335 4 0.9505390 0.0105538 0.9252465
## S gW a68 t69 336 4 0.9505390 0.0105538 0.9252465
## S gW a69 t70 337 4 0.9505390 0.0105538 0.9252465
## S gW a70 t71 338 4 0.9505390 0.0105538 0.9252465
## S gW a71 t72 339 4 0.9505390 0.0105538 0.9252465
## S gW a72 t73 340 4 0.9505390 0.0105538 0.9252465
## S gW a73 t74 341 4 0.9505390 0.0105538 0.9252465
## S gW a74 t75 342 4 0.9505390 0.0105538 0.9252465
## S gW a75 t76 343 4 0.9505390 0.0105538 0.9252465
## S gW a76 t77 344 4 0.9505390 0.0105538 0.9252465
## S gW a77 t78 345 4 0.9505390 0.0105538 0.9252465
## S gW a78 t79 346 4 0.9505390 0.0105538 0.9252465
## S gW a79 t80 347 4 0.9505390 0.0105538 0.9252465
## S gW a80 t81 348 4 0.9505390 0.0105538 0.9252465
## S gW a81 t82 349 4 0.9505390 0.0105538 0.9252465
## S gW a82 t83 350 4 0.9505390 0.0105538 0.9252465
## S gW a83 t84 351 4 0.9505390 0.0105538 0.9252465
## S gW a84 t85 352 4 0.9505390 0.0105538 0.9252465
## S gW a85 t86 353 4 0.9505390 0.0105538 0.9252465
## S gW a86 t87 354 4 0.9505390 0.0105538 0.9252465
## S gW a87 t88 355 4 0.9505390 0.0105538 0.9252465
## S gW a88 t89 356 4 0.9505390 0.0105538 0.9252465
## ucl fixed note group age time Age Time Habitat
## S gN a0 t1 0.9551051 N 0 1 0 0 N
## S gN a1 t2 0.9551051 N 1 2 1 1 N
## S gN a2 t3 0.9551051 N 2 3 2 2 N
## S gN a3 t4 0.9551051 N 3 4 3 3 N
## S gN a4 t5 0.9551051 N 4 5 4 4 N
## S gN a5 t6 0.9551051 N 5 6 5 5 N
## S gN a6 t7 0.9551051 N 6 7 6 6 N
## S gN a7 t8 0.9551051 N 7 8 7 7 N
## S gN a8 t9 0.9551051 N 8 9 8 8 N
## S gN a9 t10 0.9551051 N 9 10 9 9 N
## S gN a10 t11 0.9551051 N 10 11 10 10 N
## S gN a11 t12 0.9551051 N 11 12 11 11 N
## S gN a12 t13 0.9551051 N 12 13 12 12 N
## S gN a13 t14 0.9551051 N 13 14 13 13 N
## S gN a14 t15 0.9551051 N 14 15 14 14 N
## S gN a15 t16 0.9551051 N 15 16 15 15 N
## S gN a16 t17 0.9551051 N 16 17 16 16 N
## S gN a17 t18 0.9551051 N 17 18 17 17 N
## S gN a18 t19 0.9551051 N 18 19 18 18 N
## S gN a19 t20 0.9551051 N 19 20 19 19 N
## S gN a20 t21 0.9551051 N 20 21 20 20 N
## S gN a21 t22 0.9551051 N 21 22 21 21 N
## S gN a22 t23 0.9551051 N 22 23 22 22 N
## S gN a23 t24 0.9551051 N 23 24 23 23 N
## S gN a24 t25 0.9551051 N 24 25 24 24 N
## S gN a25 t26 0.9551051 N 25 26 25 25 N
## S gN a26 t27 0.9551051 N 26 27 26 26 N
## S gN a27 t28 0.9551051 N 27 28 27 27 N
## S gN a28 t29 0.9551051 N 28 29 28 28 N
## S gN a29 t30 0.9551051 N 29 30 29 29 N
## S gN a30 t31 0.9551051 N 30 31 30 30 N
## S gN a31 t32 0.9551051 N 31 32 31 31 N
## S gN a32 t33 0.9551051 N 32 33 32 32 N
## S gN a33 t34 0.9551051 N 33 34 33 33 N
## S gN a34 t35 0.9551051 N 34 35 34 34 N
## S gN a35 t36 0.9551051 N 35 36 35 35 N
## S gN a36 t37 0.9551051 N 36 37 36 36 N
## S gN a37 t38 0.9551051 N 37 38 37 37 N
## S gN a38 t39 0.9551051 N 38 39 38 38 N
## S gN a39 t40 0.9551051 N 39 40 39 39 N
## S gN a40 t41 0.9551051 N 40 41 40 40 N
## S gN a41 t42 0.9551051 N 41 42 41 41 N
## S gN a42 t43 0.9551051 N 42 43 42 42 N
## S gN a43 t44 0.9551051 N 43 44 43 43 N
## S gN a44 t45 0.9551051 N 44 45 44 44 N
## S gN a45 t46 0.9551051 N 45 46 45 45 N
## S gN a46 t47 0.9551051 N 46 47 46 46 N
## S gN a47 t48 0.9551051 N 47 48 47 47 N
## S gN a48 t49 0.9551051 N 48 49 48 48 N
## S gN a49 t50 0.9551051 N 49 50 49 49 N
## S gN a50 t51 0.9551051 N 50 51 50 50 N
## S gN a51 t52 0.9551051 N 51 52 51 51 N
## S gN a52 t53 0.9551051 N 52 53 52 52 N
## S gN a53 t54 0.9551051 N 53 54 53 53 N
## S gN a54 t55 0.9551051 N 54 55 54 54 N
## S gN a55 t56 0.9551051 N 55 56 55 55 N
## S gN a56 t57 0.9551051 N 56 57 56 56 N
## S gN a57 t58 0.9551051 N 57 58 57 57 N
## S gN a58 t59 0.9551051 N 58 59 58 58 N
## S gN a59 t60 0.9551051 N 59 60 59 59 N
## S gN a60 t61 0.9551051 N 60 61 60 60 N
## S gN a61 t62 0.9551051 N 61 62 61 61 N
## S gN a62 t63 0.9551051 N 62 63 62 62 N
## S gN a63 t64 0.9551051 N 63 64 63 63 N
## S gN a64 t65 0.9551051 N 64 65 64 64 N
## S gN a65 t66 0.9551051 N 65 66 65 65 N
## S gN a66 t67 0.9551051 N 66 67 66 66 N
## S gN a67 t68 0.9551051 N 67 68 67 67 N
## S gN a68 t69 0.9551051 N 68 69 68 68 N
## S gN a69 t70 0.9551051 N 69 70 69 69 N
## S gN a70 t71 0.9551051 N 70 71 70 70 N
## S gN a71 t72 0.9551051 N 71 72 71 71 N
## S gN a72 t73 0.9551051 N 72 73 72 72 N
## S gN a73 t74 0.9551051 N 73 74 73 73 N
## S gN a74 t75 0.9551051 N 74 75 74 74 N
## S gN a75 t76 0.9551051 N 75 76 75 75 N
## S gN a76 t77 0.9551051 N 76 77 76 76 N
## S gN a77 t78 0.9551051 N 77 78 77 77 N
## S gN a78 t79 0.9551051 N 78 79 78 78 N
## S gN a79 t80 0.9551051 N 79 80 79 79 N
## S gN a80 t81 0.9551051 N 80 81 80 80 N
## S gN a81 t82 0.9551051 N 81 82 81 81 N
## S gN a82 t83 0.9551051 N 82 83 82 82 N
## S gN a83 t84 0.9551051 N 83 84 83 83 N
## S gN a84 t85 0.9551051 N 84 85 84 84 N
## S gN a85 t86 0.9551051 N 85 86 85 85 N
## S gN a86 t87 0.9551051 N 86 87 86 86 N
## S gN a87 t88 0.9551051 N 87 88 87 87 N
## S gN a88 t89 0.9551051 N 88 89 88 88 N
## S gP a0 t1 0.9623833 P 0 1 0 0 P
## S gP a1 t2 0.9623833 P 1 2 1 1 P
## S gP a2 t3 0.9623833 P 2 3 2 2 P
## S gP a3 t4 0.9623833 P 3 4 3 3 P
## S gP a4 t5 0.9623833 P 4 5 4 4 P
## S gP a5 t6 0.9623833 P 5 6 5 5 P
## S gP a6 t7 0.9623833 P 6 7 6 6 P
## S gP a7 t8 0.9623833 P 7 8 7 7 P
## S gP a8 t9 0.9623833 P 8 9 8 8 P
## S gP a9 t10 0.9623833 P 9 10 9 9 P
## S gP a10 t11 0.9623833 P 10 11 10 10 P
## S gP a11 t12 0.9623833 P 11 12 11 11 P
## S gP a12 t13 0.9623833 P 12 13 12 12 P
## S gP a13 t14 0.9623833 P 13 14 13 13 P
## S gP a14 t15 0.9623833 P 14 15 14 14 P
## S gP a15 t16 0.9623833 P 15 16 15 15 P
## S gP a16 t17 0.9623833 P 16 17 16 16 P
## S gP a17 t18 0.9623833 P 17 18 17 17 P
## S gP a18 t19 0.9623833 P 18 19 18 18 P
## S gP a19 t20 0.9623833 P 19 20 19 19 P
## S gP a20 t21 0.9623833 P 20 21 20 20 P
## S gP a21 t22 0.9623833 P 21 22 21 21 P
## S gP a22 t23 0.9623833 P 22 23 22 22 P
## S gP a23 t24 0.9623833 P 23 24 23 23 P
## S gP a24 t25 0.9623833 P 24 25 24 24 P
## S gP a25 t26 0.9623833 P 25 26 25 25 P
## S gP a26 t27 0.9623833 P 26 27 26 26 P
## S gP a27 t28 0.9623833 P 27 28 27 27 P
## S gP a28 t29 0.9623833 P 28 29 28 28 P
## S gP a29 t30 0.9623833 P 29 30 29 29 P
## S gP a30 t31 0.9623833 P 30 31 30 30 P
## S gP a31 t32 0.9623833 P 31 32 31 31 P
## S gP a32 t33 0.9623833 P 32 33 32 32 P
## S gP a33 t34 0.9623833 P 33 34 33 33 P
## S gP a34 t35 0.9623833 P 34 35 34 34 P
## S gP a35 t36 0.9623833 P 35 36 35 35 P
## S gP a36 t37 0.9623833 P 36 37 36 36 P
## S gP a37 t38 0.9623833 P 37 38 37 37 P
## S gP a38 t39 0.9623833 P 38 39 38 38 P
## S gP a39 t40 0.9623833 P 39 40 39 39 P
## S gP a40 t41 0.9623833 P 40 41 40 40 P
## S gP a41 t42 0.9623833 P 41 42 41 41 P
## S gP a42 t43 0.9623833 P 42 43 42 42 P
## S gP a43 t44 0.9623833 P 43 44 43 43 P
## S gP a44 t45 0.9623833 P 44 45 44 44 P
## S gP a45 t46 0.9623833 P 45 46 45 45 P
## S gP a46 t47 0.9623833 P 46 47 46 46 P
## S gP a47 t48 0.9623833 P 47 48 47 47 P
## S gP a48 t49 0.9623833 P 48 49 48 48 P
## S gP a49 t50 0.9623833 P 49 50 49 49 P
## S gP a50 t51 0.9623833 P 50 51 50 50 P
## S gP a51 t52 0.9623833 P 51 52 51 51 P
## S gP a52 t53 0.9623833 P 52 53 52 52 P
## S gP a53 t54 0.9623833 P 53 54 53 53 P
## S gP a54 t55 0.9623833 P 54 55 54 54 P
## S gP a55 t56 0.9623833 P 55 56 55 55 P
## S gP a56 t57 0.9623833 P 56 57 56 56 P
## S gP a57 t58 0.9623833 P 57 58 57 57 P
## S gP a58 t59 0.9623833 P 58 59 58 58 P
## S gP a59 t60 0.9623833 P 59 60 59 59 P
## S gP a60 t61 0.9623833 P 60 61 60 60 P
## S gP a61 t62 0.9623833 P 61 62 61 61 P
## S gP a62 t63 0.9623833 P 62 63 62 62 P
## S gP a63 t64 0.9623833 P 63 64 63 63 P
## S gP a64 t65 0.9623833 P 64 65 64 64 P
## S gP a65 t66 0.9623833 P 65 66 65 65 P
## S gP a66 t67 0.9623833 P 66 67 66 66 P
## S gP a67 t68 0.9623833 P 67 68 67 67 P
## S gP a68 t69 0.9623833 P 68 69 68 68 P
## S gP a69 t70 0.9623833 P 69 70 69 69 P
## S gP a70 t71 0.9623833 P 70 71 70 70 P
## S gP a71 t72 0.9623833 P 71 72 71 71 P
## S gP a72 t73 0.9623833 P 72 73 72 72 P
## S gP a73 t74 0.9623833 P 73 74 73 73 P
## S gP a74 t75 0.9623833 P 74 75 74 74 P
## S gP a75 t76 0.9623833 P 75 76 75 75 P
## S gP a76 t77 0.9623833 P 76 77 76 76 P
## S gP a77 t78 0.9623833 P 77 78 77 77 P
## S gP a78 t79 0.9623833 P 78 79 78 78 P
## S gP a79 t80 0.9623833 P 79 80 79 79 P
## S gP a80 t81 0.9623833 P 80 81 80 80 P
## S gP a81 t82 0.9623833 P 81 82 81 81 P
## S gP a82 t83 0.9623833 P 82 83 82 82 P
## S gP a83 t84 0.9623833 P 83 84 83 83 P
## S gP a84 t85 0.9623833 P 84 85 84 84 P
## S gP a85 t86 0.9623833 P 85 86 85 85 P
## S gP a86 t87 0.9623833 P 86 87 86 86 P
## S gP a87 t88 0.9623833 P 87 88 87 87 P
## S gP a88 t89 0.9623833 P 88 89 88 88 P
## S gR a0 t1 0.9704853 R 0 1 0 0 R
## S gR a1 t2 0.9704853 R 1 2 1 1 R
## S gR a2 t3 0.9704853 R 2 3 2 2 R
## S gR a3 t4 0.9704853 R 3 4 3 3 R
## S gR a4 t5 0.9704853 R 4 5 4 4 R
## S gR a5 t6 0.9704853 R 5 6 5 5 R
## S gR a6 t7 0.9704853 R 6 7 6 6 R
## S gR a7 t8 0.9704853 R 7 8 7 7 R
## S gR a8 t9 0.9704853 R 8 9 8 8 R
## S gR a9 t10 0.9704853 R 9 10 9 9 R
## S gR a10 t11 0.9704853 R 10 11 10 10 R
## S gR a11 t12 0.9704853 R 11 12 11 11 R
## S gR a12 t13 0.9704853 R 12 13 12 12 R
## S gR a13 t14 0.9704853 R 13 14 13 13 R
## S gR a14 t15 0.9704853 R 14 15 14 14 R
## S gR a15 t16 0.9704853 R 15 16 15 15 R
## S gR a16 t17 0.9704853 R 16 17 16 16 R
## S gR a17 t18 0.9704853 R 17 18 17 17 R
## S gR a18 t19 0.9704853 R 18 19 18 18 R
## S gR a19 t20 0.9704853 R 19 20 19 19 R
## S gR a20 t21 0.9704853 R 20 21 20 20 R
## S gR a21 t22 0.9704853 R 21 22 21 21 R
## S gR a22 t23 0.9704853 R 22 23 22 22 R
## S gR a23 t24 0.9704853 R 23 24 23 23 R
## S gR a24 t25 0.9704853 R 24 25 24 24 R
## S gR a25 t26 0.9704853 R 25 26 25 25 R
## S gR a26 t27 0.9704853 R 26 27 26 26 R
## S gR a27 t28 0.9704853 R 27 28 27 27 R
## S gR a28 t29 0.9704853 R 28 29 28 28 R
## S gR a29 t30 0.9704853 R 29 30 29 29 R
## S gR a30 t31 0.9704853 R 30 31 30 30 R
## S gR a31 t32 0.9704853 R 31 32 31 31 R
## S gR a32 t33 0.9704853 R 32 33 32 32 R
## S gR a33 t34 0.9704853 R 33 34 33 33 R
## S gR a34 t35 0.9704853 R 34 35 34 34 R
## S gR a35 t36 0.9704853 R 35 36 35 35 R
## S gR a36 t37 0.9704853 R 36 37 36 36 R
## S gR a37 t38 0.9704853 R 37 38 37 37 R
## S gR a38 t39 0.9704853 R 38 39 38 38 R
## S gR a39 t40 0.9704853 R 39 40 39 39 R
## S gR a40 t41 0.9704853 R 40 41 40 40 R
## S gR a41 t42 0.9704853 R 41 42 41 41 R
## S gR a42 t43 0.9704853 R 42 43 42 42 R
## S gR a43 t44 0.9704853 R 43 44 43 43 R
## S gR a44 t45 0.9704853 R 44 45 44 44 R
## S gR a45 t46 0.9704853 R 45 46 45 45 R
## S gR a46 t47 0.9704853 R 46 47 46 46 R
## S gR a47 t48 0.9704853 R 47 48 47 47 R
## S gR a48 t49 0.9704853 R 48 49 48 48 R
## S gR a49 t50 0.9704853 R 49 50 49 49 R
## S gR a50 t51 0.9704853 R 50 51 50 50 R
## S gR a51 t52 0.9704853 R 51 52 51 51 R
## S gR a52 t53 0.9704853 R 52 53 52 52 R
## S gR a53 t54 0.9704853 R 53 54 53 53 R
## S gR a54 t55 0.9704853 R 54 55 54 54 R
## S gR a55 t56 0.9704853 R 55 56 55 55 R
## S gR a56 t57 0.9704853 R 56 57 56 56 R
## S gR a57 t58 0.9704853 R 57 58 57 57 R
## S gR a58 t59 0.9704853 R 58 59 58 58 R
## S gR a59 t60 0.9704853 R 59 60 59 59 R
## S gR a60 t61 0.9704853 R 60 61 60 60 R
## S gR a61 t62 0.9704853 R 61 62 61 61 R
## S gR a62 t63 0.9704853 R 62 63 62 62 R
## S gR a63 t64 0.9704853 R 63 64 63 63 R
## S gR a64 t65 0.9704853 R 64 65 64 64 R
## S gR a65 t66 0.9704853 R 65 66 65 65 R
## S gR a66 t67 0.9704853 R 66 67 66 66 R
## S gR a67 t68 0.9704853 R 67 68 67 67 R
## S gR a68 t69 0.9704853 R 68 69 68 68 R
## S gR a69 t70 0.9704853 R 69 70 69 69 R
## S gR a70 t71 0.9704853 R 70 71 70 70 R
## S gR a71 t72 0.9704853 R 71 72 71 71 R
## S gR a72 t73 0.9704853 R 72 73 72 72 R
## S gR a73 t74 0.9704853 R 73 74 73 73 R
## S gR a74 t75 0.9704853 R 74 75 74 74 R
## S gR a75 t76 0.9704853 R 75 76 75 75 R
## S gR a76 t77 0.9704853 R 76 77 76 76 R
## S gR a77 t78 0.9704853 R 77 78 77 77 R
## S gR a78 t79 0.9704853 R 78 79 78 78 R
## S gR a79 t80 0.9704853 R 79 80 79 79 R
## S gR a80 t81 0.9704853 R 80 81 80 80 R
## S gR a81 t82 0.9704853 R 81 82 81 81 R
## S gR a82 t83 0.9704853 R 82 83 82 82 R
## S gR a83 t84 0.9704853 R 83 84 83 83 R
## S gR a84 t85 0.9704853 R 84 85 84 84 R
## S gR a85 t86 0.9704853 R 85 86 85 85 R
## S gR a86 t87 0.9704853 R 86 87 86 86 R
## S gR a87 t88 0.9704853 R 87 88 87 87 R
## S gR a88 t89 0.9704853 R 88 89 88 88 R
## S gW a0 t1 0.9675738 W 0 1 0 0 W
## S gW a1 t2 0.9675738 W 1 2 1 1 W
## S gW a2 t3 0.9675738 W 2 3 2 2 W
## S gW a3 t4 0.9675738 W 3 4 3 3 W
## S gW a4 t5 0.9675738 W 4 5 4 4 W
## S gW a5 t6 0.9675738 W 5 6 5 5 W
## S gW a6 t7 0.9675738 W 6 7 6 6 W
## S gW a7 t8 0.9675738 W 7 8 7 7 W
## S gW a8 t9 0.9675738 W 8 9 8 8 W
## S gW a9 t10 0.9675738 W 9 10 9 9 W
## S gW a10 t11 0.9675738 W 10 11 10 10 W
## S gW a11 t12 0.9675738 W 11 12 11 11 W
## S gW a12 t13 0.9675738 W 12 13 12 12 W
## S gW a13 t14 0.9675738 W 13 14 13 13 W
## S gW a14 t15 0.9675738 W 14 15 14 14 W
## S gW a15 t16 0.9675738 W 15 16 15 15 W
## S gW a16 t17 0.9675738 W 16 17 16 16 W
## S gW a17 t18 0.9675738 W 17 18 17 17 W
## S gW a18 t19 0.9675738 W 18 19 18 18 W
## S gW a19 t20 0.9675738 W 19 20 19 19 W
## S gW a20 t21 0.9675738 W 20 21 20 20 W
## S gW a21 t22 0.9675738 W 21 22 21 21 W
## S gW a22 t23 0.9675738 W 22 23 22 22 W
## S gW a23 t24 0.9675738 W 23 24 23 23 W
## S gW a24 t25 0.9675738 W 24 25 24 24 W
## S gW a25 t26 0.9675738 W 25 26 25 25 W
## S gW a26 t27 0.9675738 W 26 27 26 26 W
## S gW a27 t28 0.9675738 W 27 28 27 27 W
## S gW a28 t29 0.9675738 W 28 29 28 28 W
## S gW a29 t30 0.9675738 W 29 30 29 29 W
## S gW a30 t31 0.9675738 W 30 31 30 30 W
## S gW a31 t32 0.9675738 W 31 32 31 31 W
## S gW a32 t33 0.9675738 W 32 33 32 32 W
## S gW a33 t34 0.9675738 W 33 34 33 33 W
## S gW a34 t35 0.9675738 W 34 35 34 34 W
## S gW a35 t36 0.9675738 W 35 36 35 35 W
## S gW a36 t37 0.9675738 W 36 37 36 36 W
## S gW a37 t38 0.9675738 W 37 38 37 37 W
## S gW a38 t39 0.9675738 W 38 39 38 38 W
## S gW a39 t40 0.9675738 W 39 40 39 39 W
## S gW a40 t41 0.9675738 W 40 41 40 40 W
## S gW a41 t42 0.9675738 W 41 42 41 41 W
## S gW a42 t43 0.9675738 W 42 43 42 42 W
## S gW a43 t44 0.9675738 W 43 44 43 43 W
## S gW a44 t45 0.9675738 W 44 45 44 44 W
## S gW a45 t46 0.9675738 W 45 46 45 45 W
## S gW a46 t47 0.9675738 W 46 47 46 46 W
## S gW a47 t48 0.9675738 W 47 48 47 47 W
## S gW a48 t49 0.9675738 W 48 49 48 48 W
## S gW a49 t50 0.9675738 W 49 50 49 49 W
## S gW a50 t51 0.9675738 W 50 51 50 50 W
## S gW a51 t52 0.9675738 W 51 52 51 51 W
## S gW a52 t53 0.9675738 W 52 53 52 52 W
## S gW a53 t54 0.9675738 W 53 54 53 53 W
## S gW a54 t55 0.9675738 W 54 55 54 54 W
## S gW a55 t56 0.9675738 W 55 56 55 55 W
## S gW a56 t57 0.9675738 W 56 57 56 56 W
## S gW a57 t58 0.9675738 W 57 58 57 57 W
## S gW a58 t59 0.9675738 W 58 59 58 58 W
## S gW a59 t60 0.9675738 W 59 60 59 59 W
## S gW a60 t61 0.9675738 W 60 61 60 60 W
## S gW a61 t62 0.9675738 W 61 62 61 61 W
## S gW a62 t63 0.9675738 W 62 63 62 62 W
## S gW a63 t64 0.9675738 W 63 64 63 63 W
## S gW a64 t65 0.9675738 W 64 65 64 64 W
## S gW a65 t66 0.9675738 W 65 66 65 65 W
## S gW a66 t67 0.9675738 W 66 67 66 66 W
## S gW a67 t68 0.9675738 W 67 68 67 67 W
## S gW a68 t69 0.9675738 W 68 69 68 68 W
## S gW a69 t70 0.9675738 W 69 70 69 69 W
## S gW a70 t71 0.9675738 W 70 71 70 70 W
## S gW a71 t72 0.9675738 W 71 72 71 71 W
## S gW a72 t73 0.9675738 W 72 73 72 72 W
## S gW a73 t74 0.9675738 W 73 74 73 73 W
## S gW a74 t75 0.9675738 W 74 75 74 74 W
## S gW a75 t76 0.9675738 W 75 76 75 75 W
## S gW a76 t77 0.9675738 W 76 77 76 76 W
## S gW a77 t78 0.9675738 W 77 78 77 77 W
## S gW a78 t79 0.9675738 W 78 79 78 78 W
## S gW a79 t80 0.9675738 W 79 80 79 79 W
## S gW a80 t81 0.9675738 W 80 81 80 80 W
## S gW a81 t82 0.9675738 W 81 82 81 81 W
## S gW a82 t83 0.9675738 W 82 83 82 82 W
## S gW a83 t84 0.9675738 W 83 84 83 83 W
## S gW a84 t85 0.9675738 W 84 85 84 84 W
## S gW a85 t86 0.9675738 W 85 86 85 85 W
## S gW a86 t87 0.9675738 W 86 87 86 86 W
## S gW a87 t88 0.9675738 W 87 88 87 87 W
## S gW a88 t89 0.9675738 W 88 89 88 88 W
# we plot the results
# because the DSR is the same for all days, we use the time=1 values
plotdata <- get.real(mod.hab, "S", se=TRUE)
head(plotdata)
## all.diff.index par.index estimate se lcl ucl
## S gN a0 t1 1 1 0.9459833 0.0050725 0.935134 0.9551051
## S gN a1 t2 2 1 0.9459833 0.0050725 0.935134 0.9551051
## S gN a2 t3 3 1 0.9459833 0.0050725 0.935134 0.9551051
## S gN a3 t4 4 1 0.9459833 0.0050725 0.935134 0.9551051
## S gN a4 t5 5 1 0.9459833 0.0050725 0.935134 0.9551051
## S gN a5 t6 6 1 0.9459833 0.0050725 0.935134 0.9551051
## fixed note group age time Age Time Habitat
## S gN a0 t1 N 0 1 0 0 N
## S gN a1 t2 N 1 2 1 1 N
## S gN a2 t3 N 2 3 2 2 N
## S gN a3 t4 N 3 4 3 3 N
## S gN a4 t5 N 4 5 4 4 N
## S gN a5 t6 N 5 6 5 5 N
plotdata <- plotdata[plotdata$time==1,]
plotdata
## all.diff.index par.index estimate se lcl
## S gN a0 t1 1 1 0.9459833 0.0050725 0.9351340
## S gP a0 t1 90 2 0.9562953 0.0033341 0.9492739
## S gR a0 t1 179 3 0.9558094 0.0090265 0.9343297
## S gW a0 t1 268 4 0.9505390 0.0105538 0.9252465
## ucl fixed note group age time Age Time Habitat
## S gN a0 t1 0.9551051 N 0 1 0 0 N
## S gP a0 t1 0.9623833 P 0 1 0 0 P
## S gR a0 t1 0.9704853 R 0 1 0 0 R
## S gW a0 t1 0.9675738 W 0 1 0 0 W
dsr.hab.plot <- ggplot(data=plotdata, aes(x=Habitat, y=estimate))+
ggtitle("DSR in different habitats for mallards")+
geom_point()+
geom_errorbar(aes(ymin=lcl, ymax=ucl), width=.1)+
ylab("DSR (95% ci)")
dsr.hab.plot
ggsave(dsr.hab.plot,
file=file.path("..","..","..","..","MyStuff","Images","mallard-dsr-habitat.png"), h=4, w=6, units="in", dpi=300)
# you can get the nest survival probability in the usual way by using the estimated DSRs.
# we want to do this for each habitat
ns <- plyr::ldply(unique(malldata$Habitat), function(habitat, ndays){
daily.dsr <- get.real(mod.hab, "S", se=TRUE, vcv=TRUE, expand=TRUE)
# this has all habitats together; select only habitat of interest
select <- daily.dsr$estimates$Habitat==habitat
daily.dsr$estimates <- daily.dsr$estimates[select,]
daily.dsr$vcv.real <- daily.dsr$vcv.real[select,select]
#browser()
ns <- prod(daily.dsr$estimates$estimate[1:ndays])
ns.se <- deltamethod.special("prod",
daily.dsr$estimates$estimate[1:ndays],
daily.dsr$vcv.real[1:ndays, 1:ndays])
data.frame(Habitat=habitat, ns=ns, ns.se=ns.se)
}, ndays=20)
ns
## Habitat ns ns.se
## 1 R 0.4049726 0.07648941
## 2 P 0.4091100 0.02852674
## 3 N 0.3293588 0.03532138
## 4 W 0.3625759 0.08051351
# fit a null models
# 3. Fit a particular model
# This is a model with S varying by habitat type
mod.null <- RMark::mark(mall.proc, ddl=mall.ddl,
model="Nest",
model.parameters=list(
S =list(formula=~1)
)
)
##
## Output summary for Nest model
## Name : S(~1)
##
## Npar : 1
## -2lnL: 1567.116
## AICc : 1569.117
##
## Beta
## estimate se lcl ucl
## S:(Intercept) 3.005653 0.0576808 2.892598 3.118707
##
##
## Real Parameter S
## 1 2 3 4 5 6
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 7 8 9 10 11 12
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 13 14 15 16 17 18
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 19 20 21 22 23 24
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 25 26 27 28 29 30
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 31 32 33 34 35 36
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 37 38 39 40 41 42
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 43 44 45 46 47 48
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 49 50 51 52 53 54
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 55 56 57 58 59 60
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 61 62 63 64 65 66
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 67 68 69 70 71 72
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 73 74 75 76 77 78
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 79 80 81 82 83 84
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 85 86 87 88 89
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
summary(mod.null)
## Output summary for Nest model
## Name : S(~1)
##
## Npar : 1
## -2lnL: 1567.116
## AICc : 1569.117
##
## Beta
## estimate se lcl ucl
## S:(Intercept) 3.005653 0.0576808 2.892598 3.118707
##
##
## Real Parameter S
## 1 2 3 4 5 6
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 7 8 9 10 11 12
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 13 14 15 16 17 18
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 19 20 21 22 23 24
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 25 26 27 28 29 30
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 31 32 33 34 35 36
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 37 38 39 40 41 42
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 43 44 45 46 47 48
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 49 50 51 52 53 54
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 55 56 57 58 59 60
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 61 62 63 64 65 66
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 67 68 69 70 71 72
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 73 74 75 76 77 78
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 79 80 81 82 83 84
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## 85 86 87 88 89
## Group:HabitatN 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatP 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatR 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
## Group:HabitatW 0.9528288 0.9528288 0.9528288 0.9528288 0.9528288
collect.models(type="Nest")
## model npar AICc DeltaAICc weight Deviance
## 2 S(~1) 1 1569.117 0.000000 0.805384 1567.116
## 1 S(~Habitat) 4 1571.957 2.840582 0.194616 1563.951
cleanup(ask=FALSE)