# st: Re: re:locating the natural log of a number

 From "Rodrigo A. Alfaro" To Subject st: Re: re:locating the natural log of a number Date Sat, 31 Mar 2007 17:49:14 -0400

Victor,

I don't understand what are you doing with your code.
But I think that you want

set obs 22026
gen n = _n
gen double nl = log(_n)

Rodrigo.

----- Original Message ----- From: "Victor M. Zammit" <vmz@vol.net.mt>
To: <statalist@hsphsun2.harvard.edu>
Sent: Saturday, March 31, 2007 6:24 PM
Subject: st: re:locating the natural log of a number

Dear Statalist,

I am trying to write a program to locate the natural log of any number
between 1 and 22026.

In this specific case I am trying the number 22000 but I am having trouble
with the error

sign "obs. out of range" as I try to get to the exact natural log of the
number in question.

Can any one give me any suggestions as to how to handle the problem.

The following is the best I could do :

qui {

drop _all

set obs 10000

gen a = .

gen exp = .

local n = 1

while `n'<= 9999 {

local e =.001

while `e'<= 10 {

replace exp = `e' in `n'

replace a = `n' in `n'

local e = `e'+.001

local n = `n'+1

}

}

gen nl = 20000

gen x = 2.718281828^exp

gen diff = x-nl

keep diff exp

save diff,replace

keep if diff<0

qui summ exp

local l = r(max)

}

di `l'

di 2.718281828^`l'

drop _all

use diff,clear

keep if diff>0

qui summ exp

local h = r(min)

di `h'

di 2.718281828^`h'

qui {

drop _all

set obs 10000

gen a = .

gen exp = .

local n = 1

while `n'<= 9999 {

local e = `l'

while `e'<= `h' {

replace exp = `e' in `n'

replace a = `n' in `n'

local e = `e'+.00001

local n = `n'+1

}

}

gen nl = 20000

gen x = 2.718281828^exp

gen diff = x-nl

keep diff exp

save diff,replace

keep if diff<0

qui summ exp

local l = r(max)

}

di `l'

di 2.718281828^`l'

drop _all

use diff,clear

keep if diff>0

qui summ exp

local h = r(min)

di `h'

di 2.718281828^`h'

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