Oops, the line:
res[i,1] = sqrt(variance(data[|k0,1\k1,1|]))
should have been:
res[i,1] = sqrt(variance(data[|k0,2\k1,2|]))
That is, I was computing the local standard
deviation of X (the first column of data) rather than
Y (the second column of data).
Sorry,
Maarten
-----------------------------------------
Maarten L. Buis
Department of Social Research Methodology
Vrije Universiteit Amsterdam
Boelelaan 1081
1081 HV Amsterdam
The Netherlands
visiting address:
Buitenveldertselaan 3 (Metropolitan), room N515
+31 20 5986715
http://home.fsw.vu.nl/m.buis/
-----------------------------------------
*
* 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/
