diff --git a/scripts/4_check_matching_of_epigenetic_sex.Rmd b/scripts/4_check_matching_of_epigenetic_sex.Rmd index 736dafd..bf8ac07 100644 --- a/scripts/4_check_matching_of_epigenetic_sex.Rmd +++ b/scripts/4_check_matching_of_epigenetic_sex.Rmd @@ -54,6 +54,7 @@ fail the sex match check ```{r, predict sex with DNAm data, include= FALSE} predictedSex <- getSex(gRatioSet_quality, cutoff=-2) +phenotype_data <- subset(phenotype_data, arrayid %in% rownames(predictedSex),) sex_df <- phenotype_data %>% select(personid, arrayid, reported_sex = sex) %>%