Main Page: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
No edit summary
No edit summary
Line 1: Line 1:
[[Image:Helicoid.svg|right|thumb|350px|A helicoid with α=1, -1≤ρ≤1 and -π≤θ≤π.]]
In [[statistics]], the '''Neyman-Pearson [[lemma (mathematics)|lemma]]''', named after [[Jerzy Neyman]] and [[Egon Pearson]], states that when performing a [[statistical hypothesis testing|hypothesis test]] between two point  hypotheses ''H''<sub>0</sub>:&nbsp;''θ''&nbsp;=&nbsp;''θ''<sub>0</sub> and ''H''<sub>1</sub>:&nbsp;''θ''&nbsp;=&nbsp;''θ''<sub>1</sub>, then the [[likelihood-ratio test]] which rejects ''H''<sub>0</sub> in favour of ''H''<sub>1</sub> when


The '''helicoid''', after the [[Plane (geometry)|plane]] and the [[catenoid]], is the third [[minimal surface]] to be known. It was first discovered by [[Jean Baptiste Meusnier]] in 1776. Its [[Nomenclature|name]] derives from its similarity to the [[helix]]: for every [[Point (geometry)|point]] on the helicoid there is a helix contained in the helicoid which passes through that point. Since it is considered that the planar range extends through negative and positive infinity, close observation shows the appearance of two parallel or mirror planes in the sense that if the slope of one plane is traced, the co-plane can be seen to be bypassed or skipped, though in actuality the co-plane is also traced from the opposite perspective.
:<math>\Lambda(x)=\frac{ L( \theta _{0} \mid x)}{ L (\theta _{1} \mid x)} \leq \eta</math>
where
:<math>P(\Lambda(X)\leq \eta|H_0)=\alpha </math>


The helicoid is also a [[ruled surface]] (and a [[right conoid]]), meaning that it is a trace of a line. Alternatively, for any point on the surface, there is a line on the surface passing through it.  Indeed, [[Eugène Charles Catalan|Catalan]] proved in 1842 that the helicoid and the plane were the only ruled minimal surfaces.<ref>Elements of the Geometry and Topology of Minimal Surfaces in Three-dimensional Space
is the '''most powerful test''' of [[Type I and type II errors|size ''&alpha;'']] for a threshold η. If the test is most powerful for all <math>\theta_1 \in \Theta_1</math>, it is said to be [[uniformly most powerful]] (UMP) for alternatives in the set <math>\Theta_1 \, </math>.
By [[A. T. Fomenko]], A. A. Tuzhilin
Contributor A. A. Tuzhilin
Published by AMS Bookstore, 1991
ISBN 0-8218-4552-7, ISBN 978-0-8218-4552-3, p.33</ref>


The helicoid and the [[catenoid]] are parts of a family of helicoid-catenoid minimal surfaces.
In practice, the [[likelihood ratio]] is often used directly to construct tests &mdash; see [[Likelihood-ratio test]]. However it can also be used to suggest particular test-statistics that might be of interest or to suggest simplified tests &mdash; for this one considers algebraic manipulation of the ratio to see if there are key statistics in it related to the size of the ratio (i.e. whether a large statistic corresponds to a small ratio or to a large one).


