Displaying 2 results from an estimated 2 matches for "logsigma2".
2010 Jul 31
2
Is profile.mle flexible enough?
...much or if there is room for improvement.
* Problem #1 with fixed parameters. I can't manage to get profile-based
confidence intervals when some parameters of the likelihood function are
fixed:
-----------8<--------------------------------
library(stats4)
minusLogL1 <- function(mu, logsigma2) { N*log(2*pi*exp(logsigma2))/2 +
N*(var(x)+(mean(x)-mu)^2)/(2*exp(logsigma2)) }
minusLogL2 <- function(mu) { logsigma2 <- 0;
N*log(2*pi*exp(logsigma2))/2 +
N*(var(x)+(mean(x)-mu)^2)/(2*exp(logsigma2)) }
N <- 100
x <- rnorm(N, 0, 1)
fit <- mle(minusLogL1, start=list(mu=0, logsig...
2007 Aug 14
2
State Space Modelling
Hey all,
I am trying to work under a State Space form, but I didn't get the help
exactly.
Have anyone eles used this functions?
I was used to work with S-PLUS, but I have some codes I need to adpt.
Thanks alot,
Bernardo
[[alternative HTML version deleted]]