かずのたしなみ(統計環境Rによる) @Wiki

ローレンツ順序付け

最終更新:

匿名ユーザー

- view
だれでも歓迎! 編集
##ローレンツ曲線をひく

#まず1985年の5階層別の平均所得を入力する
shotoku1985<-c(0,247897,337088,419016,510747,709482)
#1985年の累積所得を計算する
ruiseki_shotoku1985<-cumsum(shotoku1985)
ruiseki_shotoku1985
[1]       0  247897  584985 1004001 1514748 2224230
#1985年の相対累積所得を計算する
  soutai_ruiseki_shotoku1985<-ruiseki_shotoku1985/ruiseki_shotoku1985[6]
soutai_ruiseki_shotoku1985
[1] 0.0000000 0.1114530 0.2630056 0.4513926 0.6810213 1.0000000

#つぎに1951年の5階層別の平均所得を入力する
shotoku1951<-c(0,5516,11107,14830,19428,31781)
#1985年の累積所得を計算する
ruiseki_shotoku1951<-cumsum(shotoku1951)
ruiseki_shotoku1951
 [1]     0  5516 16623 31453 50881 82662
#1985年の相対累積所得を計算する
soutai_ruiseki_shotoku1951<-ruiseki_shotoku1951/ruiseki_shotoku1951[6]
soutai_ruiseki_shotoku1951
 [1] 0.00000000 0.06672957 0.20109603 0.38050132 0.61553072 1.00000000

#相対累積人数を入力する(これは1951年、1985年に共通である)。
sotai_ruiseki_ninzu<-c(0,0.2,0.4,0.6,0.8,1.0)

##
#1985年のローレンツ曲線を引き、重ねて1951年のローレンツ曲線を引く)
 plot(sotai_ruiseki_ninzu,soutai_ruiseki_shotoku1985,xlim=c(0,1),ylim=c(0,1),type="b")
par(new=T)
plot(sotai_ruiseki_ninzu,soutai_ruiseki_shotoku1951,xlim=c(0,1),ylim=c(0,1),type="b",ann=F,lty=2)
#最後に原点
(0,0)と(1,1)をつなぐ直線(完全平等線)を引く(abline(a,b)は、y=a+bxの直線を引く命令)
abline(0,1,col="red")


タグ:

+ タグ編集
  • タグ:

このサイトはreCAPTCHAによって保護されており、Googleの プライバシーポリシー利用規約 が適用されます。

目安箱バナー