R/hwe.modelE.sol.R

In HardyWeinberg: Statistical Tests and Graphics for Hardy-Weinberg Equilibrium

hwe.modelE.sol <-
function(x,parinit=c(0.50,0.50,0.00),verbose=FALSE,tracing=1) {
  hwe.out <- solnp(parinit, fun = hwe.modelE.obj, eqfun = NULL, eqB = NULL,  
                   ineqfun = eqnE, ineqLB = c(0,0,0,-Inf,0,-Inf),
                   ineqUB = c(Inf,Inf,Inf,0,Inf,0), 
                   LB = c(0,0,-1), UB = c(1,1,1), 
                   control = list(trace = tracing, tol = 1e-8, inner.iter = 1000), x)
  k <- hwe.out$pars
  H <- hwe.out$hessian
  return(list(k=k,H=H,hwe.out=hwe.out))
}