RRBS data analysis using Bismark and methylkit
Once we have read in the cytosine report files with methRead(), how to generate a table/dataframe of beta values or equivalent for all samples as columns and CpG sites as rows? Is the percMethylation() function the correct function to use? Beta values I guess are from microarray, is percent-methylation the correct equivalent?
• 1,650 views
•
link
updated
by
Papyrus
•
written
by
Arindam Ghosh
•
written
by
Arindam Ghosh
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Plots for DNA methylation vs Gene Expression
written by sswang25 •
-
How to draw histogram plot of methylation values for all CpGs, CHGs and CHHs contexts
written by anisa •
-
Generating 5mCG or 5hmCG values per gene from BS/OXBS data
written by hchetia •
-
Bismark to MethylKit
written by Arindam Ghosh
-
How to check for the presence/absence of unique rows in two different matrices?
written by mario.red8976
-
Smoothing of methylation signal within windows
written by gamma.jian •
-
Integrating methylation in survival analysis
written by Pin.Bioinf
-
Comparison of Read Count per Genomic Position in R
written by ioannis •
-
Removing batch effects from methylation beta values for molecular subtyping/cluster analysis
written by alecmchiu •
-
Is it necessary to scale CpG data while using partitioning method and hierarchical clustering as a …
written by ivivek_ngs
You can actually divide the C counts by the total counts (C+T) to get methylation percentages. I have never used that particular function
percMethylation(), but you can check the 2 approaches and see if they return you the same values.Regarding the nomenclature, yes, "beta value" is array terminology, and I would not use it for bisulfite sequencing data. For sequencing data, in the literature you usually see a variety of terms such as "methylation value", "methylation percentage", etc. And people use both the [0,100] and [0,1] scales, so it is up to you (there's not much difference).