Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.

# Re: st: documentation on iteration for a non linear regression

 From Austin Nichols To statalist@hsphsun2.harvard.edu Subject Re: st: documentation on iteration for a non linear regression Date Wed, 13 Jun 2012 11:39:35 -0400

```tashi lama <ltashi32@hotmail.com>:
You seem to be assuming additive Gaussian errors for a count outcome,
which strikes me as a bad idea. Compare the plausibility of the
various equivalent models below, and read the corresponding manual
entries for those commands.

clear
input x y
1 211
2 60
3 28
4 16
end
nl exp2: y x, nolog
di [b1]_b[_cons], [b2]_b[_cons]
nl exp2 y x, nolog
di _b[b1], _b[b2]
glm y x, link(log) fam(gaussian) nolog
di exp(_b[_cons]), exp(_b[x])
glm y x, link(log) fam(poisson) nolog
di exp(_b[_cons]), exp(_b[x])
poisson y x, nolog
di exp(_b[_cons]), exp(_b[x])

On Wed, Jun 13, 2012 at 11:04 AM, tashi lama <ltashi32@hotmail.com> wrote:
> That would help. Thanx... but it doesn't hurt to learn more although I am planning to spend a day or two to go through the manual on matrix, sth I haven't looked at before. I had a dataset as follows
>
>
>
>            day     hits   b          v |
>     |---------------------------------|
>  1. | 03jan2011    211   1   209.0752 |
>  2. | 04jan2011     60 2 69.08601 |
>  3. | 05jan2011     28 3 22.82851 |
>  4. | 06jan2011     16   4   7.543365 |
> ----------------------------------------
>
>
>
> v is my predicted value. day and hits are my original value and I generated b to use it as independant variable instead of day for simplification. Drawing a scatter plot helped me identify the nature of the curve which is exponential so i decided to use one of the built in exp models. So, I said
>
>
>
> nl exp2 :hits b
>

*
*   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/
```