swh:1:snp:0c004a03453a29b80f921a24433f7e780b9ceb53
Tip revision: c6bb0e960554a23eb712690cbd7f8e3f7d79ca0d authored by Alexey Sergushichev on 02 July 2021, 13:19:49 UTC
Update js
Update js
Tip revision: c6bb0e9
gen_toy_mouse_matrix_h5.R
library(GEOquery)
library(rhdf5)
es <- getGEO("GSE99709")[[1]]
srcfile <- "./cache/mouse_matrix.h5"
samples <- h5read(srcfile, "meta/Sample_geo_accession")
genes <- as.character(h5read(srcfile, "meta/genes"))
keepGSMs <- head(intersect(es$geo_accession, samples), 3)
es <- es[, es$geo_accession %in% keepGSMs]
sampleIndexes <- match(es$geo_accession,
samples)
expression <- h5read(srcfile,
"data/expression",
index=list(seq_along(genes),
stats::na.omit(sampleIndexes)))
rownames(expression) <- genes
colnames(expression) <- colnames(es)[!is.na(sampleIndexes)]
H5close()
destfile <- "./inst/testdata/mouse_matrix.h5"
h5createFile(destfile)
h5createGroup(destfile, "data")
# h5createDataset(file=destfile,
# dataset="data/expression",
# dims=dim(expression),
# storage.mode="integer")
h5write(expression,
file=destfile,
name="data/expression",
index=list(
NULL,
NULL))
h5createGroup(destfile, "meta")
# h5createDataset(file=destfile,
# dataset="meta/genes",
# dims=length(genes),
# storage.mode="character",
# size = max(nchar(genes)*2))
h5write(genes,
file=destfile,
name="meta/genes")
h5write(colnames(expression),
file=destfile,
name="meta/Sample_geo_accession")
H5close()