similar to: New Meta-Analysis Package (metafor)

Dear R users, In metafor 1.6-0, the Mantel-Haenszel method is implemented by the rma.mh() function. I have observed that the sum of the weights computed by weights(x) doesn't add to 100% when x is an object of class rma.mh. The consequences of this fact can be clearly seen when a forest diagram is drawn with forest(x), which calls weights(x) (or more precisely, the method weights.rma.mh()

Dear R users: In the example of meta-analysis (cochrane, package rmeta), I can not found the p-value of Test for overall effect, and some other indices (Z, I, weight and et al). How can I get the these indices listed? > library(rmeta) > data(cochrane) > cochrane name ev.trt n.trt ev.ctrl n.ctrl 1 Auckland 36 532 60 538 2 Block 1 69 5

Dear Dushanthi, Please keep your e-mails on the R-Help list, where Michael has already given you some excellent advice. As Michael already explained, metafor can handle proportions, but does not have any specific functionality for categorical variables with more than 2 levels (at the moment). So, if it is logical and possible to do so, you could collapse the levels of the categorical outcome to 2

Hello, I experimented the Metafor and Meta packages in the scope of replacing Excel for meta-analysis. I performed the first working example provided in Michael Borenstein's book "Introduction to Meta-Analysis" with Excel, Metafor and Meta. The numbers given by my spreadsheet, which I validated from Borenstein's book, conrespond quite closely to those given by Meta, but

Hello, Great to see the new metafor package for meta-analysis. I would like to perform a meta-analysis in which I initially calculate the intercept of the model with a nested random-effects structure. In lme, this would be model<- lme(v3~1, random=~1|species/study, weights = varFixed(~Wt), method = "REML") where multiple effects sizes are measured for some studies and more than

Hi - I am conducting a meta-analysis and I have a matrix of f-statistics, Ns and dfs from a series of studies that tested for an interaction in a 2x2 anova. I'd like to test whether the 2x2 interaction is significant in the aggregate. Similarly, I have a matrix of chi-square statistics that I'd like to meta-analyze. How can I input these test statistics into the "metafor"

Hello, I once used the metawin software to perform a meta-analysis (see metawinsoft, Rosenberg et al.) and produced normal qqplot to test for a potential bias in the dataset. I now want to re-use the same dataset with the package metafor by W. Viechtbauer (great package btw). I run the qqnorm.rma.uni function. I use standardized effect sizes as in metawin. QQplot generated with metafor differs

Greetings, I wish to conduct a meta-analysis for which the outcome is a continuous variable measured on the same individuals before and after an intervention. Hence, the comparison is not made between two groups, but within groups, at diffrent times. Each study reports the mean outcome and SD before the intervention and the mean outcome and SD after the intervention. While p-values for

Dear all, We are trying to do a meta regression in R using the lme function. The reason for doing this with lme function is that we have covariates and studies within references. In S-Plus this is possible by using the following command: lme(outcome ~ covars, random = ~1 | reference/study, weights = varFixed(~var.outcome), data = mydata, control = lmeControl(sigma = 1)) This means that the

Dear All, I have been checking the metafor package against R-devel. R CMD check --as-cran metafor yields one note: FOSS licence with BuildVignettes: false Yes, I have 'BuildVignettes: FALSE' in my DESCRIPTION file. I see at http://developer.r-project.org/blosxom.cgi/R-devel/NEWS: Tue, 25 Jun 2013 CHANGES IN R-devel UTILITIES 'R CMD check --as-cran' warns about a false

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

Hi, I would like to do a meta-analysis, i.e., a mixed-effects regression, but I don't seem to get what I want using both the nlme or metafor packages. My question: is there indeed no way to do it? And if so, is there another package I could use? Here are the details: In my meta-analysis I'm comparing different studies that report a measure at time zero and after a certain followup

Dear Wolfgang Viechtbauer and R users, I have few questions regarding the development of the package 'metafor. As you suggested , I post to the R-help mailing list. I read you're planning an extension of this method to the multivariate case. I think it would be a useful tool. I'm currently performing some analyses with R on multiple outcomes, using the Stata command mvmeta to get

Hello, I have some satellite tag time-at-depth (TAD) frequency data that I would like some help with. The data was transmitted via satellite as percent time spent in each of 7 depth bins (0m, 0-1m, 1-10m, 10-50m etc.), binned over 6-hour intervals. I categorized each row of data corresponding to a date and time into summer vs. winter, and day vs. night, and then summed and averaged the given

Dear R-experts! Running a meta-analysis (using the magnificent metafor-package), I use cohen's d as a main outcome measure in a random-effects model. For most of the samples cohen's d is derived form a comparison of two groups (A & B). However some studies report results from an ANOVA (one-factor with three levels: C,D,E) whereas two groups (C,D) correspond to one group in the other

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

Hi all, Does anyone know if an R function for the Mantel-Haenszel chi-square for a 2x2 table exists? I’ve also seen it called the randomization Q statistic. Note that I’m not looking for the Cochran-Mantel-Haenszel…I did see that out there as cmh.test. Thanks in advance, Jeff Internal Virus Database is out-of-date. Checked by AVG Free Edition. 11:27 AM [[alternative HTML version

Hello, Dr. Viechtbauer. I am trying to perform a meta-analyis on a group of before-after studies using Metafor. I read your webpage including your correspondence with Dr. Dewey (https://stat.ethz.ch/pipermail/r-help/2012-April/308946.html), who also conducted a similar study. These information is very hepful, but I have one additonal question which I wonder if you can give me some instruction.

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 have a print method for a set of statistical tests, vcdExtra::CMHtest, for which I'd like to have more sensible printing of pvalues, as in print.anova(). [Testing this requires the latest version of vcdExtra, from R-Forge **|install.packages("vcdExtra", repos="http://R-Forge.R-project.org")|** ] With my current print method, I get results like this, but all Prob values