The distribution of Power Law is seen though it bought for 2008/10 by delivering to home food.
"Power Law:(20)Jobless celebration?" と同様に、品物の金額に関して傾向を見る。
The tendency is seen as well as "Power Law:(20)Jobless celebration?" for the amount of money of goods.
http://humanbeing-etcman.blogspot.com/2008/11/power-law20jobless-celebration.html
~~~
(１週間目,One weeks from the beginning)
近似曲線,Approximation curve)
LN(y) = -1.0666*LN(x) + 6.9917
LN
(x, y)=(1, exp(6.9917)),(6, exp(-1.0666*LN(6) + 6.9917))
=(1, 1087.57),(6, 160.87)
Relative value)
=(1, 1087.57/160.87),(6, 1)
=(1, 6.761),(6, 1)
A=-1.0666
B=6.761
~~~
(２週間目,Two weeks from the beginnings)
Approximation curve)
LN(y) = -0.8878*LN(x) + 7.2903
LN
(x, y)=(1, exp(7.2903)),(11, exp(-0.8878*LN(11) + 7.2903))
=(1, 1466),(11, 174.42)
Relative value)
=(1, 1466/174.42),(11, 1)
=(1, 8.405),(11, 1)
A=-0.8878
B=8.405
~~~
(３週間目,Three weeks from the beginnings)
Approximation curve)
LN(y) = -1.1581*LN(x) + 8.1678
LN
(x, y)=(1, exp(8.1678)),(15, exp(-1.1581*LN(15) + 8.1678)))
=(1, 3525.58),(15, 153.18)
Relative value)
=(1, 3525.58/153.18),(15, 1)
=(1, 23.02),(15, 1)
A=-1.1582
B=23.02
~~~
(４週間目,Four weeks from the beginnings)
Approximation curve)
LN(y) = -1.1581*LN(x) + 8.1678
３週目に同じ。
It is the same as the third week.
~~~
月初めから各週までのA-B変化は以下のようになった。
The A-B change from beginning of the month to each week was as follows.
1week)
A=-1.0666
B=6.761
2week)
A=-0.8878
B=8.405
3week)
A=-1.1582
B=23.02
~~~
Rで、A-Bをプロットする。
A-B is plotted by R.
powerlaw-24-a-b.txt
===
A,B,week
-1.0666,6.761,1week
-0.8878,8.405,2week
-1.1582,23.02,"3,4week"
===
data_ab = read.csv("powerlaw-24-a-b.txt");
plot(log(abs(data_ab$A[0])), log(data_ab$B[0]), xlab="log(abs(a))", ylab="log(b)", xlim=c(-0.2, 0.2), ylim=c(1, 4),
col="orange", pch=20, main="Scatter chart:a-b, Power Law(24)")
par(new=T)
#
for (i in 1:length(data_ab$A)){
plot(log(abs(data_ab$A[i])), log(data_ab$B[i]), xlim=c(-0.2, 0.2), ylim=c(1, 4), ann=F)
text(log(abs(data_ab$A[i])), log(data_ab$B[i]), data_ab$week[i], pos=4, offset=0.5)
par(new=T)
}
abline(v = log(1), col="red")
#
abline(log(4), 5, col="gray", lty=2)
text(log(1), log(4), "L(3)=4", pos=4, offset=0.5)
par(new=T)
plot(log(1), log(4), xlim=c(-0.2, 0.2), ylim=c(1, 4), pch=20, ann=F)
#
abline(log(5), 5, col="gray", lty="dotted")
text(log(1), log(5), "5?", pos=4, offset=0.5)
#
abline(log(7), 5, col="gray", lty=2)
text(log(1), log(7), "L(4)=7", pos=4, offset=0.5)
par(new=T)
plot(log(1), log(7), xlim=c(-0.2, 0.2), ylim=c(1, 4), pch=20, ann=F)
#
abline(log(11), 5, col="gray", lty=2)
text(log(1), log(11), "L(5)=11", pos=4, offset=0.5)
par(new=T)
plot(log(1), log(11), xlim=c(-0.2, 0.2), ylim=c(1, 4), pch=20, ann=F)
#
abline(log(15), 5, col="gray", lty="dotted")
text(log(1), log(15), "15?", pos=4, offset=0.5)
#
abline(log(18), 5, col="gray", lty=2)
text(log(1), log(18), "L(6)=18", pos=4, offset=0.5)
par(new=T)
plot(log(1), log(18), xlim=c(-0.2, 0.2), ylim=c(1, 4), pch=20, ann=F)
[図に関するコメント,Comment concerning figure]
[1]
１週目と２週目は、Lucasの線上にない。３週目と４週目はLucasの線上にある。
The first week and the second week, it is not on the line of Lucas.
The third week and the fourth week are on the line of Lucas.
[2]
３週目は、１週目と２週目の中間にある？
：１週目と２週目は変動して、３週目で収束した？
The third week is in the middle of the first week and the second week?
:The first week and the second week changed, and were settled on the week third?
~~~
[Tips for R]
[1]データのラベルを出したい。
I want to put out the label of data.
例)
text(x, y, label, pos=4, offset=0.5)
http://www.is.titech.ac.jp/~mase/mase/html.jp/temp/text.jp.html
pos テキストの位置指定子。
値 1, 2, 3 そして 4 はそれぞれ指定された座標の、下側、左側、上側、右側を指示する。
~~
[2]ablineのltyの使い方は、以下を参照。
Refer to the usage of lty of abline as follows.
http://cse.naro.affrc.go.jp/takezawa/r-tips/r/53.html
抜粋,Excerpt)
lty=2 ，lty="dashed"
：線分の形式 (line type) をダッシュにする．
lty=3 ，lty="dotted"
：線分の形式 (line type) をドットにする．
~~~~
番外）
昨日、人材銀行に登録した。
職種として、SE現役か、営業SEのいずれかを選択するか決めかねていたが、
スキルに関する自己紹介の文面を直した方がいいと言われた。確かに！
とりあえず、「大学時代、、、」
~~~~
end
0 件のコメント:
コメントを投稿