Displaying 1 result from an estimated 1 matches for "sumlogob".
Did you mean:
sumlogobs
2006 Jun 02
2
Problem with mle
....23860,
34.98480, 43.73100, 52.47720, 61.22340, 69.96960, 71.71880, 73.46810,
75.21730, 76.96650, 78.71580, 80.46500 )
#Sum of log of values in vector
SumLog <- function (x)
{
sum(log(x))
}
#-ln(L) with Poisson Likelihood estimator
NeglogPoisL <- function (obs, est)
{
sel <- est != 0
SumLogObs <- SumLog(obs[sel])
-sum( obs[sel] * log(est[sel]) - est[sel]- SumLogObs )
}
#-ln(L) for Negative Exponential
NENeglogPoisL <- function(a=0.2, rate=0.04)
{
NeglogPoisL( SumSeeds, a * dexp( X, rate) )
}
#-ln(L) for Weibull
WBNeglogPoisL <- function(a=1000000, shape=0.12, scal...