swh:1:snp:a0c50030ca57fc477cf863f130d89b83459157c3
Tip revision: 5bcaa1984bd17e1409b28344ec900233a21c5364 authored by Alexey Sergushichev on 13 November 2021, 18:46:20 UTC
bump version in js
bump version in js
Tip revision: 5bcaa19
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()