请教：画一个分位图

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我现在又来源于CHNS的两年数据。
我想比较两年每个分位数上的增长率，比如1989-2009各个分位数上的工资增长率，该如何操作？

我首先自己想到的方法，就是：gen wage_2009= wage if year== 2009
                                                 gen wage_1989= wage if year== 1989
                                                 gen lnwage_2009 = ln(wage_2009)
                                                 gen lnwage_1989=  ln(wage_1989)
                                                 gen lnwage_gap = lnwage_2009 - lnwage 1989
                                           quantile lnwage_gap

画不出想要的图形，分析原因可能是：
1.每年的样本数量都不一样，无法相减
2.产生每年的工资时产生了缺漏值
3.分位数如何能全部提取出来？

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：quantile quant Wage year 如何操作 增长率 如何

需要贴出部分代表性数据。

sungmoo 发表于 2012-3-9 11:31
需要贴出部分代表性数据。

数据大概的样子是这样的，我节选了两年的一部分数据，请版主指点。想计算1989-2009年各个分位上的工资增长率并画图
。工资增长率用对数形式
不知该如何写命令

year        wage        worktype
1989        490        国有单位
1989        644        集体企业
1989        1216.8        国有单位
1989        2821.6        集体企业
1989        1226.4        其他单位
1989        1228.4        国有单位
1989        1190.4        国有单位
1989        100        国有单位
1989        1270.8        国有单位
1989        933.6        国有单位
1989        2059.2        国有单位
1989        3132        集体企业
1989        1096.8        其他单位
1989        584.4        其他单位
1989        730.8        其他单位
1989        12075.2        国有单位
1989        300        其他单位
1989        532.8        其他单位
1989        730.8        集体企业
1989        1502.8        国有单位
1989        913.2        国有单位
1989        501.6        国有单位
1989        1903.6        国有单位
1989        2252.8        国有单位
1989        1392.8        国有单位
1989        1575        其他单位
1989        1550        国有单位
1989        632.8        国有单位
1989        1188.8        国有单位
2009        9600        其他单位
2009        12800        其他单位
2009        3000        其他单位
2009        5400        
2009        14000        国有单位
2009        7000        其他单位
2009        2000        其他单位
2009        100        其他单位
2009        8600        其他单位
2009        3000        其他单位
2009        1300        其他单位
2009        18000        其他单位
2009        8800        
2009        1000        其他单位
2009        9460        其他单位
2009        221000        集体企业
2009        5000        其他单位
2009        9000        
2009        5600        国有单位
2009        86400        
2009        12000        国有单位
2009        8900        其他单位
2009        14000        国有单位
2009        10800        其他单位
2009        4200        
2009        7200        其他单位
2009        4800        其他单位
2009        22700        其他单位
2009        6000        其他单位
2009        200        其他单位
2009        5400        其他单位
2009        14400        其他单位
2009        12500        国有单位
2009        12800        国有单位
2009        4800        国有单位
2009        14400        
2009        18000        国有单位
2009        9600        其他单位
2009        9600        其他单位
2009        12000        其他单位
2009        12000        其他单位
2009        24800        国有单位
2009        36000        国有单位
2009        16600        其他单位
2009        20400        国有单位
2009        4000        
2009        24000        其他单位
2009        6000        其他单位
2009        6200        其他单位
2009        24000        其他单位
2009        28000        其他单位
2009        3000        国有单位
2009        12000        其他单位
2009        14400        其他单位
2009        18000        国有单位
2009        10100        其他单位
2009        4400        其他单位
2009        9960        其他单位
2009        16400        其他单位
2009        10800        国有单位
2009        450        
2009        4760        
2009        366796        其他单位
2009        18000        其他单位
2009        12000        其他单位
2009        20800        国有单位
2009        30000        国有单位
2009        10800        其他单位
2009        12000        其他单位
2009        4800        其他单位
2009        12100        其他单位
2009        4800        
2009        19100        国有单位
2009        13600        其他单位
2009        6000        其他单位
2009        15200        国有单位
2009        9600        国有单位
2009        29200        国有单位
2009        15600        其他单位
2009        12150        集体企业
2009        24200        国有单位
2009        12200        其他单位
2009        23900        国有单位
2009        7200        集体企业
2009        300        其他单位
2009        12860        其他单位
2009        14400        其他单位
2009        20802.829        国有单位
2009        9800        其他单位
2009        6000        其他单位
2009        9600        
2009        12000        其他单位
2009        4800        其他单位
2009        11700        其他单位
2009        9600        其他单位
2009        800        
2009        18000        国有单位
2009        7450        其他单位
2009        15300        国有单位
2009        12000        其他单位
2009        8600        其他单位
2009        9600        其他单位
2009        18000        其他单位
2009        6000        其他单位
2009        25580        国有单位
2009        32400        国有单位
2009        500        
2009        6000        国有单位
2009        28000        国有单位
2009        10200        其他单位
2009        8400        国有单位
2009        11250        其他单位
2009        12200        其他单位
2009        1000        
2009        7200        其他单位
2009        24000        其他单位
2009        19000        国有单位
2009        9600        
2009        12300        其他单位
2009        8000        其他单位
2009        14400        
2009        10800        
2009        4200        其他单位
2009        500        其他单位
2009        6400        其他单位
2009        5000        其他单位
2009        1560        其他单位
2009        15000        其他单位
2009        9000        其他单位
2009        6000        其他单位
2009        600        其他单位
2009        144000        其他单位
2009        14500        国有单位
2009        5800        其他单位
2009        9600        其他单位
2009        1400        其他单位
2009        10600        其他单位
2009        6000        其他单位
2009        80000        其他单位
2009        550        其他单位
2009        8500        其他单位
2009        300        其他单位
2009        288000        
2009        2800        其他单位
2009        12000        其他单位
2009        21000        国有单位
2009        36000        其他单位
2009        6000        其他单位
2009        4400        其他单位
2009        1300        其他单位
2009        9600        其他单位
2009        3000        其他单位
2009        9600        
2009        9600        其他单位
2009        15600        其他单位
2009        10600        国有单位
2009        1100        其他单位
2009        3400        其他单位
2009        3600        其他单位
2009        5000        其他单位
2009        10800        其他单位
2009        12400        国有单位
2009        5000        其他单位
2009        12000        集体企业
2009        36000        其他单位
2009        32400        其他单位
2009        7200        其他单位
2009        12000        其他单位
2009        2000        其他单位
2009        100        其他单位
2009        6600        其他单位
2009        12000        其他单位
2009        7200        其他单位
2009        1500        其他单位
2009        2200        其他单位
2009        6400        其他单位
2009        1000        其他单位
2009        30000        其他单位
2009        2550        其他单位
2009        1000        其他单位
2009        8000        其他单位
2009        29000        其他单位
2009        7400        其他单位
2009        2200        其他单位
2009        13900        其他单位
2009        8000        其他单位
2009        5400        其他单位
2009        17400        其他单位
2009        1800        其他单位
2009        7300        其他单位
2009        6000        其他单位
2009        19200        其他单位
2009        21000        其他单位
2009        12000        其他单位
2009        12100        其他单位
2009        18000        其他单位

