Ok makes sense, thank you very much for the great help! ... View more

proc import datafile='C:\Users\xiakeshan\Documents\Downloads\dtr.xlsx' out=have dbms=xlsx replace; run; data one; set have; if strip(Variety)='GEM' then Variety=' GEM'; DTR=input(DTR__days_to_ripe_,best.); run; proc sgpanel data=one ; styleattrs datacolors=(red blue) datacontrastcolors=(red blue) ; where Season in(2021, 2022, 2023) and Harvest__month_ in(1,3,4,5,6,8) and strip(Variety) in ('GEM','Hass'); panelby Season/rows=1 onepanel ; vline Harvest__month_/response=DTR group=Variety grouporder=ascending markers markerattrs=(symbol=circlefilled) TRANSPARENCY=0 stat=mean ; rowaxis values=(0 to 12 by 2) ; run;   ... View more

Thank you so much Tom for your great advice! You are right, I specified the length to $9 and it worked! That was all the issue ... View more

Thank you so much Ksharp, thats exactly what I wanted! ... View more

Great, thanks! ... View more

Thanks StatDave ... View more

Ok, thank you Season! ... View more

Oh ok, but at least it helped me a lot with proc glm and when having a binary outcome in proc logistic. Thanks a lot Ksharp! ... View more

Great, thank you so much Ksharp!! ... View more

Thank you so much for your great support on this! That makes a lot of sense so I will remove that outlier and switch to proc glm to treat weeks as a categorical variable.    Thanks a lot! ... View more

Thank yo so much Steve for helping me with this and for the valuable information.  Ok, I used the log statement before the "cards" statement so it worked. If I use the log-transformed Firmness in genmod with a normal dist.  I get AIC: 2.94 and BIC: 17.24. With a normal dist. and no transformation I get 738 and 752.  With a gamma dist. I get 898 and 913.  So I guess I will choose the log-transformed variable with the normal dist. in genmod.    To the second point on making predictions: What do I have to include in my model in order to fit a first degree linear model? What would be the code for that?   Thank you very much ... View more

That makes a lot of sense, I will do so. Thank you so much StatDave!     ... View more

Thanks a lot for your reply! It seems like with proc logistic is getting a bit challenging. I will try Proc FREQ as you suggested. What about using proc genmod, specifiying a multinomial distribution? I dont know much, so Im just asking.    I tried your code for using the unequalslopes test for Harvest and this is what I got: Type 3 Analysis of Effects Effect DF Wald Chi-Square Pr > ChiSq Harvest 0 . . U_Harvest 0 . . Variety 2 14.4212 0.0007   Analysis of Maximum Likelihood Estimates Parameter   PeelColor DF Estimate Standard Error Wald Chi-Square Pr > ChiSq Intercept   2 1 -2.0236 0.5433 13.8708 0.0002 Intercept   3 1 -1.7612 0.5217 11.3981 0.0007 Harvest 1   1 0.6247 0.5672 1.2128 0.2708 Harvest 55   0 0 . . . U_Harvest 1 2 1 0.5878 0.5562 1.1166 0.2906 U_Harvest 1 3 0 0 . . . U_Harvest 55 2 0 0 . . . U_Harvest 55 3 0 0 . . . Variety 465418_9   1 2.4091 0.6456 13.9251 0.0002 Variety BL516   1 1.3044 0.5592 5.4415 0.0197 Variety Hass   0 0 . . .   Which p-value do I need to check?   Thank you so much StatDave!!     ... View more

The syntax looks reasonable to me. Thanks, ... View more

Thank you StatsMan for your reply. I attached the SAS code I am using. I have four replications (blocks) per treatment. I have also attached a file of the map.   Thank you very much for your help! ... View more