2008年12月14日日曜日

Power Law:(31)STOCK,B-5A-Fn,(1)

(2008/12/15)
終値に誤り発見!Power Law:(33)を参照。再トライ,TODO)
Mistakes to the closing share price!Refer to Power Law:(33). Re-trial and TODO)

http://humanbeing-etcman.blogspot.com/2008/12/power-law33stocklast-price-incorrect.html

(2008/12/03-2008/12/06)

東証の食品、電機の株価の日別べき乗傾向から、A-Bを収集している。
A-B is collected from the Power Law tendency to stock prices of the Food and Electric Appliances of Tokyo Stock Exchange according to the day.
~~~
[1]
食品は、以下(Power Law:(29))からデータが増加した。
As for Food, data has increased as follows at Power Law (29).
http://humanbeing-etcman.blogspot.com/2008/12/power-law29stocktsefoods200811xx2.html



:図から何が見えるか?
:What do you see from figure?

~~~
[2]
電機も、以下(Power Law:(28))からデータが増加した。
As for Electric Appliances, data has increased as follows at Power Law (28).
http://humanbeing-etcman.blogspot.com/2008/12/power-law28stocktseelectric.html



~~~
[3]
A-Bの散布図では日別遷移がわかりづらいので、A、Bを別々に着目してみる。
It pays attention separately A and B because the transition according to the day is not understood easily in the scatter chart of A-B.

[3-0]
案)Aは、リズムしているか?
The idea) Does A do the rhythm?


---
[3-1]
食品のAの日別遷移, Transition of A of Food according to day
:時間軸は下に経過する, The time axis passes below.


フィボナッチ数:F11=89, F12=144, F13=233
④:AVERAGE(AVERAGE(F11, F13), F12)で振動している?

Number of Fibonacci: F11=89, F12=144, and F13=233
Does it vibrate by ④:AVERAGE(AVERAGE(F11, F13), F12)?

---
[3-2]
電機のAの日別遷移, Transition of A of Electric Appliances according to day
:時間軸は下に経過する, The time axis passes below.

フィボナッチ数:F10=55, F11=89, F12=144
⑤:AVERAGE(AVERAGE(F10, F12), F11)と③:AVERAGE(F10, F12)の間で振動している?

Number of Fibonacci: F10=55, F11=89, and F12=144
Do it vibrate between ⑤:AVERAGE(AVERAGE(F10, F12), F11) and ③:AVERAGE(F10, F12)?

---
[3-3]
食品のBの日別遷移, Transition of B of Food according to day
:時間軸は右に経過する, The time axis passes right.


:A-Bは連動するか?, Does A-B synchronize?

---
[3-4]
電機のBの日別遷移, Transition of B of Electric Appliances according to day
:時間軸は右に経過する, The time axis passes right.


~~~
案)B値が大きいと、A値の変化に対して反動が多い?
The idea)Are there a lot of reactionaries compared with the change of A value when B value is large?
~~~
end

2008年12月10日水曜日

廃棄物X+2号s:230MB MOs

(2008/11/03:廃棄済み)

メディアは尽きない、、、

経験年数が長いために、いろいろあります。

大学時代は、紙テープの残骸を研究室で見た。
実習は、FORTRANをカードにマークシート、その後、パンチ機の穴あけで、1ステートメントごとです。
電算室には、オンライン端末があったが、システム内部の人間にしかアクセスが許可されていなかった。

自宅ではBASICマシン。ソースの保存はラジカセの音楽テープ。300bps程度だったような。
ロード時、ボリュームの調整が微妙でした。
当時、パソコンショップで、コモドール、ペット?、ラジオシャークなどが羨望の的でした。
会社では、8インチFDからお世話になっております。ディスクコピーは音的になごみます?
3.5インチFD
3.5インチMO:こちらも容量が年々?アップした。
5インチMO:レコードのようにA面、B面あります。
スーパーFD:名前をよく覚えていない。FDの何十枚分?
リムーバブルディスク
CD-RW:よく使ったが、書き込み速度が遅すぎた。
USB256M:出始めはやはり高い。数万?
そして、今は、、、

モデムと同様、メディアへのお付き合いにもコストがかかる時代でした。