*试一下（以下以百分位点为例）：

bys year (wage): cumul wage, g(p) eq
replace p=int(100*p)/100
collapse (min)wage,by(year p)
bys p (year): g d=wage/wage[_n-1]-1
line d p

*缺失值可以自动忽略

我要所有分位上的增长率，我想画一条线，横轴是各个分位点，纵轴是增长率，类似于分位图那样的。

alabala 发表于 2012-3-9 13:26
我要所有分位上的增长率，我想画一条线，横轴是各个分位点，纵轴是增长率，类似于分位图那样的。

如果不拟合出一条连续的曲线，只根据样本，未必各分位点上都有对应的样本观测值。

我刚才用一个笨方法，您给看一下。
use 2010_total_2.dta, clear
keep if year== 2009 | year== 1989      // 只保留需要用的两年的样本

forvalues i= 1/99 {
by year, sort: egen w`i'= pctile(wage), p(`i')
gen lnw`i'= ln(w`i')

}
forvalues i= 1/99 {
duplicates drop year w`i', force

}

forvalues i= 1/99 {
drop w`i'

}

xpose, clear
gen ratio = v2 - v1
drop if v1 == 1989
gen n= _n

twoway line ratio n, sort lw(*1.1) msize(large)  ///
       xlabel(0(10)100)  ylabel(-1(1)3)  ///
           xtitle("分位")   ytitle(" ") lw(*2)

sungmoo 发表于 2012-3-9 13:39
如果不拟合出一条连续的曲线，只根据样本，未必各分位点上都有对应的样本观测值。

版主，关于这个分位图，看了您的讲解，我想提两个具体的问题：
1.可否画出各个位点的增长率，范围是1 （5）96  ，即范围为1-96的等差（每次递增5）数列。图形方式希望是connected的这种。如何画
2.我用的方法太笨了，不知道可否有简洁的命令。

可否画出各个位点的增长率，范围是1 （5）96  ，即范围为1-96的等差（每次递增5）数列。图形方式希望是connected的这种。如何画

前面已经给出了百分位点的增长率。

这里先有个问题：如果对于某个概率值，两组样本并不同时具有相应的分位值，怎么办？

sungmoo 发表于 2012-3-10 08:22
前面已经给出了百分位点的增长率。

这里先有个问题：如果对于某个概率值，两组样本并不同时具有相应 ...

版主，我觉得我的数据样本足够多，应该不存在“不同时具有相应分位值的情况”。
因为分位数的定义是：
第 p_th 百分位数值，记为 x_|p|，设 p = N_p/100，
      例如，若求取第 25 百分位的数值，则 N_p=25, p=0.25
      每个观察值的权重为 w(i)=1/N (N为样本数)，
     前 i 个观察值的权重之和为 W(i)=SUM_w(i) = i/N
      则第 p_th 百分位的数值定义为，第一个满足 W(i)>p 的观察值，即
     {x[i-1]+x}/2    if W(i-1)=p
      x              otherwise
不知道我理解的对吗