Changes in version 1.17 o BSseq() will no longer reorder inputs. Previously, the returned _BSseq_ object was ordered by ordering the loci, although this behaviour was not documented. BSseq() may still filter out loci if rmZeroCov = FALSE or collapse loci if strandCollapse = FALSE or duplicate loci are detected, but the relative order of loci in the output will match that of the input. o Fix bug with maxGap argument of BSmooth(). The bug meant that the 'maximum gap between two methylation loci' was incorrectly set to 2 * maxGap + 1 instead of maxGap. This likely did not affect results for users who left the default value of maxGap = 10^8 but may have affected results for small values of maxGap. o Cleaning up Imports and Suggests. Changes in version 1.13 o 1.13.6: Fix performance regression in BSmooth(). Thanks to Shan Andrews for the report (). o 1.13.5: Fix major bug in combine() and combineList(). *This bug led to bad BSseq objects with incorrect methylation estimates due to incorrect 'M', 'Cov', 'coef', and 'se.coef' assays.* To be safe, BSseq objects created with versions 1.13.0 to 1.13.4 should be re-created using a newer version. More specifically, any BSseq objects created with combine() or combineList() should be re-created. Also, BSseq objects created using read.bismark() or read.bsmooth() with multiple 'files' should be re-created. Thanks to Alejandro Reyes (@areyesq89) for the report (). o 1.13.4: Fix performance regression in getMeth() and getCoverage() when 'regions' were supplied. Thanks to Alejandro Reyes (@areyesq89) for the report (). o Moved vignettes to Rmd. Changes in version 1.11 o bsseq now uses DelayedMatrix objects from the DelayedArray package for all matrix-like data. This enables large data to be stored on disk rather than in memory. o Serialized (saved) BSseq, BSseqTstat, and BSseqStat objects will need to be updated by invoking x <- updateObject(x). Changes in version 1.7 o Fixing an error with reading multiple samples aligned with bismark in the format "cytosineReport". Changes in version 1.5 o new function strandCollapse for collapsing forward and reverse strand data to be unstranded. o Updated read.bismark() to support the cytosine report files; both formats are supported. Other minor updates (mostly internal) to read.bismark(). Greatly improved documentation of this function, paying particular attention to differences in file formats between versions of Bismark. Changes in version 0.11 o Converted to using Authors@R in the DESCRIPTION file. o plotRegion did not respect the col/lty/lwd argument if given explicitely as opposed to through pData(). Reported by Karen Conneely . o Fixed an issue introduced by the previous change (to plotRegion). Reported (with a fix) by Tim Triche Jr . o Fixed a serious bug in collapseBSseq reported by Julien Roux : it would use the Meth values instead of Cov for the final object's Cov. However, this will result in the return object having a methylation of 100 percent across all loci, so hopefully users will have seen this for themselves. o Fixed a bug in combineList which made combineList use "slow" code even in a special case where a faster shortcut was possible. This bug fix does not change the output of the function under any input, it only makes it faster. Reported by Julien Roux . o validity now checks for the presence of colnames (sampleNames) which was assumed to be set. Reported by Kevin Rue-Albrecht . o Fixed a man page issue. o Slightly changed CITATION. o Added GitHub URL to DESCRIPTION. Changes in version 0.9 o Fixed a problem with "width" in the title of bsseq plots. o plot.BSseqTstat now allows for BSseqTstat objects computed without correction. o validObject(BSseq) has been extended to also check for sampleNames consistency. o Fixed a bug related to validity checking. o Increased maxk from 10,000 to 50,000 in calls to locfit, to allowing fitting the model on genomes with unusally long chromosomes (Thanks to Brian Herb for reporting). o The class representation for class 'BSseq' has changed completely. The new class is build on 'SummarizedExperiment' from GenomicRanges instead of 'hasGRanges'. Use 'x <- updateObject(x)' to update serialized (saved) objects. o Fixing a problem in orderBSseq related to chromosome names. o Allowed user specification of maxk, with a default of 10,000 in BSmooth. o Many bugfixes made necessary by the new class representation. o Better argument checking in BSmooth.tstat. o A few undocumented functions are now documented. o Rewrote orderBSseq Changes in version 0.7 o Removed the returnRaw argument to read.bismark() as it was unnecessary (Bismark output files does not have additional information beyond M and Cov and genomic positions, unlike BSmooth). o Moved the Bismark example data from data to inst/extdata. o combineList() now deals with the case where the list of BSseq objects have different genomic locations. This speeds up read.bismark() substantially. o Exposed combineList() as a faster alternative to Reduce(combine, list). o Updated the code for the plotting routines (plotRegion). This should not have an impact on user-visible code. o Added read.bismark() function to parse output from the Bismark alignment suit [thanks to Pete Hickey]. o Refactorized plotting code. Changes in version 0.6 o Fixed a bug in getMeth, where type="raw" resulted in an error for non-smoothed data objects. o Updated CITATION and citations in the vignettes. o Now read.bsmooth supports both gzipped and non-gzipped files, whereas previously it assumed the output files to be gzipped. Thanks to Andreas Schoenegger for reporting this problem. o Fixed a bug in combine() that also resulted in a bug in read.bsmooth when multiple input directories were specified. Bug reported by Andreas Schoenegger. Changes in version 0.4 o Improved combine and fixed a bug. Also added a non-exported combineList for testing. o Bug fix to read.bsmooth; it now works correctly for the default settings (= returning a single object of class BSseq and not a list). o Getting ready for initial release on Bioconductor. o Updated the citations in the vignette(s) and the CITATION file.