eco <- read.delim("http://minato.sip21c.org/advanced-statistics/ecopx.txt") ecodata <- eco[, c(1, 3, 5, 6, 11, 13, 15, 16, 18, 20:23)+3] ecodata <- subset(ecodata, complete.cases(ecodata)) library(lavaan) M1 <- ' HealthyLife =~ Q22 + Q23 AvoidChem =~ Q15 + Q18 + Q20 + Q21 Saver =~ Q11 + Q13 + Q16 Recycle =~ Q01 + Q03 AvoidWaste =~ Q05 + Q06' S1 <- sem(model=M1, data=ecodata, estimator="ML") summary(object=S1, fit.measure=TRUE) library(semTools) reliability(S1) library(semPlot) LBL <- c("Q22","Q23","Q15","Q18","Q20","Q21","Q11","Q13","Q16", "Q01","Q03","Q05","Q06", "Healthy\n Life","Avoid\n Chemical","Saver","Recycle", "Avoid\n Waste") cairo_pdf("ecolavaan.pdf", width=12, height=12) semPaths(S1, what="stand", layout="circle", style="lisrel", shapeMan="rectangle", shapeLat="ellipse", sizeMan=3, residScale=9, posCol="black", negCol="red", fade=FALSE, edge.label.cex=0.8, nodeLabels=LBL) dev.off()