R/theta.max.R

Defines functions `theta.max`

`theta.max` <-
function(x, betas){
  opt <- function(theta) -sum(dbinom(x, 1, plogis(theta - betas), log = TRUE))
  out <- optim(log(sum(x)/(length(x)/sum(x))), opt, method = "BFGS", hessian = TRUE)
  out$par
  }
wasabi1989/MiscPsycho documentation built on Jan. 19, 2020, 12:29 a.m.