So I must change the random number seed for the second program. Thanks alot ... View more

thanks alot dr.rick i solve it . i want to ask you when i applied same equation X and Y in mvordinal i got on same results why? many thanks ... View more

thank you so much i want to apply same method on mvordinal with two samples but i got on same results for X AND Y.how can i solve this problem. many thanks ... View more

thanks alot dr. rick just i want to know you mean samplesize for two different samples 1=400,2=400 because i dont get on two sample size and i can not apply this command with loop approach and my program will be as below: %include "RandMVOrd.sas"; proc iml; load module=_all_;                     /* load the modules */     /* P1    P2     P3       p4     p5      p6       p7      p8     p9      p10 */ P = {0.35 0.50 0.20 0.10   0.10   0.30    0.20 0.25 0.15 0.20,      0.30 0.10 0.20 0.30   0.20   0.20    0.40 0.05 0.45    0.10,      0.15 0.20    0.35 0.25   0.50   0.20    0.25 0.55 0.10 0.10,      0.20 0.20 0.25 0.35   0.20   0.30    0.15 0.15 0.30    0.60}; /* expected values and variance for each ordinal variable */ Expected = Mean(P) // Var(P); varNames = "X1":"X10"; print Expected[r={"Mean" "Var"} c=varNames]; /* test the RandMVOrd function */ Delta = {1  0.45 0.4   0.35 0.3    0.25 0.2    0.15 0.1    0.05,         0.45 1 0.45 0.4    0.35 0.3    0.25 0.2    0.15 0.1,         0.4  0.45 1   0.45    0.4    0.35 0.3    0.25 0.2    0.15,         0.35  0.4 0.45 1    0.45 0.4    0.35 0.3    0.25 0.2,         0.3  0.35 0.4   0.45     1    0.45 0.4    0.35 0.3    0.25,         0.25  0.3 0.35 0.4    0.45   1    0.45 0.4    0.35 0.3,         0.2  0.25 0.3   0.35    0.4    0.45 1    0.45 0.4    0.35,         0.15  0.2 0.25 0.3   0.35 0.4    0.45   1    0.45 0.4,         0.1  0.15 0.2   0.25   0.3    0.35 0.4    0.45   1    0.45,         0.05  0.1 0.15 0.2 0.25 0.3    0.35 0.4   0.45     1}; /* loop approach */ Samplesize = 800; call randseed(54321);    X = RandMVOrdinal(samplesize, P, Delta); /* print results */ c=varNames; create MVO from X[colname=c];  append from X;  close MVO; quit; thanks ... View more

i think i must repeat the program one time for x and another time for y too get on different samples. thanks alot dr. rick ... View more

thanks alot dr. rick please where can i put this y = RandMVBinary(1000, p, Delta); because i have x = RandMVBinary(1000, p, Delta);     in my program , and after putting this command y = RandMVBinary(1000, p, Delta); the program will be %include "C:\Users\hp\Desktop\RandMVBinary.sas"; proc iml; load module=_all_;                     /* load the modules */ p = {0.25  0.75  0.5  0.25  0.30  0.10  0.30  0.55  0.25 0.65};       /* expected values of the X    */ Delta = {1     0.15  0.2   0.175  0.15  0.125 0.1   0.075 0.05  0.025,         0.15  1     0.225 0.15   0.175 0.075 0.125 0.1   0.075 0.05 ,         0.2   0.225 1     0.225  0.2   0.15  0.15  0.125 0.1   0.075,         0.175 0.15  0.225 1      0.225 0.2   0.175 0.15  0.125 0.1  ,         0.15  0.175 0.2   0.225  1     0.225 0.2   0.175 0.15  0.125,         0.125 0.075 0.15   0.2   0.225 1     0.225 0.125 0.175 0.1  ,         0.1   0.125 0.15   0.175 0.2   0.225 1     0.225 0.2   0.175,         0.075 0.1   0.125 0.15   0.175 0.125 0.225 1     0.225 0.2  ,         0.05  0.075 0.1   0.125  0.15  0.175 0.2   0.225 1     0.2  ,         0.025 0.05  0.075 0.1    0.125 0.1   0.175 0.2   0.2   1    }; X = RandMVBinary(1000, p, Delta); y = RandMVBinary(1000, p, Delta); /* compare sample estimates to parameters */ print p, Delta; create MVbinary from X;  append from X;  close MVB; quit; i get on just one sample i want 2 samples, you mean i must repeat the program one to x and another to y? thanks in advance THANOON ... View more

hi all anyone help me to simulate different sample size by using sas program above to get on two different mvordinal matrix. thanks in advance ... View more

thanks alot mr. denham binary data is unordered binary data but dichotomous data is ordered binary data so there is a difference between two types of data. thanks alot ... View more

hi i want to change my program from mvbinary to dichotomous is there anyone help me please?? thanks in advance ... View more

hi i still wait anyone tell me how can i get on a two different data by using the above command. regards ... View more

hi i have same problem when i want to simulate multi different sample in mvbinary. %include "RandMVBinary.sas"; proc iml; load module=_all_;                     /* load the modules */ p = {0.25  0.75  0.5  0.25  0.30  0.10  0.30  0.55  0.25 0.65};       /* expected values of the X    */ Delta = {1     0.15  0.2   0.175  0.15  0.125 0.1   0.075 0.05  0.025,         0.15  1     0.225 0.15   0.175 0.075 0.125 0.1   0.075 0.05 ,         0.2   0.225 1     0.225  0.2   0.15  0.15  0.125 0.1   0.075,         0.175 0.15  0.225 1      0.225 0.2   0.175 0.15  0.125 0.1  ,         0.15  0.175 0.2   0.225  1     0.225 0.2   0.175 0.15  0.125,         0.125 0.075 0.15   0.2   0.225 1     0.225 0.125 0.175 0.1  ,         0.1   0.125 0.15   0.175 0.2   0.225 1     0.225 0.2   0.175,         0.075 0.1   0.125 0.15   0.175 0.125 0.225 1     0.225 0.2  ,         0.05  0.075 0.1   0.125  0.15  0.175 0.2   0.225 1     0.2  ,         0.025 0.05  0.075 0.1    0.125 0.1   0.175 0.2   0.2   1    };   /* loop approach */ NumSamples = 2; call randseed(54321); do i = 1 to NumSamples;    X = RandMVOBinary(500, P, Delta);   /* do something with each sample  */    end;   /* compare sample estimates to parameters */ print p, Delta; create MVbinary from X;  append from X;  close MVB; quit; thanks in advance ... View more