R/trecase.sex.dom.A.R

In rxSeq: Combined Total and Allele Specific Reads Sequencing Study

`trecase.sex.dom.A` = 
function(yi, ind.lst, X, sex, ni, ni0, xs, start, iphi=1, theta=1, maxiter=100, eps=1E-3, tech.ctrl){
  triali = 0
  par0 = start
  n = 6 + ncol(X)
  repeat{
    triali = triali + 1
    tag = tryCatch({    
      res = trecase.full.A(yi=yi, ind.lst=ind.lst, X=X, twosex=TRUE, sex=sex, ni=ni, ni0=ni0, xs=xs,  
                     start=par0, iphi=iphi, theta=theta, maxiter=maxiter, eps=eps, tech.ctrl=tech.ctrl)
      0
    }, error=function(e) {
      1
    })
    if((tag == 0) | (triali >= tech.ctrl$maxtrial))break;
    par0 = rnorm(length(start), start, 1)
  }
  if(tag == 1){
    ret = NULL
  }else{
    if(!is.null(res)){
      ret = c(res[1:n], 0, res[-(1:n)])
    }
  }
  return(ret)        
}