Dear Hugh
1 - Your post is unreadable because you posted in HTML so it will be
very difficult for anyone to see your V
2 - there is a separate mailing list for meta-analysis in R where you
may get more helpful responses as it is monitored by authors of most of
the main packages. Please see
Dear Laura
If you are unable for some reason to share the data why not incorporate
the output into an e-mail (and please turn of HTML as it mangles
everything). Putting the plots from profiling somewhere we can read them
would be a useful addition.
This looks at first glance one of those sit
Very difficult to diagnose what is going on without actually seeing the data.
But as I said on CV: Depending on the data, the variance components may not be
estimated precisely, so negative values for those kinds of pseudo-R^2
statistics are quite possible. In fact, if a particular moderator is
Hi Laura,
As far as I am concerned, you account for multiple effects sizes from the same
study by setting up and fitting an appropriate model, not by fiddling with the
weights. Several examples are described here:
http://www.metafor-project.org/doku.php/analyses#multivariate_multilevel_meta-ana
Dear Kathryn
I think that the author of metafor has addressed this
http://www.metafor-project.org/doku.php/tips:comp_two_independent_estimates
The other tips on that site are well worth reading too.
On 29/07/2016 14:44, Morris, Kathryn wrote:
I am running a meta-analysis using metafor and get
I think you should consult with a local statistician. Generally
speaking, statistical questions like this tend to be OT here, and you
appear to be sufficiently confused about the statistical issues that
online posts would not be sufficient.
Cheers,
Bert
Bert Gunter
"The trouble with having an o
From: John Peterson [mailto:john.peterson@gmail.com]
> > Sent: Tuesday, December 08, 2015 14:25
> > To: Viechtbauer Wolfgang (STAT)
> > Cc: R-help@r-project.org
> > Subject: Re: [R] metafor package
> >
> > Hi Dr. Viechtbauer,
> >
> > The code
Psychiatry and
> Neuropsychology | Maastricht University | P.O. Box 616 (VIJV1) | 6200 MD
> Maastricht, The Netherlands | +31 (43) 388-4170 | http://www.wvbauer.com
>
> > -Original Message-
> > From: John Peterson [mailto:john.peterson@gmail.com]
> > Sent: Monda
m: John Peterson [mailto:john.peterson@gmail.com]
> Sent: Tuesday, December 08, 2015 14:25
> To: Viechtbauer Wolfgang (STAT)
> Cc: R-help@r-project.org
> Subject: Re: [R] metafor package
>
> Hi Dr. Viechtbauer,
>
> The code provided in the metafor projects website for subgroup i
r 07, 2015 22:14
> To: Viechtbauer Wolfgang (STAT)
> Subject: Re: [R] metafor package
>
> Hi Dr. Viechtbauer,
> Thank you very much for your reply. I tried your advice and was able to
> make a forest plot with only the estimates for each subgroup.For the
> estiamte for each subgroup,
The code you posted is totally mangled up, but it's just what can be found here:
http://www.metafor-project.org/doku.php/plots:forest_plot_with_subgroups
If you don't want an overall estimate, just pass the estimates and
corresponding sampling variances to the forest() function (and not the mode
-
Von: Viechtbauer Wolfgang (STAT)
[mailto:wolfgang.viechtba...@maastrichtuniversity.nl]
Gesendet: Freitag, 27. November 2015 16:01
An: Markus Kösters; 'Michael Dewey'; r-help@r-project.org
Betreff: RE: [R] metafor - Meta-Analysis of rare events / beta-binomial
regression
I would say the
-help-boun...@r-project.org] On Behalf Of Markus
> Kösters
> Sent: Friday, November 27, 2015 14:38
> To: 'Michael Dewey'; r-help@r-project.org
> Subject: Re: [R] metafor - Meta-Analysis of rare events / beta-binomial
> regression
>
> Dear Michael,
>
> Thank you ve
co.uk]
Gesendet: Freitag, 27. November 2015 13:32
An: Markus Kösters; r-help@r-project.org
Betreff: Re: [R] metafor - Meta-Analysis of rare events / beta-binomial
regression
Dear Markus
This is not a direct answer to your question, I will leave that to Wolfgang but
two thoughts:
1 - if all the st
Dear Markus
This is not a direct answer to your question, I will leave that to
Wolfgang but two thoughts:
1 - if all the studies have very sparse data
@article{bradburn07,
author = {Bradburn, M J and Deeks, J J and Berlin, J A and Localio,
A R},
title = {Much ado about nothing: a compa
europsychology | Maastricht University | P.O. Box 616 (VIJV1) | 6200 MD
> Maastricht, The Netherlands | +31 (43) 388-4170 | http://www.wvbauer.com
>
>> -Original Message-
>> From: Marco Colagrossi [mailto:marco.colagro...@gmail.com]
>> Sent: Monday, August 31, 2015 18:37
>&
email]
<http:///user/SendEmail.jtp?type=node&node=4711682&i=1>
Subject: Re: [R] Metafor and forest(); not showing 'ilab' and text
The solution that you proposed works perfectly, thank you very much.
I'll wait for Wolfgang answer as I'm having few doubts about the
;]
> > Sent: Monday, August 31, 2015 18:37
> > To: Michael Dewey
> > Cc: Viechtbauer Wolfgang (STAT); [hidden email]
> <http:///user/SendEmail.jtp?type=node&node=4711682&i=1>
> > Subject: Re: [R] Metafor and forest(); not showing 'ilab' and te
lfgang (STAT); r-help@r-project.org
> Subject: Re: [R] Metafor and forest(); not showing 'ilab' and text
>
> The solution that you proposed works perfectly, thank you very much.
>
> I'll wait for Wolfgang answer as I'm having few doubts about the models.
>
>
hatever
>>>>> device
>>>>> you are actually using. You can then use the 'xlim' argument to create
>>>>> appropriate spacing to the left/right of the part of the plot that
>>>>> shows the
>>>>> estimates and their CIs. Wit
os' values to create a nice
looking plot that has no overlapping text and no excessive white space. An
example is this:
http://www.metafor-project.org/doku.php/plots:forest_plot_with_subgroups
Note that it took me dozens of iterations to create that plot. You just
have to start experimentin
d 'ilab.xpos' values to create a nice
>>> looking plot that has no overlapping text and no excessive white space. An
>>> example is this:
>>>
>>> http://www.metafor-project.org/doku.php/plots:forest_plot_with_subgroups
>>>
>>> Note that it t
rco.colagro...@gmail.com]
Sent: Tuesday, August 25, 2015 17:59
To: Viechtbauer Wolfgang (STAT)
Cc: r-help@r-project.org; Michael Dewey
Subject: Re: [R] Metafor and forest(); not showing 'ilab' and text
Thanks again for your help. I'm sorry to bother you but I don't get
how to
ssage-
>> From: Michael Dewey [mailto:li...@dewey.myzen.co.uk]
>> Sent: Tuesday, August 25, 2015 13:23
>> To: Marco Colagrossi; Viechtbauer Wolfgang (STAT)
>> Cc: r-help@r-project.org
>> Subject: Re: [R] Metafor and forest(); not showing 'ilab' and text
&
ot. You just have to
start experimenting.
Best,
Wolfgang
> -Original Message-
> From: Marco Colagrossi [mailto:marco.colagro...@gmail.com]
> Sent: Tuesday, August 25, 2015 17:59
> To: Viechtbauer Wolfgang (STAT)
> Cc: r-help@r-project.org; Michael Dewey
> Subject: Re: [R] Me
Further comments in line as well.
> -Original Message-
> From: Michael Dewey [mailto:li...@dewey.myzen.co.uk]
> Sent: Tuesday, August 25, 2015 13:23
> To: Marco Colagrossi; Viechtbauer Wolfgang (STAT)
> Cc: r-help@r-project.org
> Subject: Re: [R] Metafor and forest();
Hello Marco
Comments in line again
On 24/08/2015 18:49, Marco Colagrossi wrote:
I tried to upload the file once again. I tweaked it a bit, now my code is:
forest(pc, var, ci95m, ci95p, slab = authoryear, psize=1, subset=(pub==1),
xlim = c(-16, 6),
ilab = cbind(SIMdv, SIMiv),
I tried to upload the file once again. I tweaked it a bit, now my code is:
forest(pc, var, ci95m, ci95p, slab = authoryear, psize=1, subset=(pub==1),
xlim = c(-16, 6),
ilab = cbind(SIMdv, SIMiv),
ilab.xpos = c(-7.5, -5.5), cex = 0.75)
op <- par(cex=.75, font=2)
text(c(-7
I cannot reproduce the issue with 'ilab' not being shown when using 'subset'.
My guess is that the values for 'ilab.xpos' specified are actually outside of
the plotting region. After you have drawn the forest plot, try:
par("usr")[1:2]
to see what the default limits actually are. Then use 'xlim
Dear Marco
Comments inline
On 24/08/2015 15:03, Marco Colagrossi wrote:
Hello folks,
I have a couple of issues with the metafor package, specifically with
the forest graphs.
I am currently conducting a Meta-Analysis in economics throughout the
metafor package.
My meta-analysis has the specifi
errors of the variances in
metafor (and do a Z-test)?
> To: wibbelt...@hotmail.com; r-help@r-project.org
> Subject: Re: [R] Metafor - rma.mv function - variance components
> From: li...@dewey.myzen.co.uk
> Date: Mon, 20 Apr 2015 20:24:48 +
>
> Carlijn Wibbelink wrote
April 21, 2015 10:43
> To: li...@dewey.myzen.co.uk; r-help@r-project.org
> Subject: Re: [R] Metafor - rma.mv function - variance components
>
> Thank you for your reaction, it worked.
> However, I'm wondering if this is the right way to test whether there is
> significant var
4 <- rma.mv(y, v, random = list(~ 1 | y, ~ 1 | ID), sigma2=c(0,NA),
data=dat)
anova(model2,model3)
anova(model2,model4)
Is it possible to receive the standard errors of the variances in metafor (and
do a Z-test)?
> To: wibbelt...@hotmail.com; r-help@r-project.org
> Subject: Re:
Carlijn Wibbelink wrote :
Dear Carlijn
I think that if you set sigma2 to a vector of length 2 it will be possible.
> Hi all,
>
> I have a question about metafor and the rma.mv function. I have fitted a
> multivariate model (effect sizes are nested within studies) and I've found two
> variances:
> Hi, I have a quesite on meta-analysis with 'metafor'.
> I would like to calculate the standardized mean difference (SMD), as
> Hedges' g, in pre-post design studies.
> I have data on baseline (sample size, mean and SD in both the
> experimental
> and the control group) and at end of treatment (sa
Comment below
On 13/04/2015 20:46, Antonello Preti wrote:
Hi, this is another quesite related to the use of 'metafor' for calculation
of standardized mean change in pre-post design studies.
Essentially, my aim is to compare different method to arrive at the same
conclusion: Does the treatment wo
an der Voet, H.
> (2012). Transformations of summary statistics as input in meta-analysis for
> linear dose-response models on a logarithmic scale: A methodology developed
> within EURRECA. BMC Medical Research Methodology, 12(57).
>
> Best,
> Wolfgang
>
>> -Original
Michael Dewey; Purssell, Ed; r-help@r-project.org
Subject: RE: [R] metafor - code for analysing geometric means
With "geometric mean 1 CI /3.92", I assume you mean "(upper bound - lower
bound) / 3.92". Two things:
1) That will give you the SE of the mean, not the SD of the observa
p-boun...@r-project.org]
> On Behalf Of Michael Dewey
> Sent: Thursday, November 13, 2014 12:36
> To: Purssell, Ed; r-help@r-project.org
> Subject: Re: [R] metafor - code for analysing geometric means
>
> On 13/11/2014 11:00, Purssell, Ed wrote:
> > ?Dear All
> >
On 13/11/2014 11:00, Purssell, Ed wrote:
?Dear All
I have some data expressed in geometric means and 95% confidence intervals.
Can I code them in metafor as:
rma(m1i=geometric mean 1, m2i=geometric mean 2, sd1i=geometric mean 1 CI /3.92,
sd2i=geometric mean 2 CI/3.92...etc, measure="
(STAT) [wolfgang.viechtba...@maastrichtuniversity.nl] Sent: 27
August 2014 17:30 To: Owen, Branwen; r-help@r-project.org Subject:
RE: [R] Metafor -can't calculate heterogeneity with non-positive
sampling variances The warning message pretty much says it: When one
of the variances is zero, then the I^
wen; r-help@r-project.org
Subject: RE: [R] Metafor -can't calculate heterogeneity with non-positive
sampling variances
The warning message pretty much says it: When one of the variances is zero,
then the I^2 statistic (and various other things) cannot be computed, at least
if one sticks to the u
The warning message pretty much says it: When one of the variances is zero,
then the I^2 statistic (and various other things) cannot be computed, at least
if one sticks to the usual equations/methods. So, if you think the 0 sampling
variances really make sense and you really want to get somethin
On 07/07/2014 12:02, Michael Dewey wrote:
At 09:08 07/07/2014, Dietlinde Schmidt wrote:
Thanks for that link, Wolfgang. Unfortunately, there comes the Warning
with it:
"(process:3634): Pango-WARNING **: Invalid UTF-8 string passed to
pango_layout_set_text()"
and decimal being "scrambled" in fore
At 09:08 07/07/2014, Dietlinde Schmidt wrote:
Thanks for that link, Wolfgang. Unfortunately,
there comes the Warning with it:
"(process:3634): Pango-WARNING **: Invalid UTF-8
string passed to pango_layout_set_text()"
and decimal being "scrambled" in forest plot and
not displaying the midline de
--
> > From: Dietlinde Schmidt [mailto:schmidt.dietli...@web.de]
> > Sent: Monday, July 07, 2014 10:09
> > To: Viechtbauer Wolfgang (STAT); r-help@r-project.org
> > Subject: Re: [R] metafor package: changing decimal in forest plot to
> > midline decimal
> >
> > Th
Thanks for that link, Wolfgang. Unfortunately, there comes the Warning
with it:
"(process:3634): Pango-WARNING **: Invalid UTF-8 string passed to
pango_layout_set_text()"
and decimal being "scrambled" in forest plot and not displaying the
midline decimal.
I think it has to do with the fact, th
olfgang (STAT); r-help@r-project.org
> Subject: Re: [R] metafor package: changing decimal in forest plot to
> midline decimal
>
> Thanks for that link, Wolfgang. Unfortunately, there comes the Warning
> with it:
> "(process:3634): Pango-WARNING **: Invalid UTF-8 string passe
I found this:
https://stat.ethz.ch/pipermail/r-help/2012-August/321057.html
So, use this before drawing the forest plot:
options(OutDec="\xB7")
Best,
Wolfgang
--
Wolfgang Viechtbauer, Ph.D., Statistician
Department of Psychiatry and Psychology
School for Mental Health and Neuroscience
Faculty
Pass 'slab' to the rma() function instead of the forest() function. Also,
'slab' is evaluated in the data frame passed via the 'data' argument, so no
need to use empa$... And there is a 'subset' argument for rma(). Might as well
use it. This should do it:
resultREML <- rma(yi=yi, vi=vi, method=
ope is
> positive or negative). If you think the suppression occurred on the other
> side than the one that is chosen, specify the side via this argument.
>
> Best,
> Wolfgang
>
> From: r-help-boun...@r-project.org [r-help-boun...@r-proj
Monday, May 26, 2014 1:15 PM
To: Michael Dewey
Cc: r-help
Subject: Re: [R] Metafor: Strange Trim and Fill Outcome?
thanks for your reply, Michael!
ad primary studies:
yup, I have a large set of studies: At the moment I consider 126 data sets
in my analysis.
ad interpretation:
thats an interesting
thanks for your reply, Michael!
ad primary studies:
yup, I have a large set of studies: At the moment I consider 126 data sets
in my analysis.
ad interpretation:
thats an interesting information. But usually there should be found an
effect.
ad "bit" :-):
I am irritiated, because in the plots I h
At 09:34 26/05/2014, Verena Weinbir wrote:
Hey guys,
I have tested the metafor trim and fill function (y:SD, x:SMD)on my data
set and yielded the following result:
1. missing studies on the right: 34
That seems a lot of missing studies unless you have a very large set
of primary studies.
r calculations as a bonus.
Best,
Wolfgang
> -Original Message-
> From: Verena Weinbir [mailto:vwein...@gmail.com]
> Sent: Tuesday, May 06, 2014 15:09
> To: Michael Dewey
> Cc: Viechtbauer Wolfgang (STAT);
<mailto:r-help@r-project.org>r-help@r-project.org
&
;
>
> Best,
>> Wolfgang
>>
>> > -Original Message-
>> > From: Verena Weinbir [mailto:vwein...@gmail.com]
>> > Sent: Tuesday, May 06, 2014 15:09
>> > To: Michael Dewey
>> > Cc: Viechtbauer Wolfgang (STAT); r-help@r-project.org
>&g
t; Sent: Tuesday, May 06, 2014 15:09
> To: Michael Dewey
> Cc: Viechtbauer Wolfgang (STAT); r-help@r-project.org
> Subject: Re: [R] Metafor: How to integrate effectsizes?
>
> Thank you very much for your illustration, Wolfgang! It helped me a
> lot. And also thank you for the pa
88-4170 | http://www.wvbauer.com
>>
>> > -Original Message-
>> > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org
>> ]
>> > On Behalf Of Michael Dewey
>> > Sent: Friday, April 25, 2014 16:23
>> > To: Verena Weinbir
).
Best,
Wolfgang
> -Original Message-
> From: Verena Weinbir [mailto:vwein...@gmail.com]
> Sent: Tuesday, May 06, 2014 15:09
> To: Michael Dewey
> Cc: Viechtbauer Wolfgang (STAT); r-help@r-project.org
> Subject: Re: [R] Metafor: How to integrate effectsizes?
>
>
o: Verena Weinbir
> Cc: r-help@r-project.org
> Subject: Re: [R] Metafor: How to integrate effectsizes?
>
> At 12:33 25/04/2014, you wrote:
> >Thank you very much for your reply and the book recommendation, Michael.
> >
> >Yes, I mean Cohen's d - sorry for th
ilto:r-help-boun...@r-project.org]
> On Behalf Of Michael Dewey
> Sent: Friday, April 25, 2014 16:23
> To: Verena Weinbir
> Cc: r-help@r-project.org
> Subject: Re: [R] Metafor: How to integrate effectsizes?
>
> At 12:33 25/04/2014, you wrote:
> >Thank you very much for your
At 12:33 25/04/2014, you wrote:
Thank you very much for your reply and the book recommendation, Michael.
Yes, I mean Cohen's d - sorry for the typo :-)
Just to make this sure for me: There is no
possibility to integrate stated Cohens' ds in an
R-Metaanalysis (or a MA at all), if there is no
At 13:15 24/04/2014, Verena Weinbir wrote:
Hello!
I am using the metafor package for my master's thesis as an R-newbie. While
calculating effectsizes from my dataset (mean values and
standarddeviations) using "escalc" shouldn't be a problem (I hope ;-)), I
wonder how I could at this point integr
At 11:56 22/04/2014, Dipl. Kfm Dominik Wagner MSc; MSc wrote:
Dear all,
I am quite new to R. Now my following easy question.
I use metafor and performed an outlier test with rstudent(res).
This is resulting in 1000 rows of 1578 and 578 omitted rows (starting with
row 598).
1. How can I dis
I think this may help:
http://www.metafor-project.org/doku.php/tips:handling_missing_data
I am not sure I understand your second question. All studies are shown (for
which the standardized residual can be computed), but since there are so many
studies, these plots are not really all that helpfu
Prof. Viechtbauer, thanks for the articles. I appreciate your help.
Yours,
John
John Williams
ALB Candidate, Harvard University (Expected May 2014)
johnwilli...@fas.harvard.edu
jawilliam...@gmail.com
--
View this message in context:
http://r.789695.n4.nabble.com/metafor-variance-ex
The standardized mean change using 'change score standardization' is described
in this article:
Gibbons, R. D., Hedeker, D. R., & Davis, J. M. (1993). Estimation of effect
size from a series of experiments involving paired comparisons. Journal of
Educational Statistics, 18(3), 271-279.
For a c
At 09:42 30/03/2014, Purssell, Ed wrote:
Dear All
I have a question about combining effect sizes generated by escalc
in metafor. I realise these may be stupid things to do; but they
are deliberately so to explain what I mean - I don't intend doing this!
I have 3 studies; each of which has
Often, there is a mix of information available from the various studies that
needs to be used to compute the effect sizes or outcomes to be used for the
meta-analysis. Then you have to 'build up' your dataset in multiple steps and
you cannot bypass first using escalc().
As a very basic example,
14 at 6:39 AM
To: "r-help@R-project.org" , "Nathan Leon Pace, MD,
MStat" , Michael Dewey , S
Ellison
Subject: RE: [R] metafor package
>Good advice already from Michael and S Ellison.
>
>I must apologize for the 'hack-job of a function' called forest() in
&
> -Original Message-
> I must apologize for the 'hack-job of a function' called forest() in metafor.
Well, certainly no apology needed there!
> I realized a while ago that people would prefer more fine-grained control over
> the various elements of the plot (this has come up a few times
+31 (43) 388-4170 | http://www.wvbauer.com
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> On Behalf Of Michael Dewey
> Sent: Wednesday, February 12, 2014 12:50
> To: Nathan Pace; r-help@R-project.org
> Subject: Re: [R
Like most plotting routines in R, forest has a ... argument that passes names
arguments through to other routines.
In forest's case, ... passes (at least) to abline, axis and segments, so all
the lines and axes in the plot are affected by any argument that is valid for
those.
In your case, speci
At 22:11 11/02/2014, Nathan Pace wrote:
Hi,
I have a random effects meta analysis of a proportion (logit
transformation) using rma.glmm.
I have created a forest plot of the proportion (inverse logic
transformation) using forest.rma.
I have added the credibility interval.
The forest plot i
Naples University Federico II
> Italy
>
> >
> Message: 50
> Date: Tue, 19 Nov 2013 11:48:33 +0100
> From: "Viechtbauer Wolfgang (STAT)"
>
> To: "petre...@unina.it" , "r-help@r-project.org"
>
> Subject: Re: [R] me
Dear Mario,
You can always just inspect the code:
escalc.default
Best,
Wolfgang
--
Wolfgang Viechtbauer, Ph.D., Statistician
Department of Psychiatry and Psychology
School for Mental Health and Neuroscience
Faculty of Health, Medicine, and Life Sciences
Maastricht University, P.O
At 22:12 01/06/2013, emhanlon wrote:
Good afternoon,
I am trying to use a mixed-effect model to find the best predictors of
wildlife response (fledgling success, offspring growth rate, etc...) to
human disturbance. I am using a Log Response Ratio as my measure of effect
size along with study va
Dear Branwen,
This means that your design matrix is not of full rank (in a more recent
version of the metafor package, the error message is a bit more informative;
i.e., please upgrade). Since "continent" is a factor, this should imply that
one of the levels never actually occurs, leading to a
Dear Owen,
What is your definition of "multivariate analysis"? Do you mean: A
meta-regression model with more than one predictor/moderator? In that case,
yes, metafor handles that. Usually, this is referred to as "multiple
regression" (as opposed to "simple regression" with a single predictor)
-project.org
Betreff: RE: [R] Metafor "SMCR" Pre-Post Effect sizes
Absolutely correct. The documentation wasn't clear on that. Changed for the
next version of the package. And the escalc() function now no longer checks
for sd2i, since it is not needed anyway (for "SMCR"). For now,
Absolutely correct. The documentation wasn't clear on that. Changed for the
next version of the package. And the escalc() function now no longer checks for
sd2i, since it is not needed anyway (for "SMCR"). For now, you will just have
to set sd2i to something (e.g., 0).
Thanks for bringing this
Please see my comments below.
Best,
Wolfgang
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> On Behalf Of Rachel Slatyer
> Sent: Wednesday, February 27, 2013 06:00
> To: R-help@r-project.org
> Subject: [R] metafor - interpretion of QM in m
At 16:08 10/09/2012, Jarrett Byrnes wrote:
Thanks for this. A few things
First, yes, my lmer syntax was indeed bad - I was writing this as an
example of what my data & code look like. Apologies. So, 1|studyID
indeed. Also 1/variance.
I've also been wondering - I often have more than 2 dru
Thanks for this. A few things
First, yes, my lmer syntax was indeed bad - I was writing this as an example of
what my data & code look like. Apologies. So, 1|studyID indeed. Also
1/variance.
I've also been wondering - I often have more than 2 drugs - so, T1, T2, T3,
etc. It depends on wh
As usual, Michael was faster than I in responding. Let me add a few thoughts of
my own. See comments below in text.
Best,
Wolfgang
--
Wolfgang Viechtbauer, Ph.D., Statistician
Department of Psychiatry and Psychology
School for Mental Health and Neuroscience
Faculty of Health, Medicin
At 15:02 07/09/2012, Jarrett Byrnes wrote:
Hello. A quick question about incorporating variation due to study
in the metafor package. I'm working with a particular data set for
meta-analysis where some studies have multiple measurements. Others
do not. So, let's say the effect I'm looking a
, so I don't know for sure.
Best,
Wolfgang
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
cpanderson [christopher.p.ander...@healthpartners.com]
Sent: Friday, August 03, 2012 6:03 PM
To: r-help@r-project.org
Subject: Re: [R] m
At 17:03 03/08/2012, cpanderson wrote:
Wolfgang,
Thanks for your quick response. You are correct- indeed I had inadvertently
left out that i had set intercept = FALSE in rma.
Dear Christopher,
It does help if you paste in the exact code you used, see also below
Following your suggestions, I
Wolfgang,
Thanks for your quick response. You are correct- indeed I had inadvertently
left out that i had set intercept = FALSE in rma.
Following your suggestions, I get the following results:
Test of Moderators (coefficient(s) 2,3):
QM(df = 2) = 0.2207, p-val = 0.8955
Model Results:
Are you sure that output was produced by: rma(yi, vi, data=dat,
mods=cbind(LateralWedge,Dome,Complex))?
Because your model does not have an intercept, which suggests that you actually
used:
rma(yi, vi, data=dat, mods=cbind(LateralWedge,Dome,Complex), intercept=FALSE)
In that case, the coeffici
1)
6200 MD Maastricht, The Netherlands
+31 (43) 388-4170 | http://www.wvbauer.com
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> On Behalf Of Jeremy Miles
> Sent: Friday, August 03, 2012 04:01
> To: Bexkens, Anika
> Cc: r
The test of moderator coefficients (QM) is chi-square distributed.You
can use the change in this value when you add a predictor to the model
as a chi-square test, with df equal to the change in df.
Jeremy
On 2 August 2012 05:54, Bexkens, Anika wrote:
> Dear Metafor users,
>
> I'd like to test a
Thanks very much for your response. Sorry for writing to you directly.
Hereafter, I'll email to the R-Help list. I agree with all your comments and
suggestions. My analysis will be descriptive type; purely observational as the
studies are non-randomized studies. I have sample sizes and counts
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). S
help@r-project.org
> From: i...@aghmed.fsnet.co.uk
> Subject: RE: [R] "metafor" package, proportions: single groups wrt to a
> categorical dependent variable
>
> At 01:18 31/07/2012, Dushanthi Pinnaduwage wrote:
>
> >Thanks very much for the reply Michael.
>
) it is a bit hard to offer
more advice. I would regard collapsing categories as a last resort
myself as you lose information.
Dushanthi
> Date: Sat, 28 Jul 2012 16:38:24 +0100
> To: dushan...@bell.net; r-help@r-project.org
> From: i...@aghmed.fsnet.co.uk
> Subject: RE: [R] "
Thanks very much for the reply Michael.I guess I have to combine some levels to
make categorical variable to a binary variable. Dushanthi > Date: Sat, 28 Jul
2012 16:38:24 +0100
> To: dushan...@bell.net; r-help@r-project.org
> From: i...@aghmed.fsnet.co.uk
> Subject: RE:
ou got (if any) or told us the
> discrepancy between the output you obtained and that which you expected.
> Date: Sat, 28 Jul 2012 14:25:08 +0100
> To: dushan...@bell.net; r-help@r-project.org
> From: i...@aghmed.fsnet.co.uk
> Subject: Re: [R] "metafor" package, proporti
you obtained and that which you expected.
> Date: Sat, 28 Jul 2012 14:25:08 +0100
> To: dushan...@bell.net; r-help@r-project.org
> From: i...@aghmed.fsnet.co.uk
> Subject: Re: [R] "metafor" package, proportions: single groups
wrt to a categorical dependent variable
>
At 01:44 28/07/2012, Dushanthi Pinnaduwage wrote:
Dear all, I am using R version 2.15.0 and 'metafor' package version
1.6-0. Can this version of the package handle proportions from
a categorical dependent variable for single studies?If so how do I
set up my dataframe for the raw data from
1 - 100 of 127 matches
Mail list logo
