  Re: st: graph mean and sd over/by category

 From vwiggins@stata.com (Vince Wiggins, StataCorp) To statalist@hsphsun2.harvard.edu Subject Re: st: graph mean and sd over/by category Date Mon, 10 Feb 2003 16:55:20 -0600

```Among other things, George Hoffman <ghoffman@mcw.edu> asks,

> [...] fitted curves under scatter plots look beautiful - can the
> regression coefficients from fplotci or qplotci be captured somehow,
> as poor-man's curve fit?

I think George is referring to the -fpfitci- and -qfitci- plot types of
-graph twoway-.  If so, he can readily perform the regressions that produced
the graphs.

-qfitci- just performs a quadratic regression.  If we use the auto data,
-sysuse auto-, the lines for the graph command,

. twoway qfitci mpg weight

are the predictions of the quadratic fit,

. gen weight2 = weight^2
. regress mpg weight weight2

The coefficients can be seen in the output of -regress-, or manipulated in the
usual way through the saved results.

If George wants to add the predictions, and their CIs to his dataset, he can
type,

. predictnl mpg_hat2 = predict() , ci(ci_low ci_high)

This is a very simple application of -predictnl-, Bobby Gutierrez
<rgutierrez@stata.com> said more in a prior post, but it lets us get both the
predictions and their CIs with one command.

We could then get a graph similar to our earlier -twoway qfitci-, by typing,

. twoway rarea ci_low ci_high weight, sort || line mpg_hat2 weight, sort

which we will immediately think looks ugly and decide to relabel the CI in the
legend, option -legend(label())-, and change the fill color of the CI to be
the standard for our scheme, option -p(ci)-.

.  twoway rarea ci_low ci_high weight , sort p(ci) ||
line mpg_hat2 weight , sort legend(label(1 "CI"))

The -fpfitci- plot type just uses -fracpoly- as the engine to produce the
fits, much like -regress- is used for the quadratic fit.  For our example, the
corresponding -fracpoly- estimation command is,

. fracpoly regress mpg weight

and we can repeat the rest of the story, or just use -fracplot-, to plot the
fit and CI.

-- Vince
vwiggins@stata.com

*
*   For searches and help try:
*   http://www.stata.com/support/faqs/res/findit.html
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/
```