eco <- read.delim("ecopx.txt") ecodata <- eco[, c(4, 6, 8, 9, 14, 16, 19, 21, 23:26)] C1 <- cor(ecodata) library(sem) M1 <- specifyModel(text=" Q22 <- HealthyLife, a1 Q23 <- HealthyLife, a2 Q21 <- HealthyLife, a3 Q06 <- HealthyLife, a4 Q18 <- AvoidChem, b1 Q20 <- AvoidChem, b2 Q21 <- AvoidChem, b3 Q16 <- AvoidChem, b4 Q11 <- Saver, c1 Q13 <- Saver, c2 Q16 <- Saver, c3 Q01 <- Recycle, d1 Q03 <- Recycle, d2 Q16 <- Recycle, d3 Q05 <- AvoidWaste, e1 Q06 <- AvoidWaste, e2 Q13 <- AvoidWaste, e3 HealthyLife <-> HealthyLife, NA, 1 AvoidChem <-> AvoidChem, NA, 1 Saver <-> Saver, NA, 1 Recycle <-> Recycle, NA, 1 AvoidWaste <-> AvoidWaste, NA, 1 ") S1 <- sem(M1, C1, N=length(ecodata[, 1])) print(S1) summary(S1, fit.indices=c("GFI","AGFI","CFI","RMSEA")) library(semPlot) LBL <- c("Q22","Q23","Q21","Q18","Q20","Q06","Q16","Q11","Q01","Q03","Q13","Q05", "Healthy\n Life","Avoid\n Chemical","Saver","Recycler","Avoid\n Waste") semPaths(S1, what="stand", layout="spring", style="lisrel", shapeMan="rectangle", shapeLat="ellipse", sizeMan=3, residScale=9, posCol="black", negCol="red", fade=FALSE, edge.label.cex=0.8, nodeLabels=LBL) # For submission, negCol also shoule be "black"