Thank you @Rick_SAS , for your patient answer and advice! Yes, the main problem is making sure the Skew and Kurtosis match up with each other , withing a legitimate distribution.  Very kind of you to help! ... View more

Thank you, @Rick_SAS ! Yes, I suppose I have to find a nonlinear solver to get me the shape values a and b, given Skew and Kurtosis. Not sure if I need a different post but what if I choose to go with a mixture of Normalsan like in the image below, but truncated on the left? I am going for something like below, but on a bounded interval (on one side only). I am looking to try cross different Skew and Kurtosis values : Skew in 0 to 3 and Kurt in (3,5,7) (From Allison J. Ames, Brian C. Leventhal & Nnamdi C. Ezike (2020) Monte Carlo Simulation in Item Response Theory Applications Using SAS, Measurement: Interdisciplinary Research and Perspectives, 18:2, 55-74, DOI: 10.1080/15366367.2019.1689762  https://doi.org/10.1080/15366367.2019.1689762 )   I found this from your blog:  https://blogs.sas.com/content/iml/2019/04/29/normal-mixture-distribution-sas.html and Implement the truncated normal distribution in SAS - The DO Loop ? Sorry, too many questions. Any advice would help get me started. Thank you again! ... View more

Thanks, this comes closest to my requirements. I originally was working with datasets with common column names, but later realised I may have to generalize this to when there are no common column names. I can, however, create a common variable - a  row number column. Or, it would seem, from your reply and @PaigeMiller , I can make the BY variable  in PROC Excel just be the datasetname (where datasetname is the data set that combines all the original datasets). Thanks again.     ... View more

Thanks for the example. Quick followup - what if I have multiple datasets that need to output to different Excel sheets? So, my data is in different datasets - dat1,dat2,...,datN and each of these need to output to their own separate tab in one single Excel file. I can't think of a single common BY variable (except perhaps the Row Number). How could I modify your example code - any suggestions would help! ... View more

What I meant is that there are several different datasets to export. How can I manage this with a single BY group ? The column names are the same but the name of each dataset is different. ... View more

Thanks! Do you have a small example to share? I would like to avoid using a Macro if possible! ... View more

Yes, thanks! I was using "PUT" instead of "LET"! ... View more

Here is the Log: 2505 dm "log; clear; "; 2506 dm 'odsresults; clear'; 2507 options mprint; 2508 /* */ 2509 %macro getBetas(); 2510 %let varnamelist = SLSS3 PANASn3 PANASp3 /*PANASpP3 PANASnP3 CES3 CHSa3 CHSp3 CHStot3 CASSSC3 2510! CASSSP3 CASSST3*/ ; 2511 %put &varnamelist; 2512 %let end=%sysfunc(countw(&varnamelist)); 2513 %put &end; 2514 %do i = 1 %to &end; 2515 %put out= %scan(&varnamelist., &i, " "); 2516 proc print data=EstCovWithBL_&out.; 2517 run; 2518 %end; 2519 %mend; 2520 2521 ods tagsets.excelxp style=sasweb 2522 file="C:\Users\radhi\OneDrive - University of Delaware - o365\Desktop\Dr Suldo 2522! IES\estbetas.xls"; NOTE: Writing TAGSETS.EXCELXP Body file: C:\Users\radhi\OneDrive - University of Delaware - o365\Desktop\Dr Suldo IES\estbetas.xls NOTE: This is the Excel XP tagset (Compatible with SAS 9.1.3 and above, v1.131, 04/23/2015). Add options(doc='help') to the ods statement for more information. 2523 %getBetas; SLSS3 PANASn3 PANASp3 3 out= SLSS3 NOTE: Writing HTML Body file: sashtml4.htm WARNING: Apparent symbolic reference OUT not resolved. NOTE 137-205: Line generated by the invoked macro "GETBETAS". 1 proc print data=EstCovWithBL_&out.; run; - 22 ERROR 22-322: Syntax error, expecting one of the following: ;, (, BLANKLINE, CONTENTS, DATA, DOUBLE, GRANDTOTAL_LABEL, GRANDTOT_LABEL, GRAND_LABEL, GTOTAL_LABEL, GTOT_LABEL, HEADING, LABEL, N, NOOBS, NOSUMLABEL, OBS, ROUND, ROWS, SPLIT, STYLE, SUMLABEL, UNIFORM, WIDTH. NOTE: Line generated by the invoked macro "GETBETAS". 1 proc print data=EstCovWithBL_&out.; run; - 200 ERROR 200-322: The symbol is not recognized and will be ignored. WARNING: Apparent symbolic reference OUT not resolved. ERROR: File WORK.ESTCOVWITHBL_.DATA does not exist. MPRINT(GETBETAS): proc print data=EstCovWithBL_&out. run; NOTE: The SAS System stopped processing this step because of errors. NOTE: PROCEDURE PRINT used (Total process time): real time 0.43 seconds cpu time 0.15 seconds out= PANASn3 WARNING: Apparent symbolic reference OUT not resolved. NOTE 137-205: Line generated by the invoked macro "GETBETAS". 3 proc print data=EstCovWithBL_&out.; run; - 22 ERROR 22-322: Syntax error, expecting one of the following: ;, (, BLANKLINE, CONTENTS, DATA, DOUBLE, GRANDTOTAL_LABEL, GRANDTOT_LABEL, GRAND_LABEL, GTOTAL_LABEL, GTOT_LABEL, HEADING, LABEL, N, NOOBS, NOSUMLABEL, OBS, ROUND, ROWS, SPLIT, STYLE, SUMLABEL, UNIFORM, WIDTH. NOTE: Line generated by the invoked macro "GETBETAS". 3 proc print data=EstCovWithBL_&out.; run; - 200 ERROR 200-322: The symbol is not recognized and will be ignored. WARNING: Apparent symbolic reference OUT not resolved. ERROR: File WORK.ESTCOVWITHBL_.DATA does not exist. MPRINT(GETBETAS): proc print data=EstCovWithBL_&out. run; NOTE: The SAS System stopped processing this step because of errors. NOTE: PROCEDURE PRINT used (Total process time): real time 0.01 seconds cpu time 0.01 seconds out= PANASp3 WARNING: Apparent symbolic reference OUT not resolved. NOTE 137-205: Line generated by the invoked macro "GETBETAS". 5 proc print data=EstCovWithBL_&out.; run; - 22 ERROR 22-322: Syntax error, expecting one of the following: ;, (, BLANKLINE, CONTENTS, DATA, DOUBLE, GRANDTOTAL_LABEL, GRANDTOT_LABEL, GRAND_LABEL, GTOTAL_LABEL, GTOT_LABEL, HEADING, LABEL, N, NOOBS, NOSUMLABEL, OBS, ROUND, ROWS, SPLIT, STYLE, SUMLABEL, UNIFORM, WIDTH. NOTE: Line generated by the invoked macro "GETBETAS". 5 proc print data=EstCovWithBL_&out.; run; - 200 ERROR 200-322: The symbol is not recognized and will be ignored. WARNING: Apparent symbolic reference OUT not resolved. ERROR: File WORK.ESTCOVWITHBL_.DATA does not exist. MPRINT(GETBETAS): proc print data=EstCovWithBL_&out. run; NOTE: The SAS System stopped processing this step because of errors. NOTE: PROCEDURE PRINT used (Total process time): real time 0.01 seconds cpu time 0.01 seconds 2524 ods tagsets.excelxp close; ... View more

Thanks so much! The values per se of the 'time' variable aren;t important - I'm just interested in the correlation over time within a Phase for a specific subject. ... View more

I have an ABAB design with repeated time observations within a single phase (A - baseline phase, B- treatment phase. Each subject gets both phases, repeated many times, about 10 to 30 per subject per phase). There are 9 subjects in total. What is the best wyay to model this data using a mixed model structure - I was trying to see if each subject is a 'unit', and it makes sense to have random effects for the Intercept and maybe the Phase for each subject. But there is also the correlation amongst 'outcome' scores for a specific subject within a Phase - surely that has to be something like AR(1) or CS? Here is data from one single subject 'case A1': case phase time outcome case A1 1 1 7 case A1 1 2 9 case A1 1 3 8 case A1 1 4 6 case A1 1 5 7 case A1 1 6 4 case A1 1 7 5 case A1 1 8 10 case A1 1 9 2 case A1 1 10 0 case A1 2 15 3 case A1 2 16 8 case A1 2 17 8 case A1 2 18 6 case A1 2 19 10 case A1 2 20 10 case A1 2 21 10 case A1 2 22 8 case A1 2 23 3 case A1 2 24 4 ... View more

Thanks so much for the reply! I will try increasing the interations. Meanwhile, I want to know how I can get G-side effects for the Intercept in the model shown below(I have specified a random Intercept - but number of columns in Z=0 per the SAS output - what is wrong with my syntax???): title "Model with only R-side CS errors"; proc glimmix data=model_dat2; class case phase; model y = flag_2 flag_3 flag_4/solution dist=binomial link=logit cl; random Intercept/sub=case type=cs residual; run; OUTPUT: Dimensions R-side Cov. Parameters 2 Columns in X 4 Columns in Z per Subject 0 Subjects (Blocks in V) 9 Max Obs per Subject 59 ... View more