Dear R People: In Feller's article, "The Asymptotic Distribution of the Range of Sums of Independent Random Variables", he demonstrates how the range and normalized range of the sum of n independent variables can be derived via Brownian motion. ok. fine. He has a section in which he uses a Bernoulli random variable with options of -1 and 1 to get exact values and then demonstrates his Brown motion approx. He then shows the expected value and variance. I was using R to duplicate the exact process. When I get my values, they are somewhat different than Feller's. For n=6, I'm ok, but for values larger than 6, my function is off slightly: Example Feller Me 4.1523 4.1679 2.0872 2.0850 Is this just due to rounding, or (probably) a mistake in my function? Thanks in advance for your help! Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown 1 Main Street Houston, TX 77002 mailto: hodgess at uhddx01.dt.uh.edu -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._