if (FALSE) {
# Load acrylic data set
data(acrylic)
# Build a netSEMp2 model
ans2 <- netSEMp2(acrylic, criterion = "AIC")
ans2_BIC <- netSEMp2(acrylic, criterion = "BIC")
# Plot the network model
plot(ans2, cutoff = c(0.3,0.6,0.8))
plot(ans2_BIC , cutoff = c(0.3,0.6,0.8))
# Drop Relationship lower than minimum cutoff
plot(ans2, cutoff = c(0.3,0.6,0.8), style = FALSE)
# plot network model with latent argument labels
plot(ans2, cutoff = c(0.3, 0.6, 0.8),
latent = c('IAD1' = 'FundAbsEdge',
'IAD2' = 'UVStab',
'IAD2p' = 'UVStab',
'IAD3' = 'YelMet'))
# plot the network diagram and save file
#plot(ans, cutoff = c(0.3,0.6,0.8), plot.save = TRUE, filename = "acrylic-netSEMp2"))
}
Run the code above in your browser using DataLab