{"id":265,"date":"2017-04-18T03:35:34","date_gmt":"2017-04-18T03:35:34","guid":{"rendered":"http:\/\/sites.warnercnr.colostate.edu\/gwhite\/?page_id=265"},"modified":"2017-04-18T03:35:34","modified_gmt":"2017-04-18T03:35:34","slug":"robust-design-multi-state-state-uncertainty","status":"publish","type":"page","link":"https:\/\/sites.warnercnr.colostate.edu\/gwhite\/robust-design-multi-state-state-uncertainty\/","title":{"rendered":"Robust Design Multi-State with State Uncertainty"},"content":{"rendered":"

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

\n

Robust Design Multi-State with State Uncertainty<\/b><\/span><\/p>\n

The basic <\/span>robust-design multi-state data types<\/a><\/span>,<\/span> as well as the <\/span>open robust-design multi-state data type<\/a><\/span>, assume that all encounters are correctly classified to their state.<\/span><\/p>\n

Four additional robust design multi-state with state uncertainty data types (2 closed, 2 open) are also available where encounters cannot be always classified into one of the available states.\u00a0 The 2nd open data type includes parameters for when an attribute for determining the state becomes available.\u00a0 This open model is discussed in more detail below.\u00a0 What follows pertains to the 2 closed data types and the open data type that includes the phi parameter.\u00a0 These 3 data types are labeled as:<\/span>
\n119. ORDMSState: Open Robust Design Multi-state with State Probabilities<\/span>
\n\u00a0\u00a0\u00a0 121. RDMSMisClass: Closed Robust Design Multi-state with State Uncertainty<\/span>
\n122. RDMS2MisClass: Closed Robust Design Multi-state with 2 States Uncertain<\/span><\/p>\n

An additional parameter, delta, is included in these models to estimate the probability of correct classification given that the animal is in a specific state.\u00a0 Two of these data types only allow a single uncertain category.\u00a0 For example, Kendall et al. (2004) were interested in the survival of manatees conditional on their reproductive status: cows with calves versus cows without calves.\u00a0 When a calf was observed with a cow, the state is known unambiguously.\u00a0 However, a cow observed without a calf may have a calf, but the calf just isn’t visible, e.g., because of water clarity.\u00a0 On the other hand, in perfectly clear water, maybe the cow can be classified unambiguously as not having a calf.\u00a0 So, 2 states can be defined: C is with calf, and N is without calf.\u00a0 For these data types with a single uncertain category, a u (lower case u) is used to indicate that a cow was observed, but not sure about whether a calf is present or not.\u00a0 So, the encounter history consists of C, N, u, 0, and dots can be included.<\/span><\/p>\n

For three of the data types and using the example above to illustrate, the parameters would be S<\/span>C<\/span>, S<\/span>N<\/span>, psi C to N, psi N to C, pi, omega, p<\/span>C<\/span>, p<\/span>N<\/span>, delta<\/span>C<\/span>, and delta<\/span>N<\/span>.\u00a0 Most of these parameters are obvious from other robust design multi-state data types, but some are new because of the uncertain states. The parameter pi is the proportion of the population released in a specific state, with 1 less pi parameters than their are states.\u00a0 So, pi(1) is the proportion of the population released as state C, pi(2) is obtained by subtraction as 1 – pi(1), the proportion of the population released in state N.\u00a0 Because animals can be released during any primary period, there is a value of pi for all primary periods except the last.\u00a0 In addition, in April, 2013, a change was made to the 3 data types. The mixture parameters for the first primary period, pi1 was reparameterized, so that pi1^s = w1^s*p1^*s \/ sum[w1^s*p1^*s] (see <\/span>Kendall et al. 2012 Ecology<\/a><\/span>). Therefore pi1 no longer exists as a parameter in the likelihood, and there are now K-2 parameters in the pi PIMs, where K is the number of primary periods. The first parameter listed is for primary period 2, and the last pi is for primary period K-1.\u00a0 There is a pi estimate only for the first S – 1 states, where S is the number of states.\u00a0 The pi for the last state is obtained by subtraction. We made this change because for the common case where a given state is never known with certainty, pi1 and therefore the survival and transition probabilities for primary period 1 for that state were not estimable.\u00a0<\/span><\/p>\n

The omega parameter is proportion of the population in each state d at each primary period.\u00a0 So again the omega values have to sum to 1, and the omega for the last state is obtained by subtraction.\u00a0 The parameters deltaC and deltaN are the probabilities of correctly classifying the state of an encountered animal, given the true state of either C or N.\u00a0 So, the probability of encountering a cow with calf and correctly determining that the cow had a calf would be pC<\/span> deltaC, whereas encountering a cow with a calf but not seeing the calf would be p C(1 – deltaC).<\/span><\/p>\n

For the data type where there is only one uncertain state in the encounter history, which must be the lower-case letter u, the above parameters are all that are needed to parametrize the model.\u00a0 For the robust-design multi-state open data type with an uncertain state, 2 additional parameters are required for each state: the\u00a0 of entry to the study area (pent) and the probability of remaining on the study area (phi).\u00a0 Both are as defined for the robust-design multi-state open model without uncertain states.\u00a0 Note that the pent estimates within a primary period must sum to <= 1, so if there are more than 2 secondary occasions within a primary period, you should either use the <\/span>mlogit link<\/a><\/span> function, or else initiate these estimates to a small value that meets the criterion of the sum <= 1.\u00a0 The last pent parameter is obtained by subtraction.\u00a0 So, if you want all the animals to be present on occasion 1, i.e., already in the population when sampling starts, <\/span>fix<\/a><\/span> pent(1) to 1 and the rest of the pent values to zero.<\/span><\/p>\n

There are no additional parameters for the data type that assumes closed primary sessions with more than 1 uncertain state.\u00a0 This data type assumes that you correctly classify an encounter into a state and specify the upper-case state identifier, e.g., C and N.\u00a0 For observations where you are uncertain, the lower-case state identifier is used, e.g., c and n (although the example breaks down here).\u00a0\u00a0<\/span>
\n\u00a0<\/span>
\nAt this time, there is no open data type equivalent for the >1 uncertain states closed data type.<\/span><\/p>\n

Open Data Type with Parameter for when the state becomes Detectable.<\/b><\/span><\/p>\n

<\/b>This data type is labeled in the output as 142. RDMSOpenMCSeas: Robust Design Multi-state Open with State Uncertainty and Seasonal Effects.<\/span>
\nThe parameters between primary sessions are identical to the 3 data types described above: S for each state, appropriate transition parameters (psi), and pi (described above).\u00a0 A <\/span>PIM<\/a><\/span> of pent parameters appear for each primary session and each state.\u00a0 Each set must sum to one, so the last pent value for each primary session obtained by subtraction.\u00a0 The <\/span>mlogit<\/a><\/span> link function should be used for pent.\u00a0 There is also a parameter for remaining on the study area, except that this data type defines this parameter opposite of the open model described above.\u00a0 In this data type, instead of phi, the parameter d is used, with d = 1 – phi.\u00a0 Think of the d parameter as departure (rather than remain as for phi).\u00a0 The estimates of d correspond to 1 – phi, with identical SE.\u00a0 There is another difference between the 2 open data types.\u00a0 The d PIM is not upper-triangular as for phi.\u00a0 That is, the phi parameter can be made age-specific.\u00a0 However, the d parameter is only time-specific.\u00a0 The reason for this difference is because of the number of states that would be possible if both the d and c (defined below) were potentially age-specific.<\/span><\/p>\n

Likewise, there is a set of p PIMs for each primary period for each state.\u00a0 There are also PIMs for the delta parameters (probability that the state is determined, given that the attribute is present to distinguish the state), again for each primary period and each state.<\/span><\/p>\n

The 2 additional parameters for this more complex open model are alpha and c.\u00a0 PIMs for each appear for each primary occasion and each state.\u00a0 Alpha is the probability that the attribute to assign the state has appeared.\u00a0 In the example above for manatees, this might be the birth of the calf.\u00a0 So, the adult manatee is present on the study area, but until she gives birth, she cannot be classified as with calf (C), or no calf (N).\u00a0 The alpha parameters in each PIM must sum to <=1, with the last alpha obtained by subtraction.\u00a0 So, if sampling starts after all pregnant females have given birth, alpha(1) is <\/span>fixed<\/a><\/span> to 1 and the rest <\/span>fixed<\/a><\/span> to zero.\u00a0\u00a0<\/span><\/p>\n

The parameter c is the probability that the attribute allowing the state to be determined still exists.\u00a0 For example, manatee cows might wean their calves towards the end of the sampling period, so that the attribute of a calf present is no long available.\u00a0 The attribute “ceases”, so that the state can no longer be determined.\u00a0 The c parameter has upper-triangular PIMs for each primary period and each state.\u00a0 These PIMs can be time-specific or age-specific if the attribute disappears as a function of the length of time it has been present.<\/span><\/p>\n

To make the 2 open models exactly equivalent, fix alpha(1) = 1 for each primary session and each state, and the rest of the alpha parameters fixed to zero.\u00a0 Also fix all the c parameters to zero, because the attribute can never disappear.\u00a0 The phi PIMs must be time-specific (not age-specific) to match the time-specific d parameters of the more complex open model.<\/span><\/p>\n

Residence and duration times for this model are descibed in <\/span>Duration Times<\/a><\/span>.<\/span><\/p>\n","protected":false},"excerpt":{"rendered":"

Contents – Index Robust Design Multi-State with State Uncertainty The basic robust-design multi-state data types, as well as the open robust-design multi-state data type, assume that all encounters are correctly classified to their state. Four additional robust design multi-state with …<\/p>\n

Robust Design Multi-State with State Uncertainty<\/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-265","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/sites.warnercnr.colostate.edu\/gwhite\/wp-json\/wp\/v2\/pages\/265","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=265"}],"version-history":[{"count":1,"href":"https:\/\/sites.warnercnr.colostate.edu\/gwhite\/wp-json\/wp\/v2\/pages\/265\/revisions"}],"predecessor-version":[{"id":266,"href":"https:\/\/sites.warnercnr.colostate.edu\/gwhite\/wp-json\/wp\/v2\/pages\/265\/revisions\/266"}],"wp:attachment":[{"href":"https:\/\/sites.warnercnr.colostate.edu\/gwhite\/wp-json\/wp\/v2\/media?parent=265"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}