# TESTS FOR ASSIGNMENT 2.

d1 <- scan("/u/radford/a2data1")
d2 <- scan("/u/radford/a2data2")
d3 <- scan("/u/radford/a2data3")

cat("\n\nPROBLEM 2\n\n")

cat("\nDATA SET 1\n\n")

mle.mu(d1,0.0,log(0.5),debug=TRUE)
cat("\n")
mle.mu(d1,2.5,log(0.5),debug=TRUE)
cat("\n")
mle.mu(d1,5.0,log(0.5),debug=TRUE)
cat("\n")
mle.mu(d1,4.5,log(0.5),debug=TRUE)

cat("\nDATA SET 2\n\n")

mle.mu(d2,0.0,log(0.5),debug=TRUE)
cat("\n")
mle.mu(d2,2.5,log(0.5),debug=TRUE)
cat("\n")
mle.mu(d2,5.0,log(0.5),debug=TRUE)

cat("\nDATA SET 3\n\n")

mle.mu(d3,0.0,log(0.5),debug=TRUE)
cat("\n")
mle.mu(d3,2.5,log(0.5),debug=TRUE)
cat("\n")
mle.mu(d3,5.0,log(0.5),debug=TRUE)

cat("\n\nPROBLEM 3\n\n")

cat("\nDATA SET 1\n\n")

print (sqrt(diag( mle.joint(d1,0.0,log(0.5),debug=TRUE)$cov )))
cat("\n")
print (sqrt(diag( mle.joint(d1,2.5,log(0.5),debug=TRUE)$cov )))
cat("\n")
print (sqrt(diag( mle.joint(d1,5.0,log(0.5),debug=TRUE)$cov )))
cat("\n")
print (sqrt(diag( mle.joint(d1,1.3,log(0.1),debug=TRUE)$cov )))
cat("\n")
print (sqrt(diag( mle.joint(d1,1.35,log(0.1),debug=TRUE)$cov )))

cat("\nDATA SET 2\n\n")

print (sqrt(diag( mle.joint(d2,0.0,log(0.5),debug=TRUE)$cov )))
cat("\n")
print (sqrt(diag( mle.joint(d2,2.5,log(0.5),debug=TRUE)$cov )))
cat("\n")
print (sqrt(diag( mle.joint(d2,5.0,log(0.5),debug=TRUE)$cov )))
cat("\n")
print (sqrt(diag( mle.joint(d2,2.0,log(1.0),debug=TRUE)$cov )))

cat("\nDATA SET 3\n\n")

print (sqrt(diag( mle.joint(d3,0.0,log(0.5),debug=TRUE)$cov )))
cat("\n")
print (sqrt(diag( mle.joint(d3,2.5,log(0.5),debug=TRUE)$cov )))
cat("\n")
print (sqrt(diag( mle.joint(d3,5.0,log(0.5),debug=TRUE)$cov )))
cat("\n")
print (sqrt(diag( mle.joint(d3,4.78,log(1.0),debug=TRUE)$cov )))