Dear All, I have a dataset with two variables (simplified). The basic statistics
Code:
. sum indep ir2d1

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
       indep |     35,432    .3323226    .1119013          0   .5555556
       ir2d1 |     34,059    .5667652    .1640421   .2095384   .9267846
When I use, say the first 1,000 observations, the command marginscontplot works and gives me the figure. However, when I increase to say, 10,000, the estimation result is
Code:
. keep in 1/10000
(47,372 observations deleted)

. reg indep ir2d1 c.ir2d1#c.ir2d1, 

      Source |       SS           df       MS      Number of obs   =     8,888
-------------+----------------------------------   F(2, 8885)      =      5.93
       Model |  .218049191         2  .109024596   Prob > F        =    0.0027
    Residual |  163.360103     8,885  .018386056   R-squared       =    0.0013
-------------+----------------------------------   Adj R-squared   =    0.0011
       Total |  163.578153     8,887  .018406454   Root MSE        =     .1356

---------------------------------------------------------------------------------
          indep |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
----------------+----------------------------------------------------------------
          ir2d1 |   .1730636   .0502702     3.44   0.001     .0745223    .2716048
                |
c.ir2d1#c.ir2d1 |  -.1495413   .0442941    -3.38   0.001    -.2363681   -.0627146
                |
          _cons |   .2565559   .0135907    18.88   0.000      .229915    .2831968
---------------------------------------------------------------------------------
and error message is
Code:
. marginscontplot ir2d1, ci
invalid numlist has too many elements
r(123);
It seems to me that this mcp/marginscontplot can not be applied to a large dataset? Any suggestions? Thanks.