R/mpCart.R

In kStatistics: Unbiased Estimators for Cumulant Products and Faa Di Bruno's Formula

mpCart <-
function (M1=NULL, M2=NULL) {
     if ( is.null(M1) ) stop("The first parameter is missing")
     if ( is.null(M2) ) stop("The second parameter is missing")
     M12<-list();k<-0;
     for (i in 1:length(M1) ) {
          for (j in 1:length(M2) ){
              k<-k+1;M12[[k]]<-list(); 
              M12[[k]][[1]] <- c( M1[[i]][[1]], M2[[j]][[1]] );
              M12[[k]][2]   <- M1[[i]][[2]]*M2[[j]][[2]] ;
    }}
  return(M12);
}