[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

From |
Steven Samuels <sjhsamuels@earthlink.net> |

To |
statalist@hsphsun2.harvard.edu |

Subject |
Re: st: RE: RE: RE: RE: RE: re: problem using predictnl after obtaining non-linear estimates |

Date |
Sun, 26 Oct 2008 14:16:11 -0400 |

-Steve On Oct 26, 2008, at 1:43 PM, Lopa Chakraborti wrote:

I do have data on upstream water quality which explains why the R-square is so high.________________________________________From: owner-statalist@hsphsun2.harvard.edu [owner-statalist@hsphsun2.harvard.edu] On Behalf Of Nick Cox[n.j.cox@durham.ac.uk]Sent: Sunday, October 26, 2008 1:41 PM To: statalist@hsphsun2.harvard.eduSubject: st: RE: RE: RE: RE: re: problem using predictnl afterobtaining non-linear estimatesI've never used this command. Others should have better advice on it.But if your model is not based on the flow data you would like,doesn'ta R-square of 98% worry you? How can you get a model fit that appears that good without a scientifically sound model? I don't have your data and I've not studied your model, but that implies overfitting to me. Nick n.j.cox@durham.ac.uk Lopa ChakrabortiThe problem I am having is figuring out what goes inside theparenthesisin predict() since it is not the linear predictor xb(). Help on this would be appreciated. It is a very rough and approximate model for Streeter-Phelps. I do not have any stream flow etc data and hence I am using the monitoring station dummies to control for them. Nick Cox [n.j.cox@durham.ac.uk]No programming is needed; this is something available to you aspart ofyour -predictnl- command.Kit is pointing you to special functions -predict()- and -xb()-that youcan use within -predictnl-. I too no longer have Stata 8.2 on mymachinebut I do have the Stata 8 manuals, which indicate to me that these functions were present in -predictnl- in Stata 8. Therefore study the on-line help first and then the manual entry.Incidentally, I feel very queasy about an 91-parameter model, evenwiththis kind of data (which is not alien to me). What's the engineering hydrology behind that? Lopa Chakraborti Thanks you have pointed me to the right direction. However, I am not using either predict() or xb() functions. Instead, I have: predictnl bod1 =($K2/($K1-$K2))*exp(-$K1*aggthirdfourthadd)*(exp(-$K2*seconddistadd)-exp(-$K1*seconddistadd)). That said my limited knowledge in programming might also explain why I cannot make the connection. Kit Baum [baum@bc.edu] Lopa said I have a general question about predictnl and that is does it evaluate any function as in any NL function of the parameters and some explanatory variables or does it only predict dependent variables as is usually done. [R] predictnl -- Obtain nonlinear predictions, standard errors, etc., after estimation That help indicates that " pnl_exp is any valid Stata expression" so that you may calculate any function for each observation, not just the standard predicted value. In Stata 10 help, When calculating inference-related quantities such as standard errors, pnl_exp is evaluated repeatedly for different values of the model parameters. Therefore, think of predict() and xb() as a means of substituting for the formula of the calculation and not a means of substituting the value of the calculation that is obtained whenthe model parameters are set to any specific values. Forexample,. predict double pred_var, predict_options . predictnl newvar = pred_var, se(newvar_se) will give standard errors (newvar_se) equal to zero, since once evaluated, pred_varwill contain values that are fixed with respect to e(b).Instead,. predictnl newvar = predict(predict_options), se(newvar_se) will produce what is intended. which states that zero standard errors will result. I have no idea whether it works this way in Stata 8.2. * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/ * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

* * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

**References**:**st: re: problem using predictnl after obtaining non-linear estimates***From:*Kit Baum <baum@bc.edu>

**st: RE: re: problem using predictnl after obtaining non-linear estimates***From:*Lopa Chakraborti <lchakraborti@arec.umd.edu>

**st: RE: RE: re: problem using predictnl after obtaining non-linear estimates***From:*"Nick Cox" <n.j.cox@durham.ac.uk>

**st: RE: RE: RE: re: problem using predictnl after obtaining non-linear estimates***From:*Lopa Chakraborti <lchakraborti@arec.umd.edu>

**st: RE: RE: RE: RE: re: problem using predictnl after obtaining non-linear estimates***From:*"Nick Cox" <n.j.cox@durham.ac.uk>

**st: RE: RE: RE: RE: RE: re: problem using predictnl after obtaining non-linear estimates***From:*Lopa Chakraborti <lchakraborti@arec.umd.edu>

- Prev by Date:
**st: RE: RE: RE: RE: RE: RE: re: problem using predictnl after obtaining non-linear estimates** - Next by Date:
**st: RE: RE: RE: RE: RE: RE: RE: re: problem using predictnl after obtaining non-linear estimates** - Previous by thread:
**st: RE: RE: RE: RE: RE: RE: RE: re: problem using predictnl after obtaining non-linear estimates** - Next by thread:
**st: prcounts to test poisson predictive ability?** - Index(es):

© Copyright 1996–2015 StataCorp LP | Terms of use | Privacy | Contact us | What's new | Site index |