R/gen2listBN.R

In monotonicity/stacmr: Conjoint Monotonic Regression approach to State-Trace Analysis

gen2listBN <- function (data) {
  # converts a data frame into a list of lists of binomial data
  # suitable for staCMRBN.R
  #
groups = data[,1]; ugroup = unique(groups); ngroup = length(ugroup)
conds = data[,2]; ucond = unique(conds); ncond = length(ucond)
vars = data[,3]; uvar = unique(vars); nvar = length(uvar)
counts = data[,4:5]
output = vector("list", ngroup)

for (igroup in 1:ngroup) {
  temp = vector("list",nvar)
  for (ivar in 1:nvar) {
    k = which(groups==ugroup[igroup] & vars==uvar[ivar])
    temp[[ivar]] = counts[k,]
    names(temp[[ivar]]) = c("hits", "misses")
    rownames(temp[[ivar]]) = c(1:length(k))
    }
  output[[igroup]] = temp
} 
output$ngroup = ngroup; output$nvar = nvar; output$ncond = ncond
return (output)
}