added male vs female test on pds, age

nmuncy · Nov 30, 2021 · ae3a53a · ae3a53a
1 parent 1d4cb3d
commit ae3a53a
Showing 1 changed file with 12 additions and 0 deletions.
diff --git a/afq_step4_manuscript.R b/afq_step4_manuscript.R
@@ -32,3 +32,15 @@ etaSquared(stat_age)
 stat_pds <- aov(PDS ~ as.factor(Group), data = df_subset)
 summary(stat_pds)
 etaSquared(stat_pds)
+
+
+# Male vs Female PDS, age -----
+pds_female <- df_subset[which(df_subset$sex == 0),]$pds
+pds_male <- df_subset[which(df_subset$sex == 1),]$pds
+t.test(pds_female, pds_male, paired=FALSE)  
+
+age_female <- df_subset[which(df_subset$sex == 0),]$age
+age_male <- df_subset[which(df_subset$sex == 1),]$age
+t.test(age_female, age_male, paired=FALSE) 
+
+