Pessoal,
estou aqui tentando, mas sem sucesso colocar a estatística Ztab, Zcal,  mas sem sucesso.

media <- 170
stder <- 8
x <- seq(media-3*stder, media+3*stder, 0.01)
y <- dnorm(x, media, stder)
plot(x, y, type='l',axes=F,xlab="",ylab="")
abline(h=0)
hx <- seq(media+2*stder, media+3*stder, 0.01)
hy <- dnorm(hx, media, stder)
n <- length(hy)
polygon(c(hx, rev(hx)), c(hy, rep(0, n)), col=2)

# Tentei 

text(locator(n=1),"Zcal=3.53")
text(locator(n=1),"Ztab=1,64")

Mas, preciso aumentar a área abaixo do gráfico para plotar estes valores e não consegui! Alguém tem alguma sugestão? Outra situação é que Z_{cal} para ficar mais elegante não funcionou com locator.



obrigado