fill.db = function(d,tbl,exper,m=3,..., reduce=list(PKC=c(4,8),pmek=6,PKA=9,PIP2=5,pakts473=c(3,7))) { reduce = lapply(reduce,function(x) !(exper%in%x)) var.list = (1:ncol(tbl)) n = if(is.matrix(tbl)) colnames(tbl) else names(tbl) d[0,] = sapply(var.list,function(i) { r = reduce[[n[i]]] BDeu.metric(tbl,i,c(),11,r) }) fast.sync(d) for(k in 1:m) { subapply(var.list,k,function(j) { j = unlist(j) d[j,] = sapply(var.list,function(i) { if(i%in%j) { -Inf } else { r = reduce[[n[i]]] BDeu.metric(tbl,i,j,11,r) } }) }) fast.sync(d) } }