additional variables

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)
 @