R/total_effect_mult_or.R

In efflog: The Causal Effects for a Causal Loglinear Model

total_effect_mult_or <-
function(x,y,z,w,t,q)
{ ((1+x*y)^(-1)*x*y*(1+w*t)^(-1)+
(1+x*y*z*q)^(-1)*x*y*z*q* (1+w*t)^(-1)*w*t)/ (1-((1+x*y)^(-1)*x*y*(1+w*t)^(-1)+
(1+x*y*z*q)^(-1)*x*y*z*q* (1+w*t)^(-1)*w*t)) *
((1+x)^(-1)*(1+w)^(-1)+ (1+x*z)^(-1)*(1+w)^(-1)*w)*(1-(1+x)^(-1)*(1+w)^(-1)- (1+x*z)^(-1)*(1+w)^(-1)*w)^(-1)}