Earnings inequality

inequality.R

n <- 100000
pe <- .2 #patriarchy benefit to male ability
ba <- rnorm(n) #baseline ability
m <- rbinom(n,1,.5) #male
h <- ba + pe * m + rnorm(n)#h index
e1 <- ba + pe *m + rnorm(n)#earnings assuming caused by ability
e2 <- h + rnorm(n)#earnings assuming caused by h-index

lm(e1~h+m)
lm(e2~h+m)