SAS/R code of MSD method for comparing dissolution [Dissolution / BCS / IVIVC]
❝ ❝ 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:
- SAS/R code of MSD method for comparing dissolution Shuanghe 2014-02-10 17:05
- Wow, what a gem! Helmut 2014-02-10 17:15
- SAS/R code of MSD method for comparing dissolution arjunroy 2014-12-19 10:56
- SAS/R code of MSD method for comparing dissolution Shuanghe 2014-12-19 17:13
- SAS/R code of MSD method for comparing dissolutionarjunroy 2014-12-22 05:42
- SAS/R code of MSD method for comparing dissolution Shuanghe 2014-12-24 11:09
- SAS/R code of MSD method for comparing dissolution arjunroy 2014-12-24 11:44
- SAS/R code of MSD method for comparing dissolution Shuanghe 2015-01-08 11:16
- SAS/R code of MSD method for comparing dissolution arjunroy 2014-12-24 11:44
- SAS/R code of MSD method for comparing dissolution Shuanghe 2014-12-24 11:09
- SAS/R code of MSD method for comparing dissolutionarjunroy 2014-12-22 05:42
- SAS/R code of MSD method for comparing dissolution Shuanghe 2014-12-19 17:13
- SAS/R code of MSD method for comparing dissolution wligtenberg 2016-12-12 14:09
- SAS/R code of MSD method for comparing dissolution Shuanghe 2016-12-19 19:03
- SAS/R code of MSD method for comparing dissolution wligtenberg 2016-12-20 08:14
- SAS/R code of MSD method for comparing dissolution Shuanghe 2016-12-22 09:23
- SAS/R code of MSD method for comparing dissolution wligtenberg 2016-12-22 11:20
- SAS/R code of MSD method for comparing dissolution Shuanghe 2016-12-22 09:23
- SAS/R code of MSD method for comparing dissolution Olga 2018-04-21 07:01
- SAS/R code of MSD method for comparing dissolution wligtenberg 2016-12-20 08:14
- SAS/R code of MSD method for comparing dissolution Shuanghe 2016-12-19 19:03