Pythonを使って方程式のゼロを見つける

Question

私はこの方程式でnを見つけるコードを書こうとしています。ユーザー定義の変数として残りの部分と

。

from scipy.optimize import fsolve 
from scipy.stats import t 

def f(alpha, beta, sigma, delta, eps): 
    n = ((t.ppf(1-alpha,2*n-2) + t.ppf((1-beta)/2,2*n-2))**2*sigma**2)/(2*  (delta-abs(eps))**2) 

私もそれを出力パラメータと結果のテーブル（例えば、入力アルファ1、アルファ2、ベータ1、ベータ2などがあり、その後のパラメータの異なるシナリオを設定し、できるようにして取得したいのですが

out [alpha1、beta1、...、n]、[alpha1、beta2、...、n]）。誰かが世俗的に私を正しい方向に向けることができれば、それをする最善の方法は何だろうか、私は確信していません。

Answer 1

統計検定方程式を満たす観測数（n）を求めようとしています。その場合、nは自然数（0,1,2など）であり、容易に反復可能です。

あなたは、ソルバーあなたが式の「結果」として、反復可能な結果と式のようにn個持っている自分を、設定することができ

：

for n in range(0, 1000): 
    result = your_function(n, other_parameters) 

その後、あなたは、単に式があるかどうかを確認する必要があります

if n >= result: 
     print "result:", n 
     break # This will exit the loop 

異なるユーザー指定のパラメータをテストすると、アルファ、ベータなどの異なる値を繰り返す別のループを設定できます。