散点图+回归线图

1. data a;
   
2. input y x;
   
3. datalines;
   
4. 0.69 8.48
   
5. 1.10 8.29
   
6. 1.39 8.28
   
7. 1.61 8.27
   
8. 1.79 8.25
   
9. 1.95 8.06
   
10. 2.08 8.06
    
11. 2.20 7.85
    
12. 2.30 7.75
    
13. 2.40 7.73
    
14. ;
    
15. 
    
16. proc reg data=a outest=est noprint;
    
17. model y=x;
    
18. run;
    
19. quit;
    
20. 
    
21. data _null_;
    
22. set est;
    
23. call symputx('a',put(Intercept,7.4));
    
24. call symputx('b',put(x,7.4));
    
25. run;
    
26. 
    
27. 
    
28. proc template;
    
29.   define statgraph reg;
    
30.   mvar a b;
    
31.   begingraph;
    
32.     entrytitle 'Linear Regression';
    
33.     layout overlay;
    
34.         entry 'y=' a b 'x' /valign=top;
    
35.         scatterplot y=y x=x;
    
36.         regressionplot y=y x=x;
    
37.         endlayout;
    
38.   endgraph;
    
39.   end;
    
40. run;
    
41. 
    
42. proc sgrender data=a template=reg;
    
43. run;

复制代码

谢谢，很高兴能帮上忙。

11# guoluo 很感谢，这两个proc都不太熟悉，学习了。

现在又有个新问题。
我在散点图上发现有些点我不想加入回归中。
也就是说，散点图要画所有的点，以防遗漏，但是回归线图，我只用符合特定条件的点。
比如我加入第三列z，回归线图只拟合z>5的x和y。

这样好像proc reg就不管用了，因为他貌似不能在model命令后面加入where语句。
不知道这时该用什么proc才能解决？

data：

y           x          z
0.69     8.48     10
1.10     8.29     9
1.39     8.28     8
1.61     8.27     7
1.79     8.25     6
1.95     8.06     5
2.08     8.06     4
2.20     7.85     3
2.30     7.75     2
2.40     7.73     1

1. data a;
   
2. input y x z;
   
3. if z > 5 then do;
   
4.   y1 = y;
   
5.   x1 = x;
   
6. end;
   
7. datalines;
   
8. 0.69 8.48 10
   
9. 1.10 8.29 9
   
10. 1.39 8.28 8
    
11. 1.61 8.27 7
    
12. 1.79 8.25 6
    
13. 1.95 8.06 5
    
14. 2.08 8.06 4
    
15. 2.20 7.85 3
    
16. 2.30 7.75 2
    
17. 2.40 7.73 1
    
18. ;
    
19. 
    
20. proc reg data=a outest=est noprint;
    
21. model y1=x1;
    
22. run;
    
23. quit;
    
24. 
    
25. data _null_;
    
26. set est;
    
27. call symputx('a',put(Intercept,7.4));
    
28. call symputx('b',put(x1,7.4));
    
29. run;
    
30. 
    
31. 
    
32. proc template;
    
33.   define statgraph reg1;
    
34.   mvar a b;
    
35.   begingraph;
    
36.     entrytitle 'Linear Regression';
    
37.     layout overlay;
    
38.           if(b<0)
    
39.         entry 'y=' a b 'x' /valign=top;
    
40.           else
    
41.             entry 'y=' a '+' b 'x' /valign=top;
    
42.           endif;
    
43.       scatterplot y=y x=x;
    
44.       regressionplot y=y1 x=x1;
    
45.     endlayout;
    
46.   endgraph;
    
47.   end;
    
48. run;
    
49. 
    
50. proc sgrender data=a template=reg1;
    
51. run;

复制代码

受益匪浅~~

学习了。谢谢

thanks for the info~~helps a lot!!!

感谢！学到了！

要將散布圖跟迴歸畫在一起
可以用SGPLOT Procedure
proc sgplot;
scatter x=x y=y;
reg x=x y=y;
run;

如果要選擇z>5的
可以先利用REG Procedure將預測值output
再利用SGPLOT畫散步圖跟線圖
data test;
input y  x z;
k=_n_;
cards;
0.69     8.48 10
1.10     8.29 9
1.39     8.28 8
1.61     8.27 7
1.79     8.25 6
1.95     8.06 5
2.08     8.06 4
2.20     7.85 3
2.30     7.75 2
2.40     7.73 6
;
proc reg;
where z>5;
model y=x;
output out=xx p=pred;
proc print data=xx;
run;
data all;
merge test xx;
by k;
proc  print;
run;
proc sgplot;
scatter x=x y=y;
series x=x y=pred;
run;

thanks for sharing