#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)