続報・正月の記録

&gt; data <span class="synStatement">&lt;-</span> read.<a class="keyword" href="http://d.hatena.ne.jp/keyword/csv">csv</a><span class="synSpecial">(</span><span class="synConstant">&quot;c:\\temp\\nenga.<a class="keyword" href="http://d.hatena.ne.jp/keyword/csv">csv</a>&quot;</span><span class="synSpecial">,</span>header=<span class="synConstant">TRUE</span><span class="synSpecial">)</span>
&gt; subdata <span class="synStatement">&lt;-</span> subset<span class="synSpecial">(</span>data<span class="synSpecial">,</span>rel==<span class="synConstant">&quot;考司&quot;</span><span class="synSpecial">)</span>
&gt; library<span class="synSpecial">(</span>polycor<span class="synSpecial">)</span>
要求されたパッケージ mvtnorm をロード中です
要求されたパッケージ sfsmisc をロード中です
&gt; polychor<span class="synSpecial">(</span>data$send<span class="synSpecial">,</span>data$recv<span class="synSpecial">)</span>
<span class="synSpecial">[</span><span class="synConstant">1</span><span class="synSpecial">]</span> <span class="synConstant">0.7943974</span>
&gt; polychor<span class="synSpecial">(</span>subdata$send<span class="synSpecial">,</span>subdata$recv<span class="synSpecial">)</span>
<span class="synSpecial">[</span><span class="synConstant">1</span><span class="synSpecial">]</span> <span class="synConstant">0.6417118</span>

全体の相関は0.8近くになったが、自分のデータだけを抽出すると0.64に下がる(笑)