[R-meta] A question about meta-regression Analysis
Viechtbauer, Wolfgang (SP)
wolfg@ng@viechtb@uer @ending from m@@@trichtuniver@ity@nl
Fri Jul 13 18:35:55 CEST 2018
Dear Hamed,
Please post in plain text format. I cannot see any plot. Either attach it as a separate image file or post it somewhere else and provide a link.
As for 2 - If you use poly(), try:
preds <- predict(res, transf = transf.ilogit, newmods = unname(poly(years, degree=2)))
Best,
Wolfgang
-----Original Message-----
From: R-sig-meta-analysis [mailto:r-sig-meta-analysis-bounces using r-project.org] On Behalf Of Hamed Kianmehr
Sent: Friday, 13 July, 2018 16:03
To: r-sig-meta-analysis using r-project.org
Subject: [R-meta] A question about meta-regression Analysis
I have some questions regarding the meta-regression analysis. I have the database on the proportion (Below bubble plot). My advisor asked me to do the polynomial regression. I have two questions:
1- Is it any polynomial regression pattern in this data?
2- I wrote below code but I couldn't get the line, Could you help me?
----------------
res<-rma(yi=y, vi=s2, data=data, method="DL",mods = ~ poly(year, degree=2))
years <- 1989:2016
preds <- predict(res, transf = transf.ilogit, newmods =cbind(years,years^2))
wi <- 0.5/sqrt(data$s2)
wi <- 0.5 + 3.0 * (wi - min(wi))/(max(wi) - min(wi))
plot(data$year, transf.ilogit(y), cex=wi,xlab="Year",ylab="Patient Expectation (%)")
lines(years, preds$pred)
More information about the R-sig-meta-analysis
mailing list