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
mydata=data.frame(OR=c(2.04545454545, 1.10434782609, 1.22588104401,
1.14102564103,
1.20579245527, 1.375, 1.16535433071),
L95=c(1.22839621997, 0.858106819302, 1.0964802088, 0.841934120955,
0.969786886818, 1.01498537023, 0.919391492382),
U95=c(3.40546755139, 1.42122051928, 1.37055308613, 1.54632513827,
1.49917372998, 1.86258857302, 1.47707220868)
)
rownames(mydata)=c(paste("Study",1:6,sep=""),'Summary')
mydata
My problem is that I don't have the raw data as rmeta _requires_ and, even
when I have my data set in the _same_ (?) format that summary(a), when I
tried plot(mydata) it doesn't work. Another approach I used was to change
the class of my object but it didn't work either. I'm running XP SP2 on
a
2.4 GHz Intel-Core 2 Duo processor and my R-session info is the following:
R version 2.7.2 RC (2008-08-18 r46388)
i386-pc-mingw32
locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
States.1252;LC_MONETARY=English_United
States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252
attached base packages:
[1] grid stats graphics grDevices utils datasets methods
base
other attached packages:
[1] rmeta_2.14 RODBC_1.2-3
loaded via a namespace (and not attached):
[1] tools_2.7.2
I would greatly appreciate any ideas about how should I proceed.
Thanks in advance,
Jorge
-------------- next part --------------
A non-text attachment was scrubbed...
Name: example MH.pdf
Type: application/pdf
Size: 4022 bytes
Desc: not available
URL:
<https://stat.ethz.ch/pipermail/r-help/attachments/20080825/6b7256ee/attachment.pdf>
See ?metaplot On Mon, Aug 25, 2008 at 2:44 PM, Jorge Ivan Velez <jorgeivanvelez at gmail.com> wrote:> 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 > > mydata=data.frame(OR=c(2.04545454545, 1.10434782609, 1.22588104401, > 1.14102564103, > 1.20579245527, 1.375, 1.16535433071), > L95=c(1.22839621997, 0.858106819302, 1.0964802088, 0.841934120955, > 0.969786886818, 1.01498537023, 0.919391492382), > U95=c(3.40546755139, 1.42122051928, 1.37055308613, 1.54632513827, > 1.49917372998, 1.86258857302, 1.47707220868) > ) > rownames(mydata)=c(paste("Study",1:6,sep=""),'Summary') > mydata > > > My problem is that I don't have the raw data as rmeta _requires_ and, even > when I have my data set in the _same_ (?) format that summary(a), when I > tried plot(mydata) it doesn't work. Another approach I used was to change > the class of my object but it didn't work either. I'm running XP SP2 on a > 2.4 GHz Intel-Core 2 Duo processor and my R-session info is the following: > > R version 2.7.2 RC (2008-08-18 r46388) > i386-pc-mingw32 > > locale: > LC_COLLATE=English_United States.1252;LC_CTYPE=English_United > States.1252;LC_MONETARY=English_United > States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 > > attached base packages: > [1] grid stats graphics grDevices utils datasets methods > base > > other attached packages: > [1] rmeta_2.14 RODBC_1.2-3 > > loaded via a namespace (and not attached): > [1] tools_2.7.2 > > > I would greatly appreciate any ideas about how should I proceed. > > Thanks in advance, > > > Jorge > > ______________________________________________ > R-help at r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > >
If I have an file like an JPG, gif, png or other, exist any kind of function
in the R to see de dimension in pixels of this image?
Thanks
Atenciosamente,
Leandro Lins Marino
Centro de Avaliacao
Fundacao CESGRANRIO
Rua Santa Alexandrina, 1011 - 2o andar
Rio de Janeiro, RJ - CEP: 20261-903
* (21) 2103-9600 R.:236
* (21) 8777-7907
* leandro em cesgranrio.org.br <mailto:leandro em cesgranrio.org.br>
"Aquele que suporta o peso da sociedade
e precisamente aquele que obtem
as menores vantagens". (SMITH, Adam)
P Antes de imprimir pense em sua responsabilidade e compromisso com o MEIO
AMBIENTE
At 19:44 25/08/2008, Jorge Ivan Velez wrote:>Dear R-list, > >I'd like to do a meta-analysis plot similar toSince these plots are known as forest plots ?forestplot might help you. As it says you have to do a bit more work but you do get much more flexibility>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 > >mydata=data.frame(OR=c(2.04545454545, 1.10434782609, 1.22588104401, >1.14102564103, >1.20579245527, 1.375, 1.16535433071), >L95=c(1.22839621997, 0.858106819302, 1.0964802088, 0.841934120955, >0.969786886818, 1.01498537023, 0.919391492382), >U95=c(3.40546755139, 1.42122051928, 1.37055308613, 1.54632513827, >1.49917372998, 1.86258857302, 1.47707220868) >) >rownames(mydata)=c(paste("Study",1:6,sep=""),'Summary') >mydata > > >My problem is that I don't have the raw data as rmeta _requires_ and, even >when I have my data set in the _same_ (?) format that summary(a), when I >tried plot(mydata) it doesn't work. Another approach I used was to change >the class of my object but it didn't work either. I'm running XP SP2 on a >2.4 GHz Intel-Core 2 Duo processor and my R-session info is the following: > >R version 2.7.2 RC (2008-08-18 r46388) >i386-pc-mingw32 > >locale: >LC_COLLATE=English_United States.1252;LC_CTYPE=English_United >States.1252;LC_MONETARY=English_United >States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 > >attached base packages: >[1] grid stats graphics grDevices utils datasets methods >base > >other attached packages: >[1] rmeta_2.14 RODBC_1.2-3 > >loaded via a namespace (and not attached): >[1] tools_2.7.2 > > >I would greatly appreciate any ideas about how should I proceed. > >Thanks in advance, > > >Jorge > >Content-Type: application/pdf; name="example MH.pdf" >X-Attachment-Id: f_fkbfhzhb0 >Content-Disposition: attachment; filename="example MH.pdf"Michael Dewey http://www.aghmed.fsnet.co.uk
Another way to do it (based on summary data): http://finzi.psych.upenn.edu/R/Rhelp02a/archive/131342.html -- Jonathan Baron, Professor of Psychology, University of Pennsylvania Home page: http://www.sas.upenn.edu/~baron
On Mon, 25 Aug 2008, Jorge Ivan Velez wrote:> My problem is that I don't have the raw data as rmeta _requires_ and, even > when I have my data set in the _same_ (?) format that summary(a), when I > tried plot(mydata) it doesn't work.No, it doesn't _require_ that. If you just want a forest plot then use forestplot() or metaplot(). If you want to do a meta-analysis from summary data use meta.summaries(). -thomas