Update 18 January 2013: This example works for Zelig versions < 4. For versions from 4 you will likely have to use Zelig's simulation.matrix command to extract the matrix of expected values. In my previous post I briefly mentioned an early draft of a working paper ( HERE ) I've written that looks into the possible causes of violence between legislators (like the violence shown in this picture from the Turkish Parliament). From The Guardian In this post I'm going to briefly discuss how I used Zelig 's rare events logistic regression ( relogit) and ggplot2 in R to simulate and plot the legislative violence probabilities that are in the paper. In this example I am plotting simulated probabilities at fitted values on three variables: Age of democracy ( Polity IV > 5) A dichotomous electoral proportionality variable where 1 is high proportionality, 0 otherwise (see here for more details) Governing parties' majority (as a %