Puberty Measures

tanner

7T

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)
...