datagen_xx<-function(n,seed=sample(1:100000,size=1)){ set.seed(seed) In0=rnorm(n, 0, 1) L0=rnorm(n, 0, 1) R0=rnorm(n, 0, 1) In1=rnorm(n, 0, 1) #A1 p_A1=plogis(In0+L0) A1=rbinom(n, 1, p_A1) #R1 R1=rnorm(n,0.5*R0,1) #L1 m_L1=0.5*L0+0.5*A1 L1=rnorm(n,mean=m_L1,sd=1) #A2 p_A2=plogis(A1+In1+In0+L1) A2=rbinom(n,1,p_A2) #Y p_Y=plogis(-2+0.5*A1+0.5*A2-L0-L1+0.5*R0+0.5*R1) Y=rbinom(n, 1, p_Y) truth<- 0.2858638 # X=data.frame(In0, L0, In1, R0, A1, R1, L1 ,A2,Y) return(X) }