Dear Benno,
THX, you made my day! Case closed.
Helmut
Puetz, Benno wrote on 2020-09-29 13:14:
I would assume the following snippet does what you want - note the use
of outer with anonymous function wrapping powerTOST:
z <- outer(xs, ys, function(x, y)power.TOST(CV = y, theta0 = x, design
= "2x2x4", method = "central", n = res[1]))
contour(xs, ys, z)
--
Ing. Helmut Schütz
BEBAC – Consultancy Services for
Bioequivalence and Bioavailability Studies
Neubaugasse 36/11
1070 Vienna, Austria
E helmut.schu...@bebac.at
W https://bebac.at/
F https://forum.bebac.at/
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
