source("AltFindPi.R")
source("Gamma.R")
source("Confusion.R")

testp <- rbind(
   c(0.2,0.2,0.6),
   c(0.0,0.2,0.8),
   c(0.8, 0.2,0.0)) 

lgt <- function(t) 1 / (1 + exp(-t))
n <- 100
x <- rnorm(n)
half <- 1:(n/2)
x[half] <- x[half] + 1
probs <- lgt(x)
u <- runif(n)
y <- ifelse(probs < u,1,0)

print("correct answers:")
print(nmaltfindpi(testp,4))
print(nmlambbias(2,1,10,25000))
testg <- glm(y ~ x,family=binomial,x=T,y=T)
print(nmconfmat(testg))