• Contenuto: grafico della funzione di densità di probabilità della variabile casuale di Birnbaum-Saunders
• Licenza: Gnu FDL
• Note: grafico prodotto con script in R

Codice R

png(file="c:/temp/Grafico_v.c._Birnbaum-Saunders.png", width=900, height=800);
BS <- function(x,a,b) { 
    return( 1/sqrt(2*pi) * (sqrt(x/b)+sqrt(b/x))/(2*a*x)*exp(-((sqrt(x/b)-sqrt(b/x))^2/(2*a^2))))
}
dx=0.001;x <- seq(dx,3,dx);

alfas <- c(0.1,0.2,0.5,1,2,3,5);
colori <- c("red","green","black","purple","navy","orange","darkgreen")

ylim <- c(0,max(BS(x,min(alfas),1),BS(x,max(alfas),1)))

plot(x,BS(x,1,1),
     main="v.c. di Birnbaum-Saunders",
     type='l',
     ylab=expression(f(x,alpha,beta==1)),
     ylim=ylim)
i<-0;for (alfa in alfas) {i<-i+1;lines(x,BS(x,alfa,1),col=colori[i],lwd=1.2)}

legenda <- legend(max(x)/2,max(ylim)*4/5,alfas,col=colori,lty=1,bty='n')
text(legenda$rect$left+legenda$rect$w*2/3,max(ylim)*4/5,expression(alpha))
dev.off()