Takuya Miyashita

Julia演習問題：線形回帰分析 †

問題 †

1. 北太平洋平均風速データから最小2乗法を用い, 1次および2次回帰の係数を求めなさい．

   DelimitedFiles.readdlm(filename) # v0.7.0以降, Base.readdlm → DelimitedFiles.readdlmに
   coef1 = Polynomials.polyfit(x, y, 1)
   coef2 = Polynomials.polyfit(x, y, 2)

2. 求めた回帰式より，2100年まで平均風速を外挿し，生データとともにまとめてプロットしなさい．

   plt = Plots.scatter(years, winddata,lab="Raw data")
   plt = Plots.plot!(plt, years, polyval(...), lab="Linear")

3. プロットした図のx, y軸の座標の名前と線の種類を明記しなさい．

   plt = Plots.plot!(plt, label="XXX", legend=:topleft, legendfont=font(14,"sans-serif"),
                    ylabel="Y-Label", guidefont=("sans-serif",16), ... )

データ †

• 入手先
  
  過去５０年間の北太平洋平均風速
• フォーマット
  • テキスト形式，空白区切り，３カラム（年 平均風速 標準偏差）

解答例（参考） †