my $reg = Statistics::Regression->new("Pain", ["C", "X", "X**2", "X**3", "Y"]);
$reg->include(1.0, [1.0, $x, $x**2, $x**3, $y]);