The helicoid is shaped like [[Archimedes' screw]], but extends infinitely in all directions.  It can be described by the following [[parametric equation]]s in [[Cartesian coordinates]]:
==Proof==
:<math> x = \rho \cos (\alpha \theta), \ </math>
Define the rejection region of the null hypothesis for the NP test as
:<math> y = \rho \sin (\alpha \theta), \ </math>
:<math>R_{NP}=\left\{ x: \frac{L(\theta_{0}|x)}{L(\theta_{1}|x)} \leq \eta\right\} .</math>
:<math> z = \theta, \ </math>
where ''ρ'' and ''θ'' range from negative [[infinity]] to [[positive number|positive]] infinity, while ''α'' is a constant.  If ''α'' is positive then the helicoid is right-handed as shown in the figure; if negative then left-handed.


The helicoid has [[principal curvature]]s <math>\pm 1/(1+ \rho ^2) \ </math>. The sum of these quantities gives the [[mean curvature]] (zero since the helicoid is a [[minimal surface]]) and the product gives the [[Gaussian curvature]].
Any other test will have a different rejection region that we define as <math>R_A</math>. Furthermore, define the probability of the data falling in region R, given parameter <math>\theta</math> as


The helicoid is [[homeomorphism|homeomorphic]] to the plane <math> \mathbb{R}^2 </math>.  To see this, let alpha decrease [[continuous function|continuous]]ly from its given value down to [[0 (number)|zero]].  Each intermediate value of ''α'' will describe a different helicoid, until ''α = 0'' is reached and the helicoid becomes a vertical [[plane (mathematics)|plane]].
:<math>P(R,\theta)=\int_R L(\theta|x)\, dx, </math>


Conversely, a plane can be turned into a helicoid by choosing a line, or ''axis'', on the plane then twisting the plane around that axis.
For both tests to have size <math>\alpha</math>, it must be true that


For example, if we take h as the maxium value at z and R the radium, the area of the surface is π[R√(R²+h²)+h2*ln(R +√(R²+h²)/h)]
:<math>\alpha= P(R_{NP}, \theta_0)=P(R_A, \theta_0) \,.</math>


==Helicoid and catenoid==
It will be useful to break these down into integrals over distinct regions:
[[File:Helicatenoid.gif|thumb|256px|Animation showing the transformation of a helicoid into a catenoid.]]
 
The helicoid and the [[catenoid]] are locally isometric surfaces, see discussion there.
:<math>P(R_{NP},\theta) = P(R_{NP} \cap R_A, \theta) + P(R_{NP} \cap R_A^c, \theta),</math>
 
and
 
:<math>P(R_A,\theta) = P(R_{NP} \cap R_A, \theta) + P(R_{NP}^c \cap R_A, \theta).</math>
 
Setting <math>\theta=\theta_0</math> and equating the above two expression yields that
 
:<math>P(R_{NP} \cap R_A^c, \theta_0) =  P(R_{NP}^c \cap R_A, \theta_0).</math>
 
Comparing the powers of the two tests, <math>P(R_{NP},\theta_1)</math> and <math>P(R_A,\theta_1)</math>, one can see that
 
:<math>P(R_{NP},\theta_1) \geq P(R_A,\theta_1) \iff
P(R_{NP} \cap R_A^c, \theta_1) \geq P(R_{NP}^c \cap R_A, \theta_1). </math>
 
Now by the definition of <math>R_{NP}</math>,
 
:<math> P(R_{NP} \cap R_A^c, \theta_1)= \int_{R_{NP}\cap R_A^c} L(\theta_{1}|x)\,dx \geq \frac{1}{\eta} \int_{R_{NP}\cap R_A^c} L(\theta_0|x)\,dx = \frac{1}{\eta}P(R_{NP} \cap R_A^c, \theta_0)</math>
:<math> = \frac{1}{\eta}P(R_{NP}^c \cap R_A, \theta_0) = \frac{1}{\eta}\int_{R_{NP}^c \cap R_A} L(\theta_{0}|x)\,dx \geq \int_{R_{NP}^c\cap R_A} L(\theta_{1}|x)dx  = P(R_{NP}^c \cap R_A, \theta_1).</math>
 
Hence the [[Inequality (mathematics)|inequality]] holds.
 
==Example==
 
Let <math>X_1,\dots,X_n</math> be a random sample from the <math>\mathcal{N}(\mu,\sigma^2)</math> distribution where the mean <math>\mu</math> is known, and suppose that we wish to test for <math>H_0:\sigma^2=\sigma_0^2</math> against <math>H_1:\sigma^2=\sigma_1^2</math>. The likelihood for this set of [[normally distributed]] data is
 
:<math>L\left(\sigma^2;\mathbf{x}\right)\propto \left(\sigma^2\right)^{-n/2} \exp\left\{-\frac{\sum_{i=1}^n \left(x_i-\mu\right)^2}{2\sigma^2}\right\}.</math>
 
We can compute the [[likelihood ratio]] to find the key statistic in this test and its effect on the test's outcome:
 
:<math>\Lambda(\mathbf{x}) = \frac{L\left(\sigma_1^2;\mathbf{x}\right)}{L\left(\sigma_0^2;\mathbf{x}\right)} =
\left(\frac{\sigma_1^2}{\sigma_0^2}\right)^{-n/2}\exp\left\{-\frac{1}{2}(\sigma_1^{-2}-\sigma_0^{-2})\sum_{i=1}^n \left(x_i-\mu\right)^2\right\}.</math>
 
This ratio only depends on the data through <math>\sum_{i=1}^n \left(x_i-\mu\right)^2</math>. Therefore, by the Neyman–Pearson lemma, the most [[Statistical power|powerful]] test of this type of [[Statistical hypothesis testing|hypothesis]] for this data will depend only on <math>\sum_{i=1}^n \left(x_i-\mu\right)^2</math>. Also, by inspection, we can see that if <math>\sigma_1^2>\sigma_0^2</math>, then <math>\Lambda(\mathbf{x})</math> is an [[increasing function]] of <math>\sum_{i=1}^n \left(x_i-\mu\right)^2</math>. So we should reject <math>H_0</math> if <math>\sum_{i=1}^n \left(x_i-\mu\right)^2</math> is sufficiently large. The rejection threshold depends on the [[Type I and type II errors|size]]  of the test. In this example, the test statistic can be shown to be a scaled Chi-square distributed random variable and an exact critical value can be obtained.


==See also==
==See also==
*[[Dini's surface]]
* [[Statistical power]]
*[[Right conoid]]
*[[Ruled surface]]


==Notes==
==References==
<references />
* {{cite journal
|title=On the Problem of the Most Efficient Tests of Statistical Hypotheses
|author=[[Jerzy Neyman]], [[Egon Pearson]]
|journal=[[Philosophical Transactions of the Royal Society of London]]. Series A, Containing Papers of a Mathematical or Physical Character
|volume=231
|year=1933
|pages=289–337
|doi=10.1098/rsta.1933.0009
|jstor=91247
}}
*[http://cnx.org/content/m11548/latest/ cnx.org: Neyman-Pearson criterion]


==External links==
==External links==
*[http://chamicewicz.com/p5/helicoid/ Interactive 3D Helicoid plotter using Processing (with code)]
* Cosma Shalizi, a professor of statistics at Carnegie Mellon University, gives an intuitive derivation of the Neyman-Pearson Lemma [http://cscs.umich.edu/~crshalizi/weblog/630.html using ideas from economics]


[[Category:Geometric shapes]]
{{DEFAULTSORT:Neyman–Pearson Lemma}}
[[Category:Minimal surfaces]]
[[Category:Statistical theorems]]
[[Category:Surfaces]]
[[Category:Statistical tests]]
[[Category:Articles containing proofs]]


[[ar:سطح حلزوني]]
[[ca:Lema de Neyman-Pearson]]
[[de:Minimalfläche#Die Wendelfläche]]
[[de:Neyman-Pearson-Lemma]]
[[et:Helikoid]]
[[es:Lema de Neyman-Pearson]]
[[fr:Hélicoïde]]
[[fr:Lemme de Neyman-Pearson]]
[[it:Elicoide]]
[[it:Lemma fondamentale di Neyman-Pearson]]
[[hu:Csavarfelület]]
[[he:הלמה של ניימן-פירסון]]
[[nl:Helicoïde]]
[[ja:ネイマン・ピアソンの補題]]
[[pl:Helikoida]]
[[ru:Геликоид]]
[[sl:Helikoid]]
[[uk:Гелікоїд]]

Revision as of 15:17, 12 August 2014

In statistics, the Neyman-Pearson lemma, named after Jerzy Neyman and Egon Pearson, states that when performing a hypothesis test between two point hypotheses H0θ = θ0 and H1θ = θ1, then the likelihood-ratio test which rejects H0 in favour of H1 when

Λ(x)=L(θ0x)L(θ1x)η

where

P(Λ(X)η|H0)=α

is the most powerful test of size α for a threshold η. If the test is most powerful for all θ1Θ1, it is said to be uniformly most powerful (UMP) for alternatives in the set Θ1.

In practice, the likelihood ratio is often used directly to construct tests — see Likelihood-ratio test. However it can also be used to suggest particular test-statistics that might be of interest or to suggest simplified tests — for this one considers algebraic manipulation of the ratio to see if there are key statistics in it related to the size of the ratio (i.e. whether a large statistic corresponds to a small ratio or to a large one).

Proof

Define the rejection region of the null hypothesis for the NP test as

RNP={x:L(θ0|x)L(θ1|x)η}.

Any other test will have a different rejection region that we define as RA. Furthermore, define the probability of the data falling in region R, given parameter θ as

P(R,θ)=RL(θ|x)dx,

For both tests to have size α, it must be true that

α=P(RNP,θ0)=P(RA,θ0).

It will be useful to break these down into integrals over distinct regions:

P(RNP,θ)=P(RNPRA,θ)+P(RNPRAc,θ),

and

P(RA,θ)=P(RNPRA,θ)+P(RNPcRA,θ).

Setting θ=θ0 and equating the above two expression yields that

P(RNPRAc,θ0)=P(RNPcRA,θ0).

Comparing the powers of the two tests, P(RNP,θ1) and P(RA,θ1), one can see that

P(RNP,θ1)P(RA,θ1)P(RNPRAc,θ1)P(RNPcRA,θ1).

Now by the definition of RNP,

P(RNPRAc,θ1)=RNPRAcL(θ1|x)dx1ηRNPRAcL(θ0|x)dx=1ηP(RNPRAc,θ0)
=1ηP(RNPcRA,θ0)=1ηRNPcRAL(θ0|x)dxRNPcRAL(θ1|x)dx=P(RNPcRA,θ1).

Hence the inequality holds.

Example

Let X1,,Xn be a random sample from the 𝒩(μ,σ2) distribution where the mean μ is known, and suppose that we wish to test for H0:σ2=σ02 against H1:σ2=σ12. The likelihood for this set of normally distributed data is

L(σ2;x)(σ2)n/2exp{i=1n(xiμ)22σ2}.

We can compute the likelihood ratio to find the key statistic in this test and its effect on the test's outcome:

Λ(x)=L(σ12;x)L(σ02;x)=(σ12σ02)n/2exp{12(σ12σ02)i=1n(xiμ)2}.

This ratio only depends on the data through i=1n(xiμ)2. Therefore, by the Neyman–Pearson lemma, the most powerful test of this type of hypothesis for this data will depend only on i=1n(xiμ)2. Also, by inspection, we can see that if σ12>σ02, then Λ(x) is an increasing function of i=1n(xiμ)2. So we should reject H0 if i=1n(xiμ)2 is sufficiently large. The rejection threshold depends on the size of the test. In this example, the test statistic can be shown to be a scaled Chi-square distributed random variable and an exact critical value can be obtained.

See also

References

  • One of the biggest reasons investing in a Singapore new launch is an effective things is as a result of it is doable to be lent massive quantities of money at very low interest rates that you should utilize to purchase it. Then, if property values continue to go up, then you'll get a really high return on funding (ROI). Simply make sure you purchase one of the higher properties, reminiscent of the ones at Fernvale the Riverbank or any Singapore landed property Get Earnings by means of Renting

    In its statement, the singapore property listing - website link, government claimed that the majority citizens buying their first residence won't be hurt by the new measures. Some concessions can even be prolonged to chose teams of consumers, similar to married couples with a minimum of one Singaporean partner who are purchasing their second property so long as they intend to promote their first residential property. Lower the LTV limit on housing loans granted by monetary establishments regulated by MAS from 70% to 60% for property purchasers who are individuals with a number of outstanding housing loans on the time of the brand new housing purchase. Singapore Property Measures - 30 August 2010 The most popular seek for the number of bedrooms in Singapore is 4, followed by 2 and three. Lush Acres EC @ Sengkang

    Discover out more about real estate funding in the area, together with info on international funding incentives and property possession. Many Singaporeans have been investing in property across the causeway in recent years, attracted by comparatively low prices. However, those who need to exit their investments quickly are likely to face significant challenges when trying to sell their property – and could finally be stuck with a property they can't sell. Career improvement programmes, in-house valuation, auctions and administrative help, venture advertising and marketing, skilled talks and traisning are continuously planned for the sales associates to help them obtain better outcomes for his or her shoppers while at Knight Frank Singapore. No change Present Rules

    Extending the tax exemption would help. The exemption, which may be as a lot as $2 million per family, covers individuals who negotiate a principal reduction on their existing mortgage, sell their house short (i.e., for lower than the excellent loans), or take part in a foreclosure course of. An extension of theexemption would seem like a common-sense means to assist stabilize the housing market, but the political turmoil around the fiscal-cliff negotiations means widespread sense could not win out. Home Minority Chief Nancy Pelosi (D-Calif.) believes that the mortgage relief provision will be on the table during the grand-cut price talks, in response to communications director Nadeam Elshami. Buying or promoting of blue mild bulbs is unlawful.

    A vendor's stamp duty has been launched on industrial property for the primary time, at rates ranging from 5 per cent to 15 per cent. The Authorities might be trying to reassure the market that they aren't in opposition to foreigners and PRs investing in Singapore's property market. They imposed these measures because of extenuating components available in the market." The sale of new dual-key EC models will even be restricted to multi-generational households only. The models have two separate entrances, permitting grandparents, for example, to dwell separately. The vendor's stamp obligation takes effect right this moment and applies to industrial property and plots which might be offered inside three years of the date of buy. JLL named Best Performing Property Brand for second year running

    The data offered is for normal info purposes only and isn't supposed to be personalised investment or monetary advice. Motley Fool Singapore contributor Stanley Lim would not personal shares in any corporations talked about. Singapore private home costs increased by 1.eight% within the fourth quarter of 2012, up from 0.6% within the earlier quarter. Resale prices of government-built HDB residences which are usually bought by Singaporeans, elevated by 2.5%, quarter on quarter, the quickest acquire in five quarters. And industrial property, prices are actually double the levels of three years ago. No withholding tax in the event you sell your property. All your local information regarding vital HDB policies, condominium launches, land growth, commercial property and more

    There are various methods to go about discovering the precise property. Some local newspapers (together with the Straits Instances ) have categorised property sections and many local property brokers have websites. Now there are some specifics to consider when buying a 'new launch' rental. Intended use of the unit Every sale begins with 10 p.c low cost for finish of season sale; changes to 20 % discount storewide; follows by additional reduction of fiftyand ends with last discount of 70 % or extra. Typically there is even a warehouse sale or transferring out sale with huge mark-down of costs for stock clearance. Deborah Regulation from Expat Realtor shares her property market update, plus prime rental residences and houses at the moment available to lease Esparina EC @ Sengkang
  • cnx.org: Neyman-Pearson criterion

External links

  • Cosma Shalizi, a professor of statistics at Carnegie Mellon University, gives an intuitive derivation of the Neyman-Pearson Lemma using ideas from economics

ca:Lema de Neyman-Pearson de:Neyman-Pearson-Lemma es:Lema de Neyman-Pearson fr:Lemme de Neyman-Pearson it:Lemma fondamentale di Neyman-Pearson he:הלמה של ניימן-פירסון ja:ネイマン・ピアソンの補題