thanks alot dr. wicklin i got on result but the results are one sample only i need 2 samples (means two different data) and when i changed number of samples in sas program i got on 1 sample also  how can i simulate two different sample such n=500 , n1=250 , n2=250. thanks again ... View more

thanks alot mr. denham i want to use categorical data with SEM then i must use link function in SEM like logit function and probit function i dont find these types of link function in sas but in other program available how can i solve this problem?? many thanks ... View more

you means unordered variables?? ... View more

dear mr. denham i dont find where is the error in sas program. ... View more

thanks alot ... View more

thanks alot doctor rick i wrote this command SHOW MODULES in sas program and got on  same results but i think the problem in loop command. 1                                                          The SAS System                         16:57 Wednesday, February 26, 2014 1          ;*';*";*/;quit;run; 2          OPTIONS PAGENO=MIN; 3          %LET _CLIENTTASKLABEL='222222222 after intiqam multi samples'; 4          %LET _CLIENTPROJECTPATH=''; 5          %LET _CLIENTPROJECTNAME=''; 6          %LET _SASPROGRAMFILE='F:\sas simulation final\sas cods program\222222222 after intiqam multi samples.sas'; 7          8          ODS _ALL_ CLOSE; 9          OPTIONS DEV=ACTIVEX; 10         GOPTIONS XPIXELS=0 YPIXELS=0; 11         FILENAME EGSR TEMP; 12         ODS tagsets.sasreport13(ID=EGSR) FILE=EGSR STYLE=HtmlBlue 12       ! STYLESHEET=(URL="file:///C:/Program%20Files/SASHome/SASEnterpriseGuide/5.1/Styles/HtmlBlue.css") NOGTITLE NOGFOOTNOTE 12       ! GPATH=&sasworklocation ENCODING=UTF8 options(rolap="on"); NOTE: Writing TAGSETS.SASREPORT13(EGSR) Body file: EGSR 13         14         GOPTIONS ACCESSIBLE; 15         proc iml; NOTE: IML Ready 16         SHOW MODULES; NOTE: No modules are currently defined. 17         18             /* P1    P2     P3       p4     p5      p6       p7      p8     p9      p10 */ 19         P = {0.35 0.50 0.20 0.10   0.10   0.30    0.20 0.25 0.15 0.20, 20              0.30 0.10 0.20 0.30   0.20   0.20    0.40 0.05 0.45    0.10, 21              0.15 0.20    0.35 0.25   0.50   0.20    0.25 0.55 0.10 0.10, 22              0.20 0.20 0.25 0.35   0.20   0.30    0.15 0.15 0.30    0.60}; 23         24         25         /* expected values and variance for each ordinal variable */ 26         Expected = Mean(P) // Var(P); 27         varNames = "X1":"X10"; 28         print Expected[r={"Mean" "Var"} c=varNames]; 29         /* test the RandMVOrd function */ 30         Delta = {1  0.45 0.4   0.35 0.3    0.25 0.2    0.15 0.1    0.05, 31                 0.45 1 0.45 0.4    0.35 0.3    0.25 0.2    0.15 0.1, 32                 0.4  0.45 1   0.45    0.4    0.35 0.3    0.25 0.2    0.15, 33                 0.35  0.4 0.45 1    0.45 0.4    0.35 0.3    0.25 0.2, 34                 0.3  0.35 0.4   0.45     1    0.45 0.4    0.35 0.3    0.25, 35                 0.25  0.3 0.35 0.4    0.45   1    0.45 0.4    0.35 0.3, 36                 0.2  0.25 0.3   0.35    0.4    0.45 1    0.45 0.4    0.35, 37                 0.15  0.2 0.25 0.3   0.35 0.4    0.45   1    0.45 0.4, 38                 0.1  0.15 0.2   0.25   0.3    0.35 0.4    0.45   1    0.45, 39                 0.05  0.1 0.15 0.2 0.25 0.3    0.35 0.4   0.45     1}; 40         /* loop approach */ 41         NumSamples = 10; 42         call randseed(54321); 43         do i = 1 to NumSamples; 44            X = RandMVOrdinal(1000, P, Delta); 45           /* do something with each sample  */ 46            end; ERROR: Invocation of unresolved module RANDMVORDINAL. statement : ASSIGN at line 44 column 4 47         /* print results */ 48         c=varNames; 49         create MVO from X[colname=c]; ERROR: Matrix X has not been set to a value. 2                                                          The SAS System                         16:57 Wednesday, February 26, 2014 statement : CREATE at line 49 column 1 49       !                                append from X; ERROR: No data set is currently open for output. statement : APPEND at line 49 column 32 49       !                                                close MVO; NOTE: Cannot close WORK.MVO; it is not open. 50         quit; NOTE: Exiting IML. NOTE: The SAS System stopped processing this step because of errors. NOTE: PROCEDURE IML used (Total process time):       real time           0.27 seconds       cpu time            0.03 seconds       51         52         GOPTIONS NOACCESSIBLE; 53         %LET _CLIENTTASKLABEL=; 54         %LET _CLIENTPROJECTPATH=; 55         %LET _CLIENTPROJECTNAME=; 56         %LET _SASPROGRAMFILE=; 57         58         ;*';*";*/;quit;run; 59         ODS _ALL_ CLOSE; 60         61         62         QUIT; RUN; 63         ... View more

