Y the observed value
X the vector of the normal distribution
n the number of X

(Y-mean(X))/(sd(X)*sqrt(1+1/n))->ts

two-side
p-value <- (1-pt(abs(ts),n-1))*2

one-side
less (distribution less than observed value)
p-value <- (1-pt(ts,n-1))
greater
p-value <- pt(ts,n-1)

benjamini hochberg method?