nma — Oct 16, 2013, 2:50 PM
z <- seq(0,1,length=1000)
par(mfrow=c(3,3))
# Parametros: alfa=0.5 beta=0.5
plot(z,dbeta(z,0.5,0.5),xlab="theta",ylab="Beta(0.5,0.5)",type="l")
# Parametros: alfa=1 beta=0.5
plot(z,dbeta(z,1,0.5),xlab="theta",ylab="Beta(1,0.5)",type="l")
# Parametros: alfa=3 beta=0.5
plot(z,dbeta(z,3,0.5),xlab="theta",ylab="Beta(3,0.5)",type="l")
# Parametros: alfa=0.5 beta=1
plot(z,dbeta(z,0.5,1),xlab="theta",ylab="Beta(0.5,1)",type="l")
# Parametros: alfa=1 beta=1
plot(z,dbeta(z,1,1),xlab="theta",ylim=c(0,1.5),ylab="Beta(1,1)",type="l")
# Parametros: alfa=3 beta=1
plot(z,dbeta(z,3,1),xlab="theta",ylab="Beta(3,1)",type="l")
# Parametros: alfa=0.5 beta=2
plot(z,dbeta(z,0.5,3),xlab="theta",ylab="Beta(0.5,3)",type="l")
# Parametros: alfa=1 beta=2
plot(z,dbeta(z,1,3),xlab="theta",ylab="Beta(1,3)",type="l")
# Parametros: alfa=3 beta=2
plot(z,dbeta(z,3,3),xlab="theta",ylab="Beta(3,3)",type="l")