From 3a78999fe4e8b6b757b449ebfbe59640f9c24aaf Mon Sep 17 00:00:00 2001
From: SamCH93 <samuel.pawel@gmail.com>
Date: Thu, 22 Dec 2022 15:51:19 +0100
Subject: [PATCH] additional variables
---
rsAbsence.Rnw | 18 +++++++++++++++++-
1 file changed, 17 insertions(+), 1 deletion(-)
diff --git a/rsAbsence.Rnw b/rsAbsence.Rnw
index 0764e9c..a3917dd 100755
--- a/rsAbsence.Rnw
+++ b/rsAbsence.Rnw
@@ -237,11 +237,21 @@ rpcb <- rpcbRaw %>%
effect = Effect..,
internalReplication = Internal.replication..,
effectType = Effect.size.type,
+ ## effect sizes, standard errors, p-values on original scale
+ ESo = Original.effect.size,
+ seESo = Original.standard.error,
+ lowerESo = Original.lower.CI,
+ upperESo = Original.upper.CI,
po = Original.p.value,
+ ESr = Replication.effect.size,
+ seESr = Replication.standard.error,
+ lowerESr = Replication.lower.CI,
+ upperESr = Replication.upper.CI,
+ pr = Replication.p.value,
+ ## effect sizes, standard errors, p-values on SMD scale
smdo = Original.effect.size..SMD.,
so = Original.standard.error..SMD.,
no = Original.sample.size,
- pr = Replication.p.value,
smdr = Replication.effect.size..SMD.,
sr = Replication.standard.error..SMD. ,
nr = Replication.sample.size
@@ -274,6 +284,12 @@ rpcb <- rpcbRaw %>%
rpcbNull <- rpcb %>%
## filter(po1 > 0.025) #?
filter(po > 0.05) #?
+
+## ## check HR studies
+## rpcb %>%
+## filter(effectType == "Hazard ratio") %>%
+## mutate(seESo2 = (log(upperESo) - log(lowerESo))/(2*1.96)) %>%
+## select(seESo, seESo2)
@
--
GitLab