Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

preprocessSWAN fails to normalize RGChannelSet #262

Open
alexandernassar24 opened this issue Jan 22, 2024 · 0 comments
Open

preprocessSWAN fails to normalize RGChannelSet #262

alexandernassar24 opened this issue Jan 22, 2024 · 0 comments

Comments

@alexandernassar24
Copy link

Hello,

I am having issues normalizing an RGChannelSet. I initialize the object with two table that have probes as rownames and methylated/unmethylated signals for each sample as column names:

#create minfi object
RGChannelSet <- RGChannelSet(Green = green_methylated, Red = red_unmethylated)

#specify annotation
RGChannelSet@annotation <- c(array = "IlluminaHumanMethylation450k", annotation = "ilmn12.hg19")

RGChannelSet

class: RGChannelSet
dim: 485512 74
metadata(0):
assays(2): Green Red
rownames(485512): cg00050873 cg00212031 ... ch.22.47579720R ch.22.48274842R
rowData names(0):
colnames(74): 6929689068_R01C01 6929689068_R04C02 ... 9296930148_R02C01
9296930138_R06C01
colData names(0):
Annotation
array: IlluminaHumanMethylation450k
annotation: ilmn12.hg19

Running the preprocessingRAW function yields an empty matrix.
Running preprocessSWAN on the dataset yields the following errors:

Warning: no non-missing arguments to min; returning Inf
Error in sample.int(length(x), size, replace, prob) : vector size cannot be infinite

I'm not sure why the preprocessing is failing. Is it because of a mistake on my behalf in initializing the object? Guidance would be appreciated.

Thank you,
Alexander

P.S. As a starting point, I am using the non-normalized table (merged for methylated/unmethylated signals for all samples) in the supplementary files of this GEO accession page: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE52556

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant