====== Puberty Measures ======

[[:tools:tanner]] 

===== 7T =====
For [[:grants:7t]], collected during [[:grants:7t:survey]] and included in [[:grants:7t:data|7T Data: merge_7t.csv]]


''/Volumes/Hera/Projects/7TBrainMech/scripts/behave/Makefile'' has the recipe

<code=makefile>
txt/puberty.csv: svys.RData
        ./puberty_from_rdata.R
</code>

where ''puberty_from_rdata.R'' is also [[https://github.com/LabNeuroCogDevel/allscripts/blob/master/root/Volumes/Hera/Projects/7TBrainMech/scripts/behave/puberty_from_rdata.R|here]]

Orma's composite score generating script is
''/Volumes/Hera/Orma/CogAnti_Puberty/scripts/clean_puberty.R'' (also [[https://github.com/LabNeuroCogDevel/allscripts/blob/master/root/Volumes/Hera/Orma/CogAnti_Puberty/scripts/clean_puberty.R#L161|here]])

An excerpt: 
<code=R>
# 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)
...
</code>