zone<-function(alternative=c("two.sided","less","greater"),conf.level=0.95) { #Description # Visualizza le zone di rifiuto e di accettazione su una normale di media zero e scarto # quadratico medio pari ad uno #Usage # zone(alternative=c("two.sided","less","greater"),conf.level=0.95) #Arguments # alternative specifica l'ipotesi alternativa # conf.level intervallo di confidenza #Value # Grafico che visualizza in nero le zone di rifiuto x<-seq(-4,4,0.1) plot(x,dnorm(x),main="Densità di una distribuzione N(0,1)",xlab="",ylab="",axes=F,type="l",sub=paste(expression(alpha),"=",1-conf.level)) if(alternative=="two.sided") { k<-qnorm((1-conf.level)/2) si<-seq(-4,k,0.01) ss<-seq(-k,4,0.01) lines(si,dnorm(si),type="h") lines(ss,dnorm(ss),type="h") text(-3.5,0.4,paste("q=",round(k,digits=2))) text(3.5,0.4,paste("q=",round(-k,digits=2))) } else { if(alternative=="greater") { k<-qnorm(1-conf.level) ss<-seq(-k,4,0.01) lines(ss,dnorm(ss),type="h") text(3.5,0.4,paste("q=",round(-k,digits=2))) } else { k<-qnorm(1-conf.level) si<-seq(-4,k,0.01) lines(si,dnorm(si),type="h") text(-3.5,0.4,paste("q=",round(k,digits=2))) } } }