For Brain Mechanisms R01 (7T), collected during 7T BrainMech Qualtrics Survey and included in 7T Data: merge_7t.csv
/Volumes/Hera/Projects/7TBrainMech/scripts/behave/Makefile has the recipe
txt/puberty.csv: svys.RData
./puberty_from_rdata.R
where puberty_from_rdata.R is also here
Orma's composite score generating script is
/Volumes/Hera/Orma/CogAnti_Puberty/scripts/clean_puberty.R (also here)
An excerpt:
# adjust gonadal score for boys
pub$gonadm2 <- pub$gonadm
pub$gonadm2[pub$gonadm==1 & pub$mpete==1] <- 1
pub$gonadm2[pub$gonadm==1 & pub$mpete > 1] <- 2
...
print("Combine girl and boy variables to create 5-point PDS Adrenal Score (PDSA), Gonadal Score (PDSG), and Total Score (PDSS)")
pub$PDSS <- rowMeans(pub[,c('gonadm2', 'adrenm2', 'gonadf2', 'adrenf2')], na.rm=TRUE)
...