I cannot identify the DMR output file column names
Hello Everyone,
I'd like to ask if you have got more number of columns that it listed in the table when use default thresholds?
For example, I put
dmr_result="location"
echo "[3/4] Running modkit dmr pair (genome-wide)"
modkit dmr pair
-a
-b
-o ${dmr_result}
--ref
--base C
--threads
--log-filepath /location \
The output is (first line): chr1 16242 16243 . -0.5035263212843841 + h:0,m:1 1 h:0,m:5 5 h:0.00,m:100.00 h:0.00,m:100.00 1 1 1 0 0 -2.1470329724605883 2.1470329724605883
Could someone help me with identifying what does each column represent? Especially last 10 columns please?
Hello @cj528-cam, there is a description of the table schema in the documentation. You can also run modkit dmr pair with the --header flag to get a BED-compatible header line which can make identifying the values easier.
Hello @cj528-cam, there is a description of the table schema in the documentation. You can also run
modkit dmr pairwith the--headerflag to get a BED-compatible header line which can make identifying the values easier.
Thank you. I checked the table in the documentation. But it only has 17 columns where the last three were Cohen's h statistic that only shows with --regions option (which I don't have).
Theoretically, I should have 14 columns based on the documentation, but now I have 19 columns. I will do the --header flag and see if it makes things easier.