> source("ass4-sol.s") > xa <- scan("/u/radford/ass4a.data") Read 1000 items > em(xa,n=400,debug=10) 0 13.37983 16.35312 0.5 -3696.968 10 6.48944 23.31044 0.5019895 -3662.43 20 5.081083 22.26133 0.4304276 -3657.39 30 4.432508 21.24330 0.3793293 -3655.438 40 4.108623 20.69848 0.3515403 -3654.853 50 3.952688 20.42734 0.3375406 -3654.702 60 3.878179 20.2958 0.3307011 -3654.666 70 3.842604 20.23256 0.3274009 -3654.658 80 3.825614 20.20226 0.3258167 -3654.656 90 3.817497 20.18776 0.3250581 -3654.655 100 3.813619 20.18083 0.3246953 -3654.655 110 3.811765 20.17752 0.3245218 -3654.655 120 3.810880 20.17593 0.3244388 -3654.655 130 3.810456 20.17518 0.3243992 -3654.655 140 3.810254 20.17482 0.3243802 -3654.655 150 3.810157 20.17464 0.3243712 -3654.655 160 3.810111 20.17456 0.3243669 -3654.655 170 3.810089 20.17452 0.3243648 -3654.655 180 3.810079 20.1745 0.3243638 -3654.655 190 3.810074 20.17449 0.3243633 -3654.655 200 3.810071 20.17449 0.3243631 -3654.655 210 3.81007 20.17449 0.324363 -3654.655 220 3.810069 20.17449 0.3243629 -3654.655 230 3.810069 20.17448 0.3243629 -3654.655 240 3.810069 20.17448 0.3243629 -3654.655 250 3.810069 20.17448 0.3243629 -3654.655 260 3.810069 20.17448 0.3243629 -3654.655 270 3.810069 20.17448 0.3243629 -3654.655 280 3.810069 20.17448 0.3243629 -3654.655 290 3.810069 20.17448 0.3243629 -3654.655 300 3.810069 20.17448 0.3243629 -3654.655 310 3.810069 20.17448 0.3243629 -3654.655 320 3.810069 20.17448 0.3243629 -3654.655 330 3.810069 20.17448 0.3243629 -3654.655 340 3.810069 20.17448 0.3243629 -3654.655 350 3.810069 20.17448 0.3243629 -3654.655 360 3.810069 20.17448 0.3243629 -3654.655 370 3.810069 20.17448 0.3243629 -3654.655 380 3.810069 20.17448 0.3243629 -3654.655 390 3.810069 20.17448 0.3243629 -3654.655 400 3.810069 20.17448 0.3243629 -3654.655 $m1 [1] 3.810069 $m2 [1] 20.17448 $p [1] 0.3243629 > xb <- scan("/u/radford/ass4b.data") Read 1000 items > em(xb,n=4000,debug=100) 0 7.653706 9.35453 0.5 -3140.066 100 6.438168 10.49732 0.4910393 -3139.162 200 6.370806 10.44064 0.4758236 -3139.152 300 6.298352 10.38355 0.4600587 -3139.142 400 6.219858 10.32579 0.4436682 -3139.13 500 6.134112 10.26707 0.4265598 -3139.117 600 6.03954 10.20705 0.4086207 -3139.103 700 5.934049 10.14528 0.3897108 -3139.086 800 5.814791 10.08122 0.3696541 -3139.067 900 5.677783 10.01415 0.3482253 -3139.045 1000 5.517296 9.943094 0.3251338 -3139.018 1100 5.324856 9.866681 0.3000035 -3138.984 1200 5.087684 9.782926 0.2723626 -3138.94 1300 4.786829 9.688925 0.2416940 -3138.881 1400 4.398157 9.580798 0.2077475 -3138.797 1500 3.911283 9.456113 0.1716906 -3138.685 1600 3.3903 9.324982 0.1383165 -3138.566 1700 2.984077 9.218821 0.1146324 -3138.493 1800 2.758522 9.1573 0.1020791 -3138.469 1900 2.657774 9.12914 0.09658274 -3138.464 2000 2.617079 9.117644 0.09438037 -3138.463 2100 2.601291 9.113165 0.09352867 -3138.463 2200 2.595261 9.111451 0.09320374 -3138.463 2300 2.592971 9.1108 0.09308043 -3138.463 2400 2.592104 9.110553 0.09303372 -3138.463 2500 2.591776 9.11046 0.09301604 -3138.463 2600 2.591651 9.110425 0.09300934 -3138.463 2700 2.591604 9.110411 0.09300681 -3138.463 2800 2.591587 9.110406 0.09300586 -3138.463 2900 2.59158 9.110404 0.0930055 -3138.463 3000 2.591577 9.110404 0.09300536 -3138.463 3100 2.591576 9.110403 0.0930053 -3138.463 3200 2.591576 9.110403 0.09300528 -3138.463 3300 2.591576 9.110403 0.09300528 -3138.463 3400 2.591576 9.110403 0.09300527 -3138.463 3500 2.591576 9.110403 0.09300527 -3138.463 3600 2.591576 9.110403 0.09300527 -3138.463 3700 2.591576 9.110403 0.09300527 -3138.463 3800 2.591576 9.110403 0.09300527 -3138.463 3900 2.591576 9.110403 0.09300527 -3138.463 4000 2.591576 9.110403 0.09300527 -3138.463 $m1 [1] 2.591576 $m2 [1] 9.110403 $p [1] 0.09300527