!@front

#define nvar 3			! n dependent variables per individual

G1: Singleton (non-pair) data 		
Data Ninput_vars=5 NGroups=2	 	! Number of variables per family

Rectangular file=SarahData.txt		! read raw data
Labels	Id depression anxiety stress sex

	
select depression anxiety stress ;

Begin matrices; 
    M Full 1 nvar free 		! mean
    L Full nvar 1 free		! factor loadings
    R Diag nvar nvar free	! Residual Variance
End matrices;
    
Begin Algebra;
    C=L*L'+R*R' ;
End Algebra;

! start values
Start 0 M 1 1 M 1 2 M 1 3 
Start .5 L 1 1 L 2 1 L 3 1 
Start .5 R 1 1 R 2 2 R 3 3 

Means M ;	! means model 
Covariances C ;		! variance/covariance model
end

Standardize
Calc 

Begin Matrices = Group 1 ;
End Matrices ;

Begin Algebra;
 V = \sqrt(\d2v(C)') ;  	!compute the standard deviation
 S = L%V | (\d2v(R)')%V ;  	!compute the standardised factor loadings
End Algebra;
End