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) >
As I noted in my earlier post, it does - had checked that ;-) 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.duncan at 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) >
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.duncan at 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) >> >