R/fastapprox.R

In bptwin: Cumulative incidence models for twin data

fastapprox <- function(t1,t2,y1) {
  if (is.matrix(t1)) {
    y <- t1[,-1]; t1 <- t1[,1]
  }    
  arglist <- list(name="FastApprox",
                  a=t1,
                  t=y1,
                  z=t2,
                  DUP=FALSE)##,PACKAGE="bptwin")
  res <- do.call(".Call",arglist)
  return(res)
}