select_poisson.R

t.test(rnorm(10,mean=0.5), mu=0)

counters = c(2111, 1728, 2767, 2778, 2391,  708, 1413,  967, 2112, 2699,  911, 1139, 2500)
c2 = c(175, 478,  5,537, 326, 394, 551, 260, 279, 668, 261, 591, 686, 315, 695, 545, 279, 679)
c3 = c(325, 257, 649, 684, 589, 689, 427, 577, 56, 412, 659, 523, 259, 671, 512, 326, 636, 700, 635, 540, 649)
c4 = c(128, 560)
c5 = c(39, 32, 38, 33, 14, 15, 37, 35, 31)
ed = density(c5)
plot(ed)


h = density(c3)
sum(h$y[h$x<56])/(sum(h$y))
plot(h)



X = matrix(0, nrow=3, ncol=10)
X[,1] = rnorm(n=3,mean=2)
X[,4] = X[,6] = X[,1] + rnorm(3, 0,0.25)

X[,2] = rnorm(3, 78,4)
X[,3] = rnorm(3, 34,4)
X[,5] = rnorm(3, 45,4)
X[,7] = rnorm(3, 67,4)
X[,8] = rnorm(3, 89,4)
X[,9] = rnorm(3, 12,4)
X[,10] = rnorm(3, 23,4)

t(X)%*%X