rmcelreath · June 18, 2021 18:39
diff --git a/RFDT1.R b/RFDT1.R
 set.seed(1908)
 N <- 200 # number of pairs
 U <- rnorm(N) # simulate confounds
 # birth order and family sizes
 B1 <- rbinom(N,size=1,prob=0.5) # 50% first borns
 M <- rnorm( N , 2*B1 + U )
 B2 <- rbinom(N,size=1,prob=0.5)
 D <- rnorm( N , 2*B2 + U + 0*M ) # change the 0 to turn on causal influence of mom
	set.seed(1908)
	N <- 200 # number of pairs
	U <- rnorm(N) # simulate confounds
	# birth order and family sizes
	B1 <- rbinom(N,size=1,prob=0.5) # 50% first borns
	M <- rnorm( N , 2*B1 + U )
	B2 <- rbinom(N,size=1,prob=0.5)
	D <- rnorm( N , 2B2 + U + 0M ) # change the 0 to turn on causal influence of mom
No results found