[R-br] Colocar estatística Z na curva normal - TH

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