quietly {
  /*
     Calculate confidence interval usinf FIsher's z transformation.
  */
  local r=.66                                   /* Insert correlation */
  local n=30                                    /* Insert sample size */
  local rstar =.5*ln((1+`r')/(1-`r'))           /* Convert r to z     */
  local se = (`n'-3)^-.5                        /* Calculate the se   */
  local z = invnorm(.975)                       /* How big is the CI? */
  local lb = `rstar' - `z'*`se'                 /* lower bound in z   */
  local ub = `rstar' + `z'*`se'                 /* upper bound in z   */
  local lb = (exp(2*`lb')-1)/(exp(2*`lb')+1)    /* lower bound in r   */
  local ub = (exp(2*`ub')-1)/(exp(2*`ub')+1)    /* upper bound in r   */

  /* Print results to screen */
  #delimit ;
  noisily di _n _col(2) "Confidence interval for correlation based on Fisher's transformation:"
             _n _col(3) "Correlation = " %4.3f `r' " with " `n' " observations" 
             " (95% CI: " %4.3f `lb' " to " %4.3f `ub' in gr ")"
  ;
  #delimit cr
}