Bioequivalence and Bioavailability Forum

The ultimate crackpot iteration! [Two-Stage / GS Designs]

posted by ElMaestro  – Denmark, 2017-08-20 21:32 (2863 d 03:09 ago) – Posting: # 17728
Views: 29,933

Hi Hötzi,

we are not on the same page. Again
I am not sure I get it all, but I think with this code you are calling the same shit over and over again and that is completely unnecessary.
In the beginning of everything, before you initiate the simulation, literally before you make use of the nsims object, you read out the assumed GMR and target power, and on basis of that you can define an array of critical CV's which you initialise with the correct values. Once and for all. Not for every call where you need a sample size.

Note that if you initialise the array of critical values with syntax such as
vecCritCVs<<- rep(0, 400)

or however many values you need then vecCV will be public (I think), due to "<<-".
You can then access the array of critical CV's from within the SampleN.Stg2 function.

We might end up with something like:

sampleN.Stg2.EM <- function(CV)
{
  #all other variable info is already accumulated thin the critical values, so:
  Nps=1
  while (vecCritCVs[Nps]<=CV) Nps=Nps+1
  return(Nps*2) 
}

Can you see what I mean now?

Complete thread:

• Initial sample size guess for the Potvin methods ElMaestro 2017-08-19 15:04 
  • Initial sample size guess for the Potvin methods Helmut 2017-08-19 16:06
    • Initial sample size guess for the Potvin methods ElMaestro 2017-08-19 16:14
      • Initial sample size guess for the Potvin methods Helmut 2017-08-19 17:12
        • Initial sample size guess for the Potvin methods ElMaestro 2017-08-19 17:17
          • Confuse-a-Cat Helmut 2017-08-19 17:33
            • Confuse-a-Cat ElMaestro 2017-08-19 17:56
      • loop ↔ vectorized ↔ direct Helmut 2017-08-20 14:40
        • loop ↔ vectorized ↔ direct ElMaestro 2017-08-20 15:22
          • loop ↔ vectorized ↔ direct Helmut 2017-08-20 16:23
            • loop ↔ vectorized ↔ direct ElMaestro 2017-08-20 17:22
  • The n ext crackpot iteration ElMaestro 2017-08-19 20:04
    • The n ext crackpot iteration Helmut 2017-08-20 02:20
  • The ultimate crackpot iteration! ElMaestro 2017-08-20 14:35
    • The ultimate crackpot iteration! Helmut 2017-08-20 15:11
      • The ultimate crackpot iteration! ElMaestro 2017-08-20 15:28
        • The ultimate crackpot iteration! Helmut 2017-08-20 16:06
          • The ultimate crackpot iteration! ElMaestro 2017-08-20 16:15
            • The ultimate crackpot iteration! Helmut 2017-08-20 18:58
              • The ultimate crackpot iteration!ElMaestro 2017-08-20 19:32
                • Suggested code ElMaestro 2017-08-21 18:13
                  • Nitpicker! Helmut 2017-08-22 13:33
                    • Nitpicker! ElMaestro 2017-08-22 17:27
                      • Nitpicker! Helmut 2017-08-22 17:49
                        • Nitpicker! ElMaestro 2017-08-22 17:59
                          • Nitpicker! Helmut 2017-08-22 19:15
                            • Benchmark code ElMaestro 2017-08-22 22:29
                              • Benchmark code Helmut 2017-08-23 01:48