That won't work unless power.TOST is vectorized. outer() will pass it
vectors of x and y values.
Duncan Murdoch
On 29/09/2020 8:11 a.m., Helmut Schütz wrote:
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)
______________________________________________
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.
