# structural = function(w1,w2,w3,gamma){
#   return (plogis(w1,scale=gamma) * plogis(w2,scale=gamma) * plogis(w3,scale=gamma))
# }

structural = function(delta,dist){
  return (delta^dist)
}