Skip to content
master
Switch branches/tags
Code

Latest commit

 

Git stats

Files

Permalink
Failed to load latest commit information.
Type
Name
Latest commit message
Commit time
Aug 25, 2020

PharmIV

This respository supplies the tools necessary to calculate the power of an instrumental variable analysis study using a single binary instrument Z to analyse the causal effect of a binary exposure X on a continuous outcome Y in the context of pharmacoepidemiology. If you use this material, please cite:

Walker VM, Davies NM, Windmeijer F, Burgess S, Martin RM. Power calculator for instrumental variable analysis in pharmacoepidemiology. Int J Epidemiol. 2017 Oct 1;46(5):1627–32.

The simulations and plots included in this paper were produced using Stata and the code can be found in the folder 'Stata_paper'.

Shiny

The folder 'Shiny' contains the code to produce the power calculator shiny app. The app is avaliable at:

https://venexia.shinyapps.io/PharmIV/

R package

The folder 'R_package' contains the power calculator R package.

Installation

The package can be installed using the following commands:

install.packages("devtools")
devtools::install_github("venexia/PharmIV", subdir="R_package")

Syntax

The syntax for this package is as follows:

PharmIV(n = NULL, delta = NULL, alpha = 0.05, sigma = 1.00, prob_x1 = NULL, prob_z1 = NULL, cond_z1 = NULL, cond_z0 = NULL)

Example

What is the power of an instrumental variable analysis study with 10,000 participants to detect a treament effect of 0.20 for a treatment with a frequency of exposure of 0.50? The binary instrument for the study has a freqency of 0.20 and the probability of exposure given the instrument Z=0 is 0.57.

alt text

Stata package

To allow web download, the power calculator Stata package is contained in its own repository:

https://github.com/venexia/PharmIV-Stata

Installation

The package can be installed using the following commands:

net install github, from("https://haghish.github.io/github/")
github install venexia/PharmIV-Stata

Syntax

The syntax for this package is as follows:

PharmIV, n(numlist) delta(numlist) alpha(numlist) sigma(numlist) prob_x1(numlist) prob_z1(numlist) cond_z1(numlist) cond_z0(numlist) 

Example

What is the power of an instrumental variable analysis study with 10,000 participants to detect a treament effect of 0.20 for a treatment with a frequency of exposure of 0.50? The binary instrument for the study has a freqency of 0.20 and the probability of exposure given the instrument Z=0 is 0.57.

alt text

About

Power calculator for instrumental variable analysis in pharmacoepidemiology

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published