calcularIntervaloDeConfiancaPara95 = function(x) mean(x) + qt(c(0.025,0.975),df = length(x)-1)*sqrt(var(x)/length(x))

calcularIt95 = function(tam) qt(c(0.025,0.975),df = tam-1)

calcularIntervaloDeConfiancaPara99 = function(x) mean(x) + qt(c(0.005,0.995),df = length(x)-1)*sqrt(var(x)/length(x))

calcularIt99 = function(tam) qt(c(0.005,0.995),df = tam-1)

calcularGraficoNormal = function(tam) curve(dnorm(x),xlim=c(-3,3))

pintarGrafico = function(x,y,tam){
	intervalo = 0.0001
	cord.x = c(x,seq(x,y,intervalo),y)
	cord.y = c(0,dt(seq(x,y,intervalo),df=tam-1),0)
	polygon(cord.x,cord.y,col = 'skyblue')
}

Como usar :

## Calculando Intervalo de Confiança (95%)
intervaloDeConfiancaPara95 = calcularIntervaloDeConfiancaPara95(whitePerCap)
it95 = calcularIt95(n)
graficoNormal95 = calcularGraficoNormal(n)
#Agora só falta pintar as áreas que estão fora do intervalo de confiança
	pintarGrafico(-4,it95[1],n)
	pintarGrafico(it95[2],4,n)
	pintarGrafico(0,0,n)
## Fim Calculando Intervalo de Confiança (95%)

## Calculando Intervalo de Confiança (99%)
intervaloDeConfiancaPara99 = calcularIntervaloDeConfiancaPara99(whitePerCap)
it99 = calcularIt99(n)
graficoNormal99 = calcularGraficoNormal(n)
#Agora só falta pintar as áreas que estão fora do intervalo de confiança
	pintarGrafico(-4,it99[1],n)
	pintarGrafico(it99[2],4,n)
	pintarGrafico(0,0,n)
## Fim Calculando Intervalo de Confiança (99%)