Hawaii Gigas Methylation Analysis Part 28
Parameterizing methylKit
Steven ran a methylKit parameterization R script for me on klone. I’m (finally) following up on some of the comments and output!
Evaluating methylKit output
The first thing I wanted to do was determine which meth.diff and min.per.group parameters provided the best balance between a conservative DML estimate and providing DML to analyze. Steven said the output from the script could be found in this folder. However, there was only output for the ploidy DML for min.per.group 9 through 12! I’m missing information for min.per.group = 8, and I’m not sure why there isn’t any pH output. I posted a comment in the same issue.
Going forward
- Continue parameter testing for DML identification
- Revise
methylKitmethods and results methylKItrandomization test- ATAC-Seq data integration
KOG-MWUfor Crassostrea methylation comparison- Revise discussion
- Revise introduction
- Transfer scripts used to a nextflow workflow
Written on September 23, 2026