>
x=c(4381,4733,4711,4152,4810,5021,298,3904,4465,4133,4705,3964)
> x
[1] 4381 4733 4711 4152 4810 5021 298 3904 4465 4133 4705
3964
> S2=var(x) ### Sample variance
> S2
[1] 1565452
> n=length(x)##sample size
> n
[1] 12
> a=0.05##level of significance
> ql=qchisq(1-(a/2),(n-1)) ###upper chisquare value
> ql
[1] 21.92005
> qu=qchisq((a/2),(n-1)) ###lower chisquare value
> qu
[1] 3.815748
> LCL=sqrt(((n-1)*S2)/ql)##lower confidence limit
> LCL
[1] 886.3301
> UCL=sqrt(((n-1)*S2)/qu)##upper confidence limit
> UCL
[1] 2124.352
> CL=c(LCL,UCL)
> CL
[1] 886.3301 2124.3517
>