Unbias a random generator: Difference between revisions

→‎{{header|Perl}}: fixed stddev for binomial distr
(→‎{{header|GAP}}: seems it was not obvious)
(→‎{{header|Perl}}: fixed stddev for binomial distr)
Line 308:
return int(rand($n) / ($n - 1));
}
 
for my $n (3 .. 6) {
print "Bias $n: ";
Line 318:
}
print "@raw, ";
printf("%3g+-%.3g%%\tfixed: ", $raw[0]/100, 100/sqrt($raw[0]));
100 * sqrt($raw[0] * $raw[1]) / ($raw[0] + $raw[1])**1.5);
print "@fixed, ";
printf("%3g+-%.3g%%\n", 100*$fixed[0]/($fixed[0] + $fixed[1]), 100/sqrt($fixed[0]));
100 * sqrt($fixed[0] * $fixed[1]) / ($fixed[0] + $fixed[1])**1.5);
}</lang>
 
}</lang>Output:<prelang>Bias 3: 66156684 33853316, 66.1584+-10.23471% fixed: 21462188 22832228, 4849.45345471+-20.16752%
Bias 4: 75257537 24752463, 75.2537+-10.15431% fixed: 18801924 18371845, 5051.5784048+-20.31814%
Bias 5: 79897993 20112007, 79.8993+-10.12401% fixed: 15731564 15781597, 49.9207478+-20.52889%
Bias 6: 82928309 17081691, 8283.9209+-10.1375% fixed: 13781403 14041410, 49.53278756+-20.69943%</prelang>
 
=={{header|Perl 6}}==
Anonymous user