The support.bioconductor.org editor has been updated to markdown! Please see more info at: Tutorial: Updated Support Site Editor

Question: Why procset(files, procol) requires exactly one column name?
0
gravatar for klv2706w
11 months ago by
klv2706w10
klv2706w10 wrote:

Hello,

I started using ArrayExpress today to read a processed data from a recent study - https://www.ebi.ac.uk/arrayexpress/experiments/E-MTAB-6153/

I used the example provided in the documentation of getAE() to retrieve and load the data of my interest.

mexp1422 = getAE("E-MEXP-1422", type = "full")
## Build a an ExpressionSet from the processed data
cnames = getcolproc(mexp1422)
MEXP1422proc = procset(mexp1422, cnames[2])

I am very puzzled why procset() requires exactly one column name. This totally doesn't make sense. In addition, it is very slow for my dataset (19000 single cells). I would like to improve speed by rewriting the code that reads the processed data. I want to use data.table in a way that I suggest here: https://github.com/vitkl/myArrayExpress/commit/53858134241a326dc73f4f70f7fd2622e1702ee5 (line 78). For that, I need to know why procset() requires exactly one column name and whether I can remove this condition.

Thank you,

Vitalii

 

arrayexpress • 149 views
ADD COMMENTlink modified 11 months ago • written 11 months ago by klv2706w10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 247 users visited in the last hour