diff --git a/R/CytoDiv.R b/R/CytoDiv.R
index a41042f..b837eaa 100644
--- a/R/CytoDiv.R
+++ b/R/CytoDiv.R
@@ -16,7 +16,7 @@ cytoDiv <- function(df, para, Ncat, log=TRUE){
         J = log(N1)/log(N0) # Evenness
         
         # output
-        indices <- data.frame(cbind(N0, H, J))
+        indices <- data.frame(cbind(N0, N1, H, J))
         return(indices)