mersenneforum.org  

Go Back   mersenneforum.org > New To GIMPS? Start Here! > Information & Answers

Reply
 
Thread Tools
Old 2011-07-03, 17:53   #12
davieddy
 
davieddy's Avatar
 
"Lucan"
Dec 2006
England

194A16 Posts
Default Chance of finding a prime via DC v first LL

People seem to be over pessimistic on this question.

The ball park error rate of LLests is 2%.

But DCs are performed on exponents half the size, which
means they take 4 times less GHz-Days, and are twice as
likely to be prime.

So by DCing, your chance of finding a MP per year is 16%
of your chance doing first time LLtests: not as dramatic a
difference as is often portrayed.

David

PS and as Christenson has mentioned elsewhere,
6 x miniscule = miniscule ;-)

Last fiddled with by davieddy on 2011-07-03 at 18:33
davieddy is offline   Reply With Quote
Old 2011-07-04, 00:35   #13
cheesehead
 
cheesehead's Avatar
 
"Richard B. Woods"
Aug 2002
Wisconsin USA

22×3×641 Posts
Default

Quote:
Originally Posted by davieddy View Post
People seem to be over pessimistic on this question.

< snip >

But DCs are <snip > are twice as
likely to be prime.

So by DCing, your chance of finding a MP per year is 16%
of your chance doing first time LLtests:
The most common error in probability calculations is failing to include the effects and consequences of knowledge.

Your calculation seems to assume that the probability of a DC assigned exponent's being that of a Mersenne prime is just as though it were the corresponding probability for first-time tests -- without using any knowledge gained from the first-time tests. That assumption's not valid!

In actuality, GIMPS does not assign DCs for exponents which are already known to be Mersenne primes!

Thus, a correct calculation would have to take into account that GIMPS is deliberately and knowingly _not_ assigning DCs for exponents of known Mersenne primes.

Quote:
not as dramatic a difference as is often portrayed.
... only because you neglect to account for the effects of prior knowledge!
cheesehead is offline   Reply With Quote
Old 2011-07-04, 00:52   #14
davieddy
 
davieddy's Avatar
 
"Lucan"
Dec 2006
England

145128 Posts
Default

I deliberately refrained from including the necessary pedantry.

Given that a first time LLtest is erroneous, (2% say), the probability
of the a DC being prime is doubled while the test time is quartered.

What earthly difference does "already found MPs" make to this
sound and simple conditional probability argument?

David
davieddy is offline   Reply With Quote
Old 2011-07-04, 01:06   #15
cheesehead
 
cheesehead's Avatar
 
"Richard B. Woods"
Aug 2002
Wisconsin USA

1E0C16 Posts
Default

Quote:
Originally Posted by davieddy View Post
I deliberately refrained from including the necessary pedantry.
... which allowed misinterpretation, because your actual statement was not your intended statement. :-)

Quote:
What earthly difference does "already found MPs" make to this
sound and simple conditional probability argument?
None, but "this" argument is not the same as the "So by DCing, your chance of finding a MP per year" statement to which I previously replied.

Last fiddled with by cheesehead on 2011-07-04 at 01:13
cheesehead is offline   Reply With Quote
Old 2011-07-04, 02:34   #16
davieddy
 
davieddy's Avatar
 
"Lucan"
Dec 2006
England

2×3×13×83 Posts
Default

Quote:
Originally Posted by cheesehead View Post
... which allowed misinterpretation, because your actual statement was not your intended statement. :-)

None, but "this" argument is not the same as the "So by DCing, your chance of finding a MP per year" statement to which I previously replied.
The statement stands.
davieddy is offline   Reply With Quote
Old 2011-07-04, 03:54   #17
wblipp
 
wblipp's Avatar
 
"William"
May 2003
New Haven

2·7·132 Posts
Default

Quote:
Originally Posted by cheesehead View Post
... only because you neglect to account for the effects of prior knowledge!
I'm skimming, so I may have missed something important.

But didn't David start with the prior knowledge that only 2% of the first tests are faulty, so if all things were equal, you would be 2% as likely to find a prime.

But these tests, if done correctly, are twice as likely to be prime, so in the same number of tests you are 4% as likely to find a prime.

