Hi I ran BiSeq yesterday using the tutorial/manual and noticed an error kept occurring when running clusterSites() on page 13 of the manual. I found the source of the error here: groups = colData(rrbs)$group.  This should be groups = factor(colData(rrbs)$group), As the way the groups variable is handled in the filterBySharedRegions.R is it takes the length of levels(groups) but without taking the factor of this first, levels(colData(rrbs)$group)=null which throws a stop error at line 47-49 even if the intended logic holds.

Marcus