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, ...)
an object of class psFit--see Details.
an object of type psData, usually obtained from
readData.
the number of terms to compute the probability distribution for.
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.
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
Run the code above in your browser using DataLab