x= c(0,1,2,2,0,0,1)
#utiliza "cumsum" para hacer la suma acumulada 
# utiliza "admargins" para hacer la suma acumlada
#paso 1: construir la matriz de frecuencias absolutas (table), relativas, acumulada 
# paso 2: añadir los totales de la matriz
a= table(x) #frecuencias absolutas 
b= a / sum (a) #fracuencias relativas 
d= cumsum (a) #frecuencias absolutas acumuladas 
e= cumsum (b) #frecuencias relativas acumuladas
z= cbind ("F.abs" = a,"F.rel"=b,"F.abs.acum"= d, "F.rel.acum"=e)
z= addmargins (z)
#z= z [,1:4]
z= z[,-5]
z[4,c(3,4)]= NA
z [,c(3,4)]
z[,c("F.abs.acum","F.rel.acum")]
z[c(1,3),]
z[c("0","2"),]