Revision f9477d6b8382fee6926be7fd29b99afc15873fe8 authored by John Lees on 05 July 2022, 11:05:34 UTC, committed by John Lees on 05 July 2022, 11:05:34 UTC
1 parent 83396ac
roc_plot.R
library(pROC)
setwd("Documents/Imperial/Papers/carriage_age/analysis")
nl_self_predict <- read.delim("nl_self_predict.txt")
nl.nonmissing <- read.delim("nl.nonmissing.pheno")
nl_merge <- merge(nl.nonmissing, nl_self_predict, by.x = "Lane", by.y = "Sample")
maela_self_predict <- read.delim("maela_self_predict.txt")
maela_pheno <- read.delim("maela.pheno")
maela_merge <- merge(maela_pheno, maela_self_predict, by.x = "samples", by.y = "Sample")
roc(nl_merge, "Age_month_dichot", "Link", plot = TRUE, col = '#56B4E9',
smooth = TRUE,
print.auc = TRUE, print.auc.x = 0.45, print.auc.y = 0.03, print.auc.cex = 0.8,
print.auc.pattern = "Dutch cohort (AUC: %0.2f)")
roc(maela_merge, "mother", "Link", plot = TRUE, add = TRUE, col = "#E69F00",
smooth = TRUE,
print.auc = TRUE, print.auc.x = 0.45, print.auc.y = 0.1, print.auc.cex = 0.8,
print.auc.pattern = "Maela cohort (AUC: %0.2f)")
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...