Only include IDs where sex could be predicted

mpip · Mar 7, 2025 · 83cdc3c · 83cdc3c
1 parent a92d462
commit 83cdc3c
Showing 1 changed file with 1 addition and 0 deletions.
diff --git 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) %>%