#author("2018-08-28T19:08:29+09:00","default:Miyashita","Miyashita")
#author("2018-08-28T19:09:28+09:00","default:Miyashita","Miyashita")
*Julia演習問題:線形回帰分析 [#wee9972d]
***問題 [#f75263a9]
+北太平洋平均風速データから最小2乗法を用い, 1次および2次回帰の係数を求めなさい.
#codeprettify(lang-julia){{
DelimitedFiles.readdlm(filename) # v1.0.0以降, Base.readdlm → DelimitedFiles.readdlmに
DelimitedFiles.readdlm(filename) # v0.7.0以降, Base.readdlm → DelimitedFiles.readdlmに
coef1 = Polynomials.polyfit(x, y, 1)
coef2 = Polynomials.polyfit(x, y, 2)
}}
+求めた回帰式より,2100年まで平均風速を外挿し,生データとともにまとめてプロットしなさい.
#codeprettify(lang-julia){{
plt = Plots.scatter(years, winddata,lab="Raw data")
plt = Plots.plot!(plt, years, polyval(...), lab="Linear")
}}
+プロットした図のx, y軸の座標の名前と線の種類を明記しなさい.
#codeprettify(lang-julia){{
plt = Plots.plot!(plt, label="XXX", legend=:topleft, legendfont=font(14,"sans-serif"),
                 ylabel="Y-Label", guidefont=("sans-serif",16), ... )
}}
***データ [#dd7ba937]
-入手先~
#ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/image/python/windspeed.dat,過去50年間の北太平洋平均風速)
-フォーマット
--テキスト形式,空白区切り,3カラム(年 平均風速 標準偏差)
***解答例(参考) [#lfd55a07]
#ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/image/julia/regression.png,600x450)

Front page   Edit Diff Attach Copy Rename Reload   New List of pages Search Recent changes   Help   RSS of recent changes