~~~
今回、MOメディアの廃棄です。
MOドライブは、MOメディアの残骸があるかもしれないので、確認後、最後に廃棄予定。


---
在庫)3.5インチMOディスク
TEIJIN(230MB):67枚
MITSUBISHI KASEI(230MB):4枚
MITSUBISHI KASEI(128MB):1枚
SONY(128MB):1枚
TOSOH(230MB):1枚
計:74枚
ちなみに、全容量=16.816GB

---
MOドライブは、以下。
Logitec Eclace LMO-406

---
既に初期化済みであるが、SCSIカードの初期化動作を確認する。

CPU:AMD K6 200MHz
OS:RedHat 9
SCSIカード:AHA 1540CF/1542CF

SCSIカードのBIOSで初期化してみる。※SCSIカードでの処理なので、CPU速度に基本的に依存しない?

TEIJIN(230MB)で、11M13s以内。結構かかります。

---
MOの値段があったので、とりあえず、記録。

TEIJIN(230MB)
10枚セット
ヨドバシカメラ
5,280円

でした。
~~~
end

廃棄物X+1号s:Disk Changer

(2008/11/02:廃棄済み)

SCSI-DISK切替器。自作品です。
SCSI-DISKの電源ONを切り替えることで、PCに接続するDISKを切り替えるもの。

1995年当時、SCSI-DISKの容量は、1G前後。

1Gのディスクを購入するために、
狸小路界隈の九十九電機には在庫なしで、同じ通りにあった「マハーポーシャ」で購入。
うわさは聞いていたが、始めて店内に入った。地下鉄サリン事件の2カ月前(1995/01/xx)。
2年前、日比谷線で通勤していたので、転職していなかったら、事件に巻き込まれた可能性大。

ちなみに、駅前に移転した九十九電機、倒産です(2008/10/30)。
~~~
4台のDISKを内蔵。2台組で、切り替えた。
---
[Switch=A]none
---
[Switch=B]
Vendor:QUANTUM Model:EMPIRE_1080S Rev:1220 sda (1080MB)
Vendor:CONNER Model:CFP1060S 1.05GB Rev:2035 sdb (1062MB)
---
[Switch=C]
Vendor:QUANTUM Model:LIGHTNING 730S Rev:241E sda (733MB)
Vendor:QUANTUM Model:LIGHTNING 730S Rev:241E sdb (733MB)
---
[Switch=D]none

~~~
用途)
RedHatの初代リリースのインストールテスト。秋葉から通販でCDを購入。
Modula-3の動作確認。結局、ソースからコンパイルし、サンプルを試しただけに終わった。
DB:Informix(Linux版)のインストール確認。
BSD/OSのインストール
FreeBSDのインストール
:このころは、インストールオタクだった。オープンソースのコンパイルも楽しかった。
configure
make
make install
は定番でしたが。

~~~
廃棄前に中身を記念に。


:前面に穴をあけて、ロータリースイッチをつけた。ABCD切替用。


:SCSI-50pinフルのコネクタです。特に改造なし。



:SCSI-DISK4台内蔵。


:ロータリースイッチの裏のケーブル接続状況。





:切替SWを取り出した。

~~~
end

廃棄物X号s:Cable Changer(2),HUB,Modem,Telephone exchange?,Auto Power Boot

久々に整理開始。
~~~
[1]ケーブル切替器(1),Cable Changer

SANWA SUPPLY
MULTI CHANGER SW-R41N, RS-232C 4WAY TYPE

RS-232Cケーブルの切替器で、1対4ポート。
ボーレートの異なる複数のモデム(過去の資産、テスト用)の切替に使用した。

~~~
[2]ケーブル切替器(2),Cable Changer

Justy
JS-153
DATA TRANSFER SWITCH

D-sub 15pinケーブルの切替器で、1対3ポート。
こちらもモデムなどの通信機器の切替用に使用。

~~~
[3]ハブ,HUB

メーカー不明
HT-810TP, 8-PORTS HUB

ネットワークケーブルは、3種類に対応。
10BASE-5, 10BASE-2, 10BASE-T
現在は、10BASE-Tが現存。
10BASE-5は、イエローケーブルでお世話になりました。

