diff --git a/DESCRIPTION b/DESCRIPTION index 7349fed..c531ea8 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -1,8 +1,8 @@ Package: dreamlet Type: Package Title: Scalable differential expression analysis of single cell transcriptomics datasets with complex study designs -Version: 1.3.2 -Date: 2024-06-10 +Version: 1.3.3 +Date: 2024-07-19 Authors@R: person("Gabriel", "Hoffman", role=c("aut", "cre"), email="gabriel.hoffman@mssm.edu", comment=c(ORCID="0000-0002-0957-0224")) Description: Recent advances in single cell/nucleus transcriptomic technology has enabled collection of cohort-scale datasets to study cell type specific gene expression differences associated disease state, stimulus, and genetic regulation. The scale of these data, complex study designs, and low read count per cell mean that characterizing cell type specific molecular mechanisms requires a user-frieldly, purpose-build analytical framework. We have developed the dreamlet package that applies a pseudobulk approach and fits a regression model for each gene and cell cluster to test differential expression across individuals associated with a trait of interest. Use of precision-weighted linear mixed models enables accounting for repeated measures study designs, high dimensional batch effects, and varying sequencing depth or observed cells per biosample. VignetteBuilder: knitr @@ -18,6 +18,7 @@ Suggests: muscat, ExperimentHub, RUnit, + muscData, scater, scuttle biocViews: @@ -70,6 +71,8 @@ Imports: tidyr, rlang, BiocGenerics, + S4Arrays, + SparseArray, DelayedArray, gtools, reshape2, diff --git a/NAMESPACE b/NAMESPACE index c06b2e0..af1f9ca 100644 --- a/NAMESPACE +++ b/NAMESPACE @@ -66,13 +66,13 @@ exportMethods(topTable) exportMethods(zenith_gsa) import(Rcpp) import(ggplot2) +importClassesFrom(SparseArray,SparseArray) importClassesFrom(limma,EList) importFrom(BiocParallel,SerialParam) importFrom(BiocParallel,bplapply) importFrom(DelayedArray,colsum) importFrom(DelayedArray,defaultAutoGrid) importFrom(DelayedArray,getAutoBlockSize) -importFrom(DelayedArray,read_block) importFrom(DelayedArray,setAutoBlockSize) importFrom(DelayedMatrixStats,colSums2) importFrom(GSEABase,geneIds) @@ -91,6 +91,7 @@ importFrom(MatrixGenerics,rowSds) importFrom(MatrixGenerics,rowSums2) importFrom(MatrixGenerics,rowVars) importFrom(Rdpack,reprompt) +importFrom(S4Arrays,read_block) importFrom(S4Vectors,"metadata<-") importFrom(S4Vectors,DataFrame) importFrom(S4Vectors,as.data.frame) diff --git a/NEWS.md b/NEWS.md index 9fcc825..1afe0c6 100644 --- a/NEWS.md +++ b/NEWS.md @@ -1,4 +1,9 @@ +# dreamlet 1.3.3 + - July 19, 2024 + - fix in `.read_matrix_block()` + - See https://github.com/GabrielHoffman/dreamlet/pull/23/commits/c453ac98ebc0329279b4dd3ae26a674df0e9b1f2 + # dreamlet 1.3.1 - May 31, 2024 - bump Bioc version diff --git a/R/colsum_fast.R b/R/colsum_fast.R index 311f3f3..fbbd544 100644 --- a/R/colsum_fast.R +++ b/R/colsum_fast.R @@ -11,12 +11,13 @@ create_idxlist <- function(fct) { res } -# copied from DelayedArray:::.read_matrix_block +# patched: https://github.com/GabrielHoffman/dreamlet/pull/23/commits/c453ac98ebc0329279b4dd3ae26a674df0e9b1f2 #' @importFrom methods as -#' @importFrom DelayedArray read_block +#' @importFrom S4Arrays read_block +#' @importClassesFrom SparseArray SparseArray .read_matrix_block <- function(...) { block <- read_block(..., as.sparse = NA) - if (is(block, "SparseArraySeed")) { + if (is(block, "SparseArray")) { block <- as(block, "CsparseMatrix") } block diff --git a/docs/404.html b/docs/404.html index b0024d5..c9860d1 100644 --- a/docs/404.html +++ b/docs/404.html @@ -14,73 +14,55 @@ - - - + - Skip to contents - -