;;;;;;;;;;;;;
*lr_init = *lr'lr_init; #';*lr_sample = *lr'lr_sample; #';*lr_Y = *lr'lr_Y; #';*lr_X = *lr'lr_X; #';*lr_r = *lr'lr_r; #';*lr_cov = *lr'lr_cov; #';*lr_A = *lr'lr_A; #';*lr_B = *lr'lr_B; #';*lr_sigma = *lr'lr_sigma; #';*lr_mean = *lr'lr_mean; #';
package lr;
sub tagify
{
local($tag) = @_;
if (defined($tag))
{
*lr_n = eval "*${tag}_n";
*lr_sx = eval "*${tag}_sx";
*lr_sx2 = eval "*${tag}_sx2";
*lr_sxy = eval "*${tag}_sxy";
*lr_sy = eval "*${tag}_sy";
*lr_sy2 = eval "*${tag}_sy2";
}
}
sub lr_init
{
&tagify($_[$[]) if defined($_[$[]);
$lr_n = 0;
$lr_sx = 0.0;
$lr_sx2 = 0.0;
$lr_sxy = 0.0;
$lr_sy = 0.0;
$lr_sy2 = 0.0;
}
sub lr_sample
{
local($_x, $_y) = @_;
&tagify($_[$[+2]) if defined($_[$[+2]);
$lr_n++;
$lr_sx += $_x;
$lr_sy += $_y;
$lr_sxy += $_x * $_y;
$lr_sx2 += $_x**2;
$lr_sy2 += $_y**2;
}
sub lr_B
{
&tagify($_[$[]) if defined($_[$[]);
return 1 unless ($lr_n * $lr_sx2 - $lr_sx**2);
return ($lr_n * $lr_sxy - $lr_sx * $lr_sy) / ($lr_n * $lr_sx2 - $lr_sx**2);
}
sub lr_A
{
&tagify($_[$[]) if defined($_[$[]);
return ($lr_sy - &lr_B * $lr_sx) / $lr_n;
}
sub lr_Y
{
&tagify($_[$[]) if defined($_[$[]);
return &lr_A + &lr_B * $_[$[];
}
sub lr_X
{
&tagify($_[$[]) if defined($_[$[]);
return ($_[$[] - &lr_A) / &lr_B;
}
sub lr_r
{
&tagify($_[$[]) if defined($_[$[]);
local($s) = ($lr_n * $lr_sx2 - $lr_sx**2) * ($lr_n * $lr_sy2 - $lr_sy**2);
return 1 unless $s;
return ($lr_n * $lr_sxy - $lr_sx * $lr_sy) / sqrt($s);
}
sub lr_cov
{
&tagify($_[$[]) if defined($_[$[]);
return ($lr_sxy - $lr_sx * $lr_sy / $lr_n) / ($lr_n - 1);
}
sub lr_sigma
{
&tagify($_[$[]) if defined($_[$[]);
return 0 if $lr_n <= 1;
return sqrt(($lr_sy2 - ($lr_sy * $lr_sy) / $lr_n) / ($lr_n));
}
sub lr_mean
{
&tagify($_[$[]) if defined($_[$[]);
return 0 if $lr_n <= 0;
return $lr_sy / $lr_n;
}
&lr_init();
1;