I ran this code and i could see all three categorical variables in three different DIFFPLOT. Thanks   proc stdize data=WORK.IMPORT method=mean out=work._ancova_stdize; var g__RFN20; run; proc glm data=work._ancova_stdize PLOTS=all; class Meat Soda Fish; model after_g__RFN20= g__RFN20 Meat Soda Fish/solution; %makeintest; %makeslest; lsmeans Soda Meat Fish / adjust=tukey pdiff alpha=.05 plots=(diffplot meanplot(cl) ); output out=work.Ancova_stats p=predicted lclm=lclm uclm=uclm / alpha=.05; ods graphics on; run; quit; proc delete data=work._ancova_stdize; run; ... View more

Thanks for the link, this is useful,   I followed your instruction for lsmeans but now DIFFPLOT disappear, my code is in the below   proc stdize data=WORK.IMPORT method=mean out=work._ancova_stdize; var _g__Lactococcus; run; proc glm data=work._ancova_stdize PLOTS=all; class Meat Soda Fish; model _after_g__Lactococcus= _g__Lactococcus Meat Fish Soda/solution; %makeintest; %makeslest; lsmeans Meat _g__Lactococcus / adjust=tukey pdiff alpha=.05 plots=(diffplot meanplot(cl) ); output out=work.Ancova_stats p=predicted lclm=lclm uclm=uclm / alpha=.05; ods graphics on; run; quit; proc delete data=work._ancova_stdize; run; ... View more

Do you @PaigeMiller mean this attached plot? If yes, i got only for one categorical variable (meat). Please can you suggest how can i add the rest of the two other categorical variables?       ... View more

Hi @PaigeMiller , Thanks for the suggestion, i am sorry about the image, Here is my code, when i ran this code, i got diagnostic, mean and residual images, Since i am new in this analysis and running three different categorical variables (diets) and hoping to see their effect in a figure before and after bacterial (lactococcus) treatment.   Please let me know if it is possible or not. Meanwhile, i will do more research to show you an example figure (as of now i haven't found for three categorical), so that might be helpful. Thanks     proc glm data=work._ancova_stdize PLOTS=ALL; class Meat Soda Fish; model _after_g__Lactococcus= _g__Lactococcus Meat Fish Soda/solution; %makeintest; %makeslest; lsmeans Meat / adjust=tukey pdiff alpha=.05 plots=(diffplot meanplot(cl) ); output out=work.Ancova_stats p=predicted lclm=lclm uclm=uclm / alpha=.05; ods graphics on; run; quit; proc delete data=work._ancova_stdize; run;   ... View more

Hi @PaigeMiller , I dont know exactly how to represent the significant changes in a graph,   Please see my attached code   I am able to get only diagnostic plots, even after choosing in options section for all kind of plots, Please let me know if something i am doing wrong, Thanks   ... View more

Hi @Rick_SAS @PaigeMiller  I am doing ANCOVA with three categorical variable. But i am confused about how to represent the results in a graph/plot,   a box/scatter, i dont know how it works, this is my first time doing using SAS studio, Please someone can guide me on this, i would appreciate the support, Thanks,     ... View more

Hi @Rick_SAS ,   Thanks for the tips, I added ADD =  -1.681563224 but i couldn't get the same intercept estimate (please see the image). By the way, how do you calculate the value for  ADD? I asked my colleague for the details, how i can get the code, Just for curiosity, does it possible to get code from SAS enterprise guide version?   Thanks, Khem   ... View more

Hello @Rick_SAS ,   Thanks for the explanation, I tried to understand your example and compare it with my code, Looks like my code is similar to your second example (please correct me if i am wrong), If I want to get the same intercept estimate like my colleague (they used SAS enterprise and i don't have their code), do i need to remove the method part?  AS you defined X= rand...., i cant do guess for CRP_A x = rand("Normal", 20, 1); or add = -20 I am sorry, but i don't know what change to i need to make? i am still a newbie to SAS.   Please guide me, what change is needed in my code? Thanks, Khem ... View more

Hi @PaigeMiller ,   Sorry to bother you again, is it possible to reproduce your model for ANCOVA in R? If yes, can you suggest to me how? what should i add for solution (SAS studio) option in R?   Thanks, Khem ... View more

Oh, Here i attached my codes, if it may help you,   Thanks, Khem ... View more

Hi @PaigeMiller  Thanks for the suggestion, it worked, Here is the image   I got the same estimate and Pr>t values for all three diets, Only one doubt, the estimate of intercept is not matching, My collaborator had 1.26 (in previous post figure) but mine is 2.60. Does it matter? do i need to make any changes? Thanks again, Khem ... View more

Hi @PaigeMiller ,   Second analysis i did,  based on your first suggestion but when i ran your equation in SAS studio, I didn't get the intercept, estimates and Pr>t values.  do i need to add something into your previous equation?   Thanks, Khem ... View more

Hi @PaigeMiller  I am really sorry, I uploaded the wrong file of my collaborator, Here i attached the right file of my collaborator, along with my same results (after your suggestion, black&whilte figure) and data from my previous post,       Here I got similar values for TypeIII and Pr>F like my collaborator, now i need to get same intercept, estimates, Pr>t values.   Sorry again, Khem     ... View more

Hi @PaigeMiller , Okay, Here i dropped the screen images for the results of my collaborator (first one with blue font) and mine (image with black&white font). My results are based on your previous suggestion. I got similar values for TypeIII and Pr>F like my collaborator, which is great. Your suggestion is in the perfect and right direction for my data. Please kindly suggest me how can I get similar intercept and Pr>t values too? for the reference, I also put the image of my datasheet.   Note: Here In the data zero (o) means "without diet", and one (1) means "with diet"   Thanks, Khem ... View more

Hi @PaigeMiller  Thanks for the quick response, I ran again with your suggested code, I could up to typeIII values and F values, this is great, Exact same like my collaborator, I attached my results too (SAS_studio...pdf) But i am not able to get the estimate and t and Pr>t values (please check the attached file [ancova_cassie...] for results from a collaborator) like my collaborators had, Please can you help me on this too, They used SAS enterprise guide, so i am not sure, if they will have codes for them, Thanks, Khem   ... View more