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