modkit icon indicating copy to clipboard operation
modkit copied to clipboard

I cannot identify the DMR output file column names

Open cj528-cam opened this issue 4 months ago • 2 comments

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?

cj528-cam avatar Aug 18 '25 14:08 cj528-cam

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.

ArtRand avatar Aug 18 '25 15:08 ArtRand

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.

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.

cj528-cam avatar Aug 18 '25 15:08 cj528-cam