# \donttest{
library(sfa)
cs_data_trial <- data_gen_cs(N= 1000, rand = 1, sig_u = 0.3, sig_v = 0.3,
cons = 0.5, beta1 = 0.5, beta2 = 0.5, a = 4, mu = 1)
cs.nhnz <- sfm(formula = y_pcs_z ~ x1 +x2| z, model_name = "NHN",
data = cs_data_trial, PSopt = TRUE)
summary(cs.nhnz)
# }
Run the code above in your browser using DataLab