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) library(semTools) library(semPlot) M1 <- ' HealthyLife =~ Q22 + Q23 AvoidChem =~ Q15 + Q18 + Q20 + Q21 Saver =~ Q11 + Q13 + Q16 Recycle =~ Q01 + Q03 AvoidWaste =~ Q05 + Q06 EcoP =~ HealthyLife + AvoidChem + Saver + Recycle + AvoidWaste' S1 <- sem(model=M1, data=ecodata, estimator="ML") summary(object=S1, fit.measure=TRUE) reliability(S1) 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", "EcoPoint") cairo_pdf("ecolavaan2.pdf", width=12, height=12) semPaths(S1, what="stand", layout="tree", 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()