Qualifier / getQAStats: Problem to get statisric from a GatingSet
1
0
Entering edit mode
veronikag • 0
@veronikag-9110
Last seen 7.6 years ago
Germany

Hi,

I have one more question to the getQAStats command:

First, I used the following command and it worked:

test <- getQAStats(gs[[1]],isMFI=TRUE)

Second, I tried this command  :

test2 <- getQAStats(gs,isMFI=TRUE)

And I´m getting the following Error:

[1] "extracting stats..."

Error in getQAStats(gs, isMFI = TRUE) :

  Not all IDs for the current sample set are found in meta data of this GatingSet!

 

gs is my GatingSet with 6 samples. With the first command I can extract the statistic for each sample without any errors and warnings.

I would really appreciate your help.

Thank you in advance,

Veronika

qualifier • 1.7k views
ADD COMMENT
0
Entering edit mode

getQAStats method for 'GatingSet' was written in the context of QA report and thus not intended to be used without QA preprocess steps. However here is the workaround,

res <- lapply(gs, getQAStats, isMFI = TRUE)

Once you get a list of data.table, it is not difficult to combine them through 'rbindlist'

 

ADD REPLY
0
Entering edit mode

Hi Mike,

unfortunately, I´m getting the following Error with your command:

Error in as.list.default(X) : 
  no method for coercing this S4 class to a vector

and I´m not sure, what this is telling me?

Again, I would really appreciate your help!

Best,

Veronika

ADD REPLY
0
Entering edit mode

I don't know why this is happening. Either your 'gs' variable is no longer a valid 'GatingSet' object or 'lapply' somehow gets masked by some other package. Anyway, start up a fresh R session and try the following reproducible example to see if this works:

library(QUALIFIER)
gs <- load_gs(system.file("/extdata/gs_manual/", package = "flowWorkspaceData"))
res <- lapply(gs, getQAStats, isMFI = TRUE)

Alternatively, try to put namespace explicitly before the lapply:

res <- flowWorkspace::lapply(gs, getQAStats, isMFI = TRUE)

 

ADD REPLY
0
Entering edit mode

Thank you! Now, it is working. It was probably masked by some other package.

Best,

Veronika

ADD REPLY
0
Entering edit mode

Hi Mike,

If I try to do the same, I get: 

library(QUALIFIER)
gs <- load_gs(system.file("/extdata/gs_manual/", package = "flowWorkspaceData"))
res <- flowWorkspace::lapply(gs, getQAStats, isMFI = TRUE)
Error in data.table(channel = thisChnl, stain = thisStain, stats = names(statsOfNode),  : 
  column or argument 3 is NULL

Not sure why?

 

ADD REPLY
0
Entering edit mode
Jiang, Mike ★ 1.3k
@jiang-mike-4886
Last seen 2.6 years ago
(Private Address)

I pushed the bug fix to github. Try devtools::install_github("RGLab/QUALIFIER", ref = "trunk)

ADD COMMENT
0
Entering edit mode

Yes, that works now. Thanks.

ADD REPLY

Login before adding your answer.

Traffic: 765 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6