I am performing a principle component analysis with the between subject variable of fall status (non faller, single faller, recurrent faller) and within subject variable of sex to determine whether fall risk differs by fall status and sex. I would like to control for age and exposure to the exercise intervention in the principle component analysis, but I cannot find any code on how to go about doing this. The following pca code works, but does not control for covariates:

pca BMI_n OverallFallsRiskScore VisualContrast Proprio ReactionTime BestQuad_n EO_Firm_Sway_n EO_Foam_Sway_n Dual_TUG NIA_TotalScore GDS shortFES VAS FCI BarthelTotal_score IADL MMSE MOCA TrailA_n TrailB_n Stroop1_n Stroop2_n Stroop3_n DigitForward DigitBackward DSST_n VerbalFluency if NonSingleMultipleFaller > 1& Sex==0

This is a secondary analysis of RCT data.

Thanks for your help!