Thanks for your help. Ideally I would vectorize everything but my issue with doing that is I am estimating the parameters on a game by game basis as well as throughout the actual game based on when goals were actually scored. So while I might have a row of parameters for each game, I also have a list of times of goals for that specific game. With this structure I am not sure exactly how I would vectorize. I would have to really map out the model to figure it out. I am not even sure its feasible