Skip to content

Model studies in one analysis using treatment as a five level moderator in a meta-regression

5 messages · Jay Zola, vito muggeo, PIKAL Petr +1 more

Jay Zola
# 
Hello,


I am medical student, writing a meta-analysis on complication and reoperation rates after the five most common treatments of distal radius fractures. I have been busy with the statistics for months by my self, but find it quite hard since our classes were very basic. Now I want to compare the treatment modalities to see if there are significant differences. Using R I was able to synthesize the complication rates and reoperation rates for each treatment method. But I never had any R course and managed by trial and error, so the code probably doesn't look that great. Someone told me I could best model the data in one analysis using treatment as a five level moderator in a meta-regression. Can some help me with the R code to do this? Your help would be very much appreciated.


Thank you,


Jay


Study| Event Type| Treatment| Number of Events (n)| N| n/N|

Kumaravel| Complications| EF| 3| 23| 0,1304348|

Franck| Complications| EF| 2| 20| 0,1|

Schonnemann| Complications| EF| 8| 30| 0,2666667|

Aita| Complications| EF| 1| 16| 0,0625|

Hove| Complications| EF| 31| 39| 0,7948718|

Andersen| Complications| EF| 26| 75| 0,3466667|

Krughaug| Complications| EF| 22| 75| 0,2933333|

Moroni| Complications| EF| 0| 20| 0|

Plate| Complications| IMN| 3| 30| 0,1|

Chappuis| Complications| IMN| 4| 16| 0,25|

Gradl| Complications| IMN| 12| 66| 0,1818182|

Schonnemann| Complications| IMN| 6| 31| 0,1935484|

Aita| Complications| IMN| 1| 16| 0,0625|

Dremstrop| Complications| IMN| 17| 44| 0,3863636|

Wong| Complications| PC| 1| 30| 0,0333333|

Kumaravel| Complications| PC| 4| 25| 0,16|


Dataset on my dropbox: https://www.dropbox.com/s/j1urqzr99bt76ip/Basics%20excel%20file%20complication%20and%20reoperation%20rate.xlsx?dl=0

Basics excel file complication and reoperation rate.xlsx<https://www.dropbox.com/s/j1urqzr99bt76ip/Basics%20excel%20file%20complication%20and%20reoperation%20rate.xlsx?dl=0>
www.dropbox.com
Shared with Dropbox




library(meta)
library(stargazer)
library(foreign)

All <-read.spss("C:\\Users\\313635aa.STUDENT\\Desktop\\Meta-Analyse Complications and Reoperations.sav",to.data.frame = T, use.value.labels = T)
All <- na.omit(All)

Complications <- All[which(All[,"Event_Type"] == "Complications"),]
Re_operation <- All[which(All[,"Event_Type"] == "Reoperations"),]

EF <- All[which(All[,"Treatment"] == "EF"),]
IMN <- All[which(All[,"Treatment"] == "IMN"),]
pc <- All[which(All[,"Treatment"] == "PC"),]
KW <- All[which(All[,"Treatment"] == "KW"),]
VPO <- All[which(All[,"Treatment"] == "VPO"),]

EF_C <- EF[which(EF[,"Event_Type"] == "Complications"),]
EF_R <- EF[which(EF[,"Event_Type"] == "Reoperations"),]

IMN_C <- IMN[which(IMN[,"Event_Type"] == "Complications"),]
IMN_R <- IMN[which(IMN[,"Event_Type"] == "Reoperations"),]

pc_C <- pc[which(pc[,"Event_Type"] == "Complications"),]
pc_R <- pc[which(pc[,"Event_Type"] == "Reoperations"),]

KW_C <- KW[which(KW[,"Event_Type"] == "Complications"),]
KW_R <- KW[which(KW[,"Event_Type"] == "Reoperations"),]

VPO_C <- VPO[which(VPO[,"Event_Type"] == "Complications"),]
VPO_R <- VPO[which(VPO[,"Event_Type"] == "Reoperations"),]

Output <- function(x, y, k.min=10){
file <- metaprop(Events_n, N, Study_ID, data = x)

forest.meta(file, studlab = T, pooled.totals = T, bysort = F)

dev.copy2pdf(file=y, width = 11.69, height = 8.27)
print(file)
}

R code on my dropbox: https://www.dropbox.com/s/67pnfpi10qu110v/R%20code%20voor%20forrest%20en%20funnel%20plots.rtf?dl=0

