6月から7月にかけて、足立清掃工場、板橋清掃工場、光が丘清掃工場、千歳清掃工場において相次いで、一部の焼却炉を停止する事態が発生しました。これは、ごみの中に水銀を含むものが混入したことが原因です。
最近よく女性誌などで見かける「女子」の自称は、中年女の「え?歳?28よ!(相手の肩を平手打ちしながら)」的自虐表現だと思っていたのだが、先日、真顔で使っている人間を見た。
[偏差値] = (X - [平均値]) / [標準偏差] * 10 + 50
[標準偏差] = √ Σ(([個々の値]i - [平均値])^2 / [総数])
#!/usr/bin/perl use strict; use warnings; my @dat; while (<>) { chomp; push @dat, $_; } # my @dat = map {int(rand(100))} (1..10); my $num = @dat; my $ave; # 平均点 $ave += $_ for @dat; $ave /= @dat; my $sd; # 標準偏差 $sd += ($_ - $ave)**2 for @dat; $sd = sqrt($sd/$num); print "標準偏差\t$sd\n"; print "点数\t偏差値\n"; foreach my $x (@dat) { my $ss = ($x - $ave) / $sd * 10 + 50; # 偏差値 print "$x\t$ss\n"; }
% cat a.txt 12 67 25 37 52 40 13 29 5 93 % ./ss.pl a.txt 標準偏差 25.8884143971777 点数 偏差値 12 40.2272886968473 67 61.4723132689184 25 45.2488399593368 37 49.8841180477887 52 55.6782156583535 40 51.0429375699017 13 40.613561870885 29 46.7939326554875 5 37.5233764785837 93 71.5154157938974
ワンライナー(偏差値):perl -nle '$s+=$_;$v[$.]=$_;END{$a=$s/$.; for(1..$.){$r+=($v[$_]-$a)**2} print sqrt($r/$.)}' a.txt
perl -nle '$s+=$_;$v[$.]=$_;END{$a=$s/$.; for(1..$.){$r+=($v[$_]-$a)**2}$sd=sqrt($r/$.); for(1..$.){print$v[$_]," ",($v[$_]-$a)/$sd*10+50}}' a.txt
$sd += ($_ - $ave/$num)**2 for @dat; $sd = sqrt($sd);