> source("ass2-sol.s") > data<-read.table("/u/radford/ass2.data",head=T) > m<-mle(data$c,data$t,debug=T) 0 4.75359 0.2540835 45 20308.65 1 4.60046 0.2315366 40.78563 20533.82 2 4.605101 0.1794579 41.08542 20542.20 3 4.594065 0.1976949 40.42915 20545.36 4 4.592941 0.2015862 40.35751 20545.47 5 4.59287 0.2017605 40.35319 20545.47 6 4.59287 0.2017608 40.35318 20545.47 7 4.59287 0.2017608 40.35318 20545.47 > m $eta [1] 4.59287 $kappa [1] 0.2017608 $theta [1] 40.35318 $cov [,1] [,2] [,3] [1,] 2.240135e-04 -4.431895e-05 0.003347593 [2,] -4.431895e-05 8.767253e-05 -0.002660335 [3,] 3.347593e-03 -2.660335e-03 0.201222205 $log.like [1] 20545.47 > sqrt(diag(m$cov)) [1] 0.01496708 0.00936336 0.44857798 > m<-mle(data$c,data$t*100,debug=T) 0 4.75359 0.002540835 4500 20308.65 1 4.60046 0.002315366 4078.563 20533.82 2 4.605101 0.001794579 4108.542 20542.20 3 4.594065 0.001976949 4042.915 20545.36 4 4.592941 0.002015862 4035.751 20545.47 5 4.59287 0.002017605 4035.319 20545.47 6 4.59287 0.002017608 4035.318 20545.47 7 4.59287 0.002017608 4035.318 20545.47 > m $eta [1] 4.59287 $kappa [1] 0.002017608 $theta [1] 4035.318 $cov [,1] [,2] [,3] [1,] 2.240135e-04 -4.431895e-07 3.347593e-01 [2,] -4.431895e-07 8.767253e-09 -2.660335e-03 [3,] 3.347593e-01 -2.660335e-03 2.012222e+03 $log.like [1] 20545.47 > sqrt(diag(m$cov)) [1] 1.496708e-02 9.363361e-05 4.485780e+01