[https://cfl.dropboxstatic.com/static/images/icons128/page_white_word.png]<https://www.dropbox.com/s/67pnfpi10qu110v/R%20code%20voor%20forrest%20en%20funnel%20plots.rtf?dl=0>

R code voor forrest en funnel plots.rtf<https://www.dropbox.com/s/67pnfpi10qu110v/R%20code%20voor%20forrest%20en%20funnel%20plots.rtf?dl=0>
www.dropbox.com
Shared with Dropbox
vito muggeo
# 
hi Jay,
Consult a local statistician. Statistics is not you think is (namely  
simple computations, R and probably plotting..).

regards,
vito



Jay Zola <jayjay.1988 at hotmail.nl> ha scritto:
Jay Zola
# 
Dear Vito,

Thank you for your reply. I tried to contact the statistics departement numerous times, but did not receive any reply. That is why I started to look on the internet for help.

Yours sincerely,

Jay

Verstuurd vanaf mijn iPhone
PIKAL Petr
# 
Hi

You could try to ask your question on Stackexchange, maybe somebody will answer it. This list is dedicated to helping with R code not to teaching statistics (and you did not express any direct question about the code).

Anyway, you should complain about your statistics department.

Cheers
Petr
________________________________
Tento e-mail a jak?koliv k n?mu p?ipojen? dokumenty jsou d?v?rn? a jsou ur?eny pouze jeho adres?t?m.
Jestli?e jste obdr?el(a) tento e-mail omylem, informujte laskav? neprodlen? jeho odes?latele. Obsah tohoto emailu i s p??lohami a jeho kopie vyma?te ze sv?ho syst?mu.
Nejste-li zam??len?m adres?tem tohoto emailu, nejste opr?vn?ni tento email jakkoliv u??vat, roz?i?ovat, kop?rovat ?i zve?ej?ovat.
Odes?latel e-mailu neodpov?d? za eventu?ln? ?kodu zp?sobenou modifikacemi ?i zpo?d?n?m p?enosu e-mailu.

V p??pad?, ?e je tento e-mail sou??st? obchodn?ho jedn?n?:
- vyhrazuje si odes?latel pr?vo ukon?it kdykoliv jedn?n? o uzav?en? smlouvy, a to z jak?hokoliv d?vodu i bez uveden? d?vodu.
- a obsahuje-li nab?dku, je adres?t opr?vn?n nab?dku bezodkladn? p?ijmout; Odes?latel tohoto e-mailu (nab?dky) vylu?uje p?ijet? nab?dky ze strany p??jemce s dodatkem ?i odchylkou.
- trv? odes?latel na tom, ?e p??slu?n? smlouva je uzav?ena teprve v?slovn?m dosa?en?m shody na v?ech jej?ch n?le?itostech.
- odes?latel tohoto emailu informuje, ?e nen? opr?vn?n uzav?rat za spole?nost ??dn? smlouvy s v?jimkou p??pad?, kdy k tomu byl p?semn? zmocn?n nebo p?semn? pov??en a takov? pov??en? nebo pln? moc byly adres?tovi tohoto emailu p??padn? osob?, kterou adres?t zastupuje, p?edlo?eny nebo jejich existence je adres?tovi ?i osob? j?m zastoupen? zn?m?.

This e-mail and any documents attached to it may be confidential and are intended only for its intended recipients.
If you received this e-mail by mistake, please immediately inform its sender. Delete the contents of this e-mail with all attachments and its copies from your system.
If you are not the intended recipient of this e-mail, you are not authorized to use, disseminate, copy or disclose this e-mail in any manner.
The sender of this e-mail shall not be liable for any possible damage caused by modifications of the e-mail or by delay with transfer of the email.

In case that this e-mail forms part of business dealings:
- the sender reserves the right to end negotiations about entering into a contract in any time, for any reason, and without stating any reasoning.
- if the e-mail contains an offer, the recipient is entitled to immediately accept such offer; The sender of this e-mail (offer) excludes any acceptance of the offer on the part of the recipient containing any amendment or variation.
- the sender insists on that the respective contract is concluded only upon an express mutual agreement on all its aspects.
- the sender of this e-mail informs that he/she is not authorized to enter into any contracts on behalf of the company except for cases in which he/she is expressly authorized to do so in writing, and such authorization or power of attorney is submitted to the recipient or the person represented by the recipient, or the existence of such authorization is known to the recipient of the person represented by the recipient.
1 day later
Michael Dewey
# 
Dear Jay

I am not that familiar with the meta package but it looks as though it 
does not allow you to do a meta-regression within metaprop. However 
there is a function metareg which takes the object you created with 
metaprop and allows you to add a moderator so i would try that next. By 
moderator in this case I mean your variable Treatment.
On 26/06/2017 19:19, Jay Zola wrote: