ococrook / hdxstats

hdxstats: An R-package for statistical analysis of hydrogen deuterium exchange mass-spectrometry data.
Apache License 2.0
5 stars 2 forks source link

Deuterium uptake plotting - linear vs semi-log #3

Closed nlgittens closed 2 years ago

nlgittens commented 2 years ago

Regarding the peptide uptake plots, it's very difficult to interpret anything visually from the early timepoints because the time dimension is linear in all of the plots. Therefore, the uptake plots look very similar, which may be deceptive since in a semi-log plot, the differences in deuterium uptake might become more apparent. How does the same data look plotted as a semi-log plot (I didn't find an easy way to do this in R with the data in the HdxStatModel class), and is there a rationale for the linear plot (perhaps to start at t = 0, D uptake = 0)?

Perhaps it's easier to take a figure from the paper.

Figure 4a:- "Figure 4 panel d shows the forest plot corresponding to the peptide in Figure panel a. Here, the difference is visually subtle but our functional method identifies a difference between the two samples"

I think Figure 4a it's just very difficult to see what is going on at the early timepoints conclusively whether there is a visual difference or not. image

ococrook commented 2 years ago

I'll see if I can do this easily change the the time axis to log coordinate. I think it should be easy.

Though there is warning here than the log transform causes problems at 0 and can exessively exagerate changes. I usually prefer the original coordinates unless there is a good reason

ococrook commented 2 years ago

@nlgittens There is now a way to do this in the vignette

ococrook commented 2 years ago

If you're happy could you close the issue?