Question: RPMM blc() function not completing?
0
gravatar for mjt.professional
4.7 years ago by
mjt.professional0 wrote:

 

Hello,

I don't see many questions on the forum regarding the RPMM package.  In particular, I am using the BMIQ script/method to normalize the Type 2 probes on the Illumina 450K DNA methylation chip.  This procedure uses a piece-wise beta distribution fitting procedure to then recalculate methylation values for these probes.

The BMIQ script calls the blc() function from the RPMM package.  I've noticed that fairly often, this function simply hangs.  It has both a maximum iterations parameter as well as a convergence tolerance threshold.  Neither of these seem to matter in terms of this bad behavior of the function.

Has anyone encountered this problem?

 

thanks

 

ADD COMMENTlink modified 4.7 years ago by andrew.j.skelton736.0k • written 4.7 years ago by mjt.professional0
1
gravatar for andrew.j.skelton73
4.7 years ago by
London
andrew.j.skelton736.0k wrote:

Have you considered trying the wateRmelon implementation of BMIQ, it might be more stable than RPMM's

ADD COMMENTlink written 4.7 years ago by andrew.j.skelton736.0k

I'll give that a shot. This is rather maddening as it's hanging somewhere deep in the code dependencies and seems to depend on the random seeding of the core optimization routine.

ADD REPLYlink modified 8 months ago by RamRS28k • written 4.7 years ago by mjt.professional0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 666 users visited in the last hour