Maximum likelihood estimatation of a weighted probability density function is completed. is done on a weighted distribution. The weighted distribution is a typical probability density distribution multiplied by a weight function. The weight function can be used to truncate the distribution by returning zero beyond some threshold value.
estWD(fatDist, weightFun, distribution, ...)weightedDistribution(fatDist, weightFun, distribution, ...)
Data frame of the parameter estimates with fit statistics.
Vector of fatality distanes from the turbine.
R function that is multipled by the probability distribution, see details.
Character indicating the distribution for weightedDistribution
or vector for estWD
.
Additional arguments passed to weightFun
or optim
.
The function estWD
is a convient wrapper function to weightedDistribution
, for fitting multiple distributions.
The weight function should return a (relative) probability of detection at every distance.
Typically this is the proportion of area searched.
The function weightFun
is set up to take a table of proportion of area searched and return values in a function format.
Let weightFun
), distribution
, fatDist
), and
calcAC
weightFun
## load the data
data(carcassDistance)
data(proportionAreaSearched)
###############################################
## fit for fall carcasses found on road and pad (RP)
distanceFallRP <- subset(carcassDistance,plotType=='RP'&season=='fall',
select=distanceFromTurbine,drop=TRUE)
# \donttest{
fallRPFit <- estWD(fatDist=distanceFallRP,weightFun=weightFun,
distribution=c('norm','gamma','weibull'),propTable=proportionAreaSearched,type='RP',
typeCol='plotType',distanceCol='distanceFromTurbine',propCol='proportionAreaSearched',
maxDistance=100)
# }
###############################################
## fit for fall carcasses found on full plots
distanceFallFP <- subset(carcassDistance,plotType=='FULL'&season=='fall',
select=distanceFromTurbine,drop=TRUE)
# \donttest{
fallFPFit <- estWD(fatDist=distanceFallFP,weightFun=weightFun,
distribution=c('norm','gamma','weibull'),propTable=proportionAreaSearched,type='FULL',
typeCol='plotType',distanceCol='distanceFromTurbine',propCol='proportionAreaSearched',
maxDistance=100)
# }
Run the code above in your browser using DataLab