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)

> x=seq(as.Date("2000/1/1"), by = "month", length.out = 10)
> Date=rep(x,each=3)
> Date
 [1] "2000-01-01" "2000-01-01" "2000-01-01" "2000-02-01" "2000-02-01"
 [6] "2000-02-01" "2000-03-01" "2000-03-01" "2000-03-01" "2000-04-01"
[11] "2000-04-01" "2000-04-01" "2000-05-01" "2000-05-01" "2000-05-01"
[16] "2000-06-01" "2000-06-01" "2000-06-01" "2000-07-01" "2000-07-01"
[21] "2000-07-01" "2000-08-01" "2000-08-01" "2000-08-01" "2000-09-01"
[26] "2000-09-01" "2000-09-01" "2000-10-01" "2000-10-01" "2000-10-01"

> MAT=diag(seq(1:30),30)
> dim(MAT)
[1] 30 30
> colnames(MAT)=c("Date",LETTERS[1:10],LETTERS[15:24],LETTERS[9:1])
> Data=as.data.frame(MAT)
> Data[,1]=Date
> head(Data)
        Date A B C D E F G H I J O P Q R S T U V W X I H G F E D C B A
1 2000-01-01 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
2 2000-01-01 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
3 2000-01-01 0 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 2000-02-01 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
5 2000-02-01 0 0 0 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
6 2000-02-01 0 0 0 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

> A=apply(Data[,-1],1,sum)
> head(data.frame(Date,A=A))
        Date A
1 2000-01-01 0
2 2000-01-01 2
3 2000-01-01 3
4 2000-02-01 4
5 2000-02-01 5
6 2000-02-01 6