Pooling – example [Power / Sample Size]

posted by Helmut Homepage – Vienna, Austria, 2016-08-26 13:34 (3091 d 16:49 ago) – Posting: # 16583
Views: 12,376

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,380 posts in 4,914 threads, 1,662 registered users;
26 visitors (0 registered, 26 guests [including 11 identified bots]).
Forum time: 05:24 CET (Europe/Vienna)

Skill is a function of chance.
It’s an intuitive best-use of chance situations.    Philip K. Dick

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