Package MBESS and Power curiosity [🇷 for BE/BA]

posted by ElMaestro  – Denmark, 2009-10-17 14:26 (6079 d 14:42 ago) – Posting: # 4370
Views: 45,219

Ahoy d_labes,

❝ BTW: ElMaestro, can you tell me how Apfelstrudel has found his way to Gent/Belgium :confused:?


When I made port in Bremerhafen in 1953 I decided to give Apfelstrudel a try; it turned out to be a fully acceptable alternative to the famous Vlaamse Friet and with a pleasantly matching (lack of) nutritional value.

Regarding only 50000 simulations:
You can get the famous software EFG 2.01 from HS. It will do millions of sims in a few seconds. HS - put it up for download?

Here's an update to give you more decimals, if you like:

Apfelstrudel<-function(t, df, ta, tb, ncp, ACCURACY)
{
  L1 = sqrt(2*pi);
  P=0;
  step = (tb-ta)/ACCURACY;
  foo=ta;
  for (i in 1:ACCURACY)
    {
      y1 = EatMyShorts(foo, df, t, ncp);
      y2 = EatMyShorts(foo+0.5*step, df,t, ncp);
      y3 = EatMyShorts(foo+step, df,t, ncp);
      P=P+step*(y1+y2*4+y3)/6; ## homer simpson's meffud!
      foo=foo+step;
    }
  L1*P
}



Best regards
ElMaestro - sailor and simulant.


Edit: BBCoded. [Helmut]

Complete thread:

UA Flag
Activity
 Admin contact
23,653 posts in 4,991 threads, 1,570 registered users;
124 visitors (0 registered, 124 guests [including 30 identified bots]).
Forum time: 05:08 CEST (Europe/Vienna)

In theory, there is no difference between theory and practice.
But, in practice, there is.    Jan L.A. van de Snepscheut

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5