# NOT RUN {
  data = data.frame(
    a = runif(50),
    b = runif(50, 0, 1),
    c = runif(50, 0, 1),
    d = runif(50)
  )
  
  modelList = list(
    glm(b ~ a, family = "binomial", data),
    glm(c ~ a, family = "binomial", data),
    lm(d ~ b + c, data)
  )
  
  (basis.set = sem.basis.set(modelList))
  
  endogenous.reverse(basis.set, modelList)
# }
Run the code above in your browser using DataLab