search for: d1_ave

...$loadings[,ncomp] Q <- MyModel$Yloadings[,ncomp] # D1 f1 <- P["D1a"] f2 <- P["D1b"] f3 <- P["D1c"] Sp <- f1 + f2 + f3 Sp2 <- (f1 ^ 2) + (f2^ 2) + (f3^2) Sth <- (1-(f1 ^ 2)) + (1-(f2 ^ 2)) + (1-(f3^2)) D1_ICR <- (Sp^2) / ( (Sp^2) + Sth) D1_AVE <- Sp2 / ( Sp2 + Sth) but the results does not seem to give me something meaningfull. For example, while cronbach(cbind(D1a,D1b,D1c)) gives me > 0.90, the above computed D1_ICR gives me very low numbers (< .20). Also summary says % variance explained for X = 83.1 in 1st component wh...