Bookmark and Share

Notice: On March 31, it was announced that Statalist is moving from an email list to a forum. The old list will shut down at the end of May, and its replacement, statalist.org is already up and running.


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

RE: st: saving the confidence interval as a scalar


From   "Monuteaux, Michael" <Michael.Monuteaux@childrens.harvard.edu>
To   "statalist@hsphsun2.harvard.edu" <statalist@hsphsun2.harvard.edu>
Subject   RE: st: saving the confidence interval as a scalar
Date   Fri, 29 Jul 2011 09:03:15 -0400

Steven, this is perfect, thanks so much!

-----Original Message-----
From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Steven Samuels
Sent: Thursday, July 28, 2011 5:56 PM
To: statalist@hsphsun2.harvard.edu
Subject: Re: st: saving the confidence interval as a scalar


Michael, try the following:

*****************************************
sysuse auto, clear
svyset rep78 [pw = turn]
svy, subpop(if foreign==0): total trunk
matrix b = e(b)
scalar tot = b[1,1]
matrix v=e(V)
scalar v =v[1,1]
scalar bound =  ///
sqrt(v)*invttail(e(df_r),0.5*(1-c(level)/100))
scalar  ll = tot-bound
scalar  ul = tot + bound
di  c(level) "% Confidence Interval" _newline ///
"Estimate: " %8.0fc  tot ///
 " ["%7.0fc  ll "," %8.0fc  ul "]"
 ******************************************

On Jul 28, 2011, at 12:12 PM, Monuteaux, Michael wrote:

Hello everyone,

I am estimating various parameters using a survey dataset.  I am running the following command:
	svy: total myvar, subpop(adults)
I am able to save the point estimate as a scalar using the saved matrix, e(b).  But, I would like to also save the upper and lower confidence interval bounds as scalars, as well.  Is it possible to access these values as scalars or a matrix?


*
*   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/

*
*   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/


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