1

私は（私にとって）非常に洗練された問題を抱えています。特定のデータセットの補間を計算するコードを作成しました。その後、赤池情報量基準を計算して、どの補間スイートが最適かを確認します。

そのように見えます：

    polyfit=np.polyfit(x,y,5)  # x,y are data set and 5 is a root of polynomial
    poly1d=np.poly1d(polyfit)    
    print poly1d  #show final polynomial

my=[]
for i in x: 
x_=poly1d(i)
my.append(x_) #calculate list of values 

def AIC(i,j): 
for i in y:
for j in my:
    RSS=(i-j)**2

AIC=36-np.log(RSS) 
print AIC

多項式の根を変更したい場合は、コードを変更する必要があるため、このコードは好きではありません。最初に for ループを使用する必要があることはわかっています。

for i in xrange(40):
    polyfit=np.polyfit(x,y,i)

しかし、すべての多項式をリストに保存する方法がわかりません。それを使用して任意のルートの def AIC() を計算するとわかっていれば、データに最適な曲線を見つけることができます。みんなを助けてください、それはすでに2週間私を悩ませています. ご不明な点がございましたら、お気軽にお尋ねください。

4

0 に答える 0