ComputeMLE {curstatCI} | R Documentation |
The function ComputeMLE computes the Maximum Likelihood Estimator of the distribution function under current status data.
ComputeMLE(data)
data |
Dataframe with three variables:
|
In the current status model, the variable of interest X with distribution function F is not observed directly.
A censoring variable T is observed instead together with the indicator Δ = (X ≤ T).
ComputeMLE computes the MLE of F based on a sample of size n <- sum(data$freq2)
.
Dataframe with two variables :
jump locations of the MLE
MLE evaluated at the jump locations
Groeneboom, P. and Hendrickx, K. (2017). The nonparametric bootstrap for the current status model. Electronic Journal of Statistics 11(2):3446-3848.
library(Rcpp) library(curstatCI) # sample size n <- 1000 # Uniform data U(0,2) set.seed(2) y <- runif(n,0,2) t <- runif(n,0,2) delta <- as.numeric(y <= t) A<-cbind(t[order(t)], delta[order(t)], rep(1,n)) mle <-ComputeMLE(A) plot(mle$x, mle$mle,type ='s', ylim=c(0,1), main= "",ylab="",xlab="",las=1)