R function which finds confidence interval for binomial variance

I need to construct confidence intervals for the binomial variance. 
This is the usual estimate

	v = x*(n-x)/n

or its unbiased counterpart

	v' = x*(n-x)/(n-1)

where x = binomial number of successes observed in n Bernoulli trials 
from proportion p.

The usual X^2 method for variance confidence intervals will not work, 
because of the strong non-normal character of the sampling 
distribution for v (or v').

Does anyone know of an R package with R function that computes a 
reasonable confidence interval for v or v'?

Thanks.
================================================================
Robert A. LaBudde, PhD, PAS, Dpl. ACAFS  e-mail: ral at lcfltd.com
Least Cost Formulations, Ltd.            URL: http://lcfltd.com/
824 Timberlake Drive                     Tel: 757-467-0954
Virginia Beach, VA 23464-3239            Fax: 757-467-2947

"Vere scire est per causas scire"