Function to compute CVD-score made by Peter Wurtz et al. made by Deelen et al. on Nightingale metabolomics data-set.
comp.CVD_score(met, phen, betas, quiet =FALSE)
Arguments
met: numeric data-frame with Nightingale-metabolomics
phen: data-frame containing phenotypic information of the samples (specifically: sex, systolic_blood_pressure, current_smoking, diabetes, blood_pressure_lowering_med, lipidmed, totchol, and hdlchol)
betas: The betas of the linear regression composing the CVD-score
quiet: logical to suppress the messages in the console
Returns
data-frame containing the value of the CVD-score on the uploaded data-set
Examples
library(MiMIR)#load the datasetmet <- synthetic_metabolic_dataset
phen<-synthetic_phenotypic_dataset
#Prepare the metabolic features fo the mortality scoreCVDscore<-comp.CVD_score(met= met, phen=phen, betas=MiMIR::CVD_score_betas, quiet=TRUE)
References
This function is constructed to be able to apply the CVD-score as described in: Wurtz,P. et al. (2015) Metabolite profiling and cardiovascular event risk: a prospective study of 3 population-based cohorts. Circulation, 131, 774-785, doi:10.1161/CIRCULATIONAHA.114.013116