CVs (R, T) from pooled CV [RSABE / ABEL]

posted by d_labes  – Berlin, Germany, 2015-07-02 11:09 (3623 d 12:50 ago) – Posting: # 15014
Views: 31,386

Dear Helmut!

library(PowerTOST)

CV <- 0.07

sep2s <- function(CV, ratio=1) { # split CV to s-ratio

  sp  <- CV2mse(CV) # pooled s²

  # (1) (s²T + s²R)/2 = s²

  # (2) sT/sR = ratio

  # solve (1+2) for sT, sR

  swt <- sqrt(2)*sqrt(sp)*ratio/sqrt(ratio^2+1)

  swr <- sqrt(2)*sqrt(sp)/sqrt(ratio^2+1)

  r   <- se2CV(c(swt, swr))

  return(r)

} ...


Sometimes it would be helpful to RTFM :-D.
Try
library(PowerTOST)
?CVp2CV
.

Note that the ratio is here the ratio of variances.

Regards,

Detlew

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,682 registered users;
26 visitors (0 registered, 26 guests [including 17 identified bots]).
Forum time: 00:00 CEST (Europe/Vienna)

EMEA. The European Medicines Evaluation Agency.
The drug regulatory agency of the European Union.
A statistician-free zone.    Stephen Senn

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