In: Statistics and Probability
data one;
input NEA Fat;
cards;
-94 4.2
-57 3.0
-29 3.7
135 2.7
143 3.2
151 3.6
245 2.4
355 1.3
392 3.8
473 1.7
486 1.6
535 2.2
571 1.0
580 0.4
620 2.3
690 1.1
;
title "The Reg Procedure";
proc sgplot data=one;
scatter x=NEA Y=Fat;
run;
proc sgplot data=one;
reg x=NEA y=Fat;
run;
This is my data for a scatterplot. I need to know how to plot the regression line by the NEA data.
Sol:
With proc reg we build linear regression of y on x
store the parameter estimates in esti
store in a macro variables estimates and call them in proc sgplot with inset
execute below code:
SAS Code:
ods graphics off;
proc reg data=one;
model Fat = NEA;
ods output ParameterEstimates=esti;
run;
data _null_;
set esti;
if _n_ = 1 then call symput('Intercpt', put(estimate,
BEST6.));
else call symput('Slope', put(estimate, BEST6.));
run;
proc sgplot data=one noautolegend;
reg x=NEA y=Fat;
inset "Intercept = &Intercpt" "Slope = &Slope" /
border title="Parameter Estimates" position=topleft;
run;
Output:
The REG Procedure
Model: MODEL1
Dependent Variable: Fat
Number of Observations Read | 16 |
---|---|
Number of Observations Used | 16 |
Analysis of Variance | |||||
---|---|---|---|---|---|
Source | DF | Sum of Squares |
Mean Square |
F Value | Pr > F |
Model | 1 | 11.79415 | 11.79415 | 21.55 | 0.0004 |
Error | 14 | 7.66335 | 0.54738 | ||
Corrected Total | 15 | 19.45750 |
Root MSE | 0.73985 | R-Square | 0.6061 |
---|---|---|---|
Dependent Mean | 2.38750 | Adj R-Sq | 0.5780 |
Coeff Var | 30.98860 |
Parameter Estimates | |||||
---|---|---|---|---|---|
Variable | DF | Parameter Estimate |
Standard Error |
t Value | Pr > |t| |
Intercept | 1 | 3.50512 | 0.30362 | 11.54 | <.0001 |
NEA | 1 | -0.00344 | 0.00074141 | -4.64 | 0.0004 |