Statalist


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

Re: RE: st: Re: Memory


From   "Martin Weiss" <[email protected]>
To   [email protected]
Subject   Re: RE: st: Re: Memory
Date   Thu, 13 Nov 2008 16:16:31 +0100

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/



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