options ls=78; title "Profile Plot - Women's Nutrition Data"; %let p=5; data nutrient; infile "D:\Statistics\STAT 505\data\nutrient.txt"; input id calcium iron protein a c; variable="calcium"; ratio=calcium/1000; output; variable="iron"; ratio=iron/15; output; variable="protein"; ratio=protein/60; output; variable="vit a"; ratio=a/800; output; variable="vit c"; ratio=c/75; output; keep variable ratio; run; proc sort; by variable; run; proc means; by variable; var ratio; output out=a n=n mean=xbar var=s2; run; data b; set a; f=finv(0.95,&p,n-&p); ratio=xbar; output; ratio=xbar-sqrt(&p*(n-1)*f*s2/(n-&p)/n); output; ratio=xbar+sqrt(&p*(n-1)*f*s2/(n-&p)/n); output; run; proc gplot; axis1 length=4 in; axis2 length=6 in; plot ratio*variable / vaxis=axis1 haxis=axis2 vref=1 lvref=21; symbol v=none i=hilot color=black; run;