I am analyzing the effect of foreign aid on democracy and I would like to test following hypothesis (I simplified the original version):
- The amount of aid the country received and the corruption level determine the inequality level of the country (inequality is measured by Gini coefficient).
- The inequality level (affected by aid) is related to the democracy level.
Firstly, I simply did two OLS analyses. ($i$ is the index for countries, $t$ is for years)
- ${\rm Inequality}_{i,t} = \beta_0 + \beta_1 {\rm Aid}_{i,t} + \beta_2 {\rm Corruption}_{i,t} + \varepsilon_{i,t} $
- $\textrm{Democracy}_{i,t} = \beta_3 + \beta_4 \textrm{Democracy}_{i,t-1} + \beta_5 {\rm Inequality}_{i,t} + \varepsilon_{i,t}$
We can say "Aid is related to Inequality" and "Inequality is related to Democracy" (let me suppose all coefficients are positive and statistically significant), but cannot say "how much is aid related to democracy" or "how much does Aid affect the level of democracy".
What kind of technique should I use to model for my purpose?
It seems 2SLS is similar to what I want to do, but since variables in the first equation are not independent from $\textrm{Democracy}_{i,t}$, I think 2SLS is inappropriate here.
I am familiar with both R and Rstan, so either using R packages or building a model with Stan would be fine.
Thank you.