~~~
[4]モデム,Modem

SUNTAC,サン電子株式会社
MS144AVF
DATA/FAX/VOICE INTELLIGENT MODEM

14400bpsのモデムです。
モデムは、あらかた廃棄したが、まだ残っていました。
パソコン通信(NIFTY)で使用していたもの。

~~~
[5]内線電話交換?,Telephone exchange

株式会社アスキー,ユニオンブロス産業株式会社
三角電伝
3KD001

内線電話ができるもの。外線1に対して、内線2です。3回線なので、三角?
PC間のファイル転送にモデムを使用していたので、そのためにつないでいたもの。
まだ、PCのLANは普及していなかったので、ファイル移動は、FDしかなかったが、
形式が異なりFDで読めない場合もあり、1200bpsとか2400bpsの低速であるが、
圧縮してモデムで転送した。
PC間通信には、RS-232Cのクロスケーブルという話もありました、、、

~~~
[6]PC自動起動,Auto Power Boot

PC側にISAボードがあったが、既に廃棄。

SYSTEM SACOM CORP.
AC POWER CONTROLLER
ポン太,SPC-02

自宅のPCにアクセスするために、電話が着信すると、PCの電源を入れる構成で
使用したもの。
現在のように、ネットの常時接続ではなく、モデムによる時間課金なので、
電話料金対策はシビアでした。

---
MetaGEN Intternational Inc.
PC ENERGY SAVER
キーボードに接続するタイプ。
DOSコマンドで、コンセントに接続したものの電源をOFFできた。

終了時に、プログラムで、PCの電源を自動的に落とすためのものに使用。
OSは、基本的に、DOSなので、現在のように、shutdownはなかった。

~~~
end

2008年12月3日水曜日

Power Law:(30)STOCK,TSE,Tendency to B

(2008/12/15)
終値に誤り発見!Power Law:(33)を参照。再トライ,TODO)
Mistakes to the closing share price!Refer to Power Law:(33). Re-trial and TODO)

http://humanbeing-etcman.blogspot.com/2008/12/power-law33stocklast-price-incorrect.html

(2008/12/01)
~~~
[1]食品, Foods

Power Law:(29)
http://humanbeing-etcman.blogspot.com/2008/12/power-law29stocktsefoods200811xx2.html


F12(144)と、(F11 + F13)/2=161 の間に分布している?
Is it distributed F12(144) and between (F11 + F13)/2=161?


~~~
[2]電機, Electric Appliances

Power Law:(28)
http://humanbeing-etcman.blogspot.com/2008/12/power-law28stocktseelectric.html

F11(89)と、(F10 + F12)/2=99.5 の間に分布している?
Is it distributed F11(89) and between (F10 + F12)/2=99.5?


~~~
[図の考察, Consideration of figure]
B値は、特定のフィボナッチ数の値と、その前後のフィボナッチ数の平均値に挟まれた範囲に分布している?
B value is distributed within the range placed between mean values of the number of Fibonacci before and behind that the values of a specific number of Fibonacci?
~~~
end

Power Law:(29)STOCK,TSE,Foods,200811xx(2)

(2008/12/15)
終値に誤り発見!Power Law:(33)を参照。再トライ,TODO)
Mistakes to the closing share price!Refer to Power Law:(33). Re-trial and TODO)

http://humanbeing-etcman.blogspot.com/2008/12/power-law33stocklast-price-incorrect.html

(2008/11/28)
http://humanbeing-etcman.blogspot.com/2008/12/power-law27stocktsefoods200811xx.html

Power Law:(27)で、A-Bの傾き=5以外の場合を見たが、
傾き=5での傾向線を見る。
The trend line by inclination =5 is seen though the cases of except inclination =5 of A-B were seen with Power Law:(27).


~~~
end

Power Law:(28)STOCK,TSE,Electric Appliances,200811xx

(2008/12/15)
終値に誤り発見!Power Law:(33)を参照。再トライ,TODO)
Mistakes to the closing share price!Refer to Power Law:(33). Re-trial and TODO)

http://humanbeing-etcman.blogspot.com/2008/12/power-law33stocklast-price-incorrect.html

(2008/11/25-2008/12/02)

