## ----message=FALSE, warning=FALSE, eval=TRUE----------------------------------
library(GEOquery)
library(SingleCellExperiment)
library(DropletUtils)


## -----------------------------------------------------------------------------
res = searchGEO("single cell mtx matrix")
head(res)


## -----------------------------------------------------------------------------
poss= getGEOSuppFiles('GSE248214', fetch_files = FALSE)
poss


## -----------------------------------------------------------------------------
library(GEOquery)
s = getGEOSuppFiles('GSE161228', fetch_files = FALSE)
s


## -----------------------------------------------------------------------------
s = getGEOSuppFiles('GSM7908437')
s

## -----------------------------------------------------------------------------
sample_prefix = sub('barcodes.tsv.gz','', grep('barcodes.tsv.gz', s$filepath, value = TRUE))
gsm_sce <- DropletUtils::read10xCounts(samples = sample_prefix, type = 'prefix')
gsm_sce
head(colData(gsm_sce))


## -----------------------------------------------------------------------------
s = getGEOSuppFiles('GSE248214')
s


## -----------------------------------------------------------------------------
tar_filename = rownames(s)[1]
exdir = tempdir()
untar(tar_filename, exdir=exdir)
sample_prefixes = sub('matrix.mtx.gz','',dir(exdir,pattern='*matrix.mtx.gz',full.names = TRUE))
gse_sce = DropletUtils::read10xCounts(samples = sample_prefixes, type='prefix')
gse_sce


## -----------------------------------------------------------------------------
getGEOSeriesFileListing('GSE288770')


## -----------------------------------------------------------------------------
s = getGEOSuppFiles('GSM8775062')
sce = DropletUtils::read10xCounts(s$filepath)
sce


## -----------------------------------------------------------------------------
s = getGEOSuppFiles('GSM8775066')
sce2 = DropletUtils::read10xCounts(samples = '/Users/davsean/Documents/git/GEOquery/vignettes/GSM8775066/GSM8775066_vivo_day7_1_', type='prefix')
sce2