curve(dbeta(x,.75,.75),col="red",lty=2,ylab="",xlab = "Value", xlim=c(0,1),ylim=c(0,3),lwd=4,axes=FALSE,main="Beta",cex.lab=1.5,col.lab=c("mediumvioletred"),col.main="mediumvioletred")
curve(dbeta(x4,4),col="blue",lty=1,ylab="",xlab = "Value",xlim=c(0,1),ylim=c(0,3),lwd=4,main="Beta",cex.lab=1.5,col.lab=c("blue"),,col.main="mediumvioletred",add=T)
axis(1,0:1, 0:1,col.axis="mediumvioletred",col="mediumvioletred",lty=1,lwd =4,cex.axis=1.5)
axis(2,  seq(0,3,.5), seq(0,3,.5),col.axis="mediumvioletred",col="mediumvioletred",lty=1,lwd =4,cex.axis=1.5)
curve(dbeta(x,.75,.75),ylab="",xlab = "Value",
      xlim=c(0,1),ylim=c(0,3),main="Beta")
curve(dbeta(x,4,4),ylab="",xlab = "Value",
      xlim=c(0,1),ylim=c(0,3),main="Beta",add=T)

A<- c(10,20 , 50, NA, 18 ,NA,40,16,.5,.01,2.2,3.5,NA,NA,3,4.02)
length(A)
Mat=matrix(A,4,4,byrow=TRUE)
Mat[is.na(Mat)]=0
apply(Mat,2,sum)
apply(Mat,1,sum)