R/malm11.R

In pwirapati/acdx: Aggregated Cell Differential Expression

malm11 <- function( y, w = NULL, family=0,
  iter_max=30, epsilon=1e-4,u_0=0,v_0=0,verbose=0 )
{
  n <- dim(y)[2]
  m <- dim(y)[3]
  if(is.null(w)) w <- rep(1,n)
  if(length(u_0)==1) u_0 <- rep(u_0,n)
  r <- .C("malm11",
    dim=as.integer(c(n,m)),
    y = as.double(y + c(rbind(u_0,v_0))),
    w=as.double(w),
    E = as.double(rep(1,2*n*m)),
    alpha=as.double(ifelse(w==0,NA,1)),
    beta0=as.double(rep(1,m)),
    psi=as.double(rep(1,m)),
    L=double(1),
    iopt=as.integer(c(family,verbose,iter_max)),
    dopt=as.double(c(epsilon)),
    NAOK=TRUE,
    PACKAGE="acdx"
    )
  r$y <- NULL
  r$E <- NULL
  names(r$alpha) <- dimnames(y)[[2]]
  names(r$beta0) <- dimnames(y)[[3]]
  names(r$psi) <- dimnames(y)[[3]]
  r
}