Stata绘图：图形中加入公式

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

偶然看到一个2018年的帖子，求助Stata绘图中加入公式，有坛友认为Stata画不了下面的图，肯定是Excel画出来的，大家看看下面的图能用Stata画出来吗？这是https://bbs.pinggu.org/forum.php ... &page=1#pid74041454中要求画的图。

由于坛友没有给出数据，我模拟了一份数据来作图，最终效果如下：

Stata代码如下：

1.    **************************生成模拟数据******************************
   
2.     clear
   
3.     set obs 20
   
4.     set seed 123456
   
5.     gen e = rnormal()     // e~N(0,1)
   
6.     gen x = rnormal(80,30)  // x~N(0,3^2)
   
7.     gen y1 = .
   
8.     replace y1 = 1.1091*x -10.2+ e
   
9.     gen y2=.
   
10.     replace y2 = 0.7486*x +25.069+ e
    
11.         
    
12.         *****************************作图********************************
    
13.         #delimit ;
    
14.      twoway (scatter y1 x,msymbol(s)) (lfit y1 x)
    
15.                  (scatter y2 x,msymbol(d)) (lfit y2 x)
    
16.                    ,
    
17.                    xlabel(0(20)140) ylabel(0(20)140,grid angle(0))
    
18.                    legend(order(3 "PL_gfcf" 1 "PL_cons") col(1) ring(0)
    
19.                    pos(11) region(fcolor(none) lpattern(blank)))
    
20.                    text(74 98 "y = 1.1091x -10.2" "R{sup:2}=0.992" 90 65
    
21.                    "y = 0.7486x +25.069" "R{sup:2}=0.844")
    
22.                    scheme(s1mono);
    
23.         #delimit cr

复制代码

再提供一个与原图配色方案基本一致的版本：
代码如下：

1.    **************************生成模拟数据******************************
   
2.    clear
   
3.     set obs 20
   
4.     set seed 123456
   
5.     gen e = rnormal(1,5)     // e~N(0,1)
   
6.     gen x = rnormal(80,30)  // x~N(80,30^2)
   
7.     gen y1 = .
   
8.     replace y1 = 1.1091*x -10.2+ e
   
9.     gen y2=.
   
10.     replace y2 = 0.7486*x +25.069+ e
    
11. 
    
12. *****************************作图********************************
    
13.         #delimit ;
    
14.     twoway (scatter y1 x,msymbol(s) mcolor(red)) (lfit y1 x,lcolor(red) lwid(*1.8))
    
15.                 (scatter y2 x,msymbol(d) mcolor(ebblue)) (lfit y2 x,lcolor(ebblue) lwid(*1.8))
    
16.                    ,
    
17.                    ylabel(0(20)140,grid angle(0)) xlabel(0(20)140) xtitle("")
    
18.                    legend(order(3 "PL_gfcf" 1 "PL_cons") col(1) ring(0)
    
19.                    pos(11) region(fcolor(none) lpattern(blank)))
    
20.                    text(74 98 "y = 1.1091x -10.2" "R{sup:2}=0.992",color(red))
    
21.                    text(90 65 "y = 0.7486x +25.069" "R{sup:2}=0.844" ,color(ebblue));
    
22.         #delimit cr

复制代码

大家觉得这三幅图一样吗？Stata绘图还是很强大的哦！觉得对你有帮助的话，点个赞，评个分，支持一下，谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏3 回帖

关键词：stata绘图 Stata tata Scatter delimit

真心不错，支持！

调色的话，可以help palette_all,会给出多种颜色让你选择。

可以在text里加选项linegap(1)调整文字间距，可以对红色red改为red*1.2,跟原图就更像了。