src/library/stats4/tests/confint.R

## PR#14646

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
set.seed(123)
x <- rnorm(N, 0, 1)

fit <- mle(minusLogL1, start = list(mu=0, logsigma2=0))
confint(fit)

fit2 <- mle(minusLogL1, start = list(mu=0), fixed = list(logsigma2=0))
confint(fit2) # failed

fit3 <- mle(minusLogL2, start = list(mu=0))
confint(fit3) # same