IDL中更强大的统计矩

Question

我正在处理IDL中的噪声数据，所以我一直在使用STDDEV和robust_sigma 。有关于健壮的skewness和kurtosis的论文，例如[1]和[2]，但有标准偏差的实现吗？ （用IDL或C？）

Answer 1

http://idlastro.gsfc.nasa.gov/ftp/pro/robust/robust_sigma.pro州的文件：

; OPTIONAL OUPTUT KEYWORD:
;       GOODVEC = Vector of non-trimmed indices of the input vector

因此，我们使用额外的参数调用robust_sigma，以跟踪＆＃34;良好的索引＆＃34;在数据中，用于计算robust_sigma的数据，与计算中忽略的数据相反。

  good_indices = lonarr(width)
  robo_2 = robust_sigma(data[*], GOODVEC=good_indices)

然后使用（仅）那些好的指数计算其他时刻。

  robo_3 = skewness(data[good_indices])
  robo_4 = kurtosis(data[good_indices])

无需特殊实施。