drawC1C2Density()
gained argument grid
.Deprecated nbrOfFiles()
methods; please use length()
on corresponding input data sets.
Removed defunct callPeaks()
for data.frame:s.
CLEANUP: Now doSegmentByPairedPSCBS()
uses future_mapply()
of future.apply instead of deprecated dsApplyInPairs()
of R.filesets.
Utilizing subsetted calculations of matrixStats (>= 0.50.0).
CLEANUP: Now importing generic extractC1C2()
from PSCBS to avoid creating a new one.
Package now requires R (>= 3.1.1) released July 2014. This allows us to use Bioconductor (>= 3.0) (October 2014).
Bumped package dependencies.
ROBUSTNESS: Explicitly importing core R functions.
makeSmoothSplinePredict()
defunct since Aug 2013. Made callPeaks()
for data.frame:s defunct (was deprecated).ROBUSTNESS: Added the first package tests.
Bumped package dependencies.
Package now requires R (>= 3.0.3) and Bioconductor (>= 2.13) which were released March 2014 and are in fact old; it’s recommended to use a more recent version of R.
process()
for AbstractCurveNormalization
would generate an error due to read-only permissions introduced by copying the target file without resetting the file permissions.Package now requires R (>= 3.0.0) and BioC (>= 2.13), which were released April 2013 and are in fact old and it’s recommended to use a more recent version of R.
Updated package dependencies.
doSegmentByPairedPSCBS()
for AromaUnitPscnBinarySet
.Updated package dependencies.
Package requires R (>= 2.15.1) and Bioconductor (>= 2.11.0).
points()
for C1C2
passes (modified) argument x
to NextMethod()
as object = x
.:::
in calls.Minor tweaks to NAMESPACE.
Updated package dependencies.
Package requires R (>= 2.15.0) and Bioconductor (>= 2.10.0).
ROBUSTNESS: Now importing only what needs to be imported and formally declaring all S3 methods in NAMESPACE.
CLEANUP: Dropped obsolete usage of autoload()
.
**aroma.cn**
Package object is also available when the package is only loaded (but not attached).CLEANUP: Now importing only what is needed from PSCBS.
Updated package dependencies.
byPath()
, byName()
, and findByPath()
for PairedPSCBSFileSet
was also affected by the bug described in the R-devel thread ‘Do not pass’…’ to NextMethod() - it’ll do it for you; missing documentation, a bug or just me?’ on Oct 16, 2012.
getPath()
for PairedPscbsModel
would throw an error on getInputDataSet()
not defined.
makeSmoothSplinePredict()
defunct.SPEEDUP: Replaced all rm()
calls with NULL assignments. Also removed several explicit garbage collector calls.
CLEANUP: The formal package dependency on Bioconductor package aroma.light has been relaxed so the package can be installed without it.
CLEANUP: Package now only imports R.oo.
Updated package dependencies.
CRAN POLICY: Now all Rd \usage{}
lines are at most 90 characters long.
CRAN POLICY: Now all Rd example lines are at most 100 characters long.
findNeutralCopyNumberState()
is now in PSCBS.startupMessage()
of R.oo.Updated package dependencies.
CLEANUP: No longer using deprecated PSCBS methods.
ROBUSTNESS: {load,save}Cache()
from R.cache are now explicitly imported in the namespace.
PairedPSCBS
to recognize when other mean-level estimators than the sample mean have been used.process()
for PairedPscbsCaller
used the global verbose
.example()
scripts used non-defined values.Now applicable classes utilize the new ParametersInterface
.
DOCUMENTATION: Hiding more internal methods from the help indices.
cache:
field modified instead of cached:
. After correction, all clearCache()
methods could be dropped.seq_along(x)
instead of seq(along = x)
everywhere. Similarly, seq(ds)
where ds
is GenericDataFileSet
is now replaced by seq_along(ds)
. Likewise, seq_len(x)
replaces seq(length = x)
, and length(ds)
replaces nbrOfFiles(ds)
.whichVector()
with which()
, because the latter is now the fastest again.R CMD check
. The reason was that some of the internal methods call PSCBS methods, which only happens if PSCBS is loaded in the first place but R CMD check
cannot known that.Arguments$get{Read,Writ}ablePath()
instead of filePath(..., expandLinks = "any")
.Arguments$getWritablePath()
everywhere instead of mkdirs()
, because the former will do a better job in creating and asserting directories on slow shared file systems, and when it fails it gives a more informative error message.Object
methods that calls “next” methods, utilizes NextMethod()
, which became possible with R.oo v1.10.0....
to NextMethod()
, cf. R-devel thread ‘Do not pass’…’ to NextMethod() - it’ll do it for you; missing documentation, a bug or just me?’ on Oct 16, 2012.getOutputFileClass()
and getOutputFileExtension()
for TotalCnSmoothing
.Now PairedPscbsCaller()
passes ...
to the internal callers, which makes it possible to for instance specify the number of bootstrap samples done for the AB caller.
Now PairedPscbsModel()
excludes the actual gaps from the known segments it passes to segmentByPairedPSCBS()
.
PairedPscbsCaller
.callPeaks(..., flavor="all")
for PeaksAndValleys
would return an error.calculateTumorPSCNByGenotype()
.fit()
for PairedPscbsModel
generates pair names iff tumor and normal names don’t match, e.g. GSM517071_vs_GSM517072
(if match then just Patient1
). It also generated pair
tags.DOCUMENTATION: Improved help on TotalCnBinnedSmoothing
.
Bumped up the package dependencies.
PairedPscbsModel
.Adopted findAtomicAberrations()
for CBS
from ditto of PairedPSCBS
.
GENERALIZATION: Now plotTracks()
for PruneCNA
supports CBS
segmentation results in additional to PairedPSCBS
ones.
GENERALIZATION: Now pruneCNA()
is implemented for AbstractCBS
, not just PairedPSCBS
objects.
Merged updates for findAtomicAberrations()
for PairedPSCBS
and some additional internal “equality” test functions.
normalizePrincipalCurve()
.drawC1C2Density()
for PairedPSCBS
would throw an exception if there was only one segment, or less than two finite (C1,C2):s.Updated package dependencies.
Additional internal updates.
TotalCnBinnedSmoothing()
.CLEANUP: The example code for the internal PairedPSCBS methods now only runs if environment variable _R_CHECK_FULL_
is set. This makes the package easier on the CRAN servers.
ROBUSTNESS: Updated package dependencies.
ROBUSTNESS: Now process()
of TotalCnSmoothing
writes atomically.
Additional internal updates.
Updated package dependencies.
Additional internal updates.
Updated the package dependencies.
Some internal updates.
callPeaks()
for PeaksAndValleys
.Now deShearC1C2()
, translateC1C2()
, and transformC1C2()
also update C1 and C2 mean levels.
Now using getLocusData()
and getSegments()
internally for all PairedPSCBS
objects wherever applicable.
cat()
and getOption()
of R.utils (instead of base).Forgot to update some examples and test scripts which were still referring to the old psCBS package (should be PSCBS).
Updated internal code to work with the new column names in PSCBS v0.11.0.
hpaste()
internally wherever applicable.PairedPSCBSFile
and PairedPSCBSFileSet
.R CMD check
.ROBUSTNESS: Now all bootstrap methods utilize resample()
.
Added more internal utility functions for future usage.
preserveScale
to TumorBoostNormalization
for correcting for signal compression in heterozygous SNPs. The defaults is to do this correction.callXXorXY()
no longer calls gender from chr Y, when gender is estimated as XX
from chr X.Package submitted to CRAN.
Updated citation information.
Package now requires aroma.core v1.6.0.
Package pass R CMD check
on R v2.11.0 and v2.12.0 devel.
normalizeDifferencesToAverage()
, normalizeTumorBoost()
, callNaiveGenotypes()
, and internal findPeaksAndValleys()
to aroma.light v1.5.3.normalizeTumorBoost()
could introduce NaN:s if betaN
was exactly zero or exactly one.Added (for now internal) option to change the degrees of freedom of the fitted principal curves in MSCN.
Added plotSmoothedPairsOne()
to MultiSourceCopyNumberNormalization
.
h(.)
and g(.)
to AbstractCurveNormalization
, which allows us to fit say on the log scale, e.g. h(x) = log2(x)
, g(y) = 2^y
.getOutputDataSet()
of AbstractCurveNormalization
returned all files, not just the ones matching the input set.ROBUSTNESS: Using new Arguments$getInstanceOf()
were possible.
ROBUSTNESS: Now all index arguments are validated correctly using the new max
argument of Arguments$getIndices()
. Before the case where "max == 0"
was not handled correctly.
flavor = "v4"
of TumorBoostNormalization
the default, and if used then no "flavor"
tag is added.Now callXXorXY()
and callNaiveGenotypes()
handles missing values and non-finite values. They also censor outliers to become infinite/extreme values.
Added callXXorXY()
.
Added an example()
to the Rd help of callNaiveGenotypes()
.
Added Rd help to findPeaksAndValleys()
.
Now argument tol
of findPeaksAndValleys()
can be zero; before it had to be at least the smallest possible double.
CLEANUP: Removed suggested dependency on princurve, which is now indirectly suggested/requested via aroma.light.
More recent dependencies on Bioconductor packages.
Package passes R CMD check
on R v2.10.0.
Added normalizeTumorBoost()
for RawAlleleBFractions
.
Added callGenotypes()
for RawAlleleBFractions
.
Added RawGenotypeCalls
.
byPath()
instead fromFiles()
.Renamed argument alignByChromosome
for the constructor of the MultiSourceCopyNumberNormalization
class to align
in order to allow for more types of aligned.
The alignment of MultiSourceCopyNumberNormalization
is now done using normalizeDifferencesToAverage()
, which is robust against outliers, waviness, etc. The previous method which normalized toward the same overall median is no longer available.
normalizeDifferencesToAverage()
.getTags()
of MultiSourceCopyNumberNormalization
would return all asterisk tags as merged, e.g. c("mscn,align", "tagA", "tagB")
.XYCurveNormalization
and PrincipalCurveNormalization
.TumorBoostNormalization
: the srcFiles
attribute in file footer of the result files contained a duplicated default footer instead of the tumor-normal pair.callNaiveGenotype()
and normalizeTumorBoost()
.Added model flavor "v4"
which corrects heterozygots according to "v2"
and homozygotes according to "v1"
.
Added new model flavor ("v3"
) of TumorBoostNormalization
that is an extension of last weeks model flavor.
"v2"
) of TumorBoostNormalization
that avoids over correcting (especially at the heterozygotes), but adjusting the correction factor. Use argument flavor = "v2"
.TumorBoostNormalization
now only takes an AromaUnitGenotypeCallSet
for argument gcN
. It no longer takes an AromaUnitFracBCnBinarySet
object, which was only an ad hoc solution.Added argument alignByChromosomes
to MultiSourceCopyNumberNormalization
. If TRUE, the signals are shifted per chromosome such that the mean of the normalized smoothed signals is the same for all sources. This can for instance remove systematic effects on sex chromosomes added by some ad hoc preprocessing methods.
Added a clearCache()
to MultiSourceCopyNumberNormalization
.
ALPHA: Added TumorBoostNormalization
.
INTERNAL: Added foundations for TumorBoost, i.e. in memory classes such as TotalAndFracBSnpData
.
INTERNAL: Added findPeaksAndValleys()
.
ROBUSTNESS: Now all constructors report on unknown arguments.
ROBUSTNESS: Now MultiSourceCopyNumberNormalization
first write normalized data to a temporary file, which is then renamed. This lower the risk for having incomplete data in case of interrupts.
Now getOutputDataSets()
of MultiSourceCopyNumberNormalization
only returns output data files with a matching fullname in the input set.
Added missing argument verbose
in getTargetPositions()
of TotalCnSmoothing
. This caused unwanted verbose output in some cases.
process()
of TotalCnSmoothing
would not “recognize” fullname translators, that is, the output filenames were always identical to the input ones.
R CMD check
and all redundancy tests.RawGenomicSignals
.Added redundancy tests to package.
Further cleanup. Some functions are now in aroma.light.
{fit|backtransform}PrincipalCurve()
were moved to aroma.light v1.11.1.
Several classes and methods were moved to aroma.core v1.0.0.
Adopted the package to the new classes of aroma.core.
backtransformPrincipalCurve()
.MultiSourceCopyNumberNormalization
.R CMD check
on R v2.7.1 and v2.8.0.extractRawCopyNumbers()
of AromaTotalCnBinaryFile
adds annotation data fields to the returned object, e.g. platform, chipType, and the fullname of the source file.normalizePrincipalCurve()
and fitPrincipalCurve()
.ALPHA: Added extractRawCopyNumbers()
to AromaTotalCnBinaryFile
.
ALPHA: Added TotalCnSmoothing
.
Aroma{Total|FreqB}CnSignal{File|Set}
. With the more generalized aroma.core package, it is now possible retrieve the AromaUgpFile
for the above. This provides the necessary basic methods for plotting data along chromosomes.