ID      rho    lambda     value    counts    convergence
1       0.89    2.19      13.71    69         0
.        .      .         .        .          .
.        .      .         .        .          .
.        .      .         .        .          .

For (i in 1:length(ID)){
  f <- function(param, df){
   rho <- param[1]
   alpha <- param[2]
   logl <- #model

 }
 optim(...)
 # I may write a table to export the data, but that's not that organized and clean for manipulation.
}