# Solution to STA 410/2102 Assignment #3, Spring 2003 - Test runs.


# PROBLEM 1.  MODEL FOR THE ABSOLUTE VALUES OF NORMAL DATA POINTS.

cat("\n\nTests for problem 1\n\n")

data1 <- scan("/u/radford/data1")
print(data1)

cat("\n"); print(abs.norm.em(data1,50,debug=T))        # Default initial value

cat("\n"); print(abs.norm.em(data1,50,0,debug=T))      # Other initial values
cat("\n"); print(abs.norm.em(data1,50,0.01,debug=T))
cat("\n"); print(abs.norm.em(data1,50,2,debug=T))
cat("\n"); print(abs.norm.em(data1,50,10,debug=T))


# -----------------------------------------------------------------------------


# PROBLEM 2.  MIXTURE OF NORMALS WITH MEANS OF 2*THETA AND 3*THETA.

cat("\n\nTests for problem 2\n\n")

data2 <- scan("/u/radford/data2")
print(data2)

cat("\n"); print(mix.em(data2,50,debug=T))             # Default initial values

cat("\n"); print(mix.em(data2,50,p=0.5,mean(data2),debug=T))    # Other initial
cat("\n"); print(mix.em(data2,50,p=0.5,0,debug=T))              #   values
cat("\n"); print(mix.em(data2,50,p=0.9,mean(data2)/2,debug=T))
cat("\n"); print(mix.em(data2,50,p=0.1,mean(data2)/3,debug=T))