@ausrinesv detected that currently Weidner clock returns wrong results because of an error in CpG matching as well as possibly wrong coefficients.
The error is in this line:
|
Weidner <- 111.83 + rowSums(DNAm[,na.omit(match(colnames(DNAm), Weidner_CpGs))] * c(-64.57, -42.57, + 75.15)[present]) |
Where colnames(DNAm) and Weidner_CpGs should be switched places in the match() function call.
As written the line always selects first 3 cytosines from the data.
In addition the coefficients provided in the methods section of the paper are different from the ones assigned in the code.