I'm not sure if this is the problem, but wouldn't it be
post `1' (r(rho) - `y')
and not
post `1' (`r(rho)' - `y') ?
Joseph Coveney
----------------------------------------------------------------------------
Patrick Sturgis wrote:
I am using the following program to obtain the 95% confidence interval
around the difference between 2 average correlations, obtained via the alpha
command.
program define avcorwithin
version 7.0
if "`1'" == "?" {
global S_1 "avcorwithin"
exit
}
tempname y
alpha bef*
scalar `y' = r(rho)
alpha aft*
post `1' (`r(rho)'-`y')
end
When I run this, I do not receive any error messages but receive the
following output, which suggest that no bootstrapping is being done:
. bstrap avcorwithin, reps(2000)
(obs=224)
Bootstrap statistics
Variable Reps Observed Bias Std. Err. [95% Conf. Interval]
avcorwithin 0 . . . . . (N)
. . (P)
. . (BC)
N = normal, P = percentile, BC = bias-corrected
A diagnosis of my problem will be very gratefully received,
Patrick
Dr Patrick Sturgis
Department of Sociology
University of Surrey
Guildford
GU2 7XH
tel: 01483 686 974
*
* 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/