新聞の株式欄を見る。
東証第1部の電機関係の会社で、終値のべき乗傾向の推移をチェックする。

期間:20081105, 20081107, 20081114, 20081117, 20081119, 20081121, 20081126

The stocks column of the newspaper is seen.
The transition of the Power Law tendency to Last price is checked in the company related to the Electric Appliances of first section of the Tokyo Stock Exchange.

Period: 20081105, 20081107, 20081114, 20081117, 20081119, 20081121, 20081126

データは以下を使用。
Data is used as follows.
http://humanbeing-etcman.blogspot.com/2008/12/stock5tseelectric-applianceslast.html

~~~
方法,Method)
食品に同じ。以下を参照。
It is the same as food. The following references.
http://humanbeing-etcman.blogspot.com/2008/12/power-law27stocktsefoods200811xx.html


~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[20081105,Wed]
[1] 20081105,Company=order:1 to 161,ALL

LN(y) = -1.0735*LN(x) + 10.79 , R^2=0.8323


--
[2]
order:1 to 010 R^2=0.9647
order:1 to 020 R^2=0.9749
order:1 to 030 R^2=0.9837
order:1 to 040 R^2=0.9863
order:1 to 045 R^2=0.9877
order:1 to 046 R^2=0.988
order:1 to 047 R^2=0.9881
order:1 to 048 R^2=0.9881
order:1 to 049 R^2=0.9882

order:1 to 050 R^2=0.9883 <--:best

order:1 to 051 R^2=0.9882
order:1 to 052 R^2=0.988
order:1 to 053 R^2=0.9879
order:1 to 054 R^2=0.9878
order:1 to 055 R^2=0.9877
order:1 to 056 R^2=0.9876
order:1 to 057 R^2=0.9876
order:1 to 058 R^2=not checked.
order:1 to 059 R^2=not checked.
order:1 to 060 R^2=0.9874
order:1 to 065 R^2=0.988
order:1 to 070 R^2=0.988
order:1 to 080 R^2=not checked.
order:1 to 090 R^2=not checked.
order:1 to 100 R^2=not checked.
order:1 to 110 R^2=not checked.
order:1 to 120 R^2=not checked.
order:1 to 130 R^2=not checked.
order:1 to 140 R^2=not checked.
order:1 to 150 R^2=not checked.
order:1 to 161 R^2=0.8323


---
[3] 20081105,Company=order:1 to 050

LN(y) = -0.7086*LN(x) + 9.7297 , R^2=0.9883

LN
(x, y)=(1, exp(9.7297)),(50, exp(-0.7086*LN(50) + 9.7297))
=(1, 16809.51),(50, 1051.15)
Relative value)
=(1, 16809.51/1051.15),(50, 1)
=(1, 15.992),(50, 1)
A=-0.7086
B=15.992

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[20081107,Fri]
[1] 20081107,Company=order:1 to 161,ALL

LN(y) = -1.0677*LN(x) + 10.667 , R^2=0.8273

---
[2] 20081107,Company=order:1 to xxx
R^2 is checked.

order:1 to 010 R^2=0.9664
order:1 to 020 R^2=0.9763
order:1 to 030 R^2=0.9852
order:1 to 040 R^2=0.9882
order:1 to 050 R^2=0.9902
order:1 to 055 R^2=0.9906
order:1 to 056 R^2=0.9907
order:1 to 057 R^2=0.9909
order:1 to 058 R^2=0.9909
order:1 to 059 R^2=0.991
order:1 to 060 R^2=0.9911
order:1 to 061 R^2=0.9911
order:1 to 062 R^2=0.9911
order:1 to 063 R^2=0.9912
order:1 to 064 R^2=0.9912

order:1 to 065 R^2=0.9913 <--:best
order:1 to 066 R^2=0.9913 <--:best,selected.

order:1 to 067 R^2=0.9911
order:1 to 068 R^2=0.9909
order:1 to 069 R^2=0.9907
order:1 to 070 R^2=0.9906
order:1 to 080 R^2=0.9898
order:1 to 090 R^2=0.9865
order:1 to 100 R^2=0.9845
order:1 to 110 R^2=not checked.
order:1 to 120 R^2=not checked.
order:1 to 130 R^2=not checked.
order:1 to 140 R^2=not checked.
order:1 to 150 R^2=not checked.
order:1 to 161 R^2=0.8273

