BUG FIXES
- Noticed an issue with
GenomicState::gencode_genomic_state()
that ultimately was due tomakeGenomicState()
and the transition in R to have as defaultdata.frame(stringsAsFactors = FALSE)
instead ofTRUE
.
SIGNIFICANT USER-VISIBLE CHANGES
- Now
makeGenomicState()
now restores theGenomicFeatures::isActiveSeq()
on thetxdb
object before finishing to avoid issues, like runningregionReport::renderReport()
on two different sets of regions (different chrs). There's a new unit test for this.
BUG FIXES
- Removed knitrBootstrap citation code.
- Made a test less susceptible to numerical precision.
SIGNIFICANT USER-VISIBLE CHANGES
- Documentation website is now available at http://lcolladotor.github.io/derfinder/. It gets updated with every commit on the master branch (bioc-devel) using GitHub Actions and pkgdown.
SIGNIFICANT USER-VISIBLE CHANGES
- Use GenomeInfoDb::getChromInfoFromUCSC() when possible instead of data from biovizBase::hg19Ideogram for getting the hg19 chromosome lengths.
BUG FIXES
- Fixed an important bug in
findRegions()
that affected the end positions of the regions whenmaxRegionGap
was supplied with a value greater than the default of0
and the data was filtered (soposition
was not allTRUE
in thefindRegions()
call). To check this scenario now there is a new unit test undertests/testthat/test-maxRegionGap.R
. - The above bug went unnoticed in
getRegionCoverage()
and thus forregionMatrix()
for the same type of situations (filtered data with a non-zeromaxRegionGap
). The coverage values are ok, it's just the end positions of the regions returned byfindRegions()
that were incorrect and that would need to be re-computed with the fixed version. - Changed some internal tests to check
bumphunter::loessByCluster()
instead ofbumphunter::runmedByCluster()
given some issues with the second one.
- Added a
NEWS.md
file to track changes to the package.
BUG FIXES
railMatrix()
andloadCoverage()
helper functions had an issue when the input set of regions was duplicated. This could be reproduced with
sampleFile <- c('SRR387777' = 'http://duffel.rail.bio/recount/SRP009615/bw/SRR387777.bw')
regs <- GenomicRanges::GRanges('chrY', IRanges(start = c(1, 1), width = 10), strand = '-')
names(regs) <- c(1:2)
result <- rtracklayer::import(sampleFile, selection = regs, as = 'RleList')
This error affected recount and other reverse dependencies that use derfinder for processing BigWig files.
BUG FIXES
railMatrix()
andloadCoverage()
helper functions now attempt to import a BigWig file 3 times before giving up. Based on http://bioconductor.org/developers/how-to/web-query/ and https://github.com/leekgroup/recount/commit/8da982b309e2d19638166f263057d9f85bb64e3f which will make these functions more robust to occasional web access issues.
NEW FEATURES
- Add ORCID's following changes at http://bioconductor.org/developers/package-guidelines/#description
BUG FIXES
- Use R's random seeds from version 3.5.0 for the test thanks to https://twitter.com/StrictlyStat/status/1103303028751372289
SIGNIFICANT USER-VISIBLE CHANGES
- Switched from
outfile
tolog
when invokingBiocParallel::SnowParam()
. Thusdefine_cluster()
now has amc.log
argument instead ofmc.outfile
.
BUG FIXES
- Fix a message regarding the deprecated
IRanges
subset
method.
SIGNIFICANT USER-VISIBLE CHANGES
- Use
BiocManager
BUG FIXES
- Fixed a unit test that was breaking version 1.13.7.
NEW FEATURES
- Added an extra example to
regionMatrix()
in response to https://support.bioconductor.org/p/103591
BUG FIXES
- Improved the documentation regarding an error when coverageInfo$position
is NULL when running
analyzeChr()
and indirectly runningpreprocessCoverage()
. See https://support.bioconductor.org/p/99400/ for details.
SIGNIFICANT USER-VISIBLE CHANGES
- Vignette now uses the new
BiocStyle::html_document
that was recently released.
BUG FIXES
regionMatrix()
will now pass the hidden argumentsspecies
andcurrentStyle
togetRegionCoverage()
so they can be used byextendedMapSeqlevels()
. Related to https://support.bioconductor.org/p/95721/.
BUG FIXES
- Improved documentation of
extendedMapSeqlevels()
. Related to https://support.bioconductor.org/p/95521/. - Improved
filterData()
based on #38
BUG FIXES
- Fixed
define_cluster()
to match recent changes in BiocParallel and fixed an if clause inregionMatrix()
that could lead to warnings in some situations.
SIGNIFICANT USER-VISIBLE CHANGES
regionMatrix()
now has explicit argumentstotalMapped
andtargetSize
so that users will almost always normalize by library size when using this function (if they see the help page) or in the steps prior to usingregionMatrix()
.
BUG FIXES
- Clarified the documentation of
mc.cores
andmc.cores.load
infullCoverage()
thanks to feedback from Emily E Burke https://github.com/emilyburke.
SIGNIFICANT USER-VISIBLE CHANGES
- Help pages now document advanced arguments.
- Deprecated
advancedArg()
.
NEW FEATURES
- Added the function
getTotalMapped()
for calculating the total number of mapped reads for a BAM file or the area under the curve (AUC) for a BigWig file. This information can then be used with fullCoverage(),filterData()
and other functions. Note that if youtotalMapped
infullCoverage()
you should not usetotalMapped
again infilterData()
.
BUG FIXES
- Updated links to BrainSpan. Issue reported by Steve Semick https://github.com/SteveSemick.
BUG FIXES
- Now derfinder uses
DataFrame(check.names = FALSE)
to avoid naming issues.
SIGNIFICANT USER-VISIBLE CHANGES
- Dropped defunct functions.
BUG FIXES
annotateRegions()
now ignores strand by default.
SIGNIFICANT USER-VISIBLE CHANGES
- The users guide vignette now has a short section explaining how to use derfinder for differential binding analysis with ChIP-seq data.
SIGNIFICANT USER-VISIBLE CHANGES
- Added smoothing arguments for the single base-level approach based on functions from the bumphunter package. These arguments are useful for identifying differentially bounded ChIP-seq peaks.
BUG FIXES
- Fixed
railMatrix()
's flexibility for defining the cluster used for loading the BigWig files. You can now useBPPARAM.railChr
which will take priority overfile.cores
. Also, iffile.cores = 1L
, then the default will be to useSerialParam()
, which was the implementation available prior to 1.5.11.
SIGNIFICANT USER-VISIBLE CHANGES
- Now
coverageToExon()
,regionMatrix()
andrailMatrix()
can take anL
argument of length equal to the number of samples in case not all samples have the same read length. railMatrix()
has a new argument calledfile.cores
for controlling how many cores are used for loading the BigWig files. In theory this allows usingrailMatrix()
withBPPARAM.custom
equal to aBiocParallel::BatchJobsParam()
to submit 1 job per chromosome, thenfile.cores
determines the number of cores for reading the files. This is a highly experimental feature.
SIGNIFICANT USER-VISIBLE CHANGES
- Dropped the old introductory and advanced vignettes. We think that the new vignettes are clearer. In particular, they do a better job at highlighting the differences between the expressed regions-level and single base-level F-statistics implementations of the DER Finder approach to RNA-seq data.
SIGNIFICANT USER-VISIBLE CHANGES
- Added a users guide vignette which explains nearly every detail you would want to know as a user.
SIGNIFICANT USER-VISIBLE CHANGES
- Added a quick start vignette.
NEW FEATURES
- Introduced
railMatrix()
which generates similar output toregionMatrix()
but is much faster and less memory intensive. It achieves this by extracting the required information from BigWig files.
SIGNIFICANT USER-VISIBLE CHANGES
- Brought back the
mc.outfile
argument for specifying theoutfile
argument inSnowParam()
. See more details at https://stat.ethz.ch/pipermail/bioc-devel/2015-May/007531.html
SIGNIFICANT USER-VISIBLE CHANGES
- Deprecated functions with underscores in their names in favor of camelCase functions. This was done to simplify the package.
SIGNIFICANT USER-VISIBLE CHANGES
- Greatly increased the speed of the p-values calculation step. See #29 for details.
BUG FIXES
- Updated to work with qvalue 1.99.0.
SIGNIFICANT USER-VISIBLE CHANGES
- Changed citation information to reference the bioRxiv pre-print.
BUG FIXES
- Polished the interaction with bumphunter >= 1.7.3.
- Updated the default cluster option now that
BiocParallel::SnowParam()
no longer has anoutfile
argument.
SIGNIFICANT USER-VISIBLE CHANGES
analyzeChr()
now usesannotateTranscripts()
andmatchGenes()
from bumphunter version 1.7.3 (or greater). As announced at https://support.bioconductor.org/p/63568/ these changes in bumphunter allow straight forward use of non-human annotation. InanalyzeChr()
using a different organism can be used by changing thetxdb
argument: finer control can be achieved through...
. For example, by specifying theannotationPackage
argument used inannotateTranscripts()
.
BUG FIXES
makeGenomicState()
incorrectly labeled regions as intragenic. The correct name is intergenic.
BUG FIXES
- Fixed an important bug on
calculatePvalues()
! Basically, internallymaxRegionGap
was set to 300 instead of 0 in one step by default. Thus the process of mapping regions to genomic coordinates was messed up. If you have results prior to this fix you can try using https://gist.github.com/bf85e2c7d5d1f8197707 to fix the results as much as possible. Basically, regions will be correct but the p-values will be approximated with the available information from the null regions. Truly fixing the p-values can only be done by re-running derfinder.
NEW FEATURES
- Introduced function
extendedMapSeqlevels()
for usingGenomeInfoDb
when there is information regarding the species and naming style of interest. otherwise sequence names are left unchanged. If used withverbose = TRUE
, a message is printed wheneverGenomeInfoDb
could not be used or if some information had to be guessed.
BUG FIXES
NEW FEATURES
loadCoverage()
andfullCoverage()
now supportBamFile
andBigWigFile
objects.
BUG FIXES
- Fixed a bug in
loadCoverage()
when the input was aBamFileList
. Implemented tests based on the bug. Bug reported at https://support.bioconductor.org/p/62073
NEW FEATURES
- Preparing to submit to Bioconductor.
NEW FEATURES
- Added an advanced vignette.
NEW FEATURES
- Introductory vignette completed.
SIGNIFICANT USER-VISIBLE CHANGES
mergeResults()
can now calculate FWER adjusted p-values when provided withoptionsStats
. UpdatedanalyzeChr()
to supply the required information.
NEW FEATURES
- Added an introductory vignette.
NEW FEATURES
- Added
advancedArg()
and its aliasadvanced_arg()
which links to the docs for the advanced arguments by opening a browser window with the relevant information from GitHub. getRegionCoverage()
andcoverageToExon()
now have thefiles
argument which is used only whenfullCov
isNULL
. Both functions will attempt to extract the coverage data from the raw files for the regions of interest in that case.
Special care has to be taken in order to guarantee that the coverage is
the same as some reads might be discarded if the region is too narrow.
See the advanced argument protectWhich
in loadCoverage()
for more
information. Also, if totalMapped
and targetSize
were used prior to
filtering, they should be used again.
loadCoverage()
has new advanced arguments that help when reading a specific region (or regions) of the genome.
SIGNIFICANT USER-VISIBLE CHANGES
loadCoverage()
andfullCoverage()
argumentdirs
has been renamed tofiles
for greater consistency with what it represents.
SIGNIFICANT USER-VISIBLE CHANGES
regionMatrix()
now returns the output ofgetRegionCoverage()
so you don't have to run it twice if you are interested in usingderfinderPlot::plotRegionCoverage()
.regionMatrix()$regions
now guesses the seqlengths
BUG FIXES
- Fixed
.advanced_argument()
to work in nested functions - Fixed a case in
getRegionCoverage()
wherefullCov$position
was provided but it wasNULL
. - Fixed
regionMatrix(totalMapped, targetSize)
case which would previously lead to an error in thegetRegionCoverage()
step.
SIGNIFICANT USER-VISIBLE CHANGES
- Most functions had their arguments changed to increase usability. Some have advanced arguments inside the code.
SIGNIFICANT USER-VISIBLE CHANGES
- Introduced
coerceGR()
,createBwSample()
andcreateBw()
for exporting output fromfullCoverage()
into BigWig files.
SIGNIFICANT USER-VISIBLE CHANGES
- All exported functions now have aliases with underscore names for those
that prefer them over camel case names. For example,
analyze_chr()
is the new alias foranalyzeChr()
. makeBamList()
has been renamed torawFiles()
since it can be used to identify a list of BigWig files instead of BAM files.
SIGNIFICANT USER-VISIBLE CHANGES
loadCoverage()
andfullCoverage()
now have atilewidth
argument. When specifiedGenomicFiles
is used to read the coverage in chunks. In theory, this can lead to lower memory usage at the expense of time.
SIGNIFICANT USER-VISIBLE CHANGES
preprocessCoverage()
now has atoMatrix
argument which is only used whenlowMemDir
is notNULL
. It controls whether to save the chunks asDataFrame
objects ordgCMatrix
objects and the idea is that it can time by just transforming the data once instead of doing so at each permutation.
SIGNIFICANT USER-VISIBLE CHANGES
fstats.apply()
has been moved to it's own package:derfinderHelper
. This will speed up the run time when usingBiocParallel::SnowParam()
asderfinderHelper
takes much less time to load than derfinder.plotCluster()
,plotOverview()
andplotRegionCoverage()
were all moved to their new own package:derfinderPlot
. This will make maintenance easier as the dependencyggbio
is still under active development.
SIGNIFICANT USER-VISIBLE CHANGES
- Re-organized code for
fstats.apply()
. Note that improving.transformSparseMatrix()
would speed up theMatrix
method. - Note that all parallel functions have some overhead from loading derfinder
on each worker. Check
system.time(library(derfinder))
to see how long the overhead is. It only pays off to use more cores if the calculations are taking longer than the overhead.
SIGNIFICANT USER-VISIBLE CHANGES
fullCoverage()
has several new arguments and now is a full parallel implementation ofloadCoverage()
. These changes were introduced sincefullCoverage()
no longer blows up in memory since version 0.0.62 and thus the new recommended use case is to callfullCoverage()
instead of running one job withloadCoverage()
per chromosome or using alapply()
loop.
NEW FEATURES
fstats.apply()
now hasmethod
andscalefac
arguments. Themethod
argument controls which of the 3 implementations to use. The old method is calledregular
now. The new methodRle
calculates the F-statistics without de-compressing the data, which is good for memory but gets considerably slower as the number of samples increases. The default method isMatrix
which uses the Matrix package and is both faster (given that the coercion doesn't take long) and less memory intensive than theregular
method.
SIGNIFICANT USER-VISIBLE CHANGES
- Functions
analyzeChr()
,calculatePvalues()
andcalculateStats()
now have argumentsmethod
andscalefac
to match the changes infstats.apply()
.
NEW FEATURES
-
derfinder now uses
BiocParallel::blapply()
instead ofparallel::mclapply()
Whenmc.cores
is greater than 1,BiocParallel::SnowParam()
is used to construct the cluster. Otherwise,BiocParallel::SerialParam()
is used. This change reduces memory load when using the functions that have themc.cores
argument greater than 1. -
Functions
analyzeChr()
,calculatePvalues()
,calculateStats()
,coverageToExon()
,fullCoverage()
,getRegionCoverage()
,regionMatrix()
all have a new argumentmc.output
. This is passed toBiocParallel::SnowParam(outfile)
.
SIGNIFICANT USER-VISIBLE CHANGES
- You may now use
fullCoverage()
without problems and should no longer encounter errors due to longer vectors not being implemented. - Functions like
fullCoverage()
now use much less memory and do not blow up as you increasemc.cores
. Note however that the memory does increase, but now it`s close to linear. - Examples might take longer to run with
mc.cores
greater than 1, but that is due to the small setup overhead ofBiocParallel::SnowParam()
which is minimal compared to the overall speed gains with real data sets.
SIGNIFICANT USER-VISIBLE CHANGES
filterData()
andloadCoverage()
now have argumentstotalMapped
andtargetSize
getRegionCoverage()
andregionMatrix()
can now work with list output fromloadCoverage()
with a non-NULL cutoffregionMatrix()
now has an argumentrunFilter
so it can be used with previous output fromloadCoverage()
/filterData()
withreturnMean=TRUE
SIGNIFICANT USER-VISIBLE CHANGES
analyzeChr()
,annotateRegions()
,calculatePvalues()
,coverageToExon()
,findRegions()
,fullCoverage()
,getRegionCoverage()
,makeGenomicState()
,mergeResults()
,plotCluster()
,plotOverview()
now all havechrsStyle
as an argument to specify the chromosome naming convention used. Defaults to UCSC.makeGenomicState()
no longer has theaddChrPrefix
argument. It has been replaced bychrsStyle
to useGenomeInfoDb
to set the naming style.chrnums
has been renamed tochrs
infullCoverage()
andmergeResults()
chrnum
has been renamed tochr
inanalyzeChr()
NEW FEATURES
- Updated to work with BioC version 3.0
NEW FEATURES
loadCoverage()
andfullCoverage()
can now import data from BigWig files.
SIGNIFICANT USER-VISIBLE CHANGES
regionMatrix()
now relies ongetRegionCoverage()
instead ofcoverageToExon()
making it faster and less memory intensive.
NEW FEATURES
- Added
regionMatrix()
for filtering coverage data and using the resulting regions to construct a count matrix. Uses several derfinder functions.
SIGNIFICANT USER-VISIBLE CHANGES
- Made
coverageToExon()
more robust for different names infullCov
. filterData()
andloadCoverage()
have new argumentsfilter
,returnMean
, andreturnCoverage
which allow speeding upregionMatrix()
.preprocessCoverage()
was changed accordingly.getRegionCoverage()
now internally uses USCS names.
BUG FIXES
- Fixed warnings in
coverageToExon()
.
SIGNIFICANT USER-VISIBLE CHANGES
- Added a
NEWS
file with curated information from the git commit history.
SIGNIFICANT USER-VISIBLE CHANGES
- Added example data for
mergeResults()
. Now all functions have examples.
SIGNIFICANT USER-VISIBLE CHANGES
- Now derfinder includes the genomic state output for hg19 chr21. This
allowed to implement examples for
annotateRegions()
,getRegionCoverage()
, andcoverageToExon()
.
SIGNIFICANT USER-VISIBLE CHANGES
- Several examples run much faster now
SIGNIFICANT USER-VISIBLE CHANGES
makeGenomicState()
now has an example
BUG FIXES
- Fixed warnings in
makeGenomicState()
related to changes inAnnotationDbi
BUG FIXES
- Updated the help for
plotRegionCoverage()
and fixed an issue.
SIGNIFICANT USER-VISIBLE CHANGES
loadCoverage()
now allows specifying which strand you want to load. More at #16
SIGNIFICANT USER-VISIBLE CHANGES
getRegionCoverage()
now has a depth-adjustment argument
BUG FIXES
- Fixed several bugs as suggested by Andrew Jaffe
SIGNIFICANT USER-VISIBLE CHANGES
- Now requires R version 3.1
- Changed Github organization to match the Git-SVN Bioconductor bridge guidelines.
- Dropped
Rcpp
andRcppArmadillo
from F-stats calculation. More at #17 - Stored the results from some of the example data to speed up other
examples. Check
?genomeDataRaw
,?genomeFstats
,?genomeRegions
BUG FIXES
- Fixed
verbose
forgetRegionCoverage()
plotRegionCoverage()
now matches latestgetRegionCoverage()
output
SIGNIFICANT USER-VISIBLE CHANGES
- Updated
getRegionCoverage()
with a new method for sub setting the coverage matrices, allowing for coverage estimates from overlapping regions. Now also usesmclapply()
.
BUG FIXES
- Fixed
NAMESPACE
to match current bioc-devel (2.14) as suggested by Tim Triche.
SIGNIFICANT USER-VISIBLE CHANGES
- A series of performance enhancements were made to reduce the memory load (albeit a very minor time increase).
BUG FIXES
- Updated
analyzeChr()
to handle correctly the newlowMemDir
argument.
SIGNIFICANT USER-VISIBLE CHANGES
- Introduced the
lowMemDir
argument topreprocessCoverage()
,calculateStats()
,calculatePvalues()
,fstats.apply()
, andanalyzeChr()
. Reduces peak memory usage at the expense of some input-output.
SIGNIFICANT USER-VISIBLE CHANGES
mergeResults()
will not merge pre-processed data by defaultcoverageToExon()
now usesmclapply()
when possible
SIGNIFICANT USER-VISIBLE CHANGES
- Changed how the data is pre-processed and pre-splitted.
SIGNIFICANT USER-VISIBLE CHANGES
preprocessCoverage()
now usesReduce()
instead of.rowMeans()
SIGNIFICANT USER-VISIBLE CHANGES
- Updated the example in
plotCluster()
NEW FEATURES
- Added
collapseFullCoverage()
SIGNIFICANT USER-VISIBLE CHANGES
sampleDepth()
has been greatly changed. It is now based on Hector Corrada`s ideas implemented in metagenomeSeq.- Updated several man pages.
BUG FIXES
- Fixed an issue with the bumphunter dependency.
BUG FIXES
- Merged changes suggested by Michael Love
BUG FIXES
- Implemented fixes suggested by Michael Love
BUG FIXES
- Fixed an important bug in the F-stats calculation
- Implemented fixes suggested by Michael Love
BUG FIXES
loadCoverage()
now usesreadGAlignmentsFromBam()
SIGNIFICANT USER-VISIBLE CHANGES
- Added a
bai
argument tofullCoverage()
loadCoverage()
can now work with a pre-defined BamFileList object.
SIGNIFICANT USER-VISIBLE CHANGES
- Now requires R version 3.0.2
- Changed the default of
center
insampleDepth()
toFALSE
. - Added a
runAnnotation
argument toanalyzeChr()
. - Added a
bai
argument toloadCoverage()
. - Added an
adjustF
argument to all stats functions. Useful for cases when the RSS of the alternative model is very small.
BUG FIXES
- Fixed
plotRegionCoverage()
andplotCluster()
for unexpected cases.
NEW FEATURES
- Added
sampleDepth()
SIGNIFICANT USER-VISIBLE CHANGES
generateReport()
has been moved to it's own new package calledderfinderReport
. It is available at https://github.com/lcolladotor/derfinderReport- Examples and
analyzeChr()
have been updated now thatsampleDepth()
was added
SIGNIFICANT USER-VISIBLE CHANGES
- Renamed package from
derfinder2
toderfinder
to comply with Bioconductor guidelines.
SIGNIFICANT USER-VISIBLE CHANGES
- Changed how
makeModels()
deals with cases when mod and mod0 are not full rank. plotCluster()
now no longer depends on an active Internet connection forhg19 = TRUE
.
BUG FIXES
- Fixed minor graphical issues in
plotRegionCoverage()
,plotCluster()
and ingenerateReport()
BUG FIXES
- Fixed
plotRegionCoverage()
for cases when inannotateRegions(minoverlap=x)
lead to no overlaps being found between a region and annotation.
BUG FIXES
- Fixed some bugs in
calculatePvalues()
when no null regions or only some were found. - Fixed a bug using
colsubset
onanalyzeChr()
. - Fixed an issue when the F-stat cutoff used is too high and no regions are found.
- Fixed an issue when
testvars
inmakeModels()
had unused levels. - Fixed an issue when
qvalue::qvalue()
fails due to incorrect estimation ofpi0
- Fixed a bug on how the regions were being clustered.
SIGNIFICANT USER-VISIBLE CHANGES
- MA-style plots in
generateReport()
now weight the mean by the number of samples in each group. Also removed the mean coverage vs area section.generateReport()
now also has anBestClusters
argument. plotCluster()
now uses scales and has aforceLarge
argument
BUG FIXES
- Fixed an important bug in finding candidate DERs. The example
for
plotCluster()
now includes code that was used for visualizing the bug. - Fixed some NAMESPACE issues
NEW FEATURES
- Added
getRegionCoverage()
,coverageToExon()
,plotRegionCoverage()
SIGNIFICANT USER-VISIBLE CHANGES
- Renamed
plotRegion()
toplotCluster()
plus it no longer shows the exons track as it is redundant information mergeResults()
now also runsannotateRegions()
generateReport()
now usesplotRegionCoverage()
and includes MA-style plots
NEW FEATURES
- Added a
nnotateRegions()
NEW FEATURES
- Added
makeGenomicState()
NEW FEATURES
- Added
fullCoverage()
SIGNIFICANT USER-VISIBLE CHANGES
preprocessCoverage()
now uses thegroupInfo
argumentcalculatePvalues()
now calculates log2 fold changes (without scaling or adjusting for library size)- Greatly improved
generateReport()
NEW FEATURES
- Added
generateReport()
BUG FIXES
- Fixed bugs in
mergeResults()
SIGNIFICANT USER-VISIBLE CHANGES
getSegmentsRle()
was greatly simplified
BUG FIXES
- Fixed
analyzeChr()
, completedmergeResults()
NEW FEATURES
- Added
analyzeChr()
andmergeResults()
SIGNIFICANT USER-VISIBLE CHANGES
- Updated README.md
makeModels()
now usestestvars
instead ofgroup
and has a new argumentsgroupInfo
,center
andtestIntercept
calculatePvalues()
now uses area of regions instead of mean to calculate the p-values.preprocessCoverage()
now calculates the mean coverage at each base
SIGNIFICANT USER-VISIBLE CHANGES
- Users can specify significance cutoffs for
plotOverview()
andplotRegion()
SIGNIFICANT USER-VISIBLE CHANGES
calculatePvalues()
now usesqvalue::qvalue()
instead ofp.adjust()
SIGNIFICANT USER-VISIBLE CHANGES
- Added
plotOverview()
SIGNIFICANT USER-VISIBLE CHANGES
- Added
plotRegion()
BUG FIXES
makeModels()
can now handle a vector for theadjustvars
argument
SIGNIFICANT USER-VISIBLE CHANGES
calculatePvalues()
will adjust thep.values
now usingp.adjust()
makeModels()
now can handle a matrix for thegroup
argument
BUG FIXES
getSegmentsRle()
will now work properly in the case that no segments are found
SIGNIFICANT USER-VISIBLE CHANGES
- Attempted to reduce memory load in
calculateStats()
andcalculatePvalues()
SIGNIFICANT USER-VISIBLE CHANGES
fstats.apply()
now usesRcpp
andRcppArmadillo
NEW FEATURES
preprocessCoverage()
can now automatically select thechunksize
SIGNIFICANT USER-VISIBLE CHANGES
- Dropped
fstats()
and is now part offstats.apply()
NEW FEATURES
- Introduced
method
argument for getSegmetnsRle()
BUG FIXES
- No longer permuting the data matrix in
calculatePvalues()
NEW FEATURES
- Added
makeBamList()
,makeModels()
, andpreprocessCoverage()
- Added example data from the Montogemery and Pickrell studies. Check
?genomeData
and?genomeInfo
NEW FEATURES
- Added
calculatePvalues()
NEW FEATURES
- Added
clusterMakerRle()
,findRegions()
, andgetSegmentsRle()
NEW FEATURES
- Added
calculateStats()
,filterData()
,fstats()
, andfstats.apply()
SIGNIFICANT USER-VISIBLE CHANGES
- Renamed
makeCoverage()
toloadCoverage()
- Improved
NAMESPACE
NEW FEATURES
- Initialized the package (named
derfinder2
) fromderfinder
https://github.com/alyssafrazee/derfinder version 1.0.2 This version is available at https://github.com/alyssafrazee/derfinder/tree/d49f7b28c26f075da36a50ab67c9d192ab2fd63d