如何在折线图上添加概率密度函数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

在把收入对数的均值关于受教育年限做折线图之后，如何在上面添加如图所示的概率密度函数




数据格式为
受教育年限   收入对数0        1
0        2
1        3
1        2
1        1
2        3
2        4
2        5
3      
……………………

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

比较难画
你可以学习这里的例子
https://www.survey-design.com.au/tipsgraphs.html

1. // distribution plot
   
2. 
   
3. clear
   
4. sysuse auto
   
5. version 9.2
   
6. keep if foreign
   
7. sort weight
   
8. 
   
9. gen weight2 = weight^2
   
10. regress mpg weight weight2
    
11. predict fit
    
12. predict se , stdp
    
13. 
    
14. #delimit ;
    
15. twoway sc mpg weight , pstyle(p3) ms(o) ||
    
16. fn weight[3]  - 1000 * normden(x, `=fit[3]' , `=se[3]') ,
    
17. range(`=fit[3] -5' `=fit[3] +5') horiz pstyle(p1) ||
    
18. fn `=fit[3]' , range(`=weight[3]' `=weight[3]-1000*normden(0, se[3])')
    
19. pstyle(p1) ||
    
20. fn weight[17] - 1000 * normden(x, `=fit[17]', `=se[17]') ,
    
21. range(`=fit[17]-5' `=fit[17]+5') horiz pstyle(p1) ||
    
22. fn `=fit[17]', range(`=weight[17]' `=weight[17]-1000*normden(0, se[17])')
    
23. pstyle(p1) ||
    
24. fn weight[21] - 1000 * normden(x, `=fit[21]' , `=se[21]') ,
    
25. range(`=fit[21] -7' `=fit[21] +7') horiz pstyle(p1) ||
    
26. fn `=fit[21]', range(`=weight[21]' `=weight[21]-1000*normden(0, se[21])')
    
27. pstyle(p1) ||
    
28. line fit weight
    
29. ,clwidth(*2) legend(off) ytitle(Miles per gallon) xtitle(Weight)
    
30. title("Scatter with Regression Line and Confidence Interval Densities"
    
31. ,size(4.8) margin(t=0 b=1.5) span)
    
32. caption("graph46", box size(vsmall)position(5) ring(10))
    
33. ;
    
34. #delimit cr

复制代码