---
[3] 20081107,Company=order:1 to 066

LN(y) = -0.7341*LN(x) + 9.6782 , R^2=0.9913

LN
(x, y)=(1, exp(9.6782)),(66, exp(-0.7341*LN(66) + 9.6782))
=(1, 15965.73),(66, 736.99)
Relative value)
=(1, 15965.73/736.99),(66, 1)
=(1, 21.663),(66, 1)
A=-0.7341
B=21.663

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[20081114,Fri]
[1] 20081114,Company=order:1 to 161,ALL

LN(y) = -1.0724*LN(x) + 10.648 , R^2=0.8276

---
[2] 20081114,Company=order:1 to xxx
R^2 is checked.

order:1 to 010 R^2=0.9643
order:1 to 020 R^2=0.9769
order:1 to 030 R^2=0.986
order:1 to 040 R^2=0.988
order:1 to 050 R^2=0.9901
order:1 to 051 R^2=not checked.
order:1 to 052 R^2=not checked.
order:1 to 053 R^2=not checked.
order:1 to 054 R^2=0.9906
order:1 to 055 R^2=0.9908
order:1 to 056 R^2=0.9909
order:1 to 057 R^2=0.9909
order:1 to 058 R^2=0.991

order:1 to 059 R^2=0.9911 <--:best

order:1 to 060 R^2=0.991
order:1 to 061 R^2=0.991
order:1 to 062 R^2=0.9909
order:1 to 063 R^2=0.9907
order:1 to 064 R^2=0.9907
order:1 to 065 R^2=0.9907
order:1 to 070 R^2=0.9907
order:1 to 080 R^2=0.9887
order:1 to 090 R^2=not checked.
order:1 to 100 R^2=not checked.
order:1 to 110 R^2=not checked.
order:1 to 120 R^2=not checked.
order:1 to 130 R^2=not checked.
order:1 to 140 R^2=not checked.
order:1 to 150 R^2=not checked.
order:1 to 161 R^2=0.8276

---
[3] 20081114,Company=order:1 to 059

LN(y) = -0.7255*LN(x) + 9.6249 , R^2=0.9911

LN
(x, y)=(1, exp(9.6249)),(59, exp(-0.7255*LN(59) + 9.6249))
=(1, 15137.04),(59, 785.76)
Relative value)
=(1, 15137.04/785.76),(59, 1)
=(1, 19.264),(59, 1)
A=-0.7255
B=19.264

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[20081117,Mon]
[1] 20081117,Company=order:1 to 161,ALL

LN(y) = -1.0774*LN(x) + 10.673 , R^2=0.8264

---
[2] 20081117,Company=order:1 to xxx
R^2 is checked.

order:1 to 010 R^2=0.9654
order:1 to 020 R^2=0.9786
order:1 to 030 R^2=0.9869
order:1 to 040 R^2=0.9872
order:1 to 050 R^2=0.9896
order:1 to 060 R^2=0.9908
order:1 to 061 R^2=0.9909
order:1 to 062 R^2=0.9909

order:1 to 063 R^2=0.991 <--:best

order:1 to 064 R^2=0.9909
order:1 to 065 R^2=0.9909
order:1 to 066 R^2=0.9908
order:1 to 067 R^2=0.9808
order:1 to 068 R^2=0.9908
order:1 to 069 R^2=0.9909
order:1 to 070 R^2=0.9909
order:1 to 071 R^2=0.9908
order:1 to 072 R^2=0.9903
order:1 to 073 R^2=0.9899
order:1 to 074 R^2=0.9896
order:1 to 075 R^2=0.9893
order:1 to 080 R^2=0.9882
order:1 to 090 R^2=0.986
order:1 to 100 R^2=0.9838
order:1 to 110 R^2=0.9796
order:1 to 120 R^2=0.9682
order:1 to 130 R^2=0.9539
order:1 to 140 R^2=0.9373
order:1 to 150 R^2=0.9198
order:1 to 161 R^2=0.8264

