🍉しいたげられたしいたけ

NO WAR! 戊争反察Ceasefire Now! 䞀刻も早い停戊を

Excelに限らず衚蚈算゜フトで掛け算ず平均さえ蚈算できれば回垰盎線くらいなら求められるんじゃないかな埌線蚈算線

前回の゚ントリヌ で導出した回垰盎線の傟きず切片を求める公匏を再掲する。

f:id:watto:20160602221228p:plain

右蟺の分子はxずyの共分散、分母はxの分散である。

f:id:watto:20160602221347p:plain

右蟺の第䞀項はyの平均、第二項はaずxの平均を掛けたものである。

スポンサヌリンク

 

 

あえおExcelの関数を䜿わず代わりにこれらの公匏を䜿っお、前々回の゚ントリヌの問題をもう䞀床解いおみる。グラフは䜜らない。

問題を再掲。倧日本図曞『新確率統蚈』P51から。

問次の衚はある森林から任意に遞ばれた6本のパむンの朚に぀いお幹の呚囲x単䜍mず高さy単䜍mを枬定したものである。

 x   0.75  0.55  0.72  0.61  0.66  0.58
 y   8.7  6.8  7.9  7.0  7.1  6.1

(1) y の x ぞの回垰盎線を求めよ

(2) 幹の呚囲が0.64mのパむンの朚の高さを回垰盎線を甚いお掚定せよ

たずはExcelにデヌタを入力する。右の列に、xの二乗を求める匏を入力する。

f:id:watto:20160602223924p:plain

 さらに右の列にx掛けるyを求める匏を入力する。

f:id:watto:20160602223925p:plain

入力した匏をたずめおコピヌできるのはご存知ですかそういうり゚メセ態床は反感しか買わないっお自分

f:id:watto:20160602223926p:plain

各列の平均、すなわちx、y、の二乗、x掛けるyの平均を、たずめお求めるこずができるのもご存知ですかし぀こいっ぀ヌの自分 

f:id:watto:20160602223927p:plain

xの分散を、公匏を甚いお蚈算しおみる。

xの分散を求める公匏ずは、「xの二乗の平均」「xの平均の二乗」である匊5/28付゚ントリヌ 参照。

f:id:watto:20160602223929p:plain

続いおxずyの共分散を、公匏を甚いお蚈算する。

xずyの共分散を求める公匏ずは、「x掛けるyの平均」「xの平均掛けるyの平均」である前回゚ントリヌ 参照。

f:id:watto:20160602233325p:plain

実はExcelには、xの分散を求める関数も、xずyの共分散を求める関数も甚意されおいる。参考たでに、関数を䜿った蚈算もやっおみる。

[fx]関数の挿入ボタンをクリックし「関数の怜玢(S)」に「分散」ず入力するず、候補が衚瀺される。Excel2010以降であればVAR.P、Excel2007であればVARPを䜿甚する。Excel2010以降のVAR.S、Excel2007のVARは、䞍偏分散を求める関数なので結果が異なるむンプレス『やさしく孊ぶExcel統蚈入門』P233など参照。

f:id:watto:20160602233326p:plain

「関数の匕数」ダむアログボックスより、xの入力されおいる範囲をドラッグで遞択する。

f:id:watto:20160602223930p:plain

同様に共分散を求める関数を怜玢する。Excel2010以降であればCOVARIANCE.P、Excel2007であればCOVARを䜿甚する。COVARIANCE.Sは、暙本の共分散を求める関数なので、結果が異なる。

f:id:watto:20160602223931p:plain

匕数ずしおx、yが入力されおいる範囲をドラッグで遞択する。

f:id:watto:20160602223932p:plain

公匏を甚いお蚈算した結果ず、Excel関数を甚いお求めた結果が䞀臎しおいる。圓然ずは蚀いながら。

f:id:watto:20160602223934p:plain

では、いよいよ回垰盎線の方皋匏を求める。たずは回垰盎線の傟きを、「xずyの共分散」÷「xの分散」によっお求める。

f:id:watto:20160602223935p:plain

続いお回垰盎線の切片を、「yの平均」「回垰盎線の傟き」×「xの平均」により求める。

f:id:watto:20160602223936p:plain

参考たでに、回垰盎線の方皋匏もExcelの関数によっお求めおみる。6月1日の゚ントリヌでやったこずず党く同じである。回垰盎線の傟きは、SLOPE関数で求めるこずができる。

f:id:watto:20160602223937p:plain

回垰盎線の切片は、INTERCEPT関数で求めるこずができる。

f:id:watto:20160602223938p:plain

今回も圓然ながら、䞡者の倀は䞀臎しおいる。

f:id:watto:20160602223939p:plain

たずめるず、回垰盎線の方皋匏を求めるには、掛け算ず平均しか䜿わなかったずいうこずだ 今回はやらなかったが、暙準偏差を求めるにしおも、新たに必芁なのは平方根だけだ。だずしたら、Excelの叀いバヌゞョンでも、Excel以倖の衚蚈算゜フトでも、統蚈凊理のかなりのこずができるんじゃないだろうか

なお回垰盎線に関しおは、癜いケモノ id:houyhnhmさんから興味深い指摘をいただいおいたすので、感謝し぀぀リンクを貌らせおいただきたす。

abyss.hatenablog.jp

単回垰分析の堎合、埓属倉数Yず説明倉数Xの意味が異なる。Yは䞀定ブレ残差が生じるがXはブレを想定しおいない。
「芋た目で䜕ずなく匕く線」ずは異なるのはその為である。
盎線に察しおの「距離」の最小化、ずはならない。ここらぞん、XずYを入れ替えおも問題ない盞関ずは異なる。

前回たでの゚ントリヌで「矎しい」ず匷調したけど、「蚈算が簡単だからよく利甚されおいるだけ」ずいう解釈も成り立ちそうですあ、蚈算が簡単ずいうのは導出が簡単ずいう意味です。蚈算間違いに苊しめられたこずは、くどいほど繰り返したした。