title: Solve best fit polynomial and plot drop-down lines
body: <p>I'm using R 3.3.1 (64-bit) on Windows 10.  I have an x-y dataset that I've fit with a 2nd order polynomial. I'd like to solve that best-fit polynomial for x at y=4, and plot drop-down lines from y=4 to the x-axis.</p>
  
  <p>This will generate the data in a dataframe v1:</p>
  
  <pre><code>v1 &lt;- structure(list(x = c(-5.2549, -3.4893, -3.5909, -2.5546, -3.7247, 
  -5.1733, -3.3451, -2.8993, -2.6835, -3.9495, -4.9649, -2.8438, 
  -4.6926, -3.4768, -3.1221, -4.8175, -4.5641, -3.549, -3.08, -2.4153, 
  -2.9882, -3.4045, -4.6394, -3.3404, -2.6728, -3.3517, -2.6098, 
  -3.7733, -4.051, -2.9385, -4.5024, -4.59, -4.5617, -4.0658, -2.4986, 
  -3.7559, -4.245, -4.8045, -4.6615, -4.0696, -4.6638, -4.6505, 
  -3.7978, -4.5649, -5.7669, -4.519, -3.8561, -3.779, -3.0549, 
  -3.1241, -2.1423, -3.2759, -4.224, -4.028, -3.3412, -2.8832, 
  -3.3866, -0.1852, -3.3763, -4.317, -5.3607, -3.3398, -1.9087, 
  -4.431, -3.7535, -3.2545, -0.806, -3.1419, -3.7269, -3.4853, 
  -4.3129, -2.8891, -3.0572, -5.3309, -2.5837, -4.1128, -4.6631, 
  -3.4695, -4.1045, -7.064, -5.1681, -6.4866, -2.7522, -4.6305, 
  -4.2957, -3.7552, -4.9482, -5.6452, -6.0302, -5.3244, -3.9819, 
  -3.8123, -5.3085, -5.6096, -6.4557), y = c(0.99, 0.56, 0.43, 
  2.31, 0.31, 0.59, 0.62, 1.65, 2.12, 0.1, 0.24, 1.68, 0.09, 0.59, 
  1.23, 0.4, 0.36, 0.49, 1.41, 3.29, 1.22, 0.56, 0.1, 0.67, 2.38, 
  0.43, 1.56, 0.07, 0.08, 1.53, -0.01, 0.12, 0.1, 0.04, 3.42, 0.23, 
  0, 0.34, 0.15, 0.03, 0.19, 0.17, 0.2, 0.09, 2.3, 0.07, 0.15, 
  0.18, 1.07, 1.21, 3.4, 0.8, -0.04, 0.02, 0.74, 1.59, 0.71, 10.64, 
  0.64, -0.01, 1.06, 0.81, 4.58, 0.01, 0.14, 0.59, 7.35, 0.63, 
  0.17, 0.38, -0.08, 1.1, 0.89, 0.94, 1.52, 0.01, 0.1, 0.38, 0.02, 
  7.76, 0.72, 4.1, 1.36, 0.13, -0.02, 0.13, 0.42, 1.49, 2.64, 1.01, 
  0.08, 0.22, 1.01, 1.53, 4.39)), .Names = c("x", "y"), class = "data.frame", row.names = c(NA, 
  -95L))
  </code></pre>
  
  <p>Here's the code to plot y vs x, plot the best fit polynomial, and draw a line at y=4.</p>
  
  <pre><code>&gt; attach(v1)
  &gt; # simple x-y plot of the data
  &gt; plot(x,y, pch=16)
  &gt; # 2nd order polynomial fit
  &gt; fit2 &lt;- lm(y~poly(x,2,raw=TRUE))
  &gt; summary(fit2)
  &gt; # generate range of numbers for plotting polynomial
  &gt; xx &lt;- seq(-8,0, length=50)
  &gt; # overlay best fit polynomial
  &gt;lines(xx, predict(fit2, data.frame(x=xx)), col="blue")
  &gt; # add horizontal line at y=4
  &gt; abline(h=4, col="red")
  &gt;
  </code></pre>
  
  <p>It's obvious from the plot that y=4 at x of around -2 and -6.5, but I'd like to actually solve the regression polynomial for those values.</p>
  
  <p>Ideally, I'd like lines that drop down from the red-blue line intersections to the x-axis (i.e plot vertical ablines that terminate at the two y=4 solutions). If that's not possible, I'd be happy with good old vertical ablines that go all the way up the plot, so long as they at the proper x solution values.</p>
  
  <p>This graph represents parts that will be out-of-spec when y>4, so I want to use the drop-down lines to highlight the range of x values that will produce in-spec parts.</p>