Description Usage Arguments Details Value Author(s) References Examples
Computes the power of sensitivity analysis, which is based on an extension of Anderson-Rubin test and allows IV be possibly invalid within a certain range.
1 2 | ARsensitivity.power(n,k,lambda,gamma,var.z,sigma1,sigma2,rho,alpha=.05, Delta=NULL,
delta=NULL)
|
n |
Sample size. |
k |
Number of exogenous variables, intercept is also counted. If there's no observed variable, k=1. |
lambda |
True causal effect minus null hypothesis causal effect |
gamma |
Regression coefficient for effect of instrument on treatment. |
var.z |
Variance of instrument. |
sigma1 |
Standard deviation of potential outcome under control (structural error for y) |
sigma2 |
Standard deviation of error from regressing treatment on instrument |
rho |
Correlation between u(potential outcome under control ) and v (error from regressing treatment on instrument) |
alpha |
Significance level of test. |
Delta |
Range of sensitivity allowance. Usually a numeric vector of length 2. If Delta=NULL, the ARsensitivity.power function will calculate the power for standard Anderson-Rubin test with valid IV. |
delta |
True value of sensitivity parameter when calculating the power. Usual take delta=0 for the favourable situation or delta=NULL for unknown delta. |
The structural equations model assumed is: D=gamma0+gamma*z+v, Y=beta0+beta1*D+delta*Z+u. This model can also be obtained by assuming the potential outcomes model. See Jiang, Zhang and Small (2013) for details.
lambda is equal to the true beta1 minus the null hypothesis beta1.
Power of sensitivity analysis for the proposed study, which extends the Anderson-Rubin (1949) test with possibly invalid IV. The power formula is derived in Jiang, Small and Zhang (2013).
Dylan Small
Anderson, T.W. and Rubin, H. (1949), Estimation of the parameters of a single equation in a complete system of stochastic equations, Annals of Mathematical Statistics, 20, 46-63.
Jiang, Y., Zhang, N and Small, D. (2013), Sensitivity analysis and power for instrumental variable studies, Working paper.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 | ### Power for a study with in which the null hypothesis causal effect is 0,
### the true causal effect is 1, the sample size is 250, the instrument is
### binary with probability .5 (so variance = .25), the standard deviation
### of potential outcome under control is 1, the effect of the instrument
### is to increase the probability of a binary treatment being 1 from .25 to
### .75. The function sigmav.func computes the SD of v for a binary instrument,
### binary treatment. The correlation between u and v is assumed to be .5. The
### significance level for the study will be alpha = .05
sigmav.func(prob.d1.given.z1=.75,prob.d1.given.z0=.25,prob.z1=.5)
# The sigmav.func finds sigmav=.4330127
### power of Anderson-Rubin test
ARsensitivity.power(n=250, k=1, lambda=1, gamma=.5, var.z=.25, sigma1=1,
sigma2=.4330127, rho=.5, alpha = 0.05)
### power of sensitivity analysis under the favourable situation, assuming
### the range of sensitivity allowance is (-0.1, 0.1)
ARsensitivity.power(n=250, k=1, lambda=1, gamma=.5, var.z=.25, sigma1=1,
sigma2=.4330127, rho=.5, alpha = 0.05, Delta=c(-0.1, 0.1), delta=0)
### power of sensitivity analysis with unknown delta, assuming the range of sensitivity
### allowance is (-0.1, 0.1)
ARsensitivity.power(n=250, k=1, lambda=1, gamma=.5, var.z=.25, sigma1=1,
sigma2=.4330127, rho=.5, alpha = 0.05, Delta=c(-0.1, 0.1))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.