fitlogDist | R Documentation |
This function uses maximum likelihood estimation (MLE) to estimate the shape parameter of a logarithmic distribution from a set of observed counts for either the number of groups/sources of forensically interesting material (mostly glass or paint) recovered from clothing, or the number of fragments/particles in each group. This, in turn, allows the estimation of the P and S probabilities, as described by Evett and Buckleton (1990), which used in computing the likelihood ratio (LR) for activity level propositions. The data itself arises from clothing surveys. The logarithmic distribution has probability mass function
p(k) = \frac{\pi^k}{k\log_e(1-\pi)},0<\pi<1.
fitlogDist(x, nterms = 10, start = 0.5, ...)
fitLogdist(x, nterms = 10, start = 0.5, ...)
fitlogdist(x, nterms = 10, start = 0.5, ...)
x |
an object of type |
nterms |
the number of terms to compute the probability distribution for. |
start |
a starting value for the optimiser. |
... |
other parameters - not currently used. |
The function returns an object of class psFit
which is a
list
contains seven elements:
psData
– an object of class psData
–see readData
,
fit
– the fitted object from optim
,
pi
- the maximum likelihood estimate of the shape parameter,
var
– the estimated variance for the shape parameter,
fitted
– a named vector
containing the first nterms of
the fitted distribution.
an object of class psFit
–see Details.
Evett, I. W. and Buckleton, J. S., "The interpretation of glass evidence. A practical approach", Journal of the Forensic Science Society 1990: 30(4): 215–223.
plot.psFit
, print.psFit
,
probfun
.
data(Psurveys)
roux = Psurveys$roux
fit = fitlogDist(roux)
fit
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.