sampleSize=21
data1=c(18,16,18,24,23,22,22,23,26,29,32,34,34,36,36,43,42,49,46,46,57)
var=var(data1); var
[1] 134.2619
s=sqrt(var); s
[1] 11.58714
mean=mean(data1); mean
[1] 32.19048
SSE=s/sqrt(sampleSize); SSE
[1] 2.528522
#95% CI using t-distribution because of small sample size
lower.boundary=mean-2.09*SSE
upper.boundary=mean+2.09*SSE
lower.boundary;upper.boundary
[1] 26.90586
[1] 37.47509