00000000
00000018
00000030
00000048
00000060
00000078
00000090
000000A8
000000C0
000000D8
000000F0
00000108
00000120
00000138
00000150
00000168
00000180
00000198
000001B0
000001C8
000001E0
000001F8
00000210
00000228
00000240
00000258
00000270
00000288
000002A0
000002B8
000002D0
000002E8
00000300
00000318
00000330
00000348
00000360
00000378
00000390
000003A8
000003C0
000003D8
000003F0
00000408
00000420
00000438
00000450
00000468
00000480
00000498
000004B0
000004C8
000004E0
000004F8
00000510
00000528
00000540
00000558
00000570
00000588
000005A0
000005B8
000005D0
000005E8
00000600
00000618
00000630
00000648
00000660
00000678
00000690
000006A8
000006C0
000006D8
000006F0
00000708
00000720
00000738
00000750
00000768
00000780
00000798
000007B0
000007C8
000007E0
000007F8
00000810
00000828
00000840
00000858
00000870
00000888
000008A0
000008B8
000008D0
000008E8
00000900
00000918
00000930
00000948
00000960
00000978
00000990
000009A8
000009C0
000009D8
000009F0
00000A08
00000A20
00000A38
00000A50
00000A68
00000A80
00000A98
00000AB0
00000AC8
00000AE0
00000AF8
00000B10
00000B28
00000B40
00000B58
00000B70
00000B88
00000BA0
00000BB8
00000BD0
00000BE8
00000C00
00000C18
00000C30
00000C48
00000C60
00000C78
00000C90
00000CA8
00000CC0
00000CD8
00000CF0
00000D08
| 3C 3F 70 68 70 20 0A 2F 2A 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D
3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D
3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D
3D 3D 3D 3D 3D 3D 3D 3D 0A 2F 2F 20 46 69 6C 65 3A 09 4A 50 47 52 41 50
48 5F 52 45 47 53 54 41 54 2E 50 48 50 0A 2F 2F 20 44 65 73 63 72 69 70
74 69 6F 6E 3A 20 52 65 67 72 65 73 73 69 6F 6E 20 61 6E 64 20 73 74 61
74 69 73 74 69 63 61 6C 20 61 6E 61 6C 79 73 69 73 20 68 65 6C 70 65 72
20 63 6C 61 73 73 65 73 0A 2F 2F 20 43 72 65 61 74 65 64 3A 20 09 32 30
30 32 2D 31 32 2D 30 31 0A 2F 2F 20 41 75 74 68 6F 72 3A 09 4A 6F 68 61
6E 20 50 65 72 73 73 6F 6E 20 28 6A 6F 68 61 6E 70 40 61 64 69 74 75 73
2E 6E 75 29 0A 2F 2F 20 56 65 72 3A 09 09 24 49 64 3A 20 6A 70 67 72 61
70 68 5F 72 65 67 73 74 61 74 2E 70 68 70 20 32 31 20 32 30 30 35 2D 30
35 2D 33 30 20 32 30 3A 33 35 3A 33 34 5A 20 6C 6A 70 20 24 0A 2F 2F 0A
2F 2F 20 43 6F 70 79 72 69 67 68 74 20 28 63 29 20 41 64 69 74 75 73 20
43 6F 6E 73 75 6C 74 69 6E 67 2E 20 41 6C 6C 20 72 69 67 68 74 73 20 72
65 73 65 72 76 65 64 2E 0A 2F 2F 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D
3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D
3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D
3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 3D 0A 2A 2F 0A 0A 2F 2F 2D 2D 2D 2D 2D 2D
2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D
2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D
2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 0A 2F 2F 20 43 4C
41 53 53 20 53 70 6C 69 6E 65 0A 2F 2F 20 43 72 65 61 74 65 20 61 20 6E
65 77 20 64 61 74 61 20 61 72 72 61 79 20 66 72 6F 6D 20 61 6E 20 65 78
69 73 74 69 6E 67 20 64 61 74 61 20 61 72 72 61 79 20 62 75 74 20 77 69
74 68 20 6D 6F 72 65 20 70 6F 69 6E 74 73 2E 0A 2F 2F 20 54 68 65 20 6E
65 77 20 70 6F 69 6E 74 73 20 61 72 65 20 69 6E 74 65 72 70 6F 6C 61 74
65 64 20 75 73 69 6E 67 20 61 20 63 75 62 69 63 20 73 70 6C 69 6E 65 20
61 6C 67 6F 72 69 74 68 6D 0A 2F 2F 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D
2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D
2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D
2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 2D 0A 63 6C 61 73 73 20 53 70 6C 69 6E
65 20 7B 0A 20 20 20 20 2F 2F 20 33 3A 72 64 20 64 65 67 72 65 65 20 70
6F 6C 79 6E 6F 6D 20 61 70 70 72 6F 78 69 6D 61 74 69 6F 6E 0A 0A 20 20
20 20 76 61 72 20 24 78 64 61 74 61 2C 24 79 64 61 74 61 3B 20 20 20 2F
2F 20 44 61 74 61 20 76 65 63 74 6F 72 73 0A 20 20 20 20 76 61 72 20 24
79 32 3B 09 09 20 2F 2F 20 32 3A 6E 64 20 64 65 72 69 76 61 74 65 20 6F
66 20 79 64 61 74 61 09 0A 20 20 20 20 76 61 72 20 24 6E 3D 30 3B 0A 0A
20 20 20 20 66 75 6E 63 74 69 6F 6E 20 53 70 6C 69 6E 65 28 24 78 64 61
74 61 2C 24 79 64 61 74 61 29 20 7B 0A 09 24 74 68 69 73 2D 3E 79 32 20
3D 20 61 72 72 61 79 28 29 3B 0A 09 24 74 68 69 73 2D 3E 78 64 61 74 61
20 3D 20 24 78 64 61 74 61 3B 0A 09 24 74 68 69 73 2D 3E 79 64 61 74 61
20 3D 20 24 79 64 61 74 61 3B 0A 0A 09 24 6E 20 3D 20 63 6F 75 6E 74 28
24 79 64 61 74 61 29 3B 0A 09 24 74 68 69 73 2D 3E 6E 20 3D 20 24 6E 3B
0A 0A 09 2F 2F 20 4E 61 74 75 72 61 6C 20 73 70 6C 69 6E 65 20 32 3A 64
65 72 69 76 61 74 65 20 3D 3D 20 30 20 61 74 20 65 6E 64 70 6F 69 6E 74
73 0A 09 24 74 68 69 73 2D 3E 79 32 5B 30 5D 20 20 20 20 3D 20 30 2E 30
3B 0A 09 24 74 68 69 73 2D 3E 79 32 5B 24 6E 2D 31 5D 20 3D 20 30 2E 30
3B 0A 09 24 64 65 6C 74 61 5B 30 5D 20 3D 20 30 2E 30 3B 0A 0A 09 2F 2F
20 43 61 6C 63 75 6C 61 74 65 20 32 3A 6E 64 20 64 65 72 69 76 61 74 65
0A 09 66 6F 72 28 24 69 3D 31 3B 20 24 69 20 3C 20 24 6E 2D 31 3B 20 2B
2B 24 69 29 20 7B 0A 09 20 20 20 20 24 64 20 3D 20 28 24 78 64 61 74 61
5B 24 69 2B 31 5D 2D 24 78 64 61 74 61 5B 24 69 2D 31 5D 29 3B 0A 09 20
20 20 20 69 66 28 20 24 64 20 3D 3D 20 30 20 20 29 20 7B 0A 09 09 4A 70
47 72 61 70 68 45 72 72 6F 72 3A 3A 52 61 69 73 65 28 27 49 6E 76 61 6C
69 64 20 69 6E 70 75 74 20 64 61 74 61 20 66 6F 72 20 73 70 6C 69 6E 65
2E 20 54 77 6F 20 6F 72 20 6D 6F 72 65 20 63 6F 6E 73 65 63 75 74 69 76
65 20 69 6E 70 75 74 20 58 2D 76 61 6C 75 65 73 20 61 72 65 20 65 71 75
61 6C 2E 20 45 61 63 68 20 69 6E 70 75 74 20 58 2D 76 61 6C 75 65 20 6D
75 73 74 20 64 69 66 66 65 72 20 73 69 6E 63 65 20 66 72 6F 6D 20 61 20
6D 61 74 68 65 6D 61 74 69 63 61 6C 20 70 6F 69 6E 74 20 6F 66 20 76 69
65 77 20 69 74 20 6D 75 73 74 20 62 65 20 61 20 6F 6E 65 2D 74 6F 2D 6F
6E 65 20 6D 61 70 70 69 6E 67 2C 20 69 2E 65 2E 20 65 61 63 68 20 58 2D
76 61 6C 75 65 20 6D 75 73 74 20 63 6F 72 72 65 73 70 6F 6E 64 20 74 6F
20 65 78 61 63 74 6C 79 20 6F 6E 65 20 59 2D 76 61 6C 75 65 2E 27 29 3B
0A 09 20 20 20 20 7D 0A 09 20 20 20 20 24 73 20 3D 20 28 24 78 64 61 74
61 5B 24 69 5D 2D 24 78 64 61 74 61 5B 24 69 2D 31 5D 29 2F 24 64 3B 0A
09 20 20 20 20 24 70 20 3D 20 24 73 2A 24 74 68 69 73 2D 3E 79 32 5B 24
69 2D 31 5D 2B 32 2E 30 3B 0A 09 20 20 20 20 24 74 68 69 73 2D 3E 79 32
5B 24 69 5D 20 3D 20 28 24 73 2D 31 2E 30 29 2F 24 70 3B 0A 09 20 20 20
20 24 64 65 6C 74 61 5B 24 69 5D 20 3D 20 28 24 79 64 61 74 61 5B 24 69
2B 31 5D 2D 24 79 64 61 74 61 5B 24 69 5D 29 2F 28 24 78 64 61 74 61 5B
24 69 2B 31 5D 2D 24 78 64 61 74 61 5B 24 69 5D 29 20 2D 20 0A 09 09 20
20 20 20 20 20 20 20 20 28 24 79 64 61 74 61 5B 24 69 5D 2D 24 79 64 61
74 61 5B 24 69 2D 31 5D 29 2F 28 24 78 64 61 74 61 5B 24 69 5D 2D 24 78
64 61 74 61 5B 24 69 2D 31 5D 29 3B 0A 09 20 20 20 20 24 64 65 6C 74 61
5B 24 69 5D 20 3D 20 28 36 2E 30 2A 24 64 65 6C 74 61 5B 24 69 5D 2F 28
24 78 64 61 74 61 5B 24 69 2B 31 5D 2D 24 78 64 61 74 61 5B 24 69 2D 31
5D 29 2D 24 73 2A 24 64 65 6C 74 61 5B 24 69 2D 31 5D 29 2F 24 70 3B 0A
09 7D 0A 0A 09 2F 2F 20 42 61 63 6B 77 61 72 64 20 73 75 62 73 74 69 74
75 74 69 6F 6E 0A 09 66 6F 72 28 20 24 6A 3D 24 6E 2D 32 3B 20 24 6A 20
3E 3D 20 30 3B 20 2D 2D 24 6A 20 29 20 7B 0A 09 20 20 20 20 24 74 68 69
73 2D 3E 79 32 5B 24 6A 5D 20 3D 20 24 74 68 69 73 2D 3E 79 32 5B 24 6A
5D 2A 24 74 68 69 73 2D 3E 79 32 5B 24 6A 2B 31 5D 20 2B 20 24 64 65 6C
74 61 5B 24 6A 5D 3B 0A 09 7D 0A 20 20 20 20 7D 0A 0A 20 20 20 20 2F 2F
20 52 65 74 75 72 6E 20 74 68 65 20 74 77 6F 20 6E 65 77 20 64 61 74 61
20 76 65 63 74 6F 72 73 0A 20 20 20 20 66 75 6E 63 74 69 6F 6E 20 47 65
74 28 24 6E 75 6D 3D 35 30 29 20 7B 0A 09 24 6E 20 3D 20 24 74 68 69 73
2D 3E 6E 20 3B 0A 09 24 73 74 65 70 20 3D 20 28 24 74 68 69 73 2D 3E 78
64 61 74 61 5B 24 6E 2D 31 5D 2D 24 74 68 69 73 2D 3E 78 64 61 74 61 5B
30 5D 29 20 2F 20 28 24 6E 75 6D 2D 31 29 3B 0A 09 24 78 6E 65 77 3D 61
72 72 61 79 28 29 3B 0A 09 24 79 6E 65 77 3D 61 72 72 61 79 28 29 3B 0A
09 24 78 6E 65 77 5B 30 5D 20 3D 20 24 74 68 69 73 2D 3E 78 64 61 74 61
5B 30 5D 3B 0A 09 24 79 6E 65 77 5B 30 5D 20 3D 20 24 74 68 69 73 2D 3E
79 64 61 74 61 5B 30 5D 3B 0A 09 66 6F 72 28 20 24 6A 3D 31 3B 20 24 6A
20 3C 20 24 6E 75 6D 3B 20 2B 2B 24 6A 20 29 20 7B 0A 09 20 20 20 20 24
78 6E 65 77 5B 24 6A 5D 20 3D 20 24 78 6E 65 77 5B 30 5D 2B 24 6A 2A 24
73 74 65 70 3B 0A 09 20 20 20 20 24 79 6E 65 77 5B 24 6A 5D 20 3D 20 24
74 68 69 73 2D 3E 49 6E 74 65 72 70 6F 6C 61 74 65 28 24 78 6E 65 77 5B
24 6A 5D 29 3B 0A 09 7D 0A 09 72 65 74 75 72 6E 20 61 72 72 61 79 28 24
78 6E 65 77 2C 24 79 6E 65 77 29 3B 0A 20 20 20 20 7D 0A 0A 20 20 20 20
2F 2F 20 52 65 74 75 72 6E 20 61 20 73 69 6E 67 6C 65 20 69 6E 74 65 72
70 6F 6C 61 74 65 64 20 59 2D 76 61 6C 75 65 20 66 72 6F 6D 20 61 6E 20
78 20 76 61 6C 75 65 0A 20 20 20 20 66 75 6E 63 74 69 6F 6E 20 49 6E 74
65 72 70 6F 6C 61 74 65 28 24 78 70 6F 69 6E 74 29 20 7B 0A 0A 09 24 6D
61 78 20 3D 20 24 74 68 69 73 2D 3E 6E 2D 31 3B 0A 09 24 6D 69 6E 20 3D
20 30 3B 0A 0A 09 2F 2F 20 42 69 6E 61 72 79 20 73 65 61 72 63 68 20 74
6F 20 66 69 6E 64 20 69 6E 74 65 72 76 61 6C 0A 09 77 68 69 6C 65 28 20
24 6D 61 78 2D 24 6D 69 6E 20 3E 20 31 20 29 20 7B 0A 09 20 20 20 20 24
6B 20 3D 20 28 24 6D 61 78 2B 24 6D 69 6E 29 20 2F 20 32 3B 0A 09 20 20
20 20 69 66 28 20 24 74 68 69 73 2D 3E 78 64 61 74 61 5B 24 6B 5D 20 3E
20 24 78 70 6F 69 6E 74 20 29 20 0A 09 09 24 6D 61 78 3D 24 6B 3B 0A 09
20 20 20 20 65 6C 73 65 20 0A 09 09 24 6D 69 6E 3D 24 6B 3B 0A 09 7D 09
0A 0A 09 2F 2F 20 45 61 63 68 20 69 6E 74 65 72 76 61 6C 20 69 73 20 69
6E 74 65 72 70 6F 6C 61 74 65 64 20 62 79 20 61 20 33 3A 64 65 67 72 65
65 20 70 6F 6C 79 6E 6F 6D 20 66 75 6E 63 74 69 6F 6E 0A 09 24 68 20 3D
20 24 74 68 69 73 2D 3E 78 64 61 74 61 5B 24 6D 61 78 5D 2D 24 74 68 69
73 2D 3E 78 64 61 74 61 5B 24 6D 69 6E 5D 3B 0A 0A 09 69 66 28 20 24 68
20 3D 3D 20 30 20 20 29 20 7B 0A 09 20 20 20 20 4A 70 47 72 61 70 68 45
72 72 6F 72 3A 3A 52 61 69 73 65 28 27 49 6E 76 61 6C 69 64 20 69 6E 70
75 74 20 64 61 74 61 20 66 6F 72 20 73 70 6C 69 6E 65 2E 20 54 77 6F 20
6F 72 20 6D 6F 72 65 20 63 6F 6E 73 65 63 75 74 69 76 65 20 69 6E 70 75
74 20 58 2D 76 61 6C 75 65 73 20 61 72 65 20 65 71 75 61 6C 2E 20 45 61
63 68 20 69 6E 70 75 74 20 58 2D 76 61 6C 75 65 20 6D 75 73 74 20 64 69
66 66 65 72 20 73 69 6E 63 65 20 66 72 6F 6D 20 61 20 6D 61 74 68 65 6D
61 74 69 63 61 6C 20 70 6F 69 6E 74 20 6F 66 20 76 69 65 77 20 69 74 20
6D 75 73 74 20 62 65 20 61 20 6F 6E 65 2D 74 6F 2D 6F 6E 65 20 6D 61 70
70 69 6E 67 2C 20 69 2E 65 2E 20 65 61 63 68 20 58 2D 76 61 6C 75 65 20
6D 75 73 74 20 63 6F 72 72 65 73 70 6F 6E 64 20 74 6F 20 65 78 61 63 74
6C 79 20 6F 6E 65 20 59 2D 76 61 6C 75 65 2E 27 29 3B 0A 09 7D 0A 0A 0A
09 24 61 20 3D 20 28 24 74 68 69 73 2D 3E 78 64 61 74 61 5B 24 6D 61 78
5D 2D 24 78 70 6F 69 6E 74 29 2F 24 68 3B 0A 09 24 62 20 3D 20 28 24 78
70 6F 69 6E 74 2D 24 74 68 69 73 2D 3E 78 64 61 74 61 5B 24 6D 69 6E 5D
29 2F 24 68 3B 0A 09 72 65 74 75 72 6E 20 24 61 2A 24 74 68 69 73 2D 3E
79 64 61 74 61 5B 24 6D 69 6E 5D 2B 24 62 2A 24 74 68 69 73 2D 3E 79 64
61 74 61 5B 24 6D 61 78 5D 2B 0A 09 20 20 20 20 20 28 28 24 61 2A 24 61
2A 24 61 2D 24 61 29 2A 24 74 68 69 73 2D 3E 79 32 5B 24 6D 69 6E 5D 2B
28 24 62 2A 24 62 2A 24 62 2D 24 62 29 2A 24 74 68 69 73 2D 3E 79 32 5B
24 6D 61 78 5D 29 2A 28 24 68 2A 24 68 29 2F 36 2E 30 3B 0A 20 20 20 20
7D 0A 7D 0A 0A 2F 2F 20 45 4F 46 0A 3F 3E | <?php /*===============
========================
========================
======== // File: JPGRAP
H_REGSTAT.PHP // Descrip
tion: Regression and sta
tistical analysis helper
classes // Created: 20
02-12-01 // Author: Joha
n Persson (johanp@aditus
.nu) // Ver: $Id: jpgra
ph_regstat.php 21 2005-0
5-30 20:35:34Z ljp $ //
// Copyright (c) Aditus
Consulting. All rights r
eserved. //=============
========================
========================
=========== */ //------
------------------------
------------------------
------------------ // CL
ASS Spline // Create a n
ew data array from an ex
isting data array but wi
th more points. // The n
ew points are interpolat
ed using a cubic spline
algorithm //------------
------------------------
------------------------
------------ class Splin
e { // 3:rd degree p
olynom approximation
var $xdata,$ydata; /
/ Data vectors var $
y2; // 2:nd derivate o
f ydata var $n=0;
function Spline($xda
ta,$ydata) { $this->y2
= array(); $this->xdata
= $xdata; $this->ydata
= $ydata; $n = count(
$ydata); $this->n = $n;
// Natural spline 2:d
erivate == 0 at endpoint
s $this->y2[0] = 0.0
; $this->y2[$n-1] = 0.0
; $delta[0] = 0.0; //
Calculate 2:nd derivate
for($i=1; $i < $n-1; +
+$i) { $d = ($xdata
[$i+1]-$xdata[$i-1]);
if( $d == 0 ) { Jp
GraphError::Raise('Inval
id input data for spline
. Two or more consecutiv
e input X-values are equ
al. Each input X-value m
ust differ since from a
mathematical point of vi
ew it must be a one-to-o
ne mapping, i.e. each X-
value must correspond to
exactly one Y-value.');
} $s = ($xdat
a[$i]-$xdata[$i-1])/$d;
$p = $s*$this->y2[$
i-1]+2.0; $this->y2
[$i] = ($s-1.0)/$p;
$delta[$i] = ($ydata[$i
+1]-$ydata[$i])/($xdata[
$i+1]-$xdata[$i]) -
($ydata[$i]-$yda
ta[$i-1])/($xdata[$i]-$x
data[$i-1]); $delta
[$i] = (6.0*$delta[$i]/(
$xdata[$i+1]-$xdata[$i-1
])-$s*$delta[$i-1])/$p;
} // Backward substit
ution for( $j=$n-2; $j
>= 0; --$j ) { $thi
s->y2[$j] = $this->y2[$j
]*$this->y2[$j+1] + $del
ta[$j]; } } //
Return the two new data
vectors function Ge
t($num=50) { $n = $this
->n ; $step = ($this->x
data[$n-1]-$this->xdata[
0]) / ($num-1); $xnew=a
rray(); $ynew=array();
$xnew[0] = $this->xdata
[0]; $ynew[0] = $this->
ydata[0]; for( $j=1; $j
< $num; ++$j ) { $
xnew[$j] = $xnew[0]+$j*$
step; $ynew[$j] = $
this->Interpolate($xnew[
$j]); } return array($
xnew,$ynew); }
// Return a single inter
polated Y-value from an
x value function Int
erpolate($xpoint) { $m
ax = $this->n-1; $min =
0; // Binary search t
o find interval while(
$max-$min > 1 ) { $
k = ($max+$min) / 2;
if( $this->xdata[$k] >
$xpoint ) $max=$k;
else $min=$k; }
// Each interval is i
nterpolated by a 3:degre
e polynom function $h =
$this->xdata[$max]-$thi
s->xdata[$min]; if( $h
== 0 ) { JpGraphE
rror::Raise('Invalid inp
ut data for spline. Two
or more consecutive inpu
t X-values are equal. Ea
ch input X-value must di
ffer since from a mathem
atical point of view it
must be a one-to-one map
ping, i.e. each X-value
must correspond to exact
ly one Y-value.'); }
$a = ($this->xdata[$max
]-$xpoint)/$h; $b = ($x
point-$this->xdata[$min]
)/$h; return $a*$this->
ydata[$min]+$b*$this->yd
ata[$max]+ (($a*$a
*$a-$a)*$this->y2[$min]+
($b*$b*$b-$b)*$this->y2[
$max])*($h*$h)/6.0;
} } // EOF ?> |