ComBat_seq returns a list rather a matrix if rm=0

[flowerdna]

# Safely assign counts for 'rm'  original: adjust_counts_whole[rm,] <- countsOri[rm,]#this line convert to a list if rm=0
if (length(rm) > 0 && all(rm > 0)) {
  adjust_counts_whole[rm, ] <- countsOri[rm, ]
}
return(adjust_counts_whole)

}