Bookmark and Share

Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

st: CR and AVE for factor analysis with 2 factors

From   "Florian Christian Esser" <>
Subject   st: CR and AVE for factor analysis with 2 factors
Date   Mon, 28 Oct 2013 15:11:44 +0100

Hi everyone,

I am trying to do factor analysis in order to measure two strategy
constructs. I have 6 indicators that load on these two constructs. Now I
want to calculate composite reliability (CR) and average variance
I use the following code (here for CR):
factor zwmdiff_1 zwpdiff_1 zwpdiff_2 zwcost_1 zwcost_2 zwcost_3, pcf
matrix list e(Psi)
matrix list e(L)
matrix psi = e(Psi)'
matrix communalities = J(rowsof(psi),1,1)
matrix communalities = communalities - psi
matrix colnames communalities = communalities
matrix list communalities
gen nvar = e(df_m)
gen sum_factors=0
local i=1
while `i' <= nvar {
replace sum_factors = sum_factors + sqrt(communalities[`i',1])
local i=`i'+1
generate sum_psi=0
local i=1
while `i' <= nvar {
replace sum_psi = sum_psi + psi[`i',1]
local i=`i'+1
gen cr=(sum_factors*sum_factors)/((sum_factors*sum_factors)+sum_psi)
drop nvar
drop sum_factors
drop sum_psi
list cr in 1

The problem is that once the local macro starts running it tells me:
(534 real changes made)
(534 real changes made)
(534 real changes made, 534 to missing)
(0 real changes made)
(0 real changes made)
(0 real changes made)

and accordingly there is an empty result for CR.
If I add "factors(1)" in the command line above, I get a result, but I
since I have two factors, I think it is not correct to do this.

Does anyone know what to do here?

Thanks a lot in advance.

*   For searches and help try:

© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index