R/REXP_joint.R

In fbati: Gene by Environment Interaction and Conditional Gene Tests for Nuclear Families

REXP_joint <- function(  data, marker, trait, env, model, RET_a, RET_b, RET_numInf ) {
  rexport_result <- .C( 'eREXP_joint',  as.double(data), as.integer(dim(data)), as.double(marker), as.integer(length(marker)), as.double(trait), as.double(env), as.double(model), as.double(RET_a), as.integer(length(RET_a)), as.double(RET_b), as.integer(dim(RET_b)), as.double(RET_numInf), as.integer(length(RET_numInf)), NAOK=TRUE )
  list(   a = rexport_result[[8]] ,   b = rexport_result[[10]] ,  numInf = rexport_result[[12]]  )
}

# 2020-07-07 -- changed dup=FALSE to DUP=FALSE, then completely removed, it shouldn't be here anymore...