This code is modified from statsmodels.stats.anova_lm

Compute the D'Agostino-Pearson K2 test statistic continuum for a set of model residuals.

Label clusters in a binary field *b*. This function yields the same output as **scipy.ndimage.measurements.label** but is much faster for 1D field

Inverse survival function. (see also the survival function: **RFTCalculator.sf**) :Parameters: *alpha* -- upper tail probability (float;

Survival function. (Equivalent to **RFTCalculator.p.upcrossing**) Probability that 1D Gaussian fields with a smoothness *FWHM* would produce

RFT inverse survival function. (see also the survival function: **rft1d.DISTFLAG.sf**) :Parameters: *alpha* -- upper tail probability (flo

Upcrossing extents (units: nodes). :Parameters: *y* --- a 1D field *u* --- threshold height *interp* --- interpolate to thres

RFT cluster-level inference. Probability that 1D Gaussian fields with a smoothness of *FWHM* would produce an upcrossing of extent *k* when

RFT set-level inference. Probability that 1D Gaussian fields with a smoothness of *FWHM* would produce at least *c* upcrossings with a minim

Random Field Theory probabilities and expectations using unified Euler Characteristic (EC) theory. This code is based on "spm_P_RF.m" and "spm_P.m"