similar to: Create patterns within a plot?

Dear R-experts, The forest.rma() function from the metafor package creates nice forest plots for presenting the results of a meta-analysis. These plots can be annotated for e.g. giving names to the columns. E.g. as in the documentation of the package: data(dat.bcg) ### meta-analysis of the log relative risks using a random-effects model res <- rma(ai=tpos, bi=tneg, ci=cpos, di=cneg,

I am using metafor to create forest plots. This code gives me the expected plot (setting x axis tick marks): forest(forest$OR, ci.lb=forest$Low, ci.ub=forest$High, at=log(c(.05, .25, 1, 10)), slab=forest$SNP, atransf=exp) As does this (setting colors): forest(forest$OR, ci.lb=forest$Low, ci.ub=forest$High, col=c(1,2,3), slab=forest$SNP, atransf=exp) But if I try to set both 'at' and

Dear R-community, I'm currently trying to assemble a forest plot using the "forest" function from package "metaphor". Works well. Even the regular "main"-argument works for adding a title to the graph. However, I would like to add one top row which explains the nature of the columns. Very much like the usual header in spreadsheet programs. For example:

Dear R User, Please, I am new to R. I want to overlay density plot for predictive interval pooled result in meta-analysis. http://addictedtor.free.fr/graphiques/graphcode.php?graph=114 Regards Frank Peter

Hi everyone, I'm using the metafor package to meta-analyze a set of proportions. This is working really well for the raw proportions, but is there a way to back-transform the arcsine transformed proportions in the rma or forest functions with the atransf option? The estimates and CIs for the transformed proportions need to be back-transformed to be the sin of the estimate squared.

Hello, Is it possible to have two meta-plots in one graph (not par(mfrow=c(2,1))? But somthing like  library(metafor)  library("igraph")    if (interactive()) {     forest(dat.Treat$RR, ci.lb=dat.Treat$lower, ci.ub=dat.Treat$upper, xlab="Relative Risk",slab=dat.Treat$ID,refline=1)     forest(dat.Control$RR, ci.lb=dat.Control$lower, ci.ub=dat.Control$upper, xlab="Relative

Hi I have tryed to do the meta-regression in metafor package, but I would like to get the standardized coefficients for each variable, however in command:   Ø  res<-rma.uni (yi, vi, method="REML", mods=~cota+DL+uso+gadiente+idade, data= turbidez)   I just have the coefficients no standardized (estimate) of the multiple regression. What I need to do? Thanks Fernanda Melo

Dear R-experts! I am working on some meta-analysis using the metafor package. I would like to extract values of the confidence intervals of the effect sizes of the single studies from an rma object. Those values are printed out when plotting a forest plot using the forest function on the rma object, however I was not able to locate them. Many thanks for your help! Jokel [[alternative HTML

Dear R-experts! In a meta-analysis (metafor) I would like to assess the effect of two categorical covariates (A & B) whereas they both have 4 levels. Is my understanding correct that this would require to dummy-code (0,1) each level of each covariate (A & B)? However I am interested in the main-effects and the interaction of these two covariates and the dummy-coding would only allow to

Dear List, Please put me on the right track. I have been searching for the R functions for measures of agreement. Where should I look? Patrik Öhagen Läkemedelsverket Box 26 751 03 Uppsala Tel. +46(0)18-174924 patrik.ohagen@mpa.se [[alternative HTML version deleted]]

Dear R community, I'm using the ,metafor'-package by Wolfgang Viechtbauer (Version: 0.5-5) to calculate random-effects meta-analyses using Correlations and Sample Sizes as the raw data. (By the way: Really a nice piece of work, Wolfgang! Thanks heaps.) I specified the "rma.uni' function so that it looks like this: MAergebnis<-rma.uni(ri=PosOutc, ni=N,

Hello metafor users, I'm using metafor to perform a single-effect summary estimate of the raw proportion of patients experiencing a post-operative complication, and I'm interested in seeing if this proportion differs between the three most commonly used surgical techniques. The software is working as expected, but I would like to double check on the interpretation of my mixed-effect model

Hello, I am performing a meta-analysis using the metafor package. My data are proportions and I used the Freeman Tukey double arcine (FT) transformation to fit the random effects model. Now I want to create a forest plot with my estimates backtransformed to the original scale of proportions. Can this be done? Regards, Patricia

If I want to find out the soltion of X1,X2 that min(3X1+2X2+X1X2) subject to 20<=X1+3X2<=50 10<=X1 which function or package can I use? Thanks.

Hi All, 1) Is it possible to set the options such that R opens a new script editor every time I start the R and 2) specify the size of windows. Thanks for the suggestion and Best regards, Krishna [[alternative HTML version deleted]]

Hi I'm trying to wrap predict.lm within a function, but I'm having problems passing arguments into it in this way. Basically I want to create a lm object, then pass it into the predict.lm function and be able to tell predict.lm which variable I want to predict for, by passing the variable name as an argument of the wrapper function. This variable will always be a factor with 3 levels,

Hi! Anybody familiar with streaming rtp from asterisk. Preferably with the xorcom asterisk patch which streams rtp from asterisk to oreka audio server. Any ideas will do just fine though! Regards / Marcus

An r hel search picks up http://tolstoy.newcastle.edu.au/R/help/01a/1162.html from 2001, which points towards the ts package. >>> ?hagen Patrik <Patrik.Ohagen at mpa.se> 19/08/2008 09:00:12 >>> Dear List, I have used all my resources (i.e. "help.search) and I still havn't been able to figure out if there is an Exponential Smoothing command in R. Thank you in

Hello again, let say we have following data: Dat1 <- structure(list(factor.sample.LETTERS.1.3...6..replace...T.. = structure(c(1L, 3L, 2L, 1L, 3L, 3L), .Label = c("A", "B", "C"), class = "factor"), factor.sample.letters.1.2...6..replace...T.. = structure(c(2L, 2L, 1L, 1L, 2L, 1L), .Label = c("a", "b"), class =

Thanks! Do you have any suggestions on what I might do next. I have the phones, I have asterisk, and I have everything setup. But i can't get the login to work with the Polycom function. Nothing happens...and I can't find any readmes' or manuals. Regards, Jan -----Ursprungligt meddelande----- Fr?n: asterisk-users-bounces@lists.digium.com