Fishpond with unbalanced dataset
1
0
Entering edit mode
samuelcjx • 0
@cdc1cdaa
Last seen 9 months ago
United States

I have two conditions: negative control and treatment. Each one has 6 replicates. Then, from PCA analysis and heatmap, I found out that one negative control replicate is very different from the rest of the samples and that it is an outlier, so I decided to remove that sample from the dataset for downstream analysis. I proceeded to use fishpond by calling scaleInfReps(), labelKeep(), and swish(). However, swish() returned the following error.

Error in infRepsArray[, cond2, ] : subscript out of bounds 

I suspect that this has to do with me having 5 samples in one condition and 6 samples in the other. Is that right? What should I do here?

EDIT: Fishpond version is 2.6.0

EDIT2: These replicates are biological, rather than technical.

EDIT3: I will include the PCA plot below

PCAplot

rna-seq salmon fishpond swish • 673 views
ADD COMMENT
0
Entering edit mode

After some debugging, I believe that there is a bug in fishpond codebase. I've opened an issue.

ADD REPLY
0
Entering edit mode
@mikelove
Last seen 59 minutes ago
United States

Thanks for the report, I will follow up on GH.

ADD COMMENT
0
Entering edit mode

Thanks again for identifying this, I've pushed what I believe to be a fix in 81a2825.

This is also now pushed to Bioc as version 2.6.1 (release) and 2.7.1 (devel), should be available in 24 hours.

ADD REPLY

Login before adding your answer.

Traffic: 516 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