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

posted by ElMaestro  – Belgium?, 2009-09-28 19:12  – Posting: # 4261
Views: 27,497

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,694 posts in 4,181 threads, 1,355 registered users;
online 11 (0 registered, 11 guests [including 6 identified bots]).
Forum time (Europe/Vienna): 13:22 CEST

A central lesson of science is that to understand complex issues
(or even simple ones), we must try to free our minds of dogma and
to guarantee the freedom to publish, to contradict, and to experiment.
Arguments from authority are unacceptable.    Carl Sagan

BEBAC Ing. Helmut Schütz