2016-06-23

2016-06-21

ベイズの公式

多めに用意していますので、講義後に差し替えます。

2016-06-17

中心極限定理

中心極限定理と、t分布など正規分布から派生する確率分布について説明しました。

相関係数

相関係数と、第二章に入って確率の定義について説明しました。

2016-06-16

2016-06-14

回帰

2変量の関係を考える回帰について説明しました。

2016-06-13

多変量解析

前回までに作った連立一次方程式の解法を多変量回帰に用います。

今日使うデータです。

都市番号 都市 気温y緯度x1経度x2標高x3
1稚内-8.0 45.42 141.68 2.8
2旭川-13.6 43.77 142.37 111.9
3札幌-9.5 43.05 141.33 17.2
4青森-5.4 40.82 140.78 3.0
5盛岡-6.7 39.70 141.17 155.2
6仙台-3.2 38.27 140.90 38.9
7金沢-0.1 36.55 136.65 26.1
8長野-5.5 36.67 138.20 418.2
9高山-7.6 36.15 137.25 560.2
10軽井沢-10.0 36.33 138.55 999.1
11名古屋-0.9 35.17 136.97 51.1
12飯田-4.7 35.52 137.83 481.8
13東京-0.4 35.68 139.77 5.3
14鳥取0.5 35.48 134.23 7.1
15京都-0.6 35.02 135.73 41.4
16広島0.2 34.37 132.43 29.3
17福岡1.5 33.58 130.38 2.5
18鹿児島2.0 31.57 130.55 4.3
19高知0.1 33.55 133.53 1.9
20那覇13.5 26.23 127.68 34.9
出典:坂元慶行、石黒真木夫、北川源四郎著「情報量統計学」共立出版株式会社

2016-06-10

標本分布

大きな標本を観測したときの、標本平均、標本分散の期待値、分散を説明しました。

分散

標本の散布度を表す数値として、分散と標準偏差について説明しました。

2016-06-07

ヒストグラムの作り方

ヒストグラムの作り方と、代表値として平均値と中央値を説明しました。

2016-06-06

call文

これまでの連立一次方程式の解法をまとめ、Subプログラムとして独立させました。

Sub 連立一次方程式(係数行列(), r)
Rem (c) Kaoru Fueda
Dim 最大値()
ReDim 最大値(r)
For i = 1 To r
  If i < r Then
    最大比 = 0
    最大行 = i
    For j = i To r
      最大値(j) = 0
      For k = i To r + 1
        If Abs(係数行列(j, k)) > 最大値(j) Then
          最大値(j) = Abs(係数行列(j, k))
        End If
      Next k
      If 最大値(j) > 0 Then
        If 最大比 < Abs(係数行列(j, i) / 最大値(j)) Then
          最大比 = Abs(係数行列(j, i) / 最大値(j))
          最大行 = j
        End If
      End If
    Next j
    If 最大行 > i Then
      For k = 1 To r + 1
        temp = 係数行列(i, k)
        係数行列(i, k) = 係数行列(最大行, k)
        係数行列(最大行, k) = temp
      Next k
    End If
Rem ここは講義では説明していません。
Rem ピボットが0の場合は最も0=0に近い行を第i行に移動して、後に無視する
    If 最大比 = 0 Then
      最小値 = 最大値(i)
      最小行 = i
      For j = i + 1 To r
        If 最小値 > 最大値(j) Then
          最小値 = 最大値(j)
          最小行 = j
        End If
      Next j
      If 最小行 > i Then
        For k = 1 To r + 1
          temp = 係数行列(i, k)
          係数行列(i, k) = 係数行列(最小行, k)
          係数行列(最小行, k) = temp
        Next k
      End If
    End If
  End If
Rem ピボットを1に
  If 係数行列(i, i) <> 0 Then
    For k = i + 1 To r + 1
      係数行列(i, k) = 係数行列(i, k) / 係数行列(i, i)
    Next k
  Else
Rem ここは講義では説明していません。
    MsgBox ("解けないので解の一つを0とします。")
    For k = i + 1 To r + 1
      係数行列(i, k) = 0
    Next k
  End If
  係数行列(i, i) = 1
Rem ピボット以外を0に
  For j = 1 To r
    If j <> i Then
      For k = i + 1 To r + 1
        係数行列(j, k) = 係数行列(j, k) - 係数行列(i, k) * 係数行列(j, i)
      Next k
      係数行列(j, i) = 0
    End If
  Next j
Next i
End Sub

2016-06-03

確率分布の独立性

確率変数の和の分散の公式の前提となる、確率変数の独立性について説明しました。

統計学I

統計学の考え方と、データの種類について説明しました。

2016-06-02