Now we can calculate the total transition potential for the first step, 1956-57
model_TTP <- (model_accessibility)*(model_suitability+1)*(model_nhood+1)*(model_random)
Note that we added 1 already to accessibility and random.
This TTP plot contains infinite values, maybe we don’t want these. To get rid of them, just recompute using “acc” instead of model_accessibility
model_TTP_2 <- (acc)*(model_suitability+1)*(model_nhood+1)*(model_random)