#!/usr/bin/Rscript
library(stats4)
data <- read.table("08_classif1.txt")
data <- read.table("08_classif2.txt")
y <- t(data[,2])

normal.lik1<-function(theta,y){
	mu<-theta[1]
	sigma2<-theta[2]
	n<-nrow(y)
	logl <- -.5*n*log(2*pi) -.5*n*log(sigma2) - (1/(2*sigma2))*sum((y-mu)**2)
	return(-logl)
}

mixResultsML <- optim(c(0,1), normal.lik1, y=y)
x0 <- mixResultsML$par[1]
p <- c(y > x0)
p2 <- data[,1] > 1
print(table(p, p2))