Problem with contour(): typo
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)