R/ndirect_effect_mult_or.R

In efflog: The Causal Effects for a Causal Loglinear Model

ndirect_effect_mult_or <-
function(x,y,z,w,t,q)
{((1+x*y*z*q)^(-1)*x*y*z*q*(1+w)^(-1)*w+ (1+x*y)^(-1)*x*y*(1+w)^(-1))*
(1-(1+x*y*z*q)^(-1)*x*y*z*q*(1+w)^(-1)*w- (1+x*y)^(-1)*x*y*(1+w)^(-1))^(-1)*
 ((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)}