SAS/R code of MSD method for comparing dis­so­lution [Dissolution / BCS / IVIVC]

posted by arjunroy – India, 2014-12-22 06:42 (3407 d 01:16 ago) – Posting: # 14135
Views: 18,442
(edited by Ohlbe on 2014-12-22 09:39)

❝ ❝ 1. R code is working fine,

❝

❝ So I guess you have nt=nr otherwise the code will give wrong result. I noticed it some time ago so I'll modify the code when I have some time.


Thank you

❝ ❝ though SAS IML v9.1.3 did not, got problem in calculating Cov matrix. Execution error noted in log despite it can be done manually.

❝

❝ Do you add the COV function before the run? As I mentioned in the post, for SAS 9.2 and older version there's no COV function so you should add it manually. The example of COV I gave only work for no missing value which in most of the case it should be. But if not, well, you have to adapt the code.


Here is SAS IMLv9.1.3 log error for your perusal.
proc iml;
NOTE: IML Ready
21   use one;
21 !          read all into one;
21 !                             close one;
22
23   mt = one[1:6,];
24   mr = one[7:12,];
25
26   START COV(m);
27       n = NROW(m);
27 !                               /* assume no missing values */
28       d = m - m[:,];
29
30       /* subtract mean to center the data */
31       RETURN((d`*d)/(n-1));
32     FINISH;
NOTE: Module COV defined.
33
34   nt=6;
34 !       nr=6;
34 !             p=8;
34 !                  dg=15;
35
36   tbar = (mt[:,])`;
36 !                               /* mean profile of T */
37   rbar = (mr[:,])`;
37 !                      /* mean profile of R */
38
39
40   mdiff = tbar-rbar;
40 !                         /* mean difference between T&R */
41
42   st = Cov(mt);
ERROR: (execution) Matrices do not conform to the operation.
 operation : - at line 28 column 11
 operands  : M, _TEM1001

M      6 rows      8 cols    (numeric)

     19.99      36.7     47.77     55.08     65.69     81.37     92.39      97.1
     22.08     39.29     49.46     56.79     67.22     82.42     89.93     95.62
     21.93     38.54     47.76     55.14     65.25     83.49     90.19     95.62
     22.44     39.46     49.72     58.67     69.21     84.93     94.12     95.51
     25.67     42.35     52.68     59.71     71.51     86.61      93.8      96.7
     26.37     41.34     51.01     57.75     69.44      85.9     94.45     98.07

_TEM1001      1 row       8 cols    (numeric)

     23.08 39.613333 49.733333     57.19 68.053333     84.12     92.48 96.436667

 statement : ASSIGN at line 28 column 5
 traceback : module COV at line 28 column 5

NOTE: Paused in module COV.
42 !                               /* variance-covariance matrix of T & R */
43   sr = Cov(mr);
ERROR: (execution) Matrix has not been set to a value.

Thank you so much.

Edit: Full quote removed. Please delete everything from the text of the original poster which is not necessary in understanding your answer; see also this post! [Ohlbe]

Complete thread:

UA Flag
Activity
 Admin contact
22,988 posts in 4,825 threads, 1,654 registered users;
78 visitors (0 registered, 78 guests [including 1 identified bots]).
Forum time: 08:59 CEST (Europe/Vienna)

The whole purpose of education is
to turn mirrors into windows.    Sydney J. Harris

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