Perl統計::説明とトリム標準偏差

Question

Statistics::Descriptiveモジュールは、調整済み平均を計算するためのFullクラスを提供します。そこからトリミングされた標準偏差を得る方法はありますか？それは既にトリミングされた平均を得るためにデータをソートしているので、難しい部分は完了です。

Answer 1

あなたのデータはすべて正、またはすべて負の数値である場合、あなたはわずか2つの$スタットオブジェクト、無地の値の1と乗値のための1つを維持し、トリミングSTDDEVを生成するtrimmed_mean結果を組み合わせることができます。

my $stat_x = Statistics::Descriptive::Full->new(); 
my $stat_xx = Statistics::Descriptive::Full->new(); 
sub add_x { 
    $stat_x->add_data(@_); 
    $stat_xx->add_data(map {$_*$_} @_); 
} 
sub trimmed_stddev { 
    my $lower = shift; 
    my $upper = @_ ? shift : $lower; 

    my $trm_x = $stat_x->trimmed_mean($lower, $upper); 
    my $trm_xx = $stat_xx->trimmed_mean($lower, $upper); 

    my $nn = int($stat_x->count() * (1.0 - $lower - $upper)); 

    return ($nn <= 1 ? 0.0 : sqrt(($trm_xx - $trm_x * $trm_x) * $nn/($nn -1))); 
} 

# add some data 
add_x(1, 5, 12, 29, 10, 11, 2); 

print trimmed_stddev(0.0); 
print trimmed_stddev(0.05); 
print trimmed_stddev(0.2); 
# output 
# 9.45163125250522 
# 9.58570364061576 
# 4.44222166638871