All the important answers are here! [R for BE/BA]

posted by ElMaestro  – Denmark, 2009-09-28 19:12  – Posting: # 4261
Views: 26,060

Dear HS,

» Oh, I took these boots on many years ago. And they got muddy almost immediatelly. I tried to get a numerical approximation of the cdf of the t-distribution which lead me to the gamma-function and headaches.
» Finally I ended up with numerical approximations.[1,2]

Let me highfive you for this. I did the same, in C, and it took me quite some time, only to find out later that C actually comes with a beautiful function called lgamma which does it and returns the logarithm. I wasted a good part of my youth reinventing the wheel it seems :angry:

» BTW, nitpicking Diletti’s Table 1[6] (power 70%, T/R=1, CV=7.5%) and nQuery 7 give n=4, whilst my R code, StudySize 2.0.1, and FARTSSIE 1.6 come up with n=6. Now what? Go with a democratic vote of 3:2 for 6?

Two solutions:
  1. Ask dlabes
  2. Get your power from a program that bases its evaluation on brute force rather than on numerical integration.

Complete thread:

 Mix view
Bioequivalence and Bioavailability Forum |  Admin contact
19,486 posts in 4,135 threads, 1,334 registered users;
online 14 (0 registered, 14 guests [including 7 identified bots]).
Forum time (Europe/Vienna): 22:58 CEST

If you don’t know anything about computers,
just remember that they are machines that do exactly what you tell them
but often surprise you in the result.    Richard Dawkins

BEBAC Ing. Helmut Schütz