But the tests are four times faster, so in the same time you are 16% as likely to find a prime.

I don't see what additional prior knowledge would be relevant to the calculation. In particular, knowledge of known Mersenne Primes is not knowledge about faulty tests - those tests are known to be correct, so not relevant to this analysis.

William
wblipp is offline   Reply With Quote
Old 2011-07-04, 05:53   #18
Brian-E
 
Brian-E's Avatar
 
"Brian"
Jul 2007
The Netherlands

7×467 Posts
Default

Quote:
Originally Posted by cheesehead View Post
Thus, a correct calculation would have to take into account that GIMPS is deliberately and knowingly _not_ assigning DCs for exponents of known Mersenne primes.
My own understanding, and it is likely to be faulty since I have studied no probability theory in more than 25 years, is that the above statement is quite correct in that it identifies a biasing factor. However surely the bias in this case is so tiny (there are so few identified Mersenne Primes) that it will have no effect on the very rough 16% estimate which David arrives at.
Brian-E is offline   Reply With Quote
Old 2011-07-04, 09:39   #19
davieddy
 
davieddy's Avatar
 
"Lucan"
Dec 2006
England

2·3·13·83 Posts
Default

Quote:
Originally Posted by Brian-E View Post
the very rough 16% estimate which David arrives at.
THX Brian (and William).

Not all that rough though. OK someone can work out a more accurate
current "error rate" (E%) by the DCs whose residues don't match the
first LL, but I made it clear that "2%" was a ball park guess at E%.

Doubling the exponent increases the GHz-Days by slightly more than 4 times.
The probability of being prime is complicated by "how far TFed" etc, but it involves 1/exponent.

If I said "8E%" instead of "16%" I would not be too far out.

David

Last fiddled with by davieddy on 2011-07-04 at 09:45
davieddy is offline   Reply With Quote
Old 2011-07-04, 17:18   #20
LiquidNitrogen
 
LiquidNitrogen's Avatar
 
Jun 2011
Henlopen Acres, Delaware

100001012 Posts
Default

Blah blah blah blah blah, you guys should run for office.
LiquidNitrogen is offline   Reply With Quote
Old 2011-07-05, 10:39   #21
davieddy
 
davieddy's Avatar
 
"Lucan"
Dec 2006
England

2×3×13×83 Posts
Default Estimating the "error rate"

Quote:
Originally Posted by cheesehead View Post
The most common error in probability calculations is failing to include the effects and consequences of knowledge.

Your calculation seems to assume that the probability of a DC assigned exponent's being that of a Mersenne prime is just as though it were the corresponding probability for first-time tests -- without using any knowledge gained from the first-time tests. That assumption's not valid!

In actuality, GIMPS does not assign DCs for exponents which are already known to be Mersenne primes!

Thus, a correct calculation would have to take into account that GIMPS is deliberately and knowingly _not_ assigning DCs for exponents of known Mersenne primes.

... only because you neglect to account for the effects of prior knowledge!
I think all your qualms boil down to obtaining an unbiassed estimate
of the probability of the non-zero residue of an unverified LL being wrong.

No doubt a determined quibbler will find this too simple,
but I would take a suitable* large sample of "verified composite via
two matching residues", and find the % where the first test proved
to be wrong.

This totally bypasses your "known prime" "difficulty".

David

* I am "covering my back" this time by allowing
you to decide what "suitable" means!

Last fiddled with by davieddy on 2011-07-05 at 11:11
davieddy is offline   Reply With Quote
Old 2011-07-05, 11:27   #22
davieddy
 
davieddy's Avatar
 
"Lucan"
Dec 2006
England

11001010010102 Posts
Default

Quote:
Originally Posted by LiquidNitrogen View Post
Blah blah blah blah blah, you guys should run for office.
More than you needed to know ATM?

"Us guys" (comparative veterans) still enjoy "getting it straight".


David
davieddy is offline   Reply With Quote
Reply



All times are UTC. The time now is 22:28.


Fri Aug 6 22:28:21 UTC 2021 up 14 days, 16:57, 1 user, load averages: 2.94, 3.19, 3.17

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2021, Jelsoft Enterprises Ltd.

This forum has received and complied with 0 (zero) government requests for information.

Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation.
A copy of the license is included in the FAQ.