search

dayoonkwon / BioAge

Biological Age Calculations Using Several Biomarker Algorithms
GNU General Public License v3.0
72 stars 24 forks source link

kdm0 and phenoage0 #24

Closed qu515 closed 6 months ago

qu515 commented 1 year ago

Hello Dayoon, very sorry to bother you, but I'm trying to replicate the kdm0 and phenoage0 values from the NHANES4 file that you provided, which I know are the original kdm/phenoage by Levine. For kdm, I've tried 

kdm = kdm_nhanes(biomarkers=c("fev","sbp","totchol","hba1c","albumin","creat","lncrp","alp","bun"))

and also

train = kdm_calc(NHANES3, biomarkers = c("fev","sbp","totchol","hba1c","albumin","creat","lncrp","alp","bun"))
kdm = kdm_calc(NHANES4, biomarkers = c("fev","sbp","totchol","hba1c","albumin","creat","lncrp","alp","bun"), fit = train$fit, s_ba2 = train$fit$s_ba2)

but the results are not the same as what you calculated. May I ask what's the code you implemented to calculate the kdm0 and phenoage0 values? Thank you again for all your replies!

dayoonkwon commented 6 months ago

Replicating the exact values using those functions might not be directly achievable due to the differences in the datasets. kdm0 and phenoage0 variables were computed based on the published coefficients by Levine and colleagues (Levine et al, 2013 and 2018).