Hi Sir, The PROC CATMOD is a good method to run multinomial logit models (response variable with more than 2 categories). The estimation is usually based on maximum lilihood and a ANOVA table is produced. I am wondering if the AVOVA table is estimated based on Type I Sums of Squares or Type III. Just cannot find any reference for this. ... View more

Hi, I used PROC LOGISTIC to run a cumulative logit model for ordinal response. The dependent variable has 3 levels:  group=1, group=2, and group=3. For the cumulative logit model here, the predicition is for the probability in the lower category (group 1 vs group 2, 3,  and group 1, 2 vs group 3). I don't undertand the meaning of the intercepts, which are produced as: _LABEL_ Intercept: group=3     2.325 Intercept: group=2     5.273 I am very confused with the meaning of intercepts here. Can you please help? ... View more

I have a sample with one dependent dummary variable y (0, 1). The data has been sorted. Now the question is how I can get easily get a summary (e.g. mean, freq) for the top 5%, 10%, 15%, 20%, 25%, ... 95% sample. I hope there is a way to do this easily, rather than use separate steps to produce the results. Any procedure to do this, PROC UNIVARIATE, PROC FREQ, PROC MEANS?  ... View more

Hi, I used the program below with the purpose of setting age group 15 as the reference group for comparison. But the result shows that age group 9 has been set as the reference group. Don't understand where is the problem?:smileyconfused: PROC LOGISTIC DATA=work.test;   CLASS age (PARAM=REF REF='15') sex /param=GLM;   MODEL score=age sex /AGGREGATE SCALE=NONE;   OUTPUT OUT=work.test02 pred=pred reschi=reschi; RUN; Result: age1 age10 age11 age2 age3 age4 age5 age6 age7 age8 age9: 0 ... View more

Hi sir, I just run a macro as below. The macro is used to merge multiple datasets. The %else statement needs two consecutive semicolons ';' here. Otherwise, SAS reported that it cannot find the data. The log message shows that: DATA work.model; MPRINT(DOMERGE):   SET work.model work.data05 RUN; So it does miss a semicolon here. /*The macro program as below*/ %let num=16; %macro domerge;     %do i=1 %to #       DATA work.model;      %if &i=1 %then %do;          SET work.data01 work.data02;          %end;      %else SET work.model work.data%sysfunc(putn(&i+1.,z2.));;      RUN;   %end; %mend domerge; %domerge ... View more

I took a few very tedious steps to first create a few random variables, then sort the data based on the random variables. After that, I used the data set to select the first 200 responses. As I need to get 100 random samples to test a model, I would have to repeat this process for 100 times. Very poor idea and labor work. I am curious to learn if there is some code to automatically create 100 random samples from the orginal dataset. Something like: %let i=100; data work.data01 work.data02 work.data03 ... work.data&i;   set work.fulldata; /*fulldata has 1000 observations*/   do j=1 to 100;   ...output to different data sets with 200 observations;   end; run; :smileyconfused: ... View more