#author("2018-03-07T05:38:52+09:00","default:Miyashita","Miyashita") #author("2018-07-25T23:22:54+09:00","default:Miyashita","Miyashita") *Python演習問題：線形回帰分析 [#wee9972d] ***インポート・エイリアス設定 [#xef27201] import numpy as np import matplotlib.pyplot as plt from numpy import genfromtxt #codeprettify(lang-python){{ import numpy as np import matplotlib.pyplot as plt from numpy import genfromtxt }} ***問題 [#f75263a9] +北太平洋平均風速データから最小2乗法を用い, 1次および2次回帰の係数を求めなさい． np.loadtxt または genfromtxt np.polyfit #codeprettify(lang-python){{ np.loadtxt または genfromtxt np.polyfit }} +求めた回帰式より，2100年まで平均風速を外挿し，生データとともにまとめてプロットしなさい．また，グラフにグリッド線を描きなさい． np.polyval ax.plot ax.grid(which='major',color='#xxxxxx',linestyle='--') #codeprettify(lang-python){{ np.polyval ax.plot ax.grid(which='major',color='#xxxxxx',linestyle='--') }} +プロットした図のx, y軸の座標の名前と線の種類を明記しなさい． ax.set_xlabel ax.set_ylabel ax.legend([xxx,xxx,xxx],loc=num) #codeprettify(lang-python){{ ax.set_xlabel ax.set_ylabel ax.legend([xxx,xxx,xxx],loc=num) }} ***データ [#dd7ba937] -入手先~ #ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/image/python/windspeed.dat,過去５０年間の北太平洋平均風速) -フォーマット --テキスト形式，空白区切り，３カラム（年 平均風速 標準偏差） ***解答例（参考） [#lfd55a07] #ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/image/python/regression_example.png)