Pooling – example [Power / Sample Size]

posted by Helmut Homepage – Vienna, Austria, 2016-08-26 13:34 (2661 d 18:06 ago) – Posting: # 16583
Views: 11,127

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
22,822 posts in 4,786 threads, 1,633 registered users;
19 visitors (0 registered, 19 guests [including 5 identified bots]).
Forum time: 06:41 CET (Europe/Vienna)

A scientist’s aim in a discussion with his colleagues
is not to persuade, but to clarify.    Leo Szilard

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