For文の中にwhile文を入れているのですが、うまくFor文が回りません。

Question

For文の中にwhile文を入れているのですが、うまくFor文が回りません。
下のコードの部分に、うまくいかないポイントという形で入れ込んでいますが、For文がうまく機能しません。
この場合であると、１つ目の要素だけが計算されているようです。
コードを見ているのですが間違っているところがわからなく質問させていただきました。

import numpy as np
import scipy.optimize
import matplotlib.pylab as plt
import pandas as pd 
import csv

data = pd.read_excel('', sheet_name = None)

csvlist = []

name = ["Furrow", "Ridge"]
initial = [0.05, 0.1]

i = 4
for o, g in zip(name, initial):

    while i < 1305: 

        f = open(str(o)+'.variable.complete.csv', 'a')
        writer = csv.writer(f, lineterminator='\n')

        one_list = []
    
        a = data[o].iloc[i:i+400,1]

        d = data[o].iloc[i:i+400,2]
        c = a.tolist()
        b = d.tolist()

        plt.plot(c, b, 'o',label='Raw data')
    
        m = data[o].iloc[5:,2].max()
        s = data[o].iloc[5:,2].min()

        def fitfunc(h, x, z, n ):
            return x+((m-x)*((1+(z*h)**n)**(1/n-1)))


        para_ini =[g, 0.5, 1.5] 
    
        para_bounds=([0,-np.inf,1],[s,1,np.inf])

        para_opt, cov = scipy.optimize.curve_fit(fitfunc, c, b, para_ini,bounds = para_bounds )#フィAッティングパラメータの最適化

        print ("Fitted x =", str(para_opt[0]))
        print ("Fitted z =", str(para_opt[1]))
        print ("Fitted n =", str(para_opt[2]))
    
        x=para_opt[0]
        z=para_opt[1]
        n=para_opt[2]

        h1 = np.arange(min(c) - 1, max(c) + 1, 0.1)
        y_fit =  x+((m-x)*((1+(z*h1)**n)**(1/n-1)))
        plt.plot(h1, y_fit, color='RED',linewidth=2.0, label='Fitted curve')

        plt.xticks(fontsize=18)
        plt.yticks(fontsize=18)
        plt.grid(True)
        plt.title(str(o)+"Curve fitting for eqaution of states")
        plt.legend(loc='upper right')
        plt.xlabel('pressure ', fontsize=18)
        plt.ylabel('volumetric water content ', fontsize=18)
    
        l = str(i) +"  "+ "から400個のパラメータ"
        one_list.append(l)
        one_list.append(x)
        one_list.append(m)
        one_list.append(z)
        one_list.append(n)
    
        writer.writerow(one_list)
        plt.show()
        i += 200
    
    else:
        f.close()
        print(str(o)+"終了")

Comments

[t_obara]

コメントのインデントがおかしいです、動作するようにきちんと整形したものをご提示ください。

[mimi_chan]

完了いたしました。申し訳ありませんでした！

[firedfly]

コードを見ているだけではわからない場合は
デバッグ実行する、最小限のコードで実行する、
print文を要所に仕込むなどするといいです。

そのコードをいちばん長く読んでいるのは自分ですから
自分で試行錯誤するのが遠回りに見えていちばんの近道です。

[t_obara]

> １つ目の要素だけが計算
とのことですが、どのようにしてそのように判断したのかもご提示ください。

Answer 1

こんにちは。
iの初期化位置を変えましょう。

for o, g in zip(name, initial):
    i = 4
    while i < 1305: 

Comments

[mimi_chan]

答えを教えてくださるだけでなく、プログラミングへの考え方まで教えて頂き本当にありがとう御座います。精進させて頂きます。

[firedfly]

「試行錯誤のやり方」がわかってくるとプログラミングはぐんと楽しくなります。
この間違いも誰もが通る道ですので、めげずに精進して下さい。