sepscatter 之应用！

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我昨天询问 Nick Cox (Editor of Stata Journal) 有关 -sepscatter- 指令可否加入配适线，感谢他的回复，已经增加此功能。

1. // https://www.statalist.org/forums/forum/general-stata-discussion/general/3803-sepscatter-available-from-ssc
   
2. *! 1.1.0 NJC 5 December 2018
   
3. *! 1.0.2 NJC 9 May 2014
   
4. *! 1.0.1 NJC 8 May 2014
   
5. *! 1.0.0 NJC 29 April 2014
   
6. program sepscatter
   
7.     version 9
   
8. 
   
9.     capture syntax anything [if] [in] [aweight fweight pweight] ///
   
10.     , seperate(varname) [ * ]
    
11. 
    
12.     if _rc == 0 {
    
13.         noisily di _n "note: sep" as err "a" as  txt "rate() is correct spelling"
    
14. 
    
15.         local 0 `anything' `if' `in' [`weight' `exp'] ///
    
16.                 , separate(`seperate') `options'
    
17.     }
    
18. 
    
19.     syntax varlist(numeric min=2 max=2) [if] [in] ///
    
20.     [aweight fweight pweight] , SEParate(varname) ///
    
21.     [MYLAbel(varname) MYNUmeric(varname) MISSing addplot(str asis) *]
    
22. 
    
23.     capture noisily {  
    
24.    
    
25.     quietly {
    
26.         if "`mylabel'" != "" & "`mynumeric'" != "" {
    
27.             di as err "choose mylabel() or mynumeric()"
    
28.             exit 198
    
29.         }
    
30. 
    
31.         tokenize `varlist'
    
32.         args y x
    
33. 
    
34.         marksample touse
    
35.         if "`missing'" == "" markout `touse' `separate', strok  
    
36.         count if `touse'
    
37.         if r(N) == 0 exit 2000
    
38. 
    
39.         tempname stub
    
40.         separate `y' if `touse', `missing' by(`separate') ///
    
41.         gen(`stub') veryshortlabel
    
42.         local Y `r(varlist)'
    
43.         local nY : word count `Y'
    
44.     }
    
45. 
    
46.     local ytitle : var label `y'
    
47.     if `"`ytitle'"' == "" local ytitle "`y'"
    
48. 
    
49.     if "`mylabel'`mynumeric'" != "" {
    
50. 
    
51.         if "`mynumeric'" != "" {
    
52.             if "`: value label `mynumeric''" != "" {
    
53.                 tempvar mylabel
    
54.                 gen `mylabel' = `mynumeric'
    
55.             }
    
56.             else local mylabel `mynumeric'
    
57.         }
    
58.                
    
59.         local mylabel : di _dup(`nY') "`mylabel' "
    
60.         local mypos : di _dup(`nY') "0 "
    
61.         local mynone : di _dup(`nY') "none "
    
62.         local mylabel ///
    
63.         ms(`mynone') mla(`mylabel') mlabpos(`mypos') legend(off)   
    
64.     }
    
65. 
    
66.     scatter `Y' `x' if `touse' [`weight' `exp'], ///
    
67.         ytitle(`"`ytitle'"') ms(Oh plus X Th Sh Dh) `mylabel' ///
    
68.         `options' || `addplot'
    
69.     }
    
70. 
    
71.     drop `Y'
    
72. end
    
73. 
    
74. sysuse auto, clear
    
75. sepscatter mpg weight , sep(foreign) mc(red blue) addplot(qfit mpg weight if !foreign, lc(red) || qfit mpg weight if foreign, lc(blue)) legend(order(1 2))

复制代码

其图形为

扫码加我 拉你入群

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

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

分享0 收藏33 回帖

关键词：word

不错，谢谢分享

命令运行到最后没法运行，黄老师可否把ado分享下

yuqr1986 发表于 2018-12-6 18:36
命令运行到最后没法运行，黄老师可否把ado分享下

上面就是 ado 档。

不错不错，谢谢分享

黃河泉 发表于 2018-12-7 07:36
上面就是 ado 档。

黄老师，在运行sepscatter mpg weight , sep(foreign) mc(red blue) addplot(qfit mpg weight if !foreign, lc(red) || qfit mpg weight if foreign, lc(blue)) legend(order(1 2))这一行命令时，出现了option addplot() not allowed这一错误？

hudabeicun 发表于 2018-12-7 10:22
黄老师，在运行sepscatter mpg weight , sep(foreign) mc(red blue) addplot(qfit mpg weight if !forei ...

因为 sepscatter 尚未更新 (不知什么时候)，请利用上面整个档来画图。

黄老师，请教您个问题。我现在做的是类似于 ：“30个国家  每个国家10个省份  GDP与 人口的散点图”，
我的命令是 ：scatter GDP 人口,by(省份)
上面这么做，没有把国家分开。我想在by后面加 省份、国家两个控制变量。
scatter GDP 人口,by(国家,省份)   这样的命令跑不出来。
希望黄老师看到之后给予指导。谢谢

黄老师，我在线等您的回复。

2865805767 发表于 2019-2-15 11:11
黄老师，请教您个问题。我现在做的是类似于 ：“30个国家  每个国家10个省份  GDP与 人口的散点图”，
我的 ...

你若要问程序，永远附上相关资料；若附上资料，永远用 dataex 印出资料。
1.        先 ssc install dataex (并见说明)，将原始 Stata 资料中具有”代表性”的一部分资料列出，以供有意回答者实验之用，并能提供具体操作指令。
2.        并请参考 http://www.jianshu.com/p/9870080fe769,  https://bbs.pinggu.org/thread-5048204-1-1.html, 与 https://bbs.pinggu.org/thread-5917273-1-1.html。