The function also returns a p-value to evaluate whether partitioning the data is statistically warranted.
partitionCountData(data, n.boots.part = 100, max.pop=500)
fitTirm
for details)
fitTirmPartition
The PART algorithm first fits the Two-Innate Rates Model (TIRM) to the entire data set. Assuming the distribution of capture counts can be approximated as either a two (null) or three (alternative) class multinomial distribution, the algorithm considers all possible ways to partition the count data under both the null and alternative distributions. The partitioning schemes which maximizes the respective likelihoods is retained. A test-statistic, the ratio of multinomial (log) likelihoods (3:2) is obtained for the observed data.
The MLE for number of individuals in class a, number of individuals in class b and the ratio of capture probabilities (alpha) are used to simulate n.boots.part
data sets
under the TIRM. The test-statistic is calculated for each data set as described above.
The observed test-statistic is then compared to the distribution of simulated test-statistics. When the observed test-statistic falls in the tail of the distribution of the test-statisitic from the simulated data, the null model can be rejected.
Stansbury C.R., D.E. Ausband, P. Zager, C.M. Mack, C.R. Miller, M.W. Pennell, and L.P. Waits. in prep. Non-invasive genetic sampling of rendezvous sites and population estimation of grey wolves in Idaho, USA
fitTirmPartition
, fitTirm
## Use dummy data set with a few individuals having very high capture counts
d <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,3,3,4,4,6,6,7,8,10,10,14,17,19,22,22,25)
data <- buildClassTable(d)
## Perform Partitioning Algorithm
part.data <- partitionCountData(data=data, n.boots.part=10, max.pop=500)
part.data
Run the code above in your browser using DataLab