Message-ID: <a733d312-f688-5253-9af8-70fa2e4d4fc8@bebac.at>
Date: 2020-09-28T19:08:49Z
From: Helmut Schütz
Subject: Problem with contour(): typo
Dear all,
sorry, my last message contained a typo. Correct:
library(PowerTOST)
x??? <- 0.90
y??? <- 0.35
res? <- as.numeric(sampleN.TOST(theta0 = x, CV = y, design = "2x2x4",
??????????????????????????????? method = "central", details = FALSE,
??????????????????????????????? print = FALSE)[7:8])
mesh <- 28
ys?? <- unique(sort(c(y, seq(y*.8, y*1.2, length.out = mesh))))
xs?? <- unique(sort(c(x, seq(x*0.95, 1, length.out = mesh))))
z??? <- matrix(nrow = length(ys), ncol = length(xs),
?????????????? dimnames = list(paste0("y.", signif(ys, 5)),
?????????????????????????????? paste0("x.", signif(xs, 5))))
for (i in seq_along(ys)) {
? for (j in seq_along(xs)) {
??? z[i, j] <- suppressMessages(
???????????????? power.TOST(CV = ys[i], theta0 = xs[j], design = "2x2x4",
??????????????????????????? method = "central", n = res[1]))
? }
}
z??? <- z[nrow(z):1, ncol(z):1]????????? # reverse rows & columns
z[paste0("y.", y), paste0("x.", x)] == res[2] # should be TRUE
contour(xs, ys, z, nlevels = 20, las = 1, labcex = 1,
??????? xlab = "x", ylab = "y", main = paste("n =", res[1]))
abline(h = y, v = x, lty = 2)
points(x, y, col = "red", cex = 1.5)
text(x, y, labels = signif(z[paste0("y.", y), paste0("x.", x)], 6),
???? col = "red", adj = c(-0.1, 1.5))
Helmut
--
Ing. Helmut Sch?tz
BEBAC?? Consultancy Services for
Bioequivalence and Bioavailability Studies
Neubaugasse 36/11
1070 Vienna, Austria
E helmut.schuetz at bebac.at
W https://bebac.at/
F https://forum.bebac.at/