Dear Rick, Thank you for your help. I was finally able to replicate your code with my data and create the biplot. Initially, I was trying to extract the scores using PROC PRINCOMP, but I couldn't manage to do it. However, I was able to get it working using PROC PRINQUAL. By any chance, do you know if it’s possible to extract the scores using PROC PRINCOMP? Thanks again!     ... View more

Thank you for your response! I tried to replicate your code with my data, but I couldn’t get it to work. So I went back to trying the approach I was using with PROC PRINCOMP. However, I’m having trouble extracting the loadings (OUTSTAT=OutStat) and combining them with the scores (OUT=PCA). Below is the code I’m trying to use. If I remove the creation of the vectors dataset and also the VECTOR statement from PROC SGPLOT, the plot is generated normally, but it doesn’t include the variable vectors (loadings). data iris;    set Sashelp.iris; run;   proc princomp data=iris plots=pattern(vector) plots=(score(ellipse)) n=2 standard  out=PCA outstat=OutStat noprint; var SepalLength                SepalWidth        PetalLength        PetalWidth; id Species; run;   data vectors;    set OutStat;    where _TYPE_='SCORE';    length Variable $12.;    Variable = _NAME_;    x = Prin1 * 2;    y = Prin2 * 2;    keep Variable x y; run;   ods graphics / reset width=800px height=600px imagename="PCA_Graph" imagefmt=png antialias=on; proc sgplot data=PCA;                 scatter x=Prin1 y=Prin2 / group=species                                                                                                                 markerattrs=(size=12 symbol=circlefilled)                                                                                                                 transparency=0.2;                                                                                                                    vector x=x y=y / datalabel=Variable                                                                               arrowattrs=(thickness=2 color=black)                     lineattrs=(pattern=solid);                                    xaxis label="CP1 (72.77%)"                                  labelattrs=(Family="Times New Roman" Size=11 Weight=bold)           valueattrs=(Family="Times New Roman" Size=10 Weight=bold);                                     yaxis label="CP2 (23.03%)"                                  labelattrs=(Family="Times New Roman" Size=11 Weight=bold)           valueattrs=(Family="Times New Roman" Size=10 Weight=bold);                  keylegend / location=outside                                                                position=right                                                                title="Species"                                                                titleattrs=(color=black family="Times New Roman" size=11)                                                                valueattrs=(color=black family="Times New Roman" size=10);                 ellipse x=Prin1 y=Prin2 / group=species                                                                                                                  lineattrs=(pattern=solid thickness=2)                                                                                                                 transparency=0.4;                   refline 0 / axis=x lineattrs=(color=gray pattern=solid thickness=1);                 refline 0 / axis=y lineattrs=(color=gray pattern=solid thickness=1);                 run;                 ods graphics / reset;     ... View more

I've already answered your question. The letter "a" means the group 300 is better the group 0 (letter c). Do you know some way to put some letters in my vbar? If you don't know, I will give up. ... View more

Ok, I ran another code with Tukey test. You can see it below. According to my statistics, I want to put the letter 'a in group 300', 'ab in groups 100, 200 and control' and put the letter 'c in group 0'. I want to do a graph like that (second graph).   data teste; input source$ rep Carbon; cards; Controle 1 10.1 Controle 2 10.63 Controle 3 10.34 Controle 4 11.51 Controle 5 12.4 0 1 9.63 0 2 12.63 0 3 10.28 0 4 10.93 0 5 10.63 100 1 10.46 100 2 11.87 100 3 9.28 100 4 14.8 100 5 13.04 200 1 13.34 200 2 12.04 200 3 9.81 200 4 9.99 200 5 13.86 300 1 11.22 300 2 12.63 300 3 14.86 300 4 17.92 300 5 15.69 ; proc glm data=teste plots=diagnostics; class source; model carbon = source; means source / hovtest=levene(type=abs) tukey; run; proc sgplot data=teste noautolegend noborder; vbar source / response=carbon stat=mean limitstat=stderr limitattrs=(color=black) fillattrs=(color=VIGB) outlineattrs=(color="black") datalabel datalabelattrs=(color=black size=10 family="Times New Roman") fillattrs=(color=VIGB) outlineattrs=(color="black") ; yaxis label="Soil carbon" LABELATTRS=(Color=black Family="Times New Roman" Size=12) valueattrs=(Color=black Family="Times New Roman" Size=10); xaxis label="Source of nutrient" LABELATTRS=(Color=black Family="Times New Roman" Size=12) valueattrs=(Color=black Family="Times New Roman" Size=10); run;     ... View more

The letters came from the Tukey test, meaning that group 0 is different from group 100, etc., so the letters help us understand which group is greater. For example, I want to put the "letter a" above stderr bar of group 300. How can I do that? If I have to modify my dataset, how can I do that? Do you understand my doubt?  ... View more

if you run my code, you will see one graph with five bars. So I want to know how I can put some letters above them. For example, "a, b, c, d and e" above each stderr bar. Do you understand?    ... View more

I change my PROC SGPLOT code with:title "pH do solo em função de fontes de adubos"; proc sgplot data=composto2 noautolegend noborder; vbar trat / response=Nsolo stat=mean limitstat=stderr limitattrs=(color=black) fillattrs=(color=VIGB) outlineattrs=(color="black") datalabel datalabelattrs=(color=black size=10 family="Times New Roman") fillattrs=(color=VIGB) outlineattrs=(color="black") ; yaxis label="pH do solo" LABELATTRS=(Color=black Family="Times New Roman" Size=12) valueattrs=(Color=black Family="Times New Roman" Size=10); xaxis label="Fonte de adubo" LABELATTRS=(Color=black Family="Times New Roman" Size=12) valueattrs=(Color=black Family="Times New Roman" Size=10); run; Now I can put stderr bar, but I couldn't put the letters to mark significant differences. Please help me. ... View more

Thanks for answer me. For "insert standard error" I wanna mean an error bar. I want also to add letters above each bar, such as "a, b, c, d and e". My current code is below:   data teste; input trat$ rep pHsolo; cards; Controle 1 7.51 Controle 2 7.1 Controle 3 7.71 Controle 4 7.6 Controle 5 7.44 0 1 7.69 0 2 6.76 0 3 6.54 0 4 7.3 0 5 7.64 100 1 7.57 100 2 7.13 100 3 7.6 100 4 7.56 100 5 7.87 200 1 7.56 200 2 7.67 200 3 7.73 200 4 7.58 200 5 7.55 300 1 7.67 300 2 7.41 300 3 7.77 300 4 7.58 300 5 7.81 ; title "pH do solo em função de fontes de adubos"; proc sgplot data=composto2 noautolegend noborder; vbar trat / response=Nsolo datalabel datalabelattrs=(color=black size=10 family=arial) fillattrs=(color=VIGB) outlineattrs=(color="black") ; yaxis label="pH do solo" LABELATTRS=(Color=black Family=Arial Size=12) valueattrs=(Color=black Family=Arial Size=10); xaxis label="Fonte de adubo" LABELATTRS=(Color=black Family=Arial Size=12) valueattrs=(Color=black Family=Arial Size=10); run; ... View more