Hi,
I am trying to generate a metabolic syndrome score for the children aged 11 to 12 years using PCA. I run the PCA with the standardised z-score of all five metabolic components But I confused on one statement mentioned in the article that I am following "Our PCA identified
2 principal components that were summed, with weights determined by the relative amount of variance, explained to generate a total MetSrisk score".
I am also missing Stata command for final analysis " we internally standardise the final met score to get mean 0 and SD 1 so that regression coefficient can be explained by SMD
I USED THE FOLLOWING COMMAND. It will be highly ap[appreciated if you can help me out with missing command lines. I am badly stuck at this point.
global xlist Z2TotChol Z2Glucose Z2HDL Z2TG fcbmizc
global ncomp 2
** PCA
pca $xlist
screeplot, yline (1)
** PCA Component
pca $xlist, comp($ncomp)
pca $xlist, comp($ncomp) blanks(0.3)
rotate, varimax
rotate, varimax blanks(0.3)
rotate, promax
rotate, promax blanks(0.3)
estat loadings
estat residuals, fitted
estat loadings, cnorm(eigen)
predict Pc1 Pc2, score
gen PCAScore=Pc1 + Pc2
Thanks in advance
Related Posts with adding summary weight in metabolic score
using * in a for-loopHello, I can't figure out what does * mean in front of a for-loop. Here you can see an example: fo…
Identifying peaks and throughsDear all, By following the method of Bry and Boschan (1971), full reference placed below, I want to…
Out-of-sample logit: applying estimated coefficients on a new series (w/o observing 0s/1s)Hello STATA community, This is both a methodological question and a programming question, for which…
Identifying peaks and throughs within a time seriesDear all, By following the method of Bry and Boschan (1971), full reference placed below, I want to…
Calculate pairwise differencesHi friends, I have data: Code: * Example generated by -dataex-. To install: ssc install dataex cl…
Subscribe to:
Post Comments (Atom)
0 Response to adding summary weight in metabolic score
Post a Comment