---
[3] 20081117,Company=order:1 to 063

LN(y) = -0.7287*LN(x) + 9.6426 , R^2=0.991

LN
(x, y)=(1, exp(9.6426)),(63, exp(-0.7287*LN(63) + 9.6426))
=(1, 15407.35),(63, 752.57)
Relative value)
=(1, 15407.35/752.57),(63, 1)
=(1, 20.473),(63, 1)
A=-0.7287
B=20.473

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[20081119,Wed]
[1] 20081119,Company=order:1 to 161,ALL

LN(y) = -1.0714*LN(x) + 10.627 , R^2=0.8346

---
[2] 20081119,Company=order:1 to xxx
R^2 is checked.

order:1 to 010 R^2=0.9604
order:1 to 020 R^2=0.9786
order:1 to 030 R^2=0.9865
order:1 to 040 R^2=0.9896
order:1 to 050 R^2=0.9914
order:1 to 055 R^2=0.9923
order:1 to 056 R^2=0.9924

order:1 to 057 R^2=0.9926 <--:best

order:1 to 058 R^2=0.9925
order:1 to 059 R^2=0.9924
order:1 to 060 R^2=0.9923

order:1 to 061 R^2=0.9923
order:1 to 062 R^2=0.9922
order:1 to 063 R^2=0.992
order:1 to 064 R^2=0.9919
order:1 to 065 R^2=0.9918
order:1 to 070 R^2=0.9915
order:1 to 080 R^2=0.9898

order:1 to 090 R^2=not checked.
order:1 to 100 R^2=not checked.
order:1 to 110 R^2=not checked.
order:1 to 120 R^2=not checked.
order:1 to 130 R^2=not checked.
order:1 to 140 R^2=not checked.
order:1 to 150 R^2=not checked.
order:1 to 161 R^2=0.8346

---
[3] 20081119,Company=order:1 to 057

LN(y) = -0.7245*LN(x) + 9.6079 , R^2=0.9926

LN
(x, y)=(1, exp(9.6079)),(57, exp(-0.7245*LN(57) + 9.6079))
=(1, 14881.89),(57, 795.29)
Relative value)
=(1, 14881.89/795.29),(57, 1)
=(1, 18.713),(57, 1)
A=-0.7245
B=18.713

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[20081121,Fri]
[1] 20081121,Company=order:1 to 161,ALL

LN(y) = -1.0657*LN(x) + 10.587 , R^2=0.8306

---
[2] 20081121,Company=order:1 to xxx
R^2 is checked.

order:1 to 010 R^2=0.9764
order:1 to 020 R^2=0.9871
order:1 to 030 R^2=0.9921
order:1 to 040 R^2=0.9939
order:1 to 050 R^2=0.9949
order:1 to 055 R^2=0.9954
order:1 to 056 R^2=0.9955
order:1 to 057 R^2=0.9956
order:1 to 058 R^2=0.9957

order:1 to 059 R^2=0.9958 <--:best

order:1 to 060 R^2=0.9956
order:1 to 061 R^2=0.9954
order:1 to 062 R^2=0.9953
order:1 to 063 R^2=0.9951
order:1 to 064 R^2=0.9949
order:1 to 065 R^2=0.9946
order:1 to 070 R^2=0.994
order:1 to 080 R^2=0.9912
order:1 to 090 R^2=0.9873
order:1 to 100 R^2=0.9842
order:1 to 110 R^2=not checked.
order:1 to 120 R^2=not checked.
order:1 to 130 R^2=not checked.
order:1 to 140 R^2=not checked.
order:1 to 150 R^2=not checked.
order:1 to 161 R^2=0.8306

---
[3] 20081121,Company=order:1 to 059

LN(y) = -0.7098*LN(x) + 9.5407 , R^2=0.9958

LN
(x, y)=(1, exp(9.5407)),(59, exp(-0.7098*LN(59) + 9.5407))
=(1, 13914.68),(59, 770.06)
Relative value)
=(1, 13914.68/770.06),(59, 1)
=(1, 18.07),(59, 1)
A=-0.7098
B=18.07

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[20081126,Wed]
[1] 20081126,Company=order:1 to 161,ALL

