[R-meta] modifying the default forest plot

Mukht Kaur theoneliberatedkaur at gmail.com
Tue Feb 13 21:13:43 CET 2018


Hi Wolfgang,
For some reason I am not able to use the previous statement, i.e.,
print(forest(random, showweights=TRUE, xlab='Estimated Effect Size',
order="obs", slab=paste(play1$author, play1$year, sep=", ")))

I believe I used it successfully earlier but now the following error shows
up:
Error in forest.rma(random, showweights = TRUE, xlab = "Estimated Effect
Size",  :
  Number of outcomes does not correspond to the length of the 'slab'
argument.

Any suggestions?

Please let me know.
Thanks


On Tue, Feb 13, 2018 at 3:10 PM, Mukht Kaur <theoneliberatedkaur at gmail.com>
wrote:

> Hi Wolfgang,
> For some reason I am not able to use the
>
> On Tue, Dec 5, 2017 at 6:06 PM, Viechtbauer Wolfgang (SP) <
> wolfgang.viechtbauer at maastrichtuniversity.nl> wrote:
>
>> Please always cc the list.
>>
>> As for adjusting the arrangement of the various elements in the plot, you
>> have two options:
>>
>> 1) Make the plotting device wider. I don't know how you are creating the
>> plot, but if you are doing this 'interactively', you can just do this by
>> increasing the width of the window for the plot. If you are using something
>> like png(), pdf(), or some other plotting device to save the plot, then you
>> adjust the 'width' argument.
>>
>> 2) You can also play around with the 'xlim' argument. First do this:
>>
>> print(forest(random13, showweights=TRUE, xlab='Estimated Effect Size',
>> order="obs", slab=paste(play1$author, play1$year, sep=", ")))
>>
>> Note the default value for 'xlim'. Then rerun the forest() command as
>> above (without the print()) but adjust 'xlim', that is, use xlim=c(lower,
>> upper). Increase/decrease the lower and upper values to see how this
>> affects the plot.
>>
>> To label elements, use text(). For example:
>>
>> text(lower, 15, "Author(s) and Year", pos=4)
>> text(upper, 15, "ES [95% CI]",        pos=2)
>>
>> where lower/upper are the values you specified for xlim and you need to
>> adjust the '15' depending on how many rows are in your forest plot (the
>> second value of text() determines the y coordinate where the text is
>> written). You can use text() to also add a label for the weights.
>>
>> You might also want to go through the examples at the end of
>> help(forest.rma) to get an idea of what can be done.
>>
>> Best,
>> Wolfgang
>>
>> -----Original Message-----
>> From: Mukht Kaur [mailto:theoneliberatedkaur at gmail.com]
>> Sent: Tuesday, 05 December, 2017 23:31
>> To: Viechtbauer Wolfgang (SP)
>> Subject: Re: [R-meta] modifying the default forest plot
>>
>> Thank you guys! So I successfully made edits to my default forest plot
>> using
>> forest (random13, showweights=TRUE, xlab='Estimated Effect Size',
>> order="obs", slab=paste(play1$author, play1$year, sep=", "))
>>
>> The weights assigned are represented in the figure but they seem too
>> close to the CIs. Is there a way to label them and CIs? I looked at other
>> arguments which can be used for the 'forest' function. None seem applicable
>> to the labelling I intend to do. Any ideas?
>>
>> On Tue, Dec 5, 2017 at 4:58 AM, Viechtbauer Wolfgang (SP) <
>> wolfgang.viechtbauer at maastrichtuniversity.nl> wrote:
>> Ken already gave some pointers.
>>
>> 1) How could I get the author and year instead?
>>
>> Use the 'slab' argument. For example:
>>
>> slab=paste(dat$author, dat$year, sep=", "))
>>
>> (assuming that there are 'author' and 'year' variables in the data frame
>> 'dat').
>>
>> 2) It would be better if I can get them arranged according to the
>> observed effect size - is that possible?
>>
>> Use argument: order="obs"
>>
>> 3) Also, I would like to change the title of my X axis to 'Estimated
>> Effect Size'
>> rather than the default title 'Observed Outcome'.
>>
>> Use argument: xlab="Estimated Effect Size"
>>
>> 4) Lastly, would it be possible to specify the exact weight assigned to
>> each study in the forest plot?
>>
>> Use argument: showweights=TRUE
>>
>> Best,
>> Wolfgang
>>
>> -----Original Message-----
>> From: R-sig-meta-analysis [mailto:r-sig-meta-analysis-bo
>> unces at r-project.org] On Behalf Of Mukht Kaur
>> Sent: Tuesday, 05 December, 2017 6:30
>> To: r-sig-meta-analysis at r-project.org
>> Subject: [R-meta] modifying the default forest plot
>>
>> Hi everyone!
>>
>> I am new to meta-analysis methodology and I am using R software to do my
>> quantitative synthesis. I am using random effects model to analyse my
>> individual studies through rma function in metafor package.
>>
>> By using 'forest(mymodelname)', I have a basic forest plot where my
>> individual studies are shown as Study 1, Study 2, and so on. How could I
>> get the author and year instead? It would be better if I can get them
>> arranged according to the observed effect size - is that possible? Also, I
>> would like to change the title of my X axis to 'Estimated Effect Size'
>> rather than the default title 'Observed Outcome'. Lastly, would it be
>> possible to specify the exact weight assigned to each study in the forest
>> plot?
>> Many thanks in advance for your help.
>>
>> MK
>>
>
>

	[[alternative HTML version deleted]]



More information about the R-sig-meta-analysis mailing list