similar to: metaplot and meta.summaries

Hello, I'm ashamed to admit that I'm an R and meta-analysis neophyte just when my colleagues and I are attempting to run R for that latter purpose. We would like to construct a graph of the 95% confidence intervals of odds ratios (inverse odds ratios, actually). We are using the metaplot() function in the "rmeta" package. While we can construct a basic plot, we cannot configure

Der useRs, I'm working on meta analysis using rmeta package. Using code below I plot the forest plot: library(rmeta) data (catheter) a<-meta.MH (n.trt, n.ctrl, col.trt, col.ctrl, data=catheter, names=Name, subset=c(13,6,5,3,7,12,4,11,1,8,10,2)) summary(a) # odds ratio values and confidence intervals metaplot(a$logOR, a$selogOR, nn=a$selogOR^-2,a$names, summn=a$logMH, sumse=a$selogMH,

Hello! I have some data stored into 2 separate csv file. 1 file (called A.csv) (12 results named Group1, Group2, Group3, etc...) odds ratios, 2 file (called B.csv) 12 corresponded errors. How to import that data into R and make forest plot like I saw inside help file Rmeta and meta with included different font colors and names trough X and Y axis. I know for meta libb ... out <-

Hi, I have a problem with the metaplot-function in rmeta. After plotting a Forest Plot with "metaplot", no y-axis is shown for normal plots anymore. Does someone know, why this is and how I can solve it? Thank you in advance! Kind regards, Roman -- View this message in context: http://www.nabble.com/metaplot-in-rmeta%3A-y-axis-disappears-tp24917401p24917401.html Sent from the R

Dear all, I have a four variable: Stuy.Name, OR, 95%LCI and 95%UCI and I would like to create a meta analysis plot. I can't use meta.MH function in metaplot because I do not have n.trt, n.ctrl, col.trt, col.ctrl are not available! Is there an alternative way to do it? Many thanks in advance, Cheba [[alternative HTML version deleted]]

Hello, We are looking to adjust the font size of the axis annotation on the graph that results from use of the metaplot() function. Metaplot seems to respond to cex and cex.lab to change those graphical parameters, but it doesn't respond to cex.axis. Is there a way to work around this by creating a customized x-axis, and if so, how? Thanks for all your help. Syntax is below. Best, Dawn

Greetings I have been trying to use metaplot to create a forest plot of my meta-analysis results. My data is in the form of ORs, so I've been inputting the ORs as the point estimate (mn) and the SEs of the logORs as the SE (se). However, whilst the point estimates appear to be displayed properly, many of the confidence intervals go below zero (not possible for an OR, as far as I know). I

Dear all, I could not draw a forest plot for meta-analysis under ramdon models using the rmeta package. The rmeta has a default function for MH (fixed-effect) model. Has the rmeta package been updated for such a function? Or someone revised it and kept a private code? I would appreciate it if you could provide some information on this question. Thanks, Andrew This email is intended

Hello I'm trying to plot hazard risk values using the function metaplot with the specifications: > metaplot(HR,SE,W,labels=row.names(lc),xlab="Hazard Ratio",ylab="Covariates", logeffect=TRUE,logticks=FALSE,colors=meta.colors(box="black",lines="dark gray",zero="darkgray"),cex=1.5,cex.lab=1.5,font=3) However, in the plot the x axis

Hi, can you tel me how can I make a Forrest Plot with R? It is possible and easy or are there a more practical free software available? Than you Mic [[alternative HTML version deleted]]

Hi, I'm trying to build a Forest Plot using the second and fourth columns in the table (test.csv) below. My code is the following: curated <- data.frame("test.csv") tmp <- curated$coef tmp1 <- curated$se_coef plt <- metaplot(tmp, tmp1, xlim = c(-.45, .45)) I keep getting the following error at the last line and am not sure why: Error in if (is.na(lower[i] + upper[i]))

Dear all, I am trying to get the emperical bayes estimates together with their standard errors out of lme4. Up to now I have used MLwiN to get these estimates. I have fitted the following - very simple - model, just to find out how this works. test<-lmer(y~(1|subject),data,REML=F) ranef(test,postVar=T) str(ranef(test,postVar=T) If I use the formulation of the emperical bayes estimates and

Dear all I have a model that looks like this: m1 <- lmer(Difference ~ 1+ (1|Examiner) + (1|Item), data=englisho.data) I know it is not possible to estimate random effects but one can obtain BLUPs of the conditional modes with re1 <- ranef(m1, postVar=T) And then dotplot(re1) for the examiner and item levels gives me a nice prediction interval. But I would like to have the prediction

Hello! I wonder how can one get upper and lower limits of a prediction interval -- exact values? They are shown on caterpillar plot using ranef() with argument postVar=TRUE, but I would like to know them. A while ago, some discussions were opened on "Confidence Intervals for Random Effect BLUP's", but the answer was never clear: http://www.mail-archive.com/r-help at

Dear R-list, I'd like to do a meta-analysis plot similar to install.packages('rmeta') require(rmeta) data(catheter) a <- meta.MH(n.trt, n.ctrl, col.trt, col.ctrl, data=catheter, names=Name, subset=c(13,6,5,3,7,12,4,11,1,8,10,2)) summary(a) plot(a) (see attached file) by using my own OR (Odds Ratio) and 95% Confidence Interval data set, which looks like

Dear Ioannis Thank you very much for pointing me to meta-analysis. Although it may not solve my problem with the normalization, it gives me some other options to display the different correlation coefficients. One possibility is the use of Funnel plots, which are even available in library(rmeta). Another possibility is the use of forest-plots, as implemented in rmeta as metaplot. Sorrowly,

Has anyone developed a function to generate a forest plot, the one used a lot in meta-analysis? Joseph F. Lucke, PhD Biostatistician Center for Clinical Research and Evidence-based Medicine Department of Pediatrics School of Medicine University of Texas Health Science Center at Houston Voice: 713-500-5651 Email: Joseph.F.Lucke@uth.tmc.edu Postal Mail: PO Box 20708, Houston, TX 77225-0708

Hi, I have the following table of odds ratios (or), lower limits(ll) and upper limits(ul), which I would like to plot as horizontal lines beginning at the lower limit, ending at the upper limits and with a dot at the odds ratio on an x-axis on a log10 scale. The y axis would be the study sites. >From what I can figure out, it looks like the plotCI function will do everything except give me an

Hello, from Verzani, simpleR (pdf), p. 80, I created the following function to test the coefficient of lm() against an arbitrary value. coeff.test <- function(lm.result, var, coeffname, value) { # null hypothesis: coeff = value # alternative hypothesis: coeff != value es <- resid(lm.result) coeff <- (coefficients(lm.result))[[coeffname]] # degrees of freedom = length(var) -

Using ranef() (package nlme, version 3.1-75) with an 'lme' object I can obtain random effects for intercept and slope of a certain level (say: 1) - this corresponds to (say level 1) "residuals" in MLWin. Maybe I'm mistaken here, but the results are identical. However, if I try to get the standardized random effects adding the paramter "standard=T" to the