Mahalanobis distance - R code [Dissolution / BCS / IVIVC]

posted by d_labes  – Berlin, Germany, 2011-04-29 17:14 (5104 d 03:43 ago) – Posting: # 6973
Views: 22,444

Dear Ace,

there is a function mahalanobis() in the stats package, i.e in the base installation of R.
But unfortunately it gives no confidence interval out of the box.
And this is the tricky part I think and for that I wish you a nice weekend :-D.

Beside the Tsong paper which is cryptic in this respect and doesn't give ready to program advice have a look at this online resource for another method:
B. Reiser
"CONFIDENCE INTERVALS FOR THE MAHALANOBIS DISTANCE"
Needs some postscript reader.

Have also a look at the function ginv() of the MASS package, which implements the Moore-Penrose inverse which is necessary in calculating the Mahalanobis distance if the estimated variance-covariance matrix is not invertible.

If you are interested, I have some code ruins.

Regards,

Detlew

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,694 registered users;
25 visitors (0 registered, 25 guests [including 8 identified bots]).
Forum time: 20:58 CEST (Europe/Vienna)

Do not put your faith in what statistics say until you have carefully
considered what they do not say.    William W. Watt

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