This repository contains code for external validation of a microbiome-based machine learning model for Rheumatoid Arthritis (RA).
The model was originally trained on a large 16S rRNA dataset and evaluated on an independent dataset to assess generalizability.
- Raw FASTQ files processed using DADA2
- Amplicon Sequence Variants (ASVs) generated
- Taxonomy assigned using SILVA (v138.1)
- Features aggregated to genus level
- Prediction performed using a pre-trained XGBoost model
- Model performance evaluated using ROC-AUC and confusion matrix
- Download FASTQ files (external dataset)
- Perform quality filtering and trimming
- Learn error rates and infer ASVs
- Merge paired reads and remove chimeras
- Assign taxonomy (SILVA database)
- Aggregate ASVs to genus level
- Align features with training dataset
- Apply trained model
- Evaluate performance (ROC-AUC)
- R (version ≥ 4.3)
- Required packages:
- dada2
- dplyr
- stringr
- tibble
- pROC
- xgboost
- FASTQ files (external validation dataset)
- SILVA database (v138.1)
- Pre-trained model file (
model_xgb.model) - Training feature list (
train_asvs.rds)
- ROC-AUC value
- Confusion matrix
Sun Y et al. (2019)
Characteristics of gut microbiota in patients with rheumatoid arthritis in Shanghai, China
Front Cell Infect Microbiol. 2019;9:369.
Sabira Dabeer
This repository is provided for research and reproducibility purposes.