This may help Victor with the -postfile- issue:
http://www.stata-journal.com/article.html?article=pr0036
HTH
Martin
-------- Original-Nachricht --------
> Datum: Thu, 13 Nov 2008 14:56:05 -0000
> Von: "Nick Cox" <[email protected]>
> An: [email protected]
> Betreff: RE: st: Re: Memory
> I still don't understand what you are trying to do. But I can comment on
> your code.
>
> You are looping round 40,000 times and writing a single result to 40,000
> data files. Then you are looping round to put all those 40,000 data
> files in one.
>
> I'd do that directly this way using just one extra file:
>
> clear
> set obs 31
> gen a = .
> tempname out
> postfile `out' t using myresults.dta
> qui forval i = 1/40000 {
> replace a = invnorm(uniform())
> ttest a = 0
> post `out' (r(t))
> }
> postclose `out'
>
> I still doubt 40,000 is anywhere big enough to get an answer.
>
> Nick
> [email protected]
>
> Victor M. Zammit
>
> * a} The data that I have is from generating random samples of whatever
> size,in this case of size 31,from a normally distributed,infinitely
> large,
> population; ie
>
> local i = 1
>
> while `i'<= 40000 {
>
> drop _all
>
> set obs 31
>
> gen a = invnorm(uniform())
>
> qui ttest a = 0
>
> replace a = r(t) in 1
>
> keep in 1
>
> save a`i',replace
>
> local i = `i'+1
>
> }
>
> * I use 40000 due to memory constraint.Appending the a[i]'s together
> gives
> me a variable of 40000 observations ,ie
>
> use a1,clear
>
> local i = 2
>
> while `i'<= 40000 {
>
> append using a`i'.dta
>
> local i = `i'+1
>
> }
>
> save ais40000,replace
>
> * b) From ais40000.dta I get the density <= 1.31, presumably to get the
> density of 90% , <= 1.697 to get the density of 95% etc etc,according to
> the
> official ttable, ie
>
> capture program drop density
>
> program define density
>
> use ais40000,clear
>
> count if a<= `1'
>
> di " density >=" "`1'" " = " r(N)/40000
>
> end
>
> density 1.31
>
> density 1.697
>
> density 2.042
>
> density 2.457
>
> density 2.75
>
> * For smaller degrees of freedom,the discrepancy is much higher.I would
> like
> to know how if it is at all possible to resolve memory constraint .
>
>
> *
> * 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/
--
Martin Weiss
Hackersteigle 3
72076 Tübingen
00497071/793535
00491784597218
Sensationsangebot nur bis 30.11: GMX FreeDSL - Telefonanschluss + DSL
für nur 16,37 Euro/mtl.!* http://dsl.gmx.de/?ac=OM.AD.PD003K11308T4569a
*
* 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/
