Pooling – example [Power / Sample Size]

posted by Helmut Homepage – Vienna, Austria, 2016-08-26 13:34 (2943 d 15:54 ago) – Posting: # 16583
Views: 11,975

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,225 posts in 4,879 threads, 1,654 registered users;
41 visitors (0 registered, 41 guests [including 7 identified bots]).
Forum time: 05:28 CEST (Europe/Vienna)

The real purpose of the scientific method is to make sure
nature hasn’t misled you into thinking you know something
you actually don’t know.    Robert M. Pirsig

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