Displaying 4 results from an estimated 4 matches for "absval".
2004 Sep 13
2
Problem with mle in stats4 (R 1.9.1)
...I ran into). If you run the mle example in
stats4
library(stats4)
x <- 0:10
y <- c(26, 17, 13, 12, 20, 5, 9, 8, 5, 4, 8)
ll <- function(ymax=15, xhalf=6)
-sum(stats::dpois(y, lambda=ymax/(1+x/xhalf), log=TRUE))
(fit <- mle(ll))
plot(profile(fit), absVal=FALSE)
everything works fine. Now run
(fit <- mle(ll, method="BFGS", control=list(ndeps=c(1e-3, 1e-3))))
plot(profile(fit), absVal=FALSE)
and you will get
Error in interpSpline.default(obj[[i]]$par.vals[, i], obj[[i]]$z,
na.action = na.omit) :
only 0's may m...
2006 Jun 23
1
How to use mle or similar with integrate?
Hi
I have the following formula (I hope it is clear - if no, I can try to
do better the next time)
h(x, a, b) =
integral(0 to pi/2)
(
(
integral(D/sin(alpha) to Inf)
(
(
f(x, a, b)
)
dx
)
dalpha
)
and I want to do an mle with it.
I know how to use mle() and I also know about integrate(). My problem is
to give the parameter values a and b to the
2004 Jun 10
1
overhaul of mle
...tually fairly easy to implement (I
think: I haven't really tested that it works on anything hard, just
that it doesn't seem to break profiling) -- requires pfit to be
assigned globally within onestep() and a few lines of code further
down.
added an AIC method for mle objects
collapsed the absVal/!absVal code cases slightly
added a "sqrVal" argument for those who want to see the value of the
log-likelihood, not the square root or signed square root (could be
collapsed into a "scale" argument for the profile plot = "sqrt",
"abssqrt", "lik")...
2004 Aug 04
2
fitting distributions
Hello,
I also try to fit a skewed distribution (like skewed student t) to data
points. Do you have an idee howto do this???
thank you
fabrice dusonchet
***********************************************************************************
This email and any files transmitted with it are confidentia...{{dropped}}