Stata The Stata listserver
[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

st: convolution-trial


From   "Andreas Aschbacher" <aa_surf@gmx.at>
To   statalist@hsphsun2.harvard.edu
Subject   st: convolution-trial
Date   Wed, 24 Mar 2004 11:30:02 +0100 (MET)

Dear fellows !

/first thanks to Neil from yesterday/

please consider the following as "Pseudocode"

clear
set obs 100
gen d =.
gen d1 =.
forvalues i = 1(1)100 {
local i = `i' * 0.1     <-  0.1 *1,2,3,4,..... *-multplikation-operator
range x 1(1)i          <-   range until i-value                       
gen double y=exp(-(`i'-x)^2/2)*0.1/(x + 0.1)    <- y as function of i-value
integ y x,gen(Sy)    <- integrate until i-value
*save new y- and Sy- value in two columns :
replace d[_n] = y   
replace d1[_n] = Sy
drop Sy y
}

please any advice that it works correctly.
      sincerely aschbacher andreas

-- 
+++ NEU bei GMX und erstmalig in Deutschland: TÜV-geprüfter Virenschutz +++
100% Virenerkennung nach Wildlist. Infos: http://www.gmx.net/virenschutz

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



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