y=fdtrc(dfn,dfd,x) returns the complemented F distribution function.
scipy.special.fdtr
scipy.special.fdtri