R/surv.exp.gt.model.R

In survSNP: Power Calculations for SNP Studies with Censored Outcomes

Defines functions surv.exp.gt.model

Documented in surv.exp.gt.model

surv.exp.gt.model <- function(pilm,lm,gtprev,GRR,zmodel,interval)
  {
    ww=GRR^zmodel
    foo=function(lam0) {
        ww=-lam0*ww*lm
        pilm-sum(exp(ww)%*%gtprev)
    }
    ww*uniroot(foo,interval)$root
  }

Try the survSNP package in your browser