Hello, I am using sas 9.4 and I appreciate you all to help me solve these issues: 1) The rtf file is not opening in word file by default. I want the rtf file opens in word file by default. 2) only first title is printed in rtf file, the second title is not printed   ods rtf on ; ods html close; ods graphics on; ods rtf file="myfilelocation\filename.rtf "style=sasweb startpage=no bodytitle; title1'title1 '; proc print data=data1;run; title2'title2'; proc print data=data2;run; title''; ods startpage=yes; ods rtf close; ods html ;   ... View more

Hello, I appreciate for any guidance that help me understand how to check the unconditional (marginal) distribution which is assumed to be known for maximum likelihood estimation in generalized linear mixed model (proc glimmix) for multilevel repeated measures data with random effect. If I'm right one assumption for proc glimmix is that the random effect variables should follow normal distribution, so does the variable that has random effect play the role of (marginal) distribution ? Is there any article or any source that explains how to check the normality of marginal distribution and assumptions in proc glimmix and what is the concept of marginal distribution in generalized linear mixed model ? ... View more

Hello all, Any help will be appreciated to help me find a solution for my problem. My question is that if we apply proc mixed to find model for multilevel repeated measures while the response variable does not have normal distribution, will the estimates be biased? if yes how can i solve the issue? How can i find out what distribution my response has other than histogram? If my response has exponential or other distribution, how can i find the right distribution? ... View more

Hello all, I appreciate for any guide that helps me understand the concept of intercept in random statement in proc mixed? ... View more

hello all, I want to get R block matrix for some specific subjects by specifying in the R section in repeated statement in proc mixed, all the matrixes are the same unstructured matrix, but the explanation on top of them are for another subjects which are for subjects19, 20, 21, 22, 23 and 24,while R matrix should be given for subjects 13, 14, 15, 16, 17 and 18,  i do not know the reason? Any help will be appreciated to understand the reason. repeated day / subject=subjects*week*trt type=UN r=13,14,15,16,17,18; run;   ... View more

Hello all, I am looking for a source, an article or a book with examples that give me interpretation of V matrix and Z matrix in V=ZGZ'+R in mixed model with repeated measure, fixed and random effects. Any help will be appreciated.  ... View more

Hello all, Any help will be appreciated, below is the macro function sas code to detect outliers :  I do not understand these sections: " Q1 = &Q1 &val._P25;"," varL = &varL &val.L; ",  what sections in sas documentation explain theses parts  ?  How can i detect outliers using this macro function? How can i change the lower and upper limit to 10percentile and 90percentile? How can i get only outliers that are larger than 90percentile? %macro outliers(input=, var=, output= ); %let Q1=; %let Q3=; %let varL=; %let varH=; %let n=%sysfunc(countw(&var)); %do i= 1 %to &n; %let val = %scan(&var,&i); %let Q1 = &Q1 &val._P25; %let Q3 = &Q3 &val._P75; %let varL = &varL &val.L; %let varH = &varH &val.H; %end; /* Calculate the quartiles and inter-quartile range using proc univariate */ proc means data=&input nway noprint; var &var; output out=temp P25= P75= / autoname; run; /* Extract the upper and lower limits into macro variables */ data temp; set temp; ID = 1; array varb(&n) &Q1; array varc(&n) &Q3; array lower(&n) &varL; array upper(&n) &varH; do i = 1 to dim(varb); lower(i) = varb(i) - 3 * (varc(i) - varb(i)); upper(i) = varc(i) + 3 * (varc(i) - varb(i)); end; drop i _type_ _freq_; run; data temp1; set &input; ID = 1; run; data &output; merge temp1 temp; by ID; array var(&n) &var; array lower(&n) &varL; array upper(&n) &varH; do i = 1 to dim(var); if not missing(var(i)) then do; if var(i) >= lower(i) and var(i) <= upper(i); end; end; drop &Q1 &Q3 &varL &varH ID i; run; %mend;   ... View more

Hello all, Appreciate you helping me find a way to get two by two frequency table apply proc freq. I have many categorical variables, it is time consuming to enter all variables two by two.   ... View more