Skip to content

Commit

Permalink
Merge pull request #136 from FertigLab/nextflow
Browse files Browse the repository at this point in the history
debug adata to dgc mat
  • Loading branch information
dimalvovs authored Dec 13, 2024
2 parents dc800bc + fbf507e commit c6ade6f
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions nextflow/main.nf
Original file line number Diff line number Diff line change
Expand Up @@ -135,12 +135,12 @@ process ADATA_2DGCMAT {
i <- hdf5r::readDataSet(f[["X/indices"]]);
p <- hdf5r::readDataSet(f[["X/indptr"]]);
x <- hdf5r::readDataSet(f[["X/data"]]);
dims <- c(hdf5r::h5attributes(f[["X/"]])[["shape"]][1],
hdf5r::h5attributes(f[["X/"]])[["shape"]][2]);
dims <- c(hdf5r::h5attributes(f[["X/"]])[["shape"]][2],
hdf5r::h5attributes(f[["X/"]])[["shape"]][1]);
res <- Matrix::sparseMatrix(i = i, p = p, x = x, dims = dims, index1=FALSE);
message("Read matrix with dimensions: ", dims[1],",",dims[2]);
colnames(res) <- hdf5r::readDataSet(f[["var/_index"]]);
rownames(res) <- hdf5r::readDataSet(f[["obs/_index"]]);
rownames(res) <- hdf5r::readDataSet(f[["var/_index"]]);
colnames(res) <- hdf5r::readDataSet(f[["obs/_index"]]);
hdf5r::h5close(f);
message("Normalizing data");
res <- Seurat::NormalizeData(res);
Expand Down

0 comments on commit c6ade6f

Please sign in to comment.