[o] GLIM 3.77 update 1 (copyright)1985 Royal Statistical Society, London [o] [i] ? $units 32 $data time $dinput 11 80 $ [i] File name? glass [i] 439 572 315 258 [i] 904 690 315 258 [i] 1092 904 439 347 [i] 1105 1090 628 588 [i] 959 216 241 241 [i] 1065 315 315 241 [i] 1065 455 332 435 [i] 1087 473 380 455 [i] ? $cal temp=%gl(2,16)*10+160 $print temp $ [o] 170.0 170.0 170.0 170.0 170.0 170.0 170.0 170.0 170.0 [o] 170.0 170.0 170.0 170.0 170.0 170.0 170.0 180.0 180.0 [o] 180.0 180.0 180.0 180.0 180.0 180.0 180.0 180.0 180.0 [o] 180.0 180.0 180.0 180.0 180.0 [i] ? $cal kv=%gl(4,1)*50+150 $print kv $ [o] 200.0 250.0 300.0 350.0 200.0 250.0 300.0 350.0 200.0 [o] 250.0 300.0 350.0 200.0 250.0 300.0 350.0 200.0 250.0 [o] 300.0 350.0 200.0 250.0 300.0 350.0 200.0 250.0 300.0 [o] 350.0 200.0 250.0 300.0 350.0 [i] ? $plot time kv $ [o] 1150.0 | [o] 1100.0 | 3 T [o] 1050.0 | 2 [o] 1000.0 | [o] 950.0 | T [o] 900.0 | T T [o] 850.0 | [o] 800.0 | [o] 750.0 | [o] 700.0 | T [o] 650.0 | T [o] 600.0 | T [o] 550.0 | T [o] 500.0 | [o] 450.0 | T 2 T 2 [o] 400.0 | T [o] 350.0 | T T [o] 300.0 | T 3 [o] 250.0 | T 4 [o] 200.0 | T [o] 150.0 | [o] ----------:---------:---------:---------:---------:---------:---------: [o] 195.0 225.0 255.0 285.0 315.0 345.0 375.0 [i] ? $cal ltime=%log(time) $plot ltime kv $ [o] 7.200 | [o] 7.100 | [o] 7.000 | 5 L [o] 6.900 | L [o] 6.800 | L L [o] 6.700 | [o] 6.600 | [o] 6.500 | L [o] 6.400 | L L [o] 6.300 | L [o] 6.200 | L [o] 6.100 | L L L 2 [o] 6.000 | [o] 5.900 | L [o] 5.800 | L 4 L [o] 5.700 | [o] 5.600 | 2 [o] 5.500 | L 2 [o] 5.400 | L [o] 5.300 | [o] 5.200 | [o] ----------:---------:---------:---------:---------:---------:---------: [o] 195.0 225.0 255.0 285.0 315.0 345.0 375.0 [i] ? $yvar time $err gamma $link log $fit temp+kv $ [o] deviance = 4.0016 at cycle 3 [o] d.f. = 29 [o] [i] ? $dis e $ [o] estimate s.e. parameter [o] 1 11.88 2.322 1 [o] 2 -0.02119 0.01313 TEMP [o] 3 -0.006938 0.001175 KV [o] scale parameter taken as 0.1380 [o] [i] ? $cal t2=temp**2 : k2=kv**2 : tk=temp*kv $ [i] ? $fit +t2+k2+tk $dis e $ [o] deviance = 3.4827 (change = -0.5189) at cycle 3 [o] d.f. = 27 (change = -2 ) [o] [o] estimate s.e. parameter [o] 1 15.24 11.31 1 [o] 2 -0.01910 0.06374 TEMP [o] 3 -0.03152 0.04215 KV [o] 4 0.000 aliased T2 [o] 5 0.00005142 0.00002540 K2 [o] 6 -2.130e-05 0.0002271 TK [o] scale parameter taken as 0.1290 [o] [i] ? $cal u=1-%exp(-time/%fv) : hatu=(%cu(1)-0.5)/32 $ [i] ? $sort u $plot u hatu hatu $ [o] 1.0000 | H [o] 0.9500 | H [o] 0.9000 | HH [o] 0.8500 | H H U [o] 0.8000 | H [o] 0.7500 | HH U U UU U [o] 0.7000 | U2 UU U [o] 0.6500 | UU UU 2 H [o] 0.6000 | U UU UU HH [o] 0.5500 | U UU H [o] 0.5000 | UU U UU U H H [o] 0.4500 | H [o] 0.4000 | HH [o] 0.3500 | U H H [o] 0.3000 |U H [o] 0.2500 | HH [o] 0.2000 | H [o] 0.1500 | H H [o] 0.1000 | HH [o] 0.0500 | H [o] 0.0000 |H [o] ----------:---------:---------:---------:---------:---------:---------: [o] 0.000 0.200 0.400 0.600 0.800 1.000 1.200 [i] ? $cal d=u-hatu $print d $ [o] 0.3091 0.3152 0.3977 0.3665 0.3402 0.3387 0.3119 0.2807 0.2791 [o] 0.2479 0.2264 0.2212 0.2158 0.1851 0.1590 0.1391 0.1106 0.0795 [o] 0.0497 0.0259 0.0210 0.0034 -0.0273 -0.0346 -0.0658 -0.0897 -0.0697 [o] -0.0956 -0.1200 -0.1510 -0.1795 -0.1535 [i] ? $stop