
R version 4.2.2 (2022-10-31 ucrt) -- "Innocent and Trusting"
Copyright (C) 2022 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64 (64-bit)

R是自由软件，不带任何担保。
在某些条件下你可以将其自由散布。
用'license()'或'licence()'来看散布的详细条件。

R是个合作计划，有许多人为之做出了贡献.
用'contributors()'来看合作者的详细情况
用'citation()'会告诉你如何在出版物中正确地引用R或R程序包。

用'demo()'来看一些示范程序，用'help()'来阅读在线帮助文件，或
用'help.start()'通过HTML浏览器来看帮助文件。
用'q()'退出R.

[原来保存的工作空间已还原]

> library(xlsx)
> library(meta)
载入需要的程辑包：metadat
Loading 'meta' package (version 7.0-0).
Type 'help(meta)' for a brief overview.
Readers of 'Meta-Analysis with R (Use R!)' should install
older version of 'meta' package: https://tinyurl.com/dt4y5drs
Warning messages:
1: 程辑包‘meta’是用R版本4.2.3 来建造的 
2: 程辑包‘metadat’是用R版本4.2.3 来建造的 
> metagen <- read.xlsx("review.xlsx",13)
> se <-(metagen $ uci - metagen $ lci)/ (2*1.96)
> lnhr<- log(metagen$ hr)
> metagen <- cbind(metagen,lnhr,se)
> m <- metagen(TE=metagen$lnhr,seTE=metagen$se,studlab=paste(metagen$study,metagen$year,sep = "/"),sm="HR",backtransf=TRUE,byvar=metagen$measurement)
> m
Number of studies: k = 16

                         HR           95%-CI     z  p-value
Common effect model  1.8818 [1.6940; 2.0905] 11.78 < 0.0001
Random effects model 1.7681 [1.4286; 2.1883]  5.24 < 0.0001

Quantifying heterogeneity:
 tau^2 = 0.0568 [0.0000; 0.1211]; tau = 0.2383 [0.0000; 0.3479]
 I^2 = 45.1% [1.3%; 69.4%]; H = 1.35 [1.01; 1.81]

Test of heterogeneity:
     Q d.f. p-value
 27.30   15  0.0264

Results for subgroups (common effect model):
                         k     HR           95%-CI    Q  I^2
measurement = Dipstick   4 1.5359 [1.2641; 1.8660] 1.36 0.0%
measurement = 3/HP       5 2.3456 [2.0209; 2.7226] 0.18 0.0%
measurement = 2/HP       2 2.7199 [1.6094; 4.5968] 0.00 0.0%
measurement = 5/HP       5 1.2828 [0.9942; 1.6552] 2.60 0.0%

Test for subgroup differences (common effect model):
                   Q d.f.  p-value
Between groups 23.15    3 < 0.0001
Within groups   4.14   12   0.9807

Results for subgroups (random effects model):
                         k     HR           95%-CI tau^2 tau
measurement = Dipstick   4 1.5359 [1.2641; 1.8660]     0   0
measurement = 3/HP       5 2.3456 [2.0209; 2.7226]     0   0
measurement = 2/HP       2 2.7199 [1.6094; 4.5968]     0   0
measurement = 5/HP       5 1.2828 [0.9942; 1.6552]     0   0

Test for subgroup differences (random effects model):
                   Q d.f.  p-value
Between groups 23.15    3 < 0.0001

Details on meta-analytical method:
- Inverse variance method
- Restricted maximum-likelihood estimator for tau^2
- Q-Profile method for confidence interval of tau^2 and tau
> 
