source: trunk/client/modules/Elezioni/grafici/jpgraph_utils.inc.php

Last change on this file was 284, checked in by roby, 6 years ago
File size: 24.8 KB
RevLine 
[2]1<?php
2/*=======================================================================
[284]3 // File: JPGRAPH_UTILS.INC
4 // Description: Collection of non-essential "nice to have" utilities
5 // Created: 2005-11-20
6 // Ver: $Id: jpgraph_utils.inc.php 1777 2009-08-23 17:34:36Z ljp $
7 //
8 // Copyright (c) Asial Corporation. All rights reserved.
9 //========================================================================
10 */
[2]11
12//===================================================
13// CLASS FuncGenerator
[284]14// Description: Utility class to help generate data for function plots.
[2]15// The class supports both parametric and regular functions.
16//===================================================
17class FuncGenerator {
18 private $iFunc='',$iXFunc='',$iMin,$iMax,$iStepSize;
[284]19
20 function __construct($aFunc,$aXFunc='') {
21 $this->iFunc = $aFunc;
22 $this->iXFunc = $aXFunc;
[2]23 }
[284]24
[2]25 function E($aXMin,$aXMax,$aSteps=50) {
[284]26 $this->iMin = $aXMin;
27 $this->iMax = $aXMax;
28 $this->iStepSize = ($aXMax-$aXMin)/$aSteps;
[2]29
[284]30 if( $this->iXFunc != '' )
31 $t = 'for($i='.$aXMin.'; $i<='.$aXMax.'; $i += '.$this->iStepSize.') {$ya[]='.$this->iFunc.';$xa[]='.$this->iXFunc.';}';
32 elseif( $this->iFunc != '' )
33 $t = 'for($x='.$aXMin.'; $x<='.$aXMax.'; $x += '.$this->iStepSize.') {$ya[]='.$this->iFunc.';$xa[]=$x;} $x='.$aXMax.';$ya[]='.$this->iFunc.';$xa[]=$x;';
34 else
35 JpGraphError::RaiseL(24001);//('FuncGenerator : No function specified. ');
[2]36
[284]37 @eval($t);
[2]38
[284]39 // If there is an error in the function specifcation this is the only
40 // way we can discover that.
41 if( empty($xa) || empty($ya) )
42 JpGraphError::RaiseL(24002);//('FuncGenerator : Syntax error in function specification ');
[2]43
[284]44 return array($xa,$ya);
[2]45 }
46}
47
48
49//=============================================================================
50// CLASS DateScaleUtils
51// Description: Help to create a manual date scale
52//=============================================================================
53define('DSUTILS_MONTH',1); // Major and minor ticks on a monthly basis
54define('DSUTILS_MONTH1',1); // Major and minor ticks on a monthly basis
55define('DSUTILS_MONTH2',2); // Major ticks on a bi-monthly basis
56define('DSUTILS_MONTH3',3); // Major icks on a tri-monthly basis
57define('DSUTILS_MONTH6',4); // Major on a six-monthly basis
58define('DSUTILS_WEEK1',5); // Major ticks on a weekly basis
59define('DSUTILS_WEEK2',6); // Major ticks on a bi-weekly basis
60define('DSUTILS_WEEK4',7); // Major ticks on a quod-weekly basis
61define('DSUTILS_DAY1',8); // Major ticks on a daily basis
62define('DSUTILS_DAY2',9); // Major ticks on a bi-daily basis
63define('DSUTILS_DAY4',10); // Major ticks on a qoud-daily basis
64define('DSUTILS_YEAR1',11); // Major ticks on a yearly basis
65define('DSUTILS_YEAR2',12); // Major ticks on a bi-yearly basis
66define('DSUTILS_YEAR5',13); // Major ticks on a five-yearly basis
67
68
69class DateScaleUtils {
70 public static $iMin=0, $iMax=0;
71
72 private static $starthour,$startmonth, $startday, $startyear;
73 private static $endmonth, $endyear, $endday;
74 private static $tickPositions=array(),$minTickPositions=array();
75 private static $iUseWeeks = true;
76
77 static function UseWeekFormat($aFlg) {
[284]78 self::$iUseWeeks = $aFlg;
[2]79 }
80
81 static function doYearly($aType,$aMinor=false) {
[284]82 $i=0; $j=0;
83 $m = self::$startmonth;
84 $y = self::$startyear;
[2]85
[284]86 if( self::$startday == 1 ) {
87 self::$tickPositions[$i++] = mktime(0,0,0,$m,1,$y);
88 }
89 ++$m;
[2]90
91
[284]92 switch( $aType ) {
93 case DSUTILS_YEAR1:
94 for($y=self::$startyear; $y <= self::$endyear; ++$y ) {
95 if( $aMinor ) {
96 while( $m <= 12 ) {
97 if( !($y == self::$endyear && $m > self::$endmonth) ) {
98 self::$minTickPositions[$j++] = mktime(0,0,0,$m,1,$y);
99 }
100 ++$m;
101 }
102 $m=1;
103 }
104 self::$tickPositions[$i++] = mktime(0,0,0,1,1,$y);
105 }
106 break;
107 case DSUTILS_YEAR2:
108 $y=self::$startyear;
109 while( $y <= self::$endyear ) {
110 self::$tickPositions[$i++] = mktime(0,0,0,1,1,$y);
111 for($k=0; $k < 1; ++$k ) {
112 ++$y;
113 if( $aMinor ) {
114 self::$minTickPositions[$j++] = mktime(0,0,0,1,1,$y);
115 }
116 }
117 ++$y;
118 }
119 break;
120 case DSUTILS_YEAR5:
121 $y=self::$startyear;
122 while( $y <= self::$endyear ) {
123 self::$tickPositions[$i++] = mktime(0,0,0,1,1,$y);
124 for($k=0; $k < 4; ++$k ) {
125 ++$y;
126 if( $aMinor ) {
127 self::$minTickPositions[$j++] = mktime(0,0,0,1,1,$y);
128 }
129 }
130 ++$y;
131 }
132 break;
133 }
[2]134 }
135
136 static function doDaily($aType,$aMinor=false) {
[284]137 $m = self::$startmonth;
138 $y = self::$startyear;
139 $d = self::$startday;
140 $h = self::$starthour;
141 $i=0;$j=0;
[2]142
[284]143 if( $h == 0 ) {
144 self::$tickPositions[$i++] = mktime(0,0,0,$m,$d,$y);
145 }
146 $t = mktime(0,0,0,$m,$d,$y);
[2]147
[284]148 switch($aType) {
149 case DSUTILS_DAY1:
150 while( $t <= self::$iMax ) {
151 $t = strtotime('+1 day',$t);
152 self::$tickPositions[$i++] = $t;
153 if( $aMinor ) {
154 self::$minTickPositions[$j++] = strtotime('+12 hours',$t);
155 }
156 }
157 break;
158 case DSUTILS_DAY2:
159 while( $t <= self::$iMax ) {
160 $t = strtotime('+1 day',$t);
161 if( $aMinor ) {
162 self::$minTickPositions[$j++] = $t;
163 }
164 $t = strtotime('+1 day',$t);
165 self::$tickPositions[$i++] = $t;
166 }
167 break;
168 case DSUTILS_DAY4:
169 while( $t <= self::$iMax ) {
170 for($k=0; $k < 3; ++$k ) {
171 $t = strtotime('+1 day',$t);
172 if( $aMinor ) {
173 self::$minTickPositions[$j++] = $t;
174 }
175 }
176 $t = strtotime('+1 day',$t);
177 self::$tickPositions[$i++] = $t;
178 }
179 break;
180 }
[2]181 }
182
183 static function doWeekly($aType,$aMinor=false) {
[284]184 $hpd = 3600*24;
185 $hpw = 3600*24*7;
186 // Find out week number of min date
187 $thursday = self::$iMin + $hpd * (3 - (date('w', self::$iMin) + 6) % 7);
188 $week = 1 + (date('z', $thursday) - (11 - date('w', mktime(0, 0, 0, 1, 1, date('Y', $thursday)))) % 7) / 7;
189 $daynumber = date('w',self::$iMin);
190 if( $daynumber == 0 ) $daynumber = 7;
191 $m = self::$startmonth;
192 $y = self::$startyear;
193 $d = self::$startday;
194 $i=0;$j=0;
195 // The assumption is that the weeks start on Monday. If the first day
196 // is later in the week then the first week tick has to be on the following
197 // week.
198 if( $daynumber == 1 ) {
199 self::$tickPositions[$i++] = mktime(0,0,0,$m,$d,$y);
200 $t = mktime(0,0,0,$m,$d,$y) + $hpw;
201 }
202 else {
203 $t = mktime(0,0,0,$m,$d,$y) + $hpd*(8-$daynumber);
204 }
[2]205
[284]206 switch($aType) {
207 case DSUTILS_WEEK1:
208 $cnt=0;
209 break;
210 case DSUTILS_WEEK2:
211 $cnt=1;
212 break;
213 case DSUTILS_WEEK4:
214 $cnt=3;
215 break;
216 }
217 while( $t <= self::$iMax ) {
218 self::$tickPositions[$i++] = $t;
219 for($k=0; $k < $cnt; ++$k ) {
220 $t += $hpw;
221 if( $aMinor ) {
222 self::$minTickPositions[$j++] = $t;
223 }
224 }
225 $t += $hpw;
226 }
[2]227 }
228
229 static function doMonthly($aType,$aMinor=false) {
[284]230 $monthcount=0;
231 $m = self::$startmonth;
232 $y = self::$startyear;
233 $i=0; $j=0;
[2]234
[284]235 // Skip the first month label if it is before the startdate
236 if( self::$startday == 1 ) {
237 self::$tickPositions[$i++] = mktime(0,0,0,$m,1,$y);
238 $monthcount=1;
239 }
240 if( $aType == 1 ) {
241 if( self::$startday < 15 ) {
242 self::$minTickPositions[$j++] = mktime(0,0,0,$m,15,$y);
243 }
244 }
245 ++$m;
[2]246
[284]247 // Loop through all the years included in the scale
248 for($y=self::$startyear; $y <= self::$endyear; ++$y ) {
249 // Loop through all the months. There are three cases to consider:
250 // 1. We are in the first year and must start with the startmonth
251 // 2. We are in the end year and we must stop at last month of the scale
252 // 3. A year in between where we run through all the 12 months
253 $stopmonth = $y == self::$endyear ? self::$endmonth : 12;
254 while( $m <= $stopmonth ) {
255 switch( $aType ) {
256 case DSUTILS_MONTH1:
257 // Set minor tick at the middle of the month
258 if( $aMinor ) {
259 if( $m <= $stopmonth ) {
260 if( !($y==self::$endyear && $m==$stopmonth && self::$endday < 15) )
261 self::$minTickPositions[$j++] = mktime(0,0,0,$m,15,$y);
262 }
263 }
264 // Major at month
265 // Get timestamp of first hour of first day in each month
266 self::$tickPositions[$i++] = mktime(0,0,0,$m,1,$y);
[2]267
[284]268 break;
269 case DSUTILS_MONTH2:
270 if( $aMinor ) {
271 // Set minor tick at start of each month
272 self::$minTickPositions[$j++] = mktime(0,0,0,$m,1,$y);
273 }
[2]274
[284]275 // Major at every second month
276 // Get timestamp of first hour of first day in each month
277 if( $monthcount % 2 == 0 ) {
278 self::$tickPositions[$i++] = mktime(0,0,0,$m,1,$y);
279 }
280 break;
281 case DSUTILS_MONTH3:
282 if( $aMinor ) {
283 // Set minor tick at start of each month
284 self::$minTickPositions[$j++] = mktime(0,0,0,$m,1,$y);
285 }
286 // Major at every third month
287 // Get timestamp of first hour of first day in each month
288 if( $monthcount % 3 == 0 ) {
289 self::$tickPositions[$i++] = mktime(0,0,0,$m,1,$y);
290 }
291 break;
292 case DSUTILS_MONTH6:
293 if( $aMinor ) {
294 // Set minor tick at start of each month
295 self::$minTickPositions[$j++] = mktime(0,0,0,$m,1,$y);
296 }
297 // Major at every third month
298 // Get timestamp of first hour of first day in each month
299 if( $monthcount % 6 == 0 ) {
300 self::$tickPositions[$i++] = mktime(0,0,0,$m,1,$y);
301 }
302 break;
303 }
304 ++$m;
305 ++$monthcount;
306 }
307 $m=1;
308 }
[2]309
[284]310 // For the case where all dates are within the same month
311 // we want to make sure we have at least two ticks on the scale
312 // since the scale want work properly otherwise
313 if(self::$startmonth == self::$endmonth && self::$startyear == self::$endyear && $aType==1 ) {
314 self::$tickPositions[$i++] = mktime(0 ,0 ,0, self::$startmonth + 1, 1, self::$startyear);
315 }
[2]316
[284]317 return array(self::$tickPositions,self::$minTickPositions);
[2]318 }
319
320 static function GetTicks($aData,$aType=1,$aMinor=false,$aEndPoints=false) {
[284]321 $n = count($aData);
322 return self::GetTicksFromMinMax($aData[0],$aData[$n-1],$aType,$aMinor,$aEndPoints);
[2]323 }
324
325 static function GetAutoTicks($aMin,$aMax,$aMaxTicks=10,$aMinor=false) {
[284]326 $diff = $aMax - $aMin;
327 $spd = 3600*24;
328 $spw = $spd*7;
329 $spm = $spd*30;
330 $spy = $spd*352;
[2]331
[284]332 if( self::$iUseWeeks )
333 $w = 'W';
334 else
335 $w = 'd M';
[2]336
[284]337 // Decision table for suitable scales
338 // First value: Main decision point
339 // Second value: Array of formatting depending on divisor for wanted max number of ticks. <divisor><formatting><format-string>,..
340 $tt = array(
341 array($spw, array(1,DSUTILS_DAY1,'d M',2,DSUTILS_DAY2,'d M',-1,DSUTILS_DAY4,'d M')),
342 array($spm, array(1,DSUTILS_DAY1,'d M',2,DSUTILS_DAY2,'d M',4,DSUTILS_DAY4,'d M',7,DSUTILS_WEEK1,$w,-1,DSUTILS_WEEK2,$w)),
343 array($spy, array(1,DSUTILS_DAY1,'d M',2,DSUTILS_DAY2,'d M',4,DSUTILS_DAY4,'d M',7,DSUTILS_WEEK1,$w,14,DSUTILS_WEEK2,$w,30,DSUTILS_MONTH1,'M',60,DSUTILS_MONTH2,'M',-1,DSUTILS_MONTH3,'M')),
344 array(-1, array(30,DSUTILS_MONTH1,'M-Y',60,DSUTILS_MONTH2,'M-Y',90,DSUTILS_MONTH3,'M-Y',180,DSUTILS_MONTH6,'M-Y',352,DSUTILS_YEAR1,'Y',704,DSUTILS_YEAR2,'Y',-1,DSUTILS_YEAR5,'Y')));
[2]345
[284]346 $ntt = count($tt);
347 $nd = floor($diff/$spd);
348 for($i=0; $i < $ntt; ++$i ) {
349 if( $diff <= $tt[$i][0] || $i==$ntt-1) {
350 $t = $tt[$i][1];
351 $n = count($t)/3;
352 for( $j=0; $j < $n; ++$j ) {
353 if( $nd/$t[3*$j] <= $aMaxTicks || $j==$n-1) {
354 $type = $t[3*$j+1];
355 $fs = $t[3*$j+2];
356 list($tickPositions,$minTickPositions) = self::GetTicksFromMinMax($aMin,$aMax,$type,$aMinor);
357 return array($fs,$tickPositions,$minTickPositions,$type);
358 }
359 }
360 }
361 }
[2]362 }
363
364 static function GetTicksFromMinMax($aMin,$aMax,$aType,$aMinor=false,$aEndPoints=false) {
[284]365 self::$starthour = date('G',$aMin);
366 self::$startmonth = date('n',$aMin);
367 self::$startday = date('j',$aMin);
368 self::$startyear = date('Y',$aMin);
369 self::$endmonth = date('n',$aMax);
370 self::$endyear = date('Y',$aMax);
371 self::$endday = date('j',$aMax);
372 self::$iMin = $aMin;
373 self::$iMax = $aMax;
[2]374
[284]375 if( $aType <= DSUTILS_MONTH6 ) {
376 self::doMonthly($aType,$aMinor);
377 }
378 elseif( $aType <= DSUTILS_WEEK4 ) {
379 self::doWeekly($aType,$aMinor);
380 }
381 elseif( $aType <= DSUTILS_DAY4 ) {
382 self::doDaily($aType,$aMinor);
383 }
384 elseif( $aType <= DSUTILS_YEAR5 ) {
385 self::doYearly($aType,$aMinor);
386 }
387 else {
388 JpGraphError::RaiseL(24003);
389 }
390 // put a label at the very left data pos
391 if( $aEndPoints ) {
392 $tickPositions[$i++] = $aData[0];
393 }
[2]394
[284]395 // put a label at the very right data pos
396 if( $aEndPoints ) {
397 $tickPositions[$i] = $aData[$n-1];
398 }
[2]399
[284]400 return array(self::$tickPositions,self::$minTickPositions);
[2]401 }
402}
403
404//=============================================================================
405// Class ReadFileData
406//=============================================================================
407Class ReadFileData {
408 //----------------------------------------------------------------------------
409 // Desciption:
[284]410 // Read numeric data from a file.
411 // Each value should be separated by either a new line or by a specified
[2]412 // separator character (default is ',').
[284]413 // Before returning the data each value is converted to a proper float
414 // value. The routine is robust in the sense that non numeric data in the
[2]415 // file will be discarded.
416 //
[284]417 // Returns:
[2]418 // The number of data values read on success, FALSE on failure
419 //----------------------------------------------------------------------------
420 static function FromCSV($aFile,&$aData,$aSepChar=',',$aMaxLineLength=1024) {
[284]421 $rh = @fopen($aFile,'r');
422 if( $rh === false ) {
423 return false;
424 }
425 $tmp = array();
426 $lineofdata = fgetcsv($rh, 1000, ',');
427 while ( $lineofdata !== FALSE) {
428 $tmp = array_merge($tmp,$lineofdata);
429 $lineofdata = fgetcsv($rh, $aMaxLineLength, $aSepChar);
430 }
431 fclose($rh);
[2]432
[284]433 // Now make sure that all data is numeric. By default
434 // all data is read as strings
435 $n = count($tmp);
436 $aData = array();
437 $cnt=0;
438 for($i=0; $i < $n; ++$i) {
439 if( $tmp[$i] !== "" ) {
440 $aData[$cnt++] = floatval($tmp[$i]);
441 }
442 }
443 return $cnt;
[2]444 }
[284]445
446 //----------------------------------------------------------------------------
447 // Desciption:
448 // Read numeric data from a file.
449 // Each value should be separated by either a new line or by a specified
450 // separator character (default is ',').
451 // Before returning the data each value is converted to a proper float
452 // value. The routine is robust in the sense that non numeric data in the
453 // file will be discarded.
454 //
455 // Options:
456 // 'separator' => ',',
457 // 'enclosure' => '"',
458 // 'readlength' => 1024,
459 // 'ignore_first' => false,
460 // 'first_as_key' => false
461 // 'escape' => '\', # PHP >= 5.3 only
462 //
463 // Returns:
464 // The number of lines read on success, FALSE on failure
465 //----------------------------------------------------------------------------
466 static function FromCSV2($aFile, &$aData, $aOptions = array()) {
467 $aDefaults = array(
468 'separator' => ',',
469 'enclosure' => chr(34),
470 'escape' => chr(92),
471 'readlength' => 1024,
472 'ignore_first' => false,
473 'first_as_key' => false
474 );
475
476 $aOptions = array_merge(
477 $aDefaults, is_array($aOptions) ? $aOptions : array());
478
479 if( $aOptions['first_as_key'] ) {
480 $aOptions['ignore_first'] = true;
481 }
482
483 $rh = @fopen($aFile, 'r');
484
485 if( $rh === false ) {
486 return false;
487 }
488
489 $aData = array();
490 $aLine = fgetcsv($rh,
491 $aOptions['readlength'],
492 $aOptions['separator'],
493 $aOptions['enclosure']
494 /*, $aOptions['escape'] # PHP >= 5.3 only */
495 );
496
497 // Use numeric array keys for the columns by default
498 // If specified use first lines values as assoc keys instead
499 $keys = array_keys($aLine);
500 if( $aOptions['first_as_key'] ) {
501 $keys = array_values($aLine);
502 }
503
504 $num_lines = 0;
505 $num_cols = count($aLine);
506
507 while ($aLine !== false) {
508 if( is_array($aLine) && count($aLine) != $num_cols ) {
509 JpGraphError::RaiseL(24004);
510 // 'ReadCSV2: Column count mismatch in %s line %d'
511 }
512
513 // fgetcsv returns NULL for empty lines
514 if( !is_null($aLine) ) {
515 $num_lines++;
516
517 if( !($aOptions['ignore_first'] && $num_lines == 1) && is_numeric($aLine[0]) ) {
518 for( $i = 0; $i < $num_cols; $i++ ) {
519 $aData[ $keys[$i] ][] = floatval($aLine[$i]);
520 }
521 }
522 }
523
524 $aLine = fgetcsv($rh,
525 $aOptions['readlength'],
526 $aOptions['separator'],
527 $aOptions['enclosure']
528 /*, $aOptions['escape'] # PHP >= 5.3 only*/
529 );
530 }
531
532 fclose($rh);
533
534 if( $aOptions['ignore_first'] ) {
535 $num_lines--;
536 }
537
538 return $num_lines;
539 }
540
541 // Read data from two columns in a plain text file
542 static function From2Col($aFile, $aCol1, $aCol2, $aSepChar=' ') {
543 $lines = @file($aFile,FILE_IGNORE_NEW_LINES|FILE_SKIP_EMPTY_LINES);
544 if( $lines === false ) {
545 return false;
546 }
547 $s = '/[\s]+/';
548 if( $aSepChar == ',' ) {
549 $s = '/[\s]*,[\s]*/';
550 }
551 elseif( $aSepChar == ';' ) {
552 $s = '/[\s]*;[\s]*/';
553 }
554 foreach( $lines as $line => $datarow ) {
555 $split = preg_split($s,$datarow);
556 $aCol1[] = floatval(trim($split[0]));
557 $aCol2[] = floatval(trim($split[1]));
558 }
559
560 return count($lines);
561 }
562
563 // Read data from one columns in a plain text file
564 static function From1Col($aFile, $aCol1) {
565 $lines = @file($aFile,FILE_IGNORE_NEW_LINES|FILE_SKIP_EMPTY_LINES);
566 if( $lines === false ) {
567 return false;
568 }
569 foreach( $lines as $line => $datarow ) {
570 $aCol1[] = floatval(trim($datarow));
571 }
572
573 return count($lines);
574 }
575
576 static function FromMatrix($aFile,$aSepChar=' ') {
577 $lines = @file($aFile,FILE_IGNORE_NEW_LINES|FILE_SKIP_EMPTY_LINES);
578 if( $lines === false ) {
579 return false;
580 }
581 $mat = array();
582 $reg = '/'.$aSepChar.'/';
583 foreach( $lines as $line => $datarow ) {
584 $row = preg_split($reg,trim($datarow));
585 foreach ($row as $key => $cell ) {
586 $row[$key] = floatval(trim($cell));
587 }
588 $mat[] = $row;
589 }
590 return $mat;
591 }
592
593
[2]594}
595
[284]596define('__LR_EPSILON', 1.0e-8);
597//=============================================================================
598// Class LinearRegression
599//=============================================================================
600class LinearRegression {
601 private $ix=array(),$iy=array();
602 private $ib=0, $ia=0;
603 private $icalculated=false;
604 public $iDet=0, $iCorr=0, $iStdErr=0;
605
606 public function __construct($aDataX,$aDataY) {
607 if( count($aDataX) !== count($aDataY) ) {
608 JpGraph::Raise('LinearRegression: X and Y data array must be of equal length.');
609 }
610 $this->ix = $aDataX;
611 $this->iy = $aDataY;
612 }
613
614 public function Calc() {
615
616 $this->icalculated = true;
617
618 $n = count($this->ix);
619 $sx2 = 0 ;
620 $sy2 = 0 ;
621 $sxy = 0 ;
622 $sx = 0 ;
623 $sy = 0 ;
624
625 for( $i=0; $i < $n; ++$i ) {
626 $sx2 += $this->ix[$i] * $this->ix[$i];
627 $sy2 += $this->iy[$i] * $this->iy[$i];
628 $sxy += $this->ix[$i] * $this->iy[$i];
629 $sx += $this->ix[$i];
630 $sy += $this->iy[$i];
631 }
632
633 if( $n*$sx2 - $sx*$sx > __LR_EPSILON ) {
634 $this->ib = ($n*$sxy - $sx*$sy) / ( $n*$sx2 - $sx*$sx );
635 $this->ia = ( $sy - $this->ib*$sx ) / $n;
636
637 $sx = $this->ib * ( $sxy - $sx*$sy/$n );
638 $sy2 = $sy2 - $sy*$sy/$n;
639 $sy = $sy2 - $sx;
640
641 $this->iDet = $sx / $sy2;
642 $this->iCorr = sqrt($this->iDet);
643 if( $n > 2 ) {
644 $this->iStdErr = sqrt( $sy / ($n-2) );
645 }
646 else {
647 $this->iStdErr = NAN ;
648 }
649 }
650 else {
651 $this->ib = 0;
652 $this->ia = 0;
653 }
654
655 }
656
657 public function GetAB() {
658 if( $this->icalculated == false )
659 $this->Calc();
660 return array($this->ia, $this->ib);
661 }
662
663 public function GetStat() {
664 if( $this->icalculated == false )
665 $this->Calc();
666 return array($this->iStdErr, $this->iCorr, $this->iDet);
667 }
668
669 public function GetY($aMinX, $aMaxX, $aStep=1) {
670 if( $this->icalculated == false )
671 $this->Calc();
672
673 $yy = array();
674 $i = 0;
675 for( $x=$aMinX; $x <= $aMaxX; $x += $aStep ) {
676 $xx[$i ] = $x;
677 $yy[$i++] = $this->ia + $this->ib * $x;
678 }
679
680 return array($xx,$yy);
681 }
682
683}
684
685?>
Note: See TracBrowser for help on using the repository browser.