On 29/09/2020 9:16 a.m., Puetz, Benno wrote:
As I noted in my earlier post, it does - had checked that ;-)

Great!

Duncan Murdoch

It works by taking corresponding pair fo the input vectors (after possible 
recycling, as eluded by Helmut in his remark on working on only one vector) as 
needed for outer.

Thanks for the reminder, though,

        Benno

On 29. Sep 2020, at 15:12, Duncan Murdoch <murdoch.dun...@gmail.com> wrote:

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.

Reply via email to