varmodel_band_k — SPC variation op

• 数据种类:无 → table(仅由参数决定的算子 —— 不接受图像或数据输入)

• 调用: import fullseye as fs; fs.ledger.varmodel_band_k(n=20, far=0.001, spread='std', df_prior=None)(要直接调用实现,import varmodel; varmodel.varmodel_band_k(n=20, far=0.001, spread='std', df_prior=None);从台账取用则 opsspc.get("varmodel_band_k"))

用法

对应所声明的每像素双侧误检率的容许带半宽 k(以离散度为单位,table)。

> 以下的详细说明为原文 —— 摘要与标题已翻译。

A new good pixel x is flagged when |x - center| > k * spread. Under per-pixel independent normal good samples:

• std: (x - mean)/(s * sqrt(1+1/n)) ~ t_{n-1} so k = t_{n-1}^-1(1 - far/2) * sqrt(1+1/n).

• pooled: one variance pooled over many pixels (degrees of freedom df_prior; None = infinite) → normal quantile.

• moderated: empirical-Bayes shrunk variance with df_prior = d0 prior degrees of freedom (Smyth 2004) → t_{n-1+d0}.

• mad: median / MAD centre and spread → no closed form, but the ratio is a pivot, so k is a Monte-Carlo quantile that depends on (n, far) only (fixed seed; mc_far_rel_se is the relative standard error of the achieved rate). Needs n >= 5 and far >= 1e-3.

The threshold comes from the distribution, never from test defects. Returns k, df (0 for mad: no t quantile is involved), far, n, spread, source, mc_far_rel_se.

Raises ValueError for n < 2 (mad: n < 5), far outside [1e-9, 0.5), an unknown spread or a non-positive df_prior.

参考(示例数据、文献)

• 示例数据目录(下载 URL / 许可证) —— 2-D 用 skimage.data(BSD/公有领域)加合成图,3-D 给出真实数据源(Stanford/PDS 等)的下载 URL。

• 算子来历与参考文献 —— 该算子族所依据的研究/方法出处。

• 算法的正典(作者、年份)与用途见上面的族使用指南。

可运行的示例(实际调用该算子并已验证的样例)

• (暂无)

类型可衔接的下一个算子(可接受 table 作为输入)

msa_anova_table · msa_gauge_rr · msa_bias_linearity · msa_attribute_agreement · gum_standard_uncertainty · gum_propagate · gum_expanded · gum_monte_carlo

同类别(variation)

varmodel_alignment · varmodel_fit · varmodel_calibrate · varmodel_score


*Provenance: varmodel.py — SPC 算子登记表。本条目由 tools/opdocs.py md 自动生成(请勿手工编辑)。*

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.