se100k {restfulSE} | R Documentation |
Convenience functions using EC2 server to extract tenx neurons full or subset data
se100k(url = "http://h5s.channingremotedata.org:5000", tag = "tenx_100k_sorted") se1.3M(url = "http://h5s.channingremotedata.org:5000", tag = "tenx_full")
url |
server URL |
tag |
string giving the internal dataset name |
RESTfulSummarizedExperiment
SummarizedExperiment instance
se1.3M provides access to the full 1.3 million neurons with features in their order as given in the original HDF5 while se100k provides access to only 100k neurons with expression features sorted by genomic location
ss = se100k() # get a set of genes from Tasic et al. 2016 Nature Neuroscience tc = tasicCortex() adultCort = tc$GENEID # subset csums = apply(assay(ss[adultCort,1:500]),1,sum) names(csums) = tc$SYMBOL csums