title: Plot GLM model in R
body: <p>I am doing a linear regression model, so right now I have</p>
  
  <pre><code>mylogit &lt;- glm(Gender ~ Drug*Environment + Ethnicity + Age, data = mydata, family = "binomial")    
  </code></pre>
  
  <p>Gender (F/M), Drug (Y/N), Environment (H/L) are all factors and age is the only numeric. </p>
  
  <p>I am trying to graph a scatter plot with proportion of males (y-axis) vs. age (x-axis). There are four conditions, so four different colors/lines:</p>
  
  <p>1) Drug Y/Env H
  2) Drug N/Env H
  3) Drug N/Env L
  4) Drug Y/Env L</p>
  
  <p>I've tried using </p>
  
  <pre><code>plot(effect("Gender*Drug*Environment*Age", mylogit, xlevels=list(Age=20:60)), + multiline=TRUE, rug=FALSE)
  </code></pre>
  
  <p>But this just separates it into two graphs with two lines in each graph and does not plot the individual data points on there. </p>