rwehrens BatchCorrMetabolomics: doBC results in NA
Hi folks,
Wondering if anyone could help:
I've formatted my data to resemble the example data so I'm not getting what should be wrong here.
Trying to use doBC but only getting a logical with NAs on return.
I checked that all the inputs were the same as the example data and there doesn't seem to be an issue.
Any ideas or someone with previously similar issues?
Thank you!
> M2c <- doBC(data[,2], ref.idx = ref.idx,
+ batch.idx = batch_labels$Batch, method = "rlm",
+ seq.idx = batch_labels$SeqNr)
>
> M2c
[1] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
[38] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
[75] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
0 answers
No answers yet.
Log in to answer this question.
More posts like this
Finally solved this by taking a closer look at the raw code. Just in case someone has the same problem in the future: It was the minBsamp argument which requires a specified number of reference samples for each batch. If that number doesn't match or is not even across the dataset, the algorithm produces NA values.