R/mle_incomplete_gsem.R

In sdtemple/incompleteSEM: What the Package Does (One Line, Title Case)

#' MLE for Removal Rate
#' 
#' Compute the maximum likelihood estimate for the removal rate.
#' 
#' @param r
#' @param i
#' 
#' @return named vector
#'  
#' @export 
mle_removal_rate <- function(r, i){
  ind <- (!is.na(r)) * (!is.na(i))
  ind <- which(ind == 1, arr.ind=T)
  r <- r[ind]
  i <- i[ind]
  ri <- r - i
  return(c(mle = length(ri) / sum(ri), n = length(ri)))
}