Error in calcDNAmTL resulting in wrong estimates for the DNAmTL clock

[ausrinesv]

There seems to be an error in calcDNAmTL().

The code subtracts the intercept from the product of beta values and coefficients:

methylCIPHER/R/calcDNAmTL.R

Line 43 in e62937e

DNAmTL <- as.numeric(rowSums(tt,na.rm=T)-7.924780053)

However, in the related Supplementary Data 1 of the Lu, A. T. et al. (2019) DNA methylation-based estimator of telomere length that intercept is positive.