Adding estimated group methylation levels and dispersion to output - #287

Merged

andrewdavidsmith merged 8 commits into

radmeth-output-format

Aug 18, 2025

Commits on Aug 18, 2025

src/radmeth/radmeth_optimize.hpp: no longer return bool from fit_regression_model as it was never checked; throwing on error instead
andrewdavidsmith
committed
src/radmeth/radmeth_optimize.cpp: fit regression model now no longer accepts initial parameter values as those weren't used. Out-params have been added to return the estimates of p for the beta-bin…
andrewdavidsmith
committed
src/radmeth/radmeth_model.cpp: updated the function to print a model so that it only prints each sample's factor levels if there are few enough samples and so it prints a summary of group factor le…
andrewdavidsmith
committed
src/radmeth/radmeth_model.hpp: added accessor for number of parameters
andrewdavidsmith
committed
src/radmeth/radmeth.cpp: Chaned what is output so that for each group, an estimate of the methylation level is provided for each site and an estimate of the shared dispersion parameter. Also remove…
andrewdavidsmith
committed
src/radmeth/radmeth.cpp: removing unused variable
andrewdavidsmith
committed
src/radmeth/radmeth_model.cpp: changed max number of samples to report factor levels for when printing a design
andrewdavidsmith
committed
src/radmeth/radmeth_optimize.cpp: not throwing an exception when the iterative estimation fails to converge; there are ok reasons for not converging
andrewdavidsmith
committed