zwdzwd / sesame

🍪 SEnsible Step-wise Analysis of DNA MEthylation BeadChips
Other
60 stars 31 forks source link

Error using predictedAge #152

Open WalterPavicic opened 7 months ago

WalterPavicic commented 7 months ago

Hi, I was trying to calculate age using the predictionAge Func, but always get the same error. Using the example on guides, works:

betas_o <- getBetas(sesameDataGet("EPIC.1.SigDF"), mask=F) model <- readRDS(url("https://github.com/zhou-lab/InfiniumAnnotationV1/raw/main/Anno/HM450/Clock_Horvath353.rds")) predictAge(betas_o, model) [1] 25.7583

But not when running on my data: Ex1:

sdfs = openSesame(idat_test1s, func = NULL) beta_1s = openSesame(idat_test1s, BPPARAM = BiocParallel::MulticoreParam(2)) predictAge(beta_1s, model) Error in predictAge(beta_1s, model) : Fewer than 10 matching probes left. Age prediction abort. sum(is.na(sdfs)) [1] 1618794

Ex2:

sdfsR2 = openSesame(idat_test1s, func = NULL) sum(is.na(sdfsR2)) [1] 1618794 beta_R2_2 = openSesame(idat_test1s, BPPARAM = BiocParallel::MulticoreParam(2), mask=FALSE) predictAge(beta_R2_2, model) Error in predictAge(beta_R2_2, model) : Fewer than 10 matching probes left. Age prediction abort.

I dont know if there is something Im missing?

thanks in adfvance, Walter!

Archie41025 commented 3 weeks ago

Hi Walter,

I seem to be encountering the same problem. Did you manage to solve it? And how?