LN(y) = -1.0705*LN(x) + 10.613 , R^2=0.8256

---
[2] 20081126,Company=order:1 to xxx
R^2 is checked.

order:1 to 010 R^2=0.9805
order:1 to 020 R^2=0.9882
order:1 to 030 R^2=0.9924
order:1 to 040 R^2=0.9922
order:1 to 050 R^2=0.9939
order:1 to 051 R^2=0.994
order:1 to 052 R^2=0.9941
order:1 to 053 R^2=0.9942
order:1 to 054 R^2=0.9943
order:1 to 055 R^2=0.9944
order:1 to 056 R^2=0.9945

order:1 to 057 R^2=0.9946 <--:best

order:1 to 058 R^2=0.9945
order:1 to 059 R^2=0.9945
order:1 to 060 R^2=0.9944
order:1 to 061 R^2=0.9943
order:1 to 062 R^2=0.9942
order:1 to 063 R^2=0.9939
order:1 to 064 R^2=0.9935
order:1 to 065 R^2=0.9932
order:1 to 070 R^2=0.9926
order:1 to 080 R^2=0.9907
order:1 to 090 R^2=not checked.
order:1 to 100 R^2=not checked.
order:1 to 110 R^2=not checked.
order:1 to 120 R^2=not checked.
order:1 to 130 R^2=not checked.
order:1 to 140 R^2=not checked.
order:1 to 150 R^2=not checked.
order:1 to 161 R^2=0.8256

---
[3] 20081126,Company=order:1 to 057

LN(y) = -0.7155*LN(x) + 9.5699 , R^2=0.9946

LN
(x, y)=(1, exp(9.5699)),(57, exp(-0.7155*LN(57) + 9.5699))
=(1, 14326.98),(57, 794.01)
Relative value)
=(1, 14326.98/794.01),(57, 1)
=(1, 18.044),(57, 1)
A=-0.7155
B=18.044

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
A-Bをグラフに出力。
A-B is output to the graph.

---
ファイル:powerlaw-xx-a-b.txt を使用。The file: Powerlaw-xx-a-b.txt is used.
===
A,B,date,order,lab_disp
-0.7086,15.992,"20081105,Wed",50,right
-0.7341,21.663,"20081107,Fri",66,right
-0.7255,19.264,"20081114,Fri",59,right
-0.7287,20.473,"20081117,Mon",6,right
-0.7245,18.713,"20081119,Wed",57,right
-0.7098,18.07,"20081121,Fri",59,left
-0.7155,18.044,"20081126,Wed",57,right
===

---
Rのコードは以下。The code of R is the following.

data_ab = read.csv("powerlaw-xx-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.36, -0.25), ylim=c(2.5, 3.5),
col="orange", pch=20, main="Scatter chart:a-b, Power Law(xx),Electric Appliances")
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.36, -0.25), ylim=c(2.5, 3.5), ann=F)

if (data_ab$lab_disp[i] == "right"){
text(log(abs(data_ab$A[i])), log(data_ab$B[i]), data_ab$date[i], cex=0.75, pos=4, offset=0.5)
}
else if (data_ab$lab_disp[i] == "left"){
text(log(abs(data_ab$A[i])), log(data_ab$B[i]), data_ab$date[i], cex=0.75, pos=2, offset=0.5)
}
else{
text(log(abs(data_ab$A[i])), log(data_ab$B[i]), data_ab$date[i], cex=0.75, pos=4, offset=0.5)
}
par(new=T)
}

#abline(v = log(0.6), col="red")
abline(v = -0.26, col="red")

###
abline(log(144), 5, col="gray", lty=2)
abline(log((55+144)/2), 5, col="gray", lty=2)
abline(log((34+144)/2), 5, col="gray", lty=2)
abline(log((34+55+144+144)/4), 5, col="gray", lty=2)
#???
abline(log(180), 7, col="gray", lty=2)
#???
abline(log(350), 9, col="gray", lty=2)



~~~
[図への考察,Consideration to figure]
ラインとしては、週初めのライン、週中のライン、週末のラインがあるのか?いろいろ妄想する。
Does the line include the line on the line at the time of beginning of the week, the line during the week, and the weekend?Various delusion.
~~~
end