R/init_theta.R

In helske/GLMMF: What the package does (short line)

init_theta<-function (y, u, distribution) 
{
  if(distribution == "gaussian")
    return(y)
  if (distribution == "poisson") {
    y <- y/u
    y[y < 0.1 | is.na(y)] <- 0.1
    return(log(y))
  }
  if (distribution == "binomial") {
    return(qlogis((ifelse(is.na(y), 0.5, y) + 0.5)/(u + 1)))
  }
  if (distribution == "gamma") {
    y[is.na(y) | y < 1] <- 1
    return(log(y))
  }
  if (distribution == "negative binomial") {    
    y[is.na(y) | y < 1/6] <- 1/6
    return(log(y))
  }
  
}