Question: Why procset(files, procol) requires exactly one column name?
0
gravatar for klv2706w
13 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 • 171 views
ADD COMMENTlink modified 13 months ago • written 13 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: 130 users visited in the last hour