数学 - R言語 - 統計学 - データ解析 - 誤差 - 誤差を適切に示す(分割したモデル、回帰式、RMS誤差(Root Mean Square(二乗平均平方)、シグマ、残差標準誤差))

Head Firstデータ解析
原著

学習環境

• Surface 3 (4G LTE)、Surface 3 タイプ カバー、Surface ペン(端末)
• Windows 10 Pro (OS)
• Nebo(Windows アプリ)
• iPad Pro + Apple Pencil
• MyScript Nebo(iPad アプリ)
• Numbers - Apple(表計算ソフト)
• Excel(Office 365 Solo)
• R言語(プログラミング言語)
• 参考書籍
  • Pythonからはじめる数学入門
  • JavaScriptリファレンス 第6版
  • インタラクティブ・データビジュアライゼーション

Head Firstデータ解析 (Michael Milton (著)、大橋 真也 (監訳)、 木下 哲也 (翻訳)、オライリージャパン)の11章(誤差 - 誤差を適切に示す)、エクササイズ(p. 353)を取り組んでみる。

エクササイズ(p. 353)

コード(Emacs)

R言語

#!/usr/bin/env Rscript

url <- "http://www.oreilly.co.jp/pub/9784873114644/hfda_ch10_employees.csv"
employees <- read.csv(url, header=TRUE)

myLmBig <- lm(received[negotiated==TRUE & requested > 10]~
                  requested[negotiated==TRUE & requested > 10],
              data=employees)
myLmSmall <- lm(received[negotiated==TRUE & requested <= 10]~
                  requested[negotiated==TRUE & requested <= 10],
                data=employees)

plot(employees$requested[employees$negotiated==TRUE],
     employees$received[employees$negotiated==TRUE])

## 積極的な交渉者の回帰式の係数とRMS誤差
summary(myLmSmall)$coefficients
summary(myLmSmall)$sigma
## 臆病な交渉者の回帰式の係数とRMS誤差
summary(myLmBig)$coefficients
summary(myLmBig)$sigma

入出力結果(Terminal)

$ ./sample1.r
                                                 Estimate Std. Error   t value
(Intercept)                                     0.7933468 0.22472009  3.530378
requested[negotiated == TRUE & requested <= 10] 0.9424946 0.03151835 29.903041
                                                     Pr(>|t|)
(Intercept)                                      4.378156e-04
requested[negotiated == TRUE & requested <= 10] 6.588020e-134
[1] 1.374526
                                               Estimate Std. Error  t value
(Intercept)                                    7.813403  1.8760371 4.164845
requested[negotiated == TRUE & requested > 10] 0.302609  0.1420151 2.130824
                                                   Pr(>|t|)
(Intercept)                                    4.997597e-05
requested[negotiated == TRUE & requested > 10] 3.457618e-02
[1] 4.544424
$