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

Re: st: bstrap


From   Joseph Coveney <jcoveney@bigplanet.com>
To   Statalist <statalist@hsphsun2.harvard.edu>
Subject   Re: st: bstrap
Date   Fri, 05 Dec 2003 19:21:24 +0900

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/



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