[med-svn] [r-bioc-deseq2] 01/06: New upstream version 1.14.1
Andreas Tille
tille at debian.org
Sun Jan 15 18:49:47 UTC 2017
This is an automated email from the git hooks/post-receive script.
tille pushed a commit to branch master
in repository r-bioc-deseq2.
commit a9c33f5cda06b75407238f9abb0f82508e72b8ac
Author: Andreas Tille <tille at debian.org>
Date: Sun Dec 4 14:10:09 2016 +0100
New upstream version 1.14.1
---
DESCRIPTION | 4 ++--
NEWS | 5 +++++
build/vignette.rds | Bin 231 -> 231 bytes
inst/doc/DESeq2.pdf | Bin 673682 -> 672224 bytes
src/DESeq2.cpp | 10 +++++-----
5 files changed, 12 insertions(+), 7 deletions(-)
diff --git a/DESCRIPTION b/DESCRIPTION
index af2bd06..fb3ae36 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -2,7 +2,7 @@ Package: DESeq2
Type: Package
Title: Differential gene expression analysis based on the negative
binomial distribution
-Version: 1.14.0
+Version: 1.14.1
Author: Michael Love, Simon Anders, Wolfgang Huber
Maintainer: Michael Love <michaelisaiahlove at gmail.com>
Description: Estimate variance-mean dependence in count data from
@@ -22,4 +22,4 @@ biocViews: Sequencing, ChIPSeq, RNASeq, SAGE, DifferentialExpression,
GeneExpression, Transcription
RoxygenNote: 5.0.1
NeedsCompilation: yes
-Packaged: 2016-10-17 23:12:25 UTC; biocbuild
+Packaged: 2016-12-01 00:14:26 UTC; biocbuild
diff --git a/NEWS b/NEWS
index 512a034..09b5455 100644
--- a/NEWS
+++ b/NEWS
@@ -1,3 +1,8 @@
+CHANGES IN VERSION 1.14.1
+-------------------------
+
+ o Removed deprecated armadillo argument in DESeq2.cpp
+
CHANGES IN VERSION 1.13.8
-------------------------
diff --git a/build/vignette.rds b/build/vignette.rds
index a30845c..13dfd64 100644
Binary files a/build/vignette.rds and b/build/vignette.rds differ
diff --git a/inst/doc/DESeq2.pdf b/inst/doc/DESeq2.pdf
index 0c61266..40d0612 100644
Binary files a/inst/doc/DESeq2.pdf and b/inst/doc/DESeq2.pdf differ
diff --git a/src/DESeq2.cpp b/src/DESeq2.cpp
index 0efa29a..8fc4508 100644
--- a/src/DESeq2.cpp
+++ b/src/DESeq2.cpp
@@ -56,7 +56,7 @@ double dlog_posterior(double log_alpha, Rcpp::NumericMatrix::Row y, Rcpp::Numeri
arma::mat dw = arma::diagmat(Rcpp::as<arma::vec>(dw_diag));
arma::mat b = x.t() * w * x;
arma::mat db = x.t() * dw * x;
- double ddetb = ( det(b) * trace(b.i(true) * db) );
+ double ddetb = ( det(b) * trace(b.i() * db) );
double cr_term = -0.5 * ddetb / det(b);
double alpha_neg1 = R_pow_di(alpha, -1);
double alpha_neg2 = R_pow_di(alpha, -2);
@@ -84,10 +84,10 @@ double d2log_posterior(double log_alpha, Rcpp::NumericMatrix::Row y, Rcpp::Numer
Rcpp::NumericVector d2w_diag = 2 * pow(pow(mu, -1) + alpha, -3);
arma::mat d2w = arma::diagmat(as<arma::vec>(d2w_diag));
arma::mat b = x.t() * w * x;
- arma::mat b_i = b.i(true);
+ arma::mat b_i = b.i();
arma::mat db = x.t() * dw * x;
arma::mat d2b = x.t() * d2w * x;
- double ddetb = ( det(b) * trace(b.i(true) * db) );
+ double ddetb = ( det(b) * trace(b.i() * db) );
double d2detb = ( det(b) * (R_pow_di(trace(b_i * db), 2) - trace(b_i * db * b_i * db) + trace(b_i * d2b)) );
double cr_term = 0.5 * R_pow_di(ddetb/det(b), 2) - 0.5 * d2detb / det(b);
double alpha_neg1 = R_pow_di(alpha, -1);
@@ -345,10 +345,10 @@ Rcpp::List fitBeta(SEXP ySEXP, SEXP xSEXP, SEXP nfSEXP, SEXP alpha_hatSEXP, SEXP
beta_mat.row(i) = beta_hat.t();
// recalculate w so that this is identical if we start with beta_hat
w = diagmat(mu_hat/(1.0 + alpha_hat[i] * mu_hat));
- hat_matrix = sqrt(w) * x * (x.t() * w * x + ridge).i(true) * x.t() * sqrt(w);
+ hat_matrix = sqrt(w) * x * (x.t() * w * x + ridge).i() * x.t() * sqrt(w);
hat_diagonals.row(i) = diagvec(hat_matrix).t();
// sigma is the covariance matrix for the betas
- sigma = (x.t() * w * x + ridge).i(true) * x.t() * w * x * (x.t() * w * x + ridge).i(true);
+ sigma = (x.t() * w * x + ridge).i() * x.t() * w * x * (x.t() * w * x + ridge).i();
contrast_num.row(i) = contrast.t() * beta_hat;
contrast_denom.row(i) = sqrt(contrast.t() * sigma * contrast);
beta_var_mat.row(i) = diagvec(sigma).t();
--
Alioth's /usr/local/bin/git-commit-notice on /srv/git.debian.org/git/debian-med/r-bioc-deseq2.git
More information about the debian-med-commit
mailing list