## ----message=FALSE, eval=FALSE------------------------------------------------
# 
# ## Load the pv module data set
# data(PVmodule)
# 
# ## Run netSEMp1 model
# ans1 <- netSEMp1(PVmodule, exogenous = "Time", endogenous = "Pmax")
# #plot network model with latent variable
# plot(ans1, style = F, cutoff = c(0.3, 0.6, 0.9),
#      latent=c('EVA_hyd', 'Hac'))
# 
# ## Run netSEMp2 model
# ans2 <- netSEMp2(PVmodule, exogenous = "Time", endogenous = "Pmax")
# ## Plot the network model for principle 2
# plot(ans2)

## ----out.width="675px", echo=FALSE, fig.cap="PVmodule netSEMp1 model"---------
knitr::include_graphics("PV1.png")

## ----out.width="675px", echo=FALSE, fig.cap="PVmodule netSEMp2 model"---------
knitr::include_graphics("PV2.png")