裏 RjpWiki

Julia ときどき R, Python によるコンピュータプログラム,コンピュータ・サイエンス,統計学

確率という概念

2011年07月05日 | 統計学

一般の人は,確率というのを実感できないのだなあと実感した例。

http://sankei.jp.msn.com/politics/news/110705/lcl11070522060009-n1.htm

関西電力は5日、電力使用ピーク時に突発的な停電を避けるため、5%の電力供給余力を持たせようとすると家庭用エアコンの3分の1強を停止する必要がある-との試算を公表した。ただ、この節電を実行するには、例えば1台しかエアコンがない家庭の場合はどうするのかといった問題があり、この回答をもとに大阪市は、節電策としてエアコン停止を強調する大阪府の橋下徹知事に対し「現実的でない」と疑問を投げかけた。

解決策は,もしそのような事態が発生したら,例えば自宅の電話番号(世帯主の誕生日などでもよいだろう)が3で割り切れる,3で割ったら1余る,2余るのいずれかに該当した場合に運が悪いと思ってエアコンを切ってもらうようにすればよい。剰余が0,1,2のいずれに該当する場合に対象になるかは,要請時に通知する(要請があるたびにローテーションさせる)。まあ,要請されても従わない人もいるだろうから,偶数か奇数かで最高1/2のエアコン停止というようにしておくと安全かも知れない。

コメント
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする

R-de-SEM

2011年07月05日 | 統計学

共分散構造分析
http://www.ic.nanzan-u.ac.jp/~kamiya/r/content/sem.html
について,Graphviz によるパス図の描画で 2 点指摘

> path.diagram(sem.cancer, out.file="cancer.dat",

のように,第 2 引数を out.file= というように名前付きで指定すると図を描いてくれない(古いバージョンは out.file= だったのだろうか)。
path.diagram の第 2 引数は名前を持たないので,
path.diagram(sem.cancer, "cancer.dat",
のようにしないとだめだ。

> 標準化解の係数を図中に示すなら、値をここで書き直しておかないといけないですね。Graphviz は sem のためのソフトではなく、汎用の描画ソフトですから仕方ないですね(ひょっとすると、オプションで設定できるのかも知れませんが)。

これも,オンラインヘルプを良く読むことで解決される。

引数 standardize で指定できる。すなわち "if TRUE, display standardized coefficients; default is FALSE."

コメント
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする

PVアクセスランキング にほんブログ村

PVアクセスランキング にほんブログ村