Gmail Calendar Documents Reader Web more »
Recently Visited Groups | Help | Sign in
Google Groups Home
Message from discussion Somerhing´s not correct

View parsed - Show only message text

Path: g2news1.google.com!news2.google.com!news.glorb.com!tr22g12.aset.psu.edu!news.mathforum.org!not-for-mail
From: "Luis A. Afonso" <lic...@hotmail.com>
Newsgroups: sci.stat.math
Subject: =?UTF-8?Q?Somerhing=C2=B4s_not_correct?=
Date: Wed, 01 Jul 2009 17:41:23 EDT
Organization: The Math Forum
Lines: 46
Message-ID: <17818846.60245.1246484513828.JavaMail.jakarta@nitrogen.mathforum.org>
NNTP-Posting-Host: nitrogen.mathforum.org
Mime-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
X-Trace: support1.mathforum.org 1246484513 14604 144.118.30.135 (1 Jul 2009 21:41:53 GMT)
X-Complaints-To: news@news.mathforum.org
NNTP-Posting-Date: Wed, 1 Jul 2009 21:41:53 +0000 (UTC)

Something´s not correct

San josé State University:
www.applet-magic.com/samplemin.htm
it can be read:
******** If p(x) is the probability density function for a random variable x, let P(x) be the cumulative probability function; i.e.,
P(x) = Intg(-infinity, x).p(z)dz
The probability that the minimum of a sample of size n is x is given by
Prob.(min size n sample) = n[1- P(x)]^(n-1) * p(x)
*********************************************
My comment
Suppose, for example,
Uniform [0, 1] Distribution, P(x) = x, p(x) =1
Prob.(min size n sample) = 21[1- x]^20 * 1
If x=0.05 then results 7.528,
This is, no doubt, a remarkable thing (!!!). The Applet is so magic that it not available yet. Please, remove it!
The writer did intend to say:
____Prob (min <=x) = 1- (1-F(x))^n
____1- (1-0.05)^21=0.659438______
Simulation (data 8E6 items) gave 0.660, 0.660.
Prob (absolute deviation > 0.0005 | N=8E6) <=
2*EXP(-2*8E6*0.25*1E-6) = 2*EXP(-4) = 0.037= 3.7%
(According to Dvorestky - Kiefer - Wolfowitz inequality, states).
REM "ckontas"
        CLS
        DEFDBL A-Z
        DIM x(22)
        INPUT " x=  "; x0
        p = 1 - (1 - x0) ^ 21
        PRINT USING "#.### "; p
        all = 8000000
        RANDOMIZE TIMER
        FOR j = 1 TO all
        LOCATE 4, 4: PRINT USING "########"; all - j
        minor = 8
        FOR i = 1 TO 21
        x = RND
        IF x < minor THEN minor = x
        NEXT i
        IF minor > x0 THEN GOTO 10
        okk = okk + 1
10      NEXT j: LOCATE 10, 10
        PRINT USING "#.### "; okk / all
        END

Luis A. Afonso

Create a group - Google Groups - Google Home - Terms of Service - Privacy Policy
©2009 Google