A=c(0,1,2)
r=5 ; c=3
data=matrix(NA,nrow = r+1,ncol=c+1)
for(i in 1:c){
data[1:r,i]=sample(A,r, replace = TRUE)}
data
> data
     [,1] [,2] [,3] [,4]
[1,]    1    0    0   NA
[2,]    1    2    1   NA
[3,]    0    2    0   NA
[4,]    1    1    0   NA
[5,]    2    0    1   NA
[6,]   NA   NA   NA   NA
#0 counts for columns 
for(i in 1:c){
 count=0
 for(j in 1:r){
  if(data[j,i]==0) count=count+1 } 
  data[r+1,i]=count}
#0 counts for rows 
for(j in 1:r){
 count=0
 for(i in 1:c){
  if(data[j,i]==0) count=count+1  } 
  data[j,c+1]=count}
data
> data
     [,1] [,2] [,3] [,4]
[1,]    1    0    0    2
[2,]    1    2    1    0
[3,]    0    2    0    2
[4,]    1    1    0    1
[5,]    2    0    1    1
[6,]    1    2    3   NA