**Statistical Reflections of a Medical Doctor » R**, and kindly contributed to R-bloggers)

I decided to start blogging the R code used for some of my statistical posts, so I will start with the meta-analysis posts and move on to more difficult stuff.

As stated previously (here and here) the problem is to convert the reported relative risks(RR, ), 95% confidence interval () and p-value () into estimates for the log-relative risk ratio and its associated standard error for down-stream use (usually meta-analysis). Medical journals are in the bad habit of exponentiating (and rounding) the output of statistical software so that one needs to manipulate the reported estimates in order to recover the output of the statistical software. The simplest approach is one based on algebraic manipulations and disregards the finite precision of the medical journal table, treating the reported figures as exact (i.e. 1.50 is *exactly *1.50 not a number between 1.495-1.504):

- First we obtain an estimate for from the reported RR:
- Secondly we transform the confidence interval to estimates for :
- Averaging over yields
- A second estimate for the is obtained from the p-value using the quantile function of the normal distribution:
- Averaging over yields

The following R function can be used to carry out these transformations:

algSTE<-function(X,L,U,p) { v1<-c(log(X),log(sqrt(U*L))) alg.X<-mean(v1[is.finite(v1)]) v2<-c(log(sqrt(U/L))/qnorm(0.975),abs(alg.X/qnorm(p/2))) alg.Y<-mean(v2[is.finite(v2)]) return(list("TE"=alg.X,"se"=alg.Y)) }

A potentially more accurate approach is to acknowledge the finite precision inherent in the medical journal table for all quantities of interest (treatment effect, 95% confidence intervals and p-values). In a previous post I gave a Bayesian solution to this problem which obeys the precision of the treatment effect and confidence interval and the precision of the p-value. The R code for the Bayesian algorithm is the following:

## function that extracts treatment effect (usually in log scale) and ## standard error from treatment effect (usually exponentiated log-odds ratios ## or log-relative risks), associated pci% CI and p-value rounded to ## d (TE, 95% CI) and dp (p-value) significant digits. ## This is based on a statistical simulation algorithm that draws ## N samples from the posterior distribution of these parameters ## followed by a rejection step. Due to the latter, the final number ## of samples can be <The implementation is based around two functions with the user only needing to call the first one. Note that for some combinations of input values, a large number of random numbers (up to Nmax*N) have to be generated in order for the Monte Carlo algorithm to retain a certain level of precision (controlled by the number of samples N).

Let’s test these two approaches from within R now in a toy example:

> d<-2 ## precision TE, 95% CI (all exponentiated) to > dp<-3 ## precision of p-value > mn<-0.2209039 ## true treatment effect (log scale) > se<-0.2761299 ## true SE > > X<-round(exp(mn),d);L<-round(exp(mn-qnorm(0.975)*se),d);U<-round(exp(mn+qnorm(0.975)*se),d) > p<-round(2*(1-pnorm(abs(mn/se))),dp) > > X ## rounded/exponentiated treatment effect [1] 1.25 > L ## lower 95% CI (rounded/exponentiated) [1] 0.73 > U ## upper 95% CI (rounded/exponentiated) [1] 2.14 > p [1] 0.424 > > moko<-extractTSE(X,L,U,p,d,dp,N=10000) > mean(moko$TE) ## Bayesian TE [1] 0.2204252 > mean(moko$se) ## Bayesian SE [1] 0.2757051 > > moko2<-algSTE(X,L,U,p) > moko2$TE ## algebraic TE (log scale) [1] 0.2230955 > moko2$se ## algebraic se [1] 0.2767075

The performance of the two approaches appears similar in this particular case!

However as we show here there will be occasions which the simple(-minded) algebraic approach will fail and the (big) Bayesian guns will be required. Posting the source code will ensure that the Bayesian procedure becomes available to everyone who needs to carry out a meta-analysis!!

**leave a comment**for the author, please follow the link and comment on their blog:

**Statistical Reflections of a Medical Doctor » R**.

R-bloggers.com offers

**daily e-mail updates**about R news and tutorials on topics such as: Data science, Big Data, R jobs, visualization (ggplot2, Boxplots, maps, animation), programming (RStudio, Sweave, LaTeX, SQL, Eclipse, git, hadoop, Web Scraping) statistics (regression, PCA, time series, trading) and more...