thanks alot Mr. Steve this is the log 1                                                          The SAS System                            17:52 Monday, February 24, 2014 1          ;*';*";*/;quit;run; 2          OPTIONS PAGENO=MIN; 3          %LET _CLIENTTASKLABEL='222222222 after intiqam multi samples'; 4          %LET _CLIENTPROJECTPATH=''; 5          %LET _CLIENTPROJECTNAME=''; 6          %LET _SASPROGRAMFILE='F:\sas simulation final\sas cods program\222222222 after intiqam multi samples.sas'; 7          8          ODS _ALL_ CLOSE; 9          OPTIONS DEV=ACTIVEX; 10         GOPTIONS XPIXELS=0 YPIXELS=0; 11         FILENAME EGSR TEMP; 12         ODS tagsets.sasreport13(ID=EGSR) FILE=EGSR STYLE=HtmlBlue 12       ! STYLESHEET=(URL="file:///C:/Program%20Files/SASHome/SASEnterpriseGuide/5.1/Styles/HtmlBlue.css") NOGTITLE NOGFOOTNOTE 12       ! GPATH=&sasworklocation ENCODING=UTF8 options(rolap="on"); NOTE: Writing TAGSETS.SASREPORT13(EGSR) Body file: EGSR 13         14         GOPTIONS ACCESSIBLE; 15         proc iml; WARNING: The SAS/IML product with which IML (2) is associated will be expiring soon, and is currently in warning mode to indicate          this upcoming expiration. Please run PROC SETINIT to obtain more information on your warning period. NOTE: IML Ready 16         load module=_all_; NOTE: Opening storage library WORK.IMLSTOR 16       !                                        /* load the modules */ 17             /* P1    P2     P3       p4     p5      p6       p7      p8     p9      p10 */ 18         P = {0.35 0.50 0.20 0.10   0.10   0.30    0.20 0.25 0.15 0.20, 19              0.30 0.10 0.20 0.30   0.20   0.20    0.40 0.05 0.45    0.10, 20              0.15 0.20    0.35 0.25   0.50   0.20    0.25 0.55 0.10 0.10, 21              0.20 0.20 0.25 0.35   0.20   0.30    0.15 0.15 0.30    0.60}; 22         23         24         /* expected values and variance for each ordinal variable */ 25         Expected = Mean(P) // Var(P); 26         varNames = "X1":"X10"; 27         print Expected[r={"Mean" "Var"} c=varNames]; 28         /* test the RandMVOrd function */ 29         Delta = {1  0.45 0.4   0.35 0.3    0.25 0.2    0.15 0.1    0.05, 30                 0.45 1 0.45 0.4    0.35 0.3    0.25 0.2    0.15 0.1, 31                 0.4  0.45 1   0.45    0.4    0.35 0.3    0.25 0.2    0.15, 32                 0.35  0.4 0.45 1    0.45 0.4    0.35 0.3    0.25 0.2, 33                 0.3  0.35 0.4   0.45     1    0.45 0.4    0.35 0.3    0.25, 34                 0.25  0.3 0.35 0.4    0.45   1    0.45 0.4    0.35 0.3, 35                 0.2  0.25 0.3   0.35    0.4    0.45 1    0.45 0.4    0.35, 36                 0.15  0.2 0.25 0.3   0.35 0.4    0.45   1    0.45 0.4, 37                 0.1  0.15 0.2   0.25   0.3    0.35 0.4    0.45   1    0.45, 38                 0.05  0.1 0.15 0.2 0.25 0.3    0.35 0.4   0.45     1}; 39         /* loop approach */ 40         NumSamples = 10; 41         call randseed(54321); 42         do i = 1 to NumSamples; 43            X = RandMVOrdinal(1000, P, Delta); 44           /* do something with each sample  */ 45            end; ERROR: Invocation of unresolved module RANDMVORDINAL. statement : ASSIGN at line 43 column 4 46         /* print results */ 47         c=varNames; 2                                                          The SAS System                            17:52 Monday, February 24, 2014 48         create MVO from X[colname=c]; ERROR: Matrix X has not been set to a value. statement : CREATE at line 48 column 1 48       !                                append from X; ERROR: No data set is currently open for output. statement : APPEND at line 48 column 32 48       !                                                close MVO; NOTE: Cannot close WORK.MVO; it is not open. 49         quit; NOTE: Exiting IML. NOTE: Storage library WORK.IMLSTOR closed. NOTE: The SAS System stopped processing this step because of errors. NOTE: PROCEDURE IML used (Total process time):       real time           0.90 seconds       cpu time            0.04 seconds       50         51         GOPTIONS NOACCESSIBLE; 52         %LET _CLIENTTASKLABEL=; 53         %LET _CLIENTPROJECTPATH=; 54         %LET _CLIENTPROJECTNAME=; 55         %LET _SASPROGRAMFILE=; 56         57         ;*';*";*/;quit;run; 58         ODS _ALL_ CLOSE; 59         60         61         QUIT; RUN; 62         ... View more

