Pooling – example [Power / Sample Size]

posted by Helmut Homepage – Vienna, Austria, 2016-08-26 13:34 (3177 d 15:23 ago) – Posting: # 16583
Views: 12,644

Hi BE-proff,

❝ Let's say I have results of 3 studies (imagined):

❝ n=24, CV=0.16

❝ n=36, CV=0.28

❝ n=18, CV=0.14


❝ How to write R-code in this case to calculate CV pooled?


Adapting my code from above:

library(PowerTOST)
CV     <- c(0.16, 0.28, 0.14)
n      <- c(24, 36, 18)
source <- c("study 1", "study 2", "study 3")
CVdata <- data.frame(CV, n, rep("2x2", length(source)), source)
names(CVdata) <- names(CVdata) <- c("CV", "n", "design", "source")
print(CVdata, row.names=FALSE); print(CVpooled(CVdata, alpha=0.2), verbose=T)

You should get:

   CV  n design  source
 0.16 24    2x2 study 1
 0.28 36    2x2 study 2
 0.14 18    2x2 study 3
Pooled CV = 0.221 with 72 degrees of freedom
Upper 80% confidence limit of CV = 0.2391


Dif-tor heh smusma 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,671 registered users;
114 visitors (0 registered, 114 guests [including 9 identified bots]).
Forum time: 04:58 CEST (Europe/Vienna)

Freedom is always and exclusively
freedom for the one
who thinks differently.    Rosa Luxemburg

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5