sub expit {
  my $x = shift;
  1.0/(1.0 + exp(-$x));
}