diff --git a/DESCRIPTION b/DESCRIPTION
index 0d233a2..de3881a 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,7 +1,7 @@
 Package: glmGamPoi
 Type: Package
 Title: Fit a Gamma-Poisson Generalized Linear Model
-Version: 1.17.3
+Version: 1.17.4
 Authors@R: c(person("Constantin", "Ahlmann-Eltze", email = "artjom31415@googlemail.com", 
                     role = c("aut", "cre"), comment = c(ORCID = "0000-0002-3762-068X")),
              person("Nathan", "Lubock", email="lubock.nathan@gmail.com", role = "ctb", 
@@ -43,6 +43,7 @@ Imports:
     DelayedMatrixStats,
     matrixStats,
     MatrixGenerics,
+    SparseArray (>= 1.5.21),
     DelayedArray,
     HDF5Array,
     SummarizedExperiment,
diff --git a/R/estimate_size_factors.R b/R/estimate_size_factors.R
index 048e164..b9e56c4 100644
--- a/R/estimate_size_factors.R
+++ b/R/estimate_size_factors.R
@@ -91,7 +91,7 @@ combine_size_factors_and_offset <- function(offset, size_factors, Y, verbose = F
     stopifnot(length(offset) == 1 || length(offset) == n_samples)
     zero_offset <- all(offset == 0)
     if(make_offset_hdf5_mat){
-      offset_matrix <- DelayedArray::DelayedArray(DelayedArray::SparseArraySeed(c(n_genes, n_samples)))
+      offset_matrix <- DelayedArray::DelayedArray(SparseArray::COO_SparseArray(c(n_genes, n_samples)))
       offset_matrix <- add_vector_to_each_row(offset_matrix, offset)
     }else{
       offset_matrix <- matrix(offset, nrow=n_genes, ncol = n_samples, byrow = TRUE)