いくつかの t 値と自由度があり、それらから p 値を見つけたいと考えています (両側です)。現実の世界では、統計の教科書の後ろにある t 検定表を使用します。ただし、Python で stdtr または stats.t.sf 関数を使用しています。どちらも自由度が小さい場合は問題なく動作しますが、自由度が大きい場合は nan を指定してください。
pf = 2*stdtr(dof, -np.abs(tf))
pval = 2*stats.t.sf(np.abs(tf), dof-1)
tf = -11.374250, dof=-2176568.362223 は pf と pval= nan を与えます。
ここで内部で何が起こっているのかを理解するのを手伝ってくれませんか。また、Python のこれらの内部関数のコードをどのように読み取ることができますか。