dear dr. i still have some errors after putting the previous code can you correct it to me please 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 */ NumSamples = 2; call randseed(54321); do i = 1 to NumSamples;    X = RandMVOrdinal(500, P, Delta);   /* do something with each sample  */    end; /* print results */ c=varNames; create MVO from X[colname=c];  append from X;  close MVO; quit; many thanks ... View more

dear dr. i still have some errors after putting the previous code can you correct it to me please 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 */ NumSamples = 2; call randseed(54321); do i = 1 to NumSamples;    X = RandMVOrdinal(500, P, Delta);   /* do something with each sample  */    end; /* print results */ c=varNames; create MVO from X[colname=c];  append from X;  close MVO; quit; many thanks ... View more

this code is ok now ? 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 */ NumSamples = 10; call randseed(54321); do i = 1 to NumSamples;    X = RandMVOrdinal(1000, P, Delta);   /* do something with each sample  */ /* print results */ c=varNames; create MVO from X[colname=c];  append from X;  close MVO; quit; because i receive this warning WARNING: DO group statements were not executed due to missing END at the time of QUIT. ... View more

thanks dr. rick i want to simulate different samples like different cultures or different countries. thanks alot ... View more

hi if i want to simulate multi sample data can i just change the rand seed for each one. thanks in advance ... View more

hi when i simulate data with multivariate binary the results will be ordered or unordered binary data. thanks in advance ... View more