{"id":183,"date":"2017-04-18T02:32:52","date_gmt":"2017-04-18T02:32:52","guid":{"rendered":"http:\/\/sites.warnercnr.colostate.edu\/gwhite\/?page_id=183"},"modified":"2017-04-18T02:32:52","modified_gmt":"2017-04-18T02:32:52","slug":"model-structure","status":"publish","type":"page","link":"https:\/\/sites.warnercnr.colostate.edu\/gwhite\/model-structure\/","title":{"rendered":"Model Structure"},"content":{"rendered":"

Contents<\/a> – Index<\/a><\/span><\/p>\n

\n

Model Structure<\/b><\/span><\/p>\n

Program Mark provides parameter estimates for 142 data types<\/a>: Cormack-Jolly-Seber models (live animal recaptures that are released alive), band (ring) recovery models (dead animal recoveries), models with both live and dead re-encounters (Burnham’s model), known fate<\/a> models, nest survival<\/a> models, closed capture<\/a> models, band (ring) recovery models where the number of animals marked is unknown, robust design models for live recaptures<\/a>, Barker’s extension to Burnham’s model<\/a>, multi-strata live recapture model<\/a>, Brownie et al.’s model of band or ring recoveries, Jolly-Seber models<\/a> that include either seniority probability, recruitment rate, rate of population change, or probability of entry,\u00a0 probability of occupancy<\/a> models (including robust design occupancy models<\/a>), and mark-resight<\/a> models.\u00a0\u00a0<\/span><\/p>\n

A complete list of the currently available data types can be generated under the Help | Data Types menu selections.<\/span><\/p>\n

Estimation in Cormack-Jolly-Seber Designs<\/b><\/span><\/p>\n

Live Recaptures.<\/b>\u00a0\u00a0 Live recaptures are the basis of the standard Cormack-Jolly-Seber model. Marked animals are released into the population, often by trapping them from the populations.\u00a0 Then, marked animals are encountered by catching them alive and re-releasing them.\u00a0 If marked animals are released into the population on occasion 1, then each succeeding capture occasion is one encounter occasion.\u00a0 Consider the following scenario:<\/span><\/p>\n

\u00a0\u00a0\u00a0 Release —-S(1)—–> Encounter 2 ——-S(2)——> Encounter 3<\/span><\/p>\n

Animals survive from initial release to the second encounter occasion with probability S(1), and from second encounter occasion to the third encounter occasion with probability S(2).\u00a0 The recapture probability at encounter occasion 2 is p(2), and p(3) is the recapture probability at encounter occasion 3.\u00a0 At least 2 re-encounter occasions are required to estimate the survival rate between the first release occasion and the first re-encounter occasion, i.e., S(1). The survival rate between the last two encounter occasions is not estimable because only the product of survival and recapture probability for this occasion is identifiable.<\/span><\/p>\n

Generally, the survival rates of the CJS model are labeled as Phi(1), Phi(2), etc., because the quantity estimated is the probability of remaining available for recapture.\u00a0 Thus, animals that emigrate from the study area are not available for recapture, so appear to have died in this model.\u00a0 Thus, Phi(i) = S(i)(1 – E(i)), where E(i) is the probability of emigrating from the study area.<\/span><\/p>\n

Estimates of population size (N) or births and immigration (B) of the Jolly-Seber model are not provided in the CJS model of Program MARK.\u00a0 See the section on Jolly-Seber models below for these models in MARK. Programs that estimate population size for each occasion (where the quantity is identifiable) are POPAN-5\u00a0 (Arnason and Schwarz 1995<\/a>, Schwarz and Arnason 1996<\/a>) or JOLLY and JOLLYAGE\u00a0 (Pollock et al. 1990<\/a>).<\/span><\/p>\n

The modification of the standard Cormack-Jolly-Seber to allow estimation of different apparent survival rates for transients (Pradel et al. 1997<\/a>) can be developed using age-structured PIMs<\/a>.\u00a0 An examples using the Lazuli Bunting data is provided in the transients<\/a> entry.<\/span><\/p>\n

Estimation in Band\/Ring Recovery Designs<\/b><\/span><\/p>\n

Dead Recoveries.<\/b>\u00a0 With dead recoveries, i.e., band, fish tag, or ring recovery models, animals are captured from the population, marked, and\u00a0 released back into the population at each occasion.\u00a0 Later, marked animals are encountered as dead animals, typically from harvest or just found dead (e.g., gulls).<\/span><\/p>\n

Marked animals are assumed to survive from one release to the next with survival probability S(i).\u00a0 If they die, the dead marked animals are reported during each period between releases with probability r(i)\u00a0 Note that r(i) is called a “reporting probability”, but this is not the probability that a hunter reports the marked animal.\u00a0 Rather, r(i) is the probability that a marked animal is reported conditional on its death.\u00a0 Animals that die of natural causes would have a very low probability of being found and the mark reported.\u00a0 In contrast, harvested animals have a greater chance of being reported, but still not probability 1.\u00a0 Don’t make the mistake of equating r(i) with the probability that an animal is harvested.\u00a0 However, recognize that shifts in the mortality process do affect the estimates of r(i).\u00a0 For example, heavily harvested populations should have higher values of r than lightly harvested populations, because the probability that an animal dies from harvest is higher, and hence a greater probability that the mark is reported.\u00a0 But even though r relates to the probability of harvest, r is not the probability that a hunter reports the marked animal.\u00a0 Otis and White (2002)<\/a> provide further discussion on the interpretation of r versus the probability that a hunter reports the band, including an equation relating the 2 variables<\/a>.<\/span><\/p>\n

The survival probability and reporting probability prior to the last release can not be estimated individually in the full time-effects model, but only as a product.\u00a0 This parameterization differs from that of\u00a0 Brownie et al. (1985)<\/a> in that their f(i) is replaced as f(i) = (1 – S(i)) r(i).\u00a0 The r(i) are equivalent to the lambda(i) of\u00a0 Seber (1970)<\/a>, where the original description of this model was developed,\u00a0 and of life table models\u00a0 Anderson et al. 1985; Catchpole et al. 1995)<\/a>.\u00a0 The reason for making this change is so that the encounter process, modeled with the r(i) parameters, can be separated from the survival process, modeled with the S(i) parameters.\u00a0 With the f(i) parameterization, the 2 processes are both part of this parameter.\u00a0 Hence, developing more advanced models with the design matrix options of MARK is difficult, if not illogical with the f(i) parameterization.\u00a0 However, the negative side of this new parameterization is that the last S(i) and r(i) are confounded in the full time-effects model, as only the product (1 – S(i)) r(i)\u00a0 is identifiable, and hence estimable.\u00a0 Secondly, all the parameters are bounded between zero and 1, which seems like a benefit.\u00a0 However, parameter estimates at the boundary do not have proper estimates of their standard errors.\u00a0 An equivalent situation occurs with the binomial distribution when either no successes occur in the data, or all successes occur in the data, and the standard error is estimated as zero.\u00a0 Because the Brownie et al. parameterization overcomes these difficulties, it is also included in Program MARK, and is described below.<\/span><\/p>\n

Brownie et al. Dead Recoveries Model.<\/b>\u00a0 This model is the band or ring recovery model of Brownie et al. (1985)<\/a> with the S<\/i> and f<\/i> coding.\u00a0 The model gives the same estimates for fully time-specific models as does the S<\/i> and r<\/i> coding, except when estimates of S<\/i> are >1.\u00a0 However, different estimates will be obtained with covariates used to model survival or the recovery process.\u00a0 The advantage of this data type over the Dead Recoveries data type is that all the parameters are estimable under the fully time-specific model, and parameters are not constrained to the interval [0, 1], so that valid standard errors can be estimated with the identity or log link functions.\u00a0 More details on the differences in the S<\/i>, r<\/i> and S<\/i>, f<\/i> band recovery models are given by Otis and White (2002)<\/a>, in a discussion of the band reporting rate<\/a>.<\/span><\/p>\n

Recoveries Without Knowing Number Marked.<\/b>\u00a0 The British Trust for Ornithology (BTO) does not have computerized databases of the numbers of birds ringed.\u00a0 Thus, BTO cannot compute the cohort size for a set of ring recoveries.\u00a0 To circumvent this problem, a ring recovery model is formulated where the recovery rate (r<\/i>(i<\/i>) is assumed constant by age class and year.\u00a0 Then, the survival rate can be estimated from the observed recoveries.\u00a0 The cell probability for the j<\/i> year of recoveries given k<\/i> years of recoveries is<\/span><\/p>\n

S(1) S(2) … S(j-1) [1 – S(j)] \/ [1 – S(1) S(2) … S(k)]<\/span><\/p>\n

where the denominator is 1 minus the probability of still being alive.\u00a0 Of the k survival rates, only k-1 are identifiable.\u00a0 Common approaches to achieve identifiability are to set S(k-1) = S(k) or to set S(k) to the mean of S(1) … S(k-1) using appropriate constraints in the Design Matrix<\/a>.<\/span><\/p>\n

This model should only be used when you do not know the number of animals marked because you cannot evaluate the assumption of constant recovery rates with this model.\u00a0 If you know the number of animals marked, use the dead recovery model described above.<\/span><\/p>\n

Estimation in Live and Dead Encounters Designs<\/b><\/span><\/p>\n

Both Live and Dead Encounters.<\/b>\u00a0 The model for the joint live and dead encounter data type was first published by Burnham (1993)<\/a>, but with a slightly different parameterization than used in Program MARK.\u00a0 In MARK, the dead encounters are not modeled with the f(i) of Burnham (1993), but rather as f<\/i>(i<\/i>) = (1 – S<\/i>(i<\/i>))r<\/i>(i<\/i>), as discussed above for the dead encounter models.\u00a0 The method is a combination of the 2 above, but allows the estimation of fidelity (F<\/i>(i<\/i>) = 1 – E<\/i>(i<\/i>)), or the probability that the animal remains on the study area and is available for capture.\u00a0 As a result, the estimates of S<\/i>(i<\/i>) are estimates of the survival probability of the marked animals, and not the apparent survival (phi(i<\/i>) = S<\/i>(i<\/i>) F<\/i>(i<\/i>)) as discussed for the live encounter model.<\/span><\/p>\n

In the models discussed so far, live captures and resightings modeled with the p<\/i>(i<\/i>) parameters are assumed to occur over a short time interval, whereas dead recoveries modeled with the r<\/i>(i<\/i>)\u00a0 parameters extend over the time interval.\u00a0 The actual time of the dead recovery is not used in the estimation of survival for 2 reasons.\u00a0 First, it is often not known.\u00a0 Second, even if the exact time of recovery is known, little information is contributed if the recovery probability (r<\/i>(i<\/i>)) is varying during the time interval.<\/span><\/p>\n

Barker’s model, discussed below, is an extension of the both live and dead encounters model that uses information on live resightings between live recapture intervals to improve estimates of survival.<\/span><\/p>\n

Barker’s Model.<\/b>\u00a0 Richard Barker (1997, 1999)<\/a> has extended the model of Burnham (1993)<\/a> by allowing resightings of marked animals during the interval between trapping occasions.\u00a0 The model was motivated by a brown trout study, where fish were marked at regular intervals, but which were then caught by fisherman.\u00a0 The marked trout was considered “resighted and released” if the fisherman released the fish alive, and “resighted and killed” if the fisherman kept the fish.\u00a0 Additional parameters of this model from Burnham’s model are R<\/i>, R<\/i>‘, and F<\/i>‘.\u00a0 This model is particularly useful for situations where live sightings of marked animals are obtained between marking periods.\u00a0 If no dead encounters are recorded, the r<\/i>(i<\/i>) parameters can be set to zero.\u00a0 More details are provided on Barker’s model here<\/a>.<\/span><\/p>\n

Estimation in Radio-tracking and Known Fates Designs<\/b><\/span><\/p>\n

Known Fates.<\/b>\u00a0 Known fate data assumes that there are no nuisance parameters involved with animal captures or resightings.\u00a0 The data derive from radio-tracking studies, although some radio-tracking studies fail to follow all the marked animals and so would not meet the assumptions of this model.\u00a0 A diagram illustrating this scenario is<\/span><\/p>\n

\u00a0\u00a0\u00a0 Release\u00a0 —–S(1)—-> Encounter 2\u00a0 —–S(2)—-> Encounter 3\u00a0 —–S(3)—-> Encounter 4 …<\/span><\/p>\n

where the probability of encounter on each occasion is 1 if the animal is alive.\u00a0\u00a0<\/span><\/p>\n

The Known Fate<\/a> data type is equivalent to the Kaplan-Meier estimation method if only one mortality occurs per interval.\u00a0 More importantly, the Known Fate data type provides important advantages over the Kaplan-Meier approach: incorporation of covariates, and selection between competing models.\u00a0 However, the Known Fates data type retains the advantages of the Kaplan-Meier method: left and right censoring, staggered entry, and no assumption about the hazard rate during the interval.\u00a0 In contrast, the Heisey-Fuller model assumes a constant hazard function within an interval, which is often a detrimental feature of the approach.<\/span><\/p>\n

Nest Survival Models.<\/b>\u00a0 This model is different than the known fate model because the exact day that the animal (nest) dies is unknown. The nest survival<\/a> data type is appropriate for known fate data where the occasions are not clearly delineated.\u00a0 As a result, nest survival models provide a means of analyzing ragged radio-tracking data.\u00a0 The data type provides a survival parameter for each day (occasion) of the study.\u00a0 Typically, not all of these parameters would be estimated, but rather models would be constructed that provided some structure across these parameters, such as trend or trend^2 models.\u00a0 More details on the format of how to enter the data and the structure of the model are given here<\/a>.\u00a0 If data are entered as encounter histories, the LDLDLDLD format is required.<\/span><\/p>\n

The nest survival data type is equivalent to the Mayfield estimator if the occasion-specific survival rates are all assumed to be constant.\u00a0 The nest survival model in MARK can emulate the Heisey-Fuller if survival rates are set constant over the same intervals as in the Heisey-Fuller method.<\/span><\/p>\n

Estimation of Population Size in Closed Populations<\/b><\/span><\/p>\n

Closed Captures.<\/b>\u00a0 Closed-capture data assume that all survival probabilities are 1.0 across the short time intervals of the study (which are assumed to have zero length).\u00a0\u00a0 Because time intervals are defined to be too small for any mortality or emigration, you are not allowed to enter Time Interval<\/a> length.\u00a0 Thus, survival is not estimated.\u00a0 Rather, the probability of first capture (p<\/i>(i<\/i>)) and the probability of recapture (c<\/i>(i<\/i>)) are estimated, along with the number of animals in the population (N<\/i>).\u00a0 This data type is the same as is analyzed with Program CAPTURE (White et al. 1982)<\/a>.\u00a0 All the likelihood models in CAPTURE can be duplicated in MARK.\u00a0 However, MARK allows additional models not available in CAPTURE, plus comparisons between groups and the incorporation of time-specific and\/or group-specific covariates into the model.\u00a0 A total of 6 different closed captures models<\/a> are available in MARK.<\/span><\/p>\n

Program MARK also allows models incorporating individual heterogeneity for closed capture-recapture<\/a> models.\u00a0 These models are developed based on a mixture distribution of first capture and recapture parameters following Norris and Pollock (1995)<\/a> and Pledger (1998, 1999)<\/a>,\u00a0 Models for Mh<\/a> and Mtbh<\/a> are available.\u00a0\u00a0<\/span><\/p>\n

Individual covariates<\/a> cannot be used with the closed captures data type because animals that were never captured (and hence, whose individual covariates could never be measured) are incorporated into the likelihood as part of the estimate of population size (N<\/i>).\u00a0 Models that can incorporate individual covariates existing in the literature (Huggins 1989, 1991; Alho 1990)<\/a> have been implemented in MARK, and are described below.\u00a0 Estimates of population size are given for the Huggins’ models, but these estimates are not quite as efficient as the closed captures data type where the statistical models are equivalent to those in Program CAPTURE.\u00a0 However, the ability to incorporate individual covariates makes the Huggins’ models more appropriate if individual heterogeneity exists in the data.<\/span><\/p>\n

Huggins’ Closed Captures.<\/b>\u00a0\u00a0 Huggins’ model (Huggins 1989, 1991; Alho 1990)<\/a> allows estimation of closed population size (N<\/i>) from initial capture probabilities (p<\/i>) and recapture probabilities (c<\/i>).\u00a0 The model conditions on the animal being captured at least once during the study, so allows individual covariates<\/a> to be used to model p<\/i> and c<\/i>.\u00a0 The approach used in Huggins’ model is equivalent to the Horvitz-Thompson sampling design, where animals have unequal probability of being included in the sample.\u00a0 Only LLLL encounter histories are required for this model.<\/span><\/p>\n

Closed Captures with Heterogeneity.<\/b>\u00a0 This model is closely linked with the closed captures models.\u00a0 Pertinent literature includes Norris and Pollock (1995)<\/a> and Pledger (1998, 1999)<\/a>.\u00a0\u00a0 More details are provided here<\/a>.\u00a0\u00a0\u00a0 Only LLLL encounter histories are required for this model.<\/span><\/p>\n

Full Closed Captures with Heterogeneity.<\/b>\u00a0 Again, this model is closely linked with the closed captures models.\u00a0 Pertinent literature includes Norris and Pollock (1995)<\/a> and Pledger (1998, 1999)<\/a>.\u00a0\u00a0 More details are provided here<\/a>.\u00a0\u00a0\u00a0 Only LLLL encounter histories are required for this model.<\/span><\/p>\n

Mark-Resight Models.<\/b>\u00a0 The mark-resight models<\/a> allow population estimation when unmarked animals are not marked when resighted, but encounter histories are formed for a set of marked animals.\u00a0 The models in MARK have been developed by McClintock(in prep.)<\/a> and extend the estimators in the NOREMARK software package (White 1996)<\/a>.<\/span><\/p>\n

Estimation of Robust Designs<\/b><\/span><\/p>\n

Robust Design.<\/b>\u00a0 This model is a combination of the CJS live recapture model and the closed capture models, and is described in detail by Kendall et al. (1997, 1995)<\/a> and Kendall and Nichols (1995)<\/a>.\u00a0 Instead of just 1 capture occasion between survival intervals, multiple (>1) capture\u00a0 occasions are used that are close together in time.\u00a0 These closely-spaced encounter occasions are termed “sessions”.\u00a0 To specify the encounter sessions, the Time Interval<\/a> lengths are used.\u00a0 The time intervals between the encounter occasions within a session\u00a0 have a length of zero, whereas the time intervals between sessions have a positive (>0) length.\u00a0 An example will make this clearer.\u00a0 Assume that animals are trapped for 15 separate times.\u00a0 The first year, animals are trapped for 2 days, the second year for 2 days, the third year for 4 days, the fourth year for 5 days, and the fifth year for 2 days.\u00a0 The number of encounter occasions would be specified as 15.\u00a0 The length of the time intervals would be specified as: 0,1,0,1,0,0,0,1,0,0,0,0,1,0.\u00a0 That is, only 14 time intervals are needed, where the value 1 means that 1 year elapsed.\u00a0 This mechanism is flexible, but can be a bit\u00a0 tricky.\u00a0 Note that all sessions must have at least 2 occasions.\u00a0 Thus, you will never have 2 consecutive time intervals of length >0.<\/span><\/p>\n

For each trapping session (j<\/i>), the probability of first capture (p<\/i>(ji<\/i>)) and the probability of recapture (c<\/i>(ji<\/i>)) are estimated (where i<\/i> indexes the number of trapping occasions within the session), along with the number of animals in the population (N<\/i>(j<\/i>)).\u00a0 For the intervals between sessions, the probability of survival (S<\/i>(j<\/i>)), the probability of emigration from the study area (gamma’ ‘ (j<\/i>)), and the probability of staying away from the study area (gamma’ (j<\/i>)) are estimated.\u00a0 Indexing of these parameters follows the notation of Kendall et al. (1997)<\/a>.\u00a0 Thus, gamma’ ‘(2) applies to the second trapping session, and gamma’ (2) is not estimated because there are no marked animals outside the study area at that time.\u00a0 To provide identifiability of the parameters for the Markovian emigration model, Kendall et al. (1997)<\/a> suggest setting gamma’ ‘ (k<\/i> – 1) = gamma’ ‘(k<\/i>) and gamma'(k<\/i> – 1) = gamma'(k<\/i>), where k<\/i> is the number of primary trapping sessions.\u00a0\u00a0\u00a0 To obtain the “No Emigration” model, set all the gamma parameters to zero.\u00a0 To obtain the “Random Emigration” model, set gamma'(i<\/i>) = gamma’ ‘(i<\/i>).<\/span><\/p>\n

The robust design models in MARK can all incorporate individual heterogeneity closed capture data type<\/a> in the estimation of population size.\u00a0 Individual covariates<\/a> can be used to model the parameters S<\/i>, gamma’ ‘, and gamma’ in the Robust Design data type.\u00a0 Individual covariates<\/a> cannot be used with the Robust Design data type for the p<\/i>‘s, c<\/i>‘s, and N<\/i>‘s because animals that were never captured (and hence, whose individual covariates could never be measured) are incorporated into the likelihood as part of the estimate of population size (N).\u00a0 Models that can incorporate individual covariates existing in the literature (Huggins 1989, 1991; Alho 1990)<\/a> have been implemented in MARK (including the heterogeneity models), and are described below for the data type Robust Design (Huggins Est.).\u00a0 Estimates of population size are given for the Huggins’ models, but these estimates are not quite as efficient as the closed captures data type where the statistical models are equivalent to those in Program CAPTURE.\u00a0 However, the ability to incorporate individual covariates makes the Huggins’ models more appropriate if individual heterogeneity exists in the data.<\/span><\/p>\n

More details are provided on the robust design model here<\/a>.<\/span><\/p>\n

Robust Design (Huggins Est.).<\/b>\u00a0 The robust design model has also been extended to include Huggins’ estimator for population size (N<\/i>) for each trapping session (Huggins 1989, 1991; Alho 1990)<\/a>.\u00a0 Again, individual covariates can be used to model the initial capture probabilities (p<\/i>) and recapture probabilities (c<\/i>) for each trapping session.\u00a0 Only LLLL encounter histories are required for this model.<\/span><\/p>\n

Barker’s Model Robust Design.<\/b>\u00a0 Bill Kendall and Richard Barker extended Barker’s model to handle the robust design.\u00a0 This model is an extension of the Lindberg et al. (2001)<\/a> model because it uses encounter information from live recaptures, dead recoveries, and resightings between the intervals of live captures.\u00a0 LDLDLDLD encounter histories are required for this model, with resighting between the capture intervals given the value 2 for D.<\/span><\/p>\n

Estimation of Multi-state Designs<\/b><\/span><\/p>\n

Multi-state Model for Live Recaptures.<\/b>\u00a0 The multi-state model of Brownie et al. (1993)<\/a> and Hestbeck et al. (1991)<\/a> allows animals to move between states with transition probabilities.\u00a0 At this time, only the movement model without memory is implemented.\u00a0 More details are provided on the multi-state model here<\/a>.\u00a0 The multi-state model has also been extended to incorporate dead recoveries, described below, and has been extended to incorporate the robust design, both the open<\/a> and closed<\/a> robust design multi-state models.\u00a0 In addition, multi-state models with <\/span>state uncertainty<\/a><\/span> are available.<\/span><\/p>\n

Live and Dead Multi-strata Model.<\/b>\u00a0 The multi-strata model that incorporates both live and dead recoveries is available and described here<\/a>.\u00a0<\/span><\/p>\n

Estimation of Jolly-Seber Designs<\/b><\/span><\/p>\n

Jolly-Seber Models.<\/b>\u00a0\u00a0 In addition to the apparent survival and recapture probabilities of the Cormack-Jolly-Seber model (recaptures only model), the Jolly-Seber model allows estimation of the population size (N<\/i>) at each trapping occasion, plus the number of new animals entering the population (B<\/i>) at each occasion.\u00a0 Multiple parameterizations of the Jolly-Seber model in MARK.\u00a0 For all of these parameterizations, only the LLLL encounter histories<\/a> are required.\u00a0 The following parameterizations of the Jolly-Seber model are available in Program MARK: Burnham, Pradel (3 parameterizations: gamma, f<\/i>, and lambda), POPAN from Schwarz and Arnason (1996)<\/a>, and Link and Barker (2003)<\/a>.\u00a0 The relationships between the parameters of these models are given here<\/a>.\u00a0 Also, for the population change rates to be meaningful, the study area size must not change during the study.\u00a0 See Population Rate of Change<\/a> for more discussion of this point.<\/span><\/p>\n

Burnham’s Jolly-Seber Model.<\/b>\u00a0 This parameterization provides the population size at the start of the study, plus the rate of population change (lambda) for each interval.\u00a0 This model can be difficult to get numerical convergence of the parameter estimates.\u00a0 Although this model has been thoroughly checked, and found to be correct, the program has difficulty obtaining numerical solutions for the parameters because of the penalty constraints<\/a> required to keep the parameters consistent with each other.<\/span><\/p>\n

Pradel Recruitment Only Model.<\/b>\u00a0\u00a0 Pradel (1996<\/a>) developed a model to estimate the proportion of the population that was previously in the population.\u00a0 Thus, this model, labeled ‘Pradel Recruitment Only’, estimates recruitment to the population.\u00a0 The parameters of this model are the seniority probability, gamma (probability that an animal present at time i<\/i> was already present at time i<\/i> – 1), and recapture probability r<\/i>.\u00a0 Only LLLL encounter histories are required for this model.\u00a0 This model can be estimated by reversing the time sequence of the live encounter histories (Pradel 1996)<\/a>, an idea suggested by Pollock et al. (1974:85-85)<\/a>, and even mentioned by R. A. Fisher in about 1939 or so (Box ????).<\/span><\/p>\n

Pradel Survival and Seniority Model.<\/b>\u00a0\u00a0 Pradel (1996)<\/a> extended his recruitment only model to include apparent survival (phi).\u00a0 In MARK, this model is labeled ‘Pradel Survival and Seniority’.\u00a0 Parameters of the model are apparent survival (phi), recapture probability (p<\/i>), and seniority probability (gamma), which is the probability that an animal in the population at time i<\/i> was also in the population at time i<\/i> – 1 (i.e., the animal did not enter the population during the interval i <\/i>– 1 to i<\/i>.\u00a0 Only LLLL encounter histories are required for this model.<\/span><\/p>\n

Pradel Survival and Lambda Model.<\/b>\u00a0\u00a0 Pradel (1996)<\/a> also parameterized his model with both recruitment and apparent survival to have the parameters apparent survival (phi), recapture probability (p), and rate of population change [lambda(i<\/i>) = N<\/i>(i <\/i>+ 1)\/N<\/i>(i<\/i>)]).\u00a0 This model converges quite readily compared to the Burnham parameterization of the Jolly-Seber model described above.\u00a0 Only LLLL encounter histories are required for this model.<\/span><\/p>\n

Pradel Survival and Recruitment Model.<\/b>\u00a0\u00a0 Pradel (1996)<\/a> also parameterized his model with both recruitment and apparent survival to have the parameters apparent survival (phi), recapture probability (p<\/i>), and fecundity rate [f<\/i>(i<\/i>) = number of adults at time i<\/i> + 1 per adult at time i<\/i>].\u00a0 This model converges quite readily.\u00a0 Only LLLL encounter histories are required for this model.<\/span><\/p>\n

POPAN Model.<\/b>\u00a0 Schwarz and Arnason (1996) parameterized the Jolly-Seber model in terms of a super population (N), and the probability of entry (pent in MARK, beta in the paper).\u00a0 The POPAN data type implements this model.\u00a0 The MLogit link function<\/a> provides a constraint that makes the sum of the pent parameters <=1, with the probability of occurring in the population on the first occasion as 1 – sum(pent(t)).\u00a0 More details are given here<\/a>.<\/span><\/p>\n

Link-Barker Model.<\/b>\u00a0 Link and Barker (2003) reparameterized the POPAN model from the probability of entry to the recruitment parameter (f).\u00a0 The reason for this reparameterization was to provide a more biologically meaningful interpretation of the parameters of the model, as part of a hierarchical modeling approach.\u00a0 More details are given here<\/a>.<\/span><\/p>\n

Estimation of Virtual Population Analysis<\/b><\/span><\/p>\n

VPA — Virtual Population Analysis.<\/b>\u00a0 A version of the virtual population analysis used by fisheries biologists has been incorporated.<\/span><\/p>\n

Estimation of Occupancy Rates<\/b><\/span><\/p>\n

Occupancy Estimation.<\/b>\u00a0 This model provides estimates of the proportion of a set of sites or plots that are occupied by the species of interest when the probability of detection is <1.\u00a0 More details are provided here<\/a> and in MacKenzie et al. (2002)<\/a>.\u00a0 The robust design occupancy model<\/a> is also available in MARK, with 3 parameterizations.\u00a0 Other recent extensions of the occupancy model include the Royle and Nichols (2003<\/a>) model to account for heterogeneity from population size<\/a>, and the multiple-state occupancy model<\/a> of Nichols et al. (2007)<\/a>.\u00a0 Occupancy models for <\/span>2 species<\/a><\/span> are also available.<\/span><\/p>\n

Parameters for all the models are specified in Parameter Index Matrices<\/a>, or PIMs.\u00a0 See Constant Matrix<\/a> for how parameters are specified constant for each occasion, Time Parameter Matrices<\/a> for parameters that are specific to each occasion, Age Matrix<\/a> for parameters that are age-specific, and Time and Age Matrix<\/a> for an example where parameters are both time and age specific.<\/span><\/p>\n

Given a set of parameter matrices, the Design Matrix<\/a> can be used to provide further constraints on the set of estimable parameters.\u00a0 In addition, covariates are specified in the Design Matrix.<\/span><\/p>\n

The model is then “Run” (see Run Window<\/a>) to obtain parameter estimates.<\/span><\/p>\n","protected":false},"excerpt":{"rendered":"

Contents – Index Model Structure Program Mark provides parameter estimates for 142 data types: Cormack-Jolly-Seber models (live animal recaptures that are released alive), band (ring) recovery models (dead animal recoveries), models with both live and dead re-encounters (Burnham’s model), known …<\/p>\n

Model Structure<\/span> Read More »<\/a><\/p>\n","protected":false},"author":117,"featured_media":0,"parent":0,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":{"footnotes":""},"class_list":["post-183","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/sites.warnercnr.colostate.edu\/gwhite\/wp-json\/wp\/v2\/pages\/183","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/sites.warnercnr.colostate.edu\/gwhite\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/sites.warnercnr.colostate.edu\/gwhite\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/sites.warnercnr.colostate.edu\/gwhite\/wp-json\/wp\/v2\/users\/117"}],"replies":[{"embeddable":true,"href":"https:\/\/sites.warnercnr.colostate.edu\/gwhite\/wp-json\/wp\/v2\/comments?post=183"}],"version-history":[{"count":1,"href":"https:\/\/sites.warnercnr.colostate.edu\/gwhite\/wp-json\/wp\/v2\/pages\/183\/revisions"}],"predecessor-version":[{"id":184,"href":"https:\/\/sites.warnercnr.colostate.edu\/gwhite\/wp-json\/wp\/v2\/pages\/183\/revisions\/184"}],"wp:attachment":[{"href":"https:\/\/sites.warnercnr.colostate.edu\/gwhite\/wp-json\/wp\/v2\/media?parent=183"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}