Portal:Physics/Selected article/Week 17, 2007: Difference between revisions
en>SmackBot m Delink dates (WP:MOSUNLINKDATES) |
en>FrescoBot m Bot: fixing section wikilinks and minor changes |
||
| Line 1: | Line 1: | ||
== | {{Probability distribution | | ||
name =normal-gamma| | |||
type =density| | |||
pdf_image =| | |||
cdf_image =| | |||
parameters =<math>\mu\,</math> [[location parameter|location]] ([[real number|real]])<br /><math>\lambda > 0\,</math> (real)<br /><math>\alpha \ge 1\,</math> (real)<br /><math>\beta \ge 0\,</math> (real)| | |||
support =<math>x \in (-\infty, \infty)\,\!, \; \tau \in (0,\infty)</math>| | |||
pdf =<math>f(x,\tau|\mu,\lambda,\alpha,\beta) = \frac{\beta^\alpha \sqrt{\lambda}}{\Gamma(\alpha)\sqrt{2\pi}} \, \tau^{\alpha-\frac{1}{2}}\,e^{-\beta\tau}\,e^{ -\frac{ \lambda \tau (x- \mu)^2}{2}}</math>| | |||
cdf =| | |||
mean =<ref name=BS434>Bernardo & Smith (1993, p.434)</ref> <math>\operatorname{E}(X)=\mu\,\! ,\quad \operatorname{E}(\Tau)= \alpha \beta^{-1}</math> | | |||
median = <!-- <math>\mu\,</math> --> | | |||
mode = <math>\left(\mu, \frac{\alpha - \frac12}{\beta}\right)</math>| | |||
variance =<ref name=BS434/> <math>\operatorname{var}(X)= \frac{\beta}{\lambda (\alpha-1)} ,\quad | |||
\operatorname{var}(\Tau)=\alpha \beta^{-2} </math> | | |||
skewness =| | |||
kurtosis =| | |||
entropy =| | |||
mgf =| | |||
char =| | |||
}} | |||
In [[probability theory]] and [[statistics]], the '''normal-gamma distribution''' (or '''Gaussian-gamma distribution''') is a bivariate four-parameter family of continuous [[probability distribution]]s. It is the [[conjugate prior]] of a [[normal distribution]] with unknown [[mean]] and [[Precision (statistics)|precision]].<ref>Bernardo & Smith (1993, pages 136, 268, 434)</ref> | |||
==Definition== | |||
For a pair of [[random variable]], (''X'',''T''), suppose that the [[conditional distribution]] of ''X'' given ''T'' is given by | |||
:<math> X|T \sim N(\mu,1 /(\lambda T)) \,\! , </math> | |||
meaning that the condition distribution is a [[normal distribution]] with [[mean]] <math> \mu</math> and [[precision (statistics)|precision]] <math> \lambda T </math> — equivalently, with [[variance]] <math> 1 / (\lambda T) . </math> | |||
Suppose also that the marginal distribution of ''T'' is given by | |||
:<math>T |\alpha, \beta \sim \mathrm{Gamma}(\alpha,\beta) \! ,</math> | |||
where this means that ''T'' has a [[gamma distribution]]. Here λ, α and β are parameters of the joint distribution. | |||
Then (''X'',''T'') has a normal-gamma distribution, and this is denoted by | |||
:<math> (X,T) \sim \mathrm{NormalGamma}(\mu,\lambda,\alpha,\beta) \! . | |||
</math> | |||
==Properties== | |||
===Probability density function=== | |||
The joint [[probability density function]] of (''X'',''T'') is{{cn|date=April 2013}} | |||
: <math>f(x,\tau|\mu,\lambda,\alpha,\beta) = \frac{\beta^\alpha \sqrt{\lambda}}{\Gamma(\alpha)\sqrt{2\pi}} \, \tau^{\alpha-\frac{1}{2}}\,e^{-\beta\tau}\,e^{ -\frac{ \lambda \tau (x- \mu)^2}{2}}</math> | |||
===Marginal distributions=== | |||
By construction, the [[marginal distribution]] over <math>\tau</math> is a [[gamma distribution]], and the [[conditional distribution]] over <math>x</math> given <math>\tau</math> is a [[Gaussian distribution]]. The [[marginal distribution]] over <math>x</math> is a three-parameter non-standardized [[Student's t-distribution]] with parameters <math>(\nu, \mu, \sigma^2)=(2\alpha, \mu, \beta/(\lambda\alpha))</math>.{{cn|date=April 2013}} | |||
===Exponential family=== | |||
The normal-gamma distribution is a four-parameter [[exponential family]] with [[natural parameters]] <math>\alpha-1/2, -\beta-\lambda\mu^2/2, \lambda\mu, -\lambda/2</math> and [[natural statistics]] <math>\ln\tau, \tau, \tau x, \tau x^2</math>.{{cn|date=April 2013}} | |||
===Moments of the natural statistics=== | |||
The following moments can be easily computed using the [[exponential family#Moment generating function of the sufficient statistic|moment generating function of the sufficient statistic]]:{{cn|date=April 2013}} | |||
:<math>\operatorname{E}(\ln T)=\psi\left(\alpha\right) - \ln\beta</math>, where <math>\psi\left(\alpha\right)</math> is the [[digamma function]],<br/> | |||
:<math>\operatorname{E}(T)=\frac{\alpha}{\beta}</math>,<br/> | |||
:<math>\operatorname{E}(TX)=\mu \frac{\alpha}{\beta}</math>,<br/> | |||
:<math>\operatorname{E}(TX^2)=\frac{1}{\lambda} + \mu^2 \frac{\alpha}{\beta}</math>. | |||
===Scaling=== | |||
If <math> (X,T) \sim \mathrm{NormalGamma}(\mu,\lambda,\alpha,\beta), </math> then for any ''b'' > 0, (''bX'',''bT'') is distributed as{{cn|date=April 2013}} <math>{\rm NormalGamma}(b\mu, \lambda, \alpha, b^2\beta).</math>{{dubious|date=April 2013}} | |||
== Posterior distribution of the parameters == | |||
Assume that ''x'' is distributed according to a normal distribution with unknown mean <math>\mu</math> and precision <math>\tau</math>. | |||
:<math> x \sim \mathcal{N}(\mu, \tau^{-1}) </math> | |||
and that the prior distribution on <math>\mu</math> and <math>\tau</math>, <math>(\mu,\tau)</math>, has a normal-gamma distribution | |||
:<math> | |||
(\mu,\tau) \sim \text{NormalGamma}(\mu_0,\lambda_0,\alpha_0,\beta_0) , | |||
</math> | |||
for which the density ''π'' satisfies | |||
:<math> | |||
\pi(\mu,\tau) \propto \tau^{\alpha_0-\frac{1}{2}}\,\exp[{-\beta_0\tau}]\,\exp[{ -\frac{\lambda_0\tau(\mu-\mu_0)^2}{2}}]. | |||
</math> | |||
Given a dataset <math> \mathbf{X} </math>, consisting of <math>n</math> [[independent and identically distributed random_variables]] (i.i.d), <math> \{x_1,...,x_n\}</math>, the posterior distribution of <math>\mu</math> and <math>\tau</math> given this dataset can be analytically determined by [[Bayes' theorem]]. Explicitly,{{cn|date=April 2013}} | |||
:<math>\mathbf{P}(\tau,\mu | \mathbf{X}) \propto \mathbf{L}(\mathbf{X} | \tau,\mu) \pi(\tau,\mu)</math>, | |||
where <math>\mathbf{L}</math> is the likelihood of the data given the parameters. | |||
Since the data are i.i.d, the likelihood of the entire dataset is equal to the product of the likelihoods of the individual data samples: | |||
:<math> | |||
\mathbf{L}(\mathbf{X} | \tau, \mu) = \prod_{i=1}^n \mathbf{L}(x_i | \tau, \mu) . | |||
</math> | |||
This expression can be simplified as follows: | |||
:<math> | |||
\begin{align} | |||
\mathbf{L}(\mathbf{X} | \tau, \mu) & \propto \prod_{i=1}^n \tau^{1/2} \exp[\frac{-\tau}{2}(x_i-\mu)^2] \\ | |||
& \propto \tau^{n/2} \exp[\frac{-\tau}{2}\sum_{i=1}^n(x_i-\mu)^2] \\ | |||
& \propto \tau^{n/2} \exp[\frac{-\tau}{2}\sum_{i=1}^n(x_i-\bar{x} +\bar{x} -\mu)^2] \\ | |||
& \propto \tau^{n/2} \exp[\frac{-\tau}{2}\sum_{i=1}^n\left((x_i-\bar{x})^2 + (\bar{x} -\mu)^2\right)] \\ | |||
& \propto \tau^{n/2} \exp[\frac{-\tau}{2}\left(n s + n(\bar{x} -\mu)^2\right)] , | |||
\end{align} | |||
</math> | |||
where <math>\bar{x}= \frac{1}{n}\sum_{i=1}^n x_i</math>, the mean of the data samples, and <math>s= \frac{1}{n} \sum_{i=1}^n(x_i-\bar{x})^2</math>, the sample variance. | |||
The posterior distribution of the parameters is proportional to the prior times the likelihood. | |||
:<math> | |||
\begin{align} | |||
\mathbf{P}(\tau, \mu | \mathbf{X}) &\propto \mathbf{L}(\mathbf{X} | \tau,\mu) \pi(\tau,\mu) \\ | |||
&\propto \tau^{n/2} \exp[\frac{-\tau}{2}\left(n s + n(\bar{x} -\mu)^2\right)] | |||
\tau^{\alpha_0-\frac{1}{2}}\,\exp[{-\beta_0\tau}]\,\exp[{ -\frac{\lambda_0\tau(\mu-\mu_0)^2}{2}}] \\ | |||
&\propto \tau^{\frac{n}{2} + \alpha_0 - \frac{1}{2}}\exp[-\tau \left( \frac{1}{2} n s + \beta_0 \right) ] \exp\left[- \frac{\tau}{2}\left(\lambda_0(\mu-\mu_0)^2 + n(\bar{x} -\mu)^2\right)\right] \\ | |||
\end{align} | |||
</math> | |||
The final exponential term is simplified by completing the square. | |||
:<math> | |||
\begin{align} | |||
\lambda_0(\mu-\mu_0)^2 + n(\bar{x} -\mu)^2&=\lambda_0 \mu^2 - 2 \lambda_0 \mu \mu_0 + \lambda_0 \mu_0^2 + n \mu^2 - 2 n \bar{x} \mu + n \bar{x}^2 \\ | |||
&= (\lambda_0 + n) \mu^2 - 2(\lambda_0 \mu_0 + n \bar{x}) \mu + \lambda_0 \mu_0^2 +n \bar{x}^2 \\ | |||
&= (\lambda_0 + n)( \mu^2 - 2 \frac{\lambda_0 \mu_0 + n \bar{x}}{\lambda_0 + n} \mu ) + \lambda_0 \mu_0^2 +n \bar{x}^2 \\ | |||
&= (\lambda_0 + n)\left(\mu - \frac{\lambda_0 \mu_0 + n \bar{x}}{\lambda_0 + n} \right) ^2 + \lambda_0 \mu_0^2 +n \bar{x}^2 - \frac{\left(\lambda_0 \mu_0 +n \bar{x}\right)^2} {\lambda_0 + n} \\ | |||
&= (\lambda_0 + n)\left(\mu - \frac{\lambda_0 \mu_0 + n \bar{x}}{\lambda_0 + n} \right) ^2 + \frac{\lambda_0 n (\bar{x} - \mu_0 )^2}{\lambda_0 +n} | |||
\end{align} | |||
</math> | |||
On inserting this back into the expression above, | |||
:<math> | |||
\begin{align} | |||
\mathbf{P}(\tau, \mu | \mathbf{X}) & \propto \tau^{\frac{n}{2} + \alpha_0 - \frac{1}{2}} \exp \left[-\tau \left( \frac{1}{2} n s + \beta_0 \right) \right] \exp \left[- \frac{\tau}{2} \left( \left(\lambda_0 + n \right) \left(\mu- \frac{\lambda_0 \mu_0 + n \bar{x}}{\lambda_0 + n} \right)^2 + \frac{\lambda_0 n (\bar{x} - \mu_0 )^2}{\lambda_0 +n} \right) \right]\\ | |||
& \propto \tau^{\frac{n}{2} + \alpha_0 - \frac{1}{2}} \exp \left[-\tau \left( \frac{1}{2} n s + \beta_0 + \frac{\lambda_0 n (x - \mu_0 )^2}{2(\lambda_0 +n)} \right) \right] \exp \left[- \frac{\tau}{2} \left(\lambda_0 + n \right) \left(\mu- \frac{\lambda_0 \mu_0 + n \bar{x}}{\lambda_0 + n} \right)^2 \right] | |||
\end{align} | |||
</math> | |||
This final expression is in exactly the same form as a Normal-Gamma distribution, i.e., | |||
:<math> | |||
\mathbf{P}(\tau, \mu | \mathbf{X}) = \text{NormalGamma}\left(\frac{\lambda_0 \mu_0 + n \bar{x}}{\lambda_0 + n}, \lambda_0 + n, \alpha_0+\frac{n}{2}, \beta_0+ \frac{1}{2}\left(n s + \frac{\lambda_0 n (\bar{x} - \mu_0 )^2}{\lambda_0 +n} \right) \right) | |||
</math> | |||
=== Interpretation of parameters === | |||
The interpretation of parameters in terms of pseudo-observations is as follows: | |||
*The new mean takes a weighted average of the old pseudo-mean and the observed mean, weighted by the number of associated (pseudo-)observations. | |||
*The precision was estimated from <math>2\alpha</math> pseudo-observations (i.e. possibly a different number of pseudo-observations, to allow the variance of the mean and precision to be controlled separately) with sample mean <math>\mu</math> and sample variance <math>\frac{\beta}{\alpha}</math> (i.e. with sum of [[squared deviations]] <math>2\beta</math>). | |||
*The posterior updates the number of pseudo-observations (<math>\lambda_{0}</math>) simply by adding up the corresponding number of new observations (<math>n</math>). | |||
*The new sum of squared deviations is computed by adding the previous respective sums of squared deviations. However, a third "interaction term" is needed because the two sets of squared deviations were computed with respect to different means, and hence the sum of the two underestimates the actual total squared deviation. | |||
As a consequence, if one has a prior mean of <math>\mu_0</math> from <math> n_\mu </math> samples and a prior precision of <math> \tau_0 </math> from <math>n_\tau</math> samples, the prior distribution over <math> \mu </math> and <math> \tau </math> is | |||
:<math> | |||
\mathbf{P}(\tau,\mu | \mathbf{X}) = \text{NormalGamma}(\mu_0, n_\mu ,\frac{n_\tau}{2}, \frac{n_\tau}{2 \tau_0}) | |||
</math> | |||
and after observing <math>n</math> samples with mean <math>\mu</math> and variance <math>s</math>, the posterior probability is | |||
:<math> | |||
\mathbf{P}(\tau,\mu | \mathbf{X}) = \text{NormalGamma}\left( \frac{n_\mu \mu_0 + n \mu}{n_\mu +n}, n_\mu +n ,\frac{1}{2}(n_\tau+n), \frac{1}{2}\left(\frac{n_\tau}{\tau_0} + n s + \frac{n_\mu n (\mu-\mu_0)^2}{n_\mu+n}\right) \right) | |||
</math> | |||
Note that in some programming languages, such as [[Matlab]], the gamma distribution is implemented with the inverse definition of <math>\beta</math>, so the fourth argument of the Normal-Gamma distribution is <math> 2 \tau_0 /n_\tau</math>. | |||
== Generating normal-gamma random variates == | |||
Generation of random variates is straightforward: | |||
# Sample <math>\tau</math> from a gamma distribution with parameters <math>\alpha</math> and <math>\beta</math> | |||
# Sample <math>x</math> from a normal distribution with mean <math>\mu</math> and variance <math>1/(\lambda \tau)</math> | |||
== Related distributions == | |||
* The [[normal-inverse-gamma distribution]] is essentially the same distribution parameterized by variance rather than precision | |||
* The [[normal-exponential-gamma distribution]] | |||
==Notes== | |||
{{reflist}} | |||
== References == | |||
* Bernardo, J.M.; Smith, A.F.M. (1993) ''Bayesian Theory'', Wiley. ISBN 0-471-49464-X | |||
* Dearden et al. [http://www.aaai.org/Papers/AAAI/1998/AAAI98-108.pdf "Bayesian Q-learning"], ''Proceedings of the Fifteenth National Conference on Artificial Intelligence (AAAI-98)'', July 26–30, 1998, Madison, Wisconsin, USA. | |||
{{ProbDistributions|multivariate}} | |||
{{DEFAULTSORT:Normal-gamma distribution}} | |||
[[Category:Multivariate continuous distributions]] | |||
[[Category:Conjugate prior distributions]] | |||
[[Category:Normal distribution]] | |||
[[Category:Probability distributions]] | |||
Latest revision as of 23:32, 20 February 2013
Template:Probability distribution In probability theory and statistics, the normal-gamma distribution (or Gaussian-gamma distribution) is a bivariate four-parameter family of continuous probability distributions. It is the conjugate prior of a normal distribution with unknown mean and precision.[1]
Definition
For a pair of random variable, (X,T), suppose that the conditional distribution of X given T is given by
meaning that the condition distribution is a normal distribution with mean and precision — equivalently, with variance
Suppose also that the marginal distribution of T is given by
where this means that T has a gamma distribution. Here λ, α and β are parameters of the joint distribution.
Then (X,T) has a normal-gamma distribution, and this is denoted by
Properties
Probability density function
The joint probability density function of (X,T) isTemplate:Cn
Marginal distributions
By construction, the marginal distribution over is a gamma distribution, and the conditional distribution over given is a Gaussian distribution. The marginal distribution over is a three-parameter non-standardized Student's t-distribution with parameters .Template:Cn
Exponential family
The normal-gamma distribution is a four-parameter exponential family with natural parameters and natural statistics .Template:Cn
Moments of the natural statistics
The following moments can be easily computed using the moment generating function of the sufficient statistic:Template:Cn
- , where is the digamma function,
- ,
- ,
- .
Scaling
If then for any b > 0, (bX,bT) is distributed asTemplate:Cn To succeed in selling a home, it is advisable be competent in real estate advertising and marketing, authorized, monetary, operational aspects, and other information and skills. This is essential as a result of you want to negotiate with more and more sophisticated buyers. You could outperform rivals, use latest technologies, and stay ahead of the fast altering market.
Home is where the center is, and choosing the right house is a part of guaranteeing a contented expertise in Singapore. Most expats sign up for a two-year lease with the option to resume, so it is value taking the time to choose a neighbourhood that has the services you want. The experts at Expat Realtor have compiled the next data that will help you negotiate your means by way of the property minefield. Some government state properties for rent. Over 2000 units available for lease however occupancy is often excessive. Some properties come under a bidding system. Their property brokers embody DTZ and United Premas. Up to date serviced residences located just off Orchard Highway. one hundred sixty Orchard Highway, #06-01 Orchard Level, Singapore 238842. Institute Of Property Agents
There is no such thing as a deal too small. Property agents who're willing to find time for any deal even when the commission is small are those you want in your side. They also show humbleness and might relate with the average Singaporean higher. Relentlessly pursuing any deal, calling prospects even without being prompted. Even when they get rejected a hundred times, they still come back for more. These are the property brokers who will find consumers what they want finally, and who would be the most profitable in what they do. four. Honesty and Integrity
As a realtor, you're our own business. Due to this fact, it is imperative that you handle yours prices and spend money correctly in order to market your property successfully. Also, beware of mentors who always ask you to pay for pointless costs. Such mentors typically are recruiting to develop a staff and see you as a option to defray advertising and marketing prices. For foreigners who want to register with CEA as salespersons, they might want to have a valid Employment Cross (EP) issued by the Ministry of Manpower (MOM). They should consult an property agent that is ready to assist their future registration software, who would then examine with CEA. Thereafter, after they register for the RES Course, they might want to produce a letter of assist from the property agent."
Main Real Property Brokers with in depth local knowledge, Carole Ann, Elizabeth and their group of extremely skilled property consultants provide a personalised service, for those looking to buy, lease or promote in Singapore. Relocation companies out there. Properties for the aesthete. Boutique real property agency for architecturally distinguished, unique properties for rent and on the market. Caters to the niche market of design-savvy people. Sale, letting and property management and taxation services. three Shenton Means, #10-08 Shenton Home, Singapore 068805. Buy property, promote or leasing estate company. 430 Lorong 6 Toa Payoh, #08-01 OrangeTee Constructing, Singapore 319402. HIGH Date / Age of property Estate Agents and Home Search Services Property Information Highlights Prime Achievers
From the above info, you may see that saving on agent's commission will not cover the expenses wanted to market your home efficiently. As well as, it's essential make investments a whole lot of time, vitality and effort. By taking yourself away from your work and other endeavors, additionally, you will incur unnecessary opportunity prices. There may be additionally no assurance you could beat the market and get the outcomes you need. That is why you want an agent - not simply an ordinary agent - you want knowledgeable and competent specialist, geared up with the best instruments and knowledge to serve you and lead you to success! Within the midst of this ‘uniquely Singapore' Property GSS, our most needed foreign customers are nowhere to be seen. Different types of Public Residential properties
Based on Kelvin, other agents may also make use of your agent's listings. "If your pricing is on the excessive aspect, these brokers may use your house to persuade their patrons why Http://Trafficstooges.Com/Singapore-Property-Condominium they should purchase another residence." To counter this, Kelvin says it is crucial for your agent to supply a current market analysis before putting up your private home for sale. "This helps you worth your property appropriately and realistically." When property is made accessible (HIGH is issued) to the client. Becoming a successful property agent is a distinct story altogether! Hi, I would like to ask how I might be a property agent and whether there are courses I might take. And if I need to be at a certain age. www. Property BUYER com.sg (your impartial Mortgage Advisor) In private properties in
Posterior distribution of the parameters
Assume that x is distributed according to a normal distribution with unknown mean and precision .
and that the prior distribution on and , , has a normal-gamma distribution
for which the density π satisfies
Given a dataset , consisting of independent and identically distributed random_variables (i.i.d), , the posterior distribution of and given this dataset can be analytically determined by Bayes' theorem. Explicitly,Template:Cn
where is the likelihood of the data given the parameters.
Since the data are i.i.d, the likelihood of the entire dataset is equal to the product of the likelihoods of the individual data samples:
This expression can be simplified as follows:
where , the mean of the data samples, and , the sample variance.
The posterior distribution of the parameters is proportional to the prior times the likelihood.
The final exponential term is simplified by completing the square.
On inserting this back into the expression above,
This final expression is in exactly the same form as a Normal-Gamma distribution, i.e.,
Interpretation of parameters
The interpretation of parameters in terms of pseudo-observations is as follows:
- The new mean takes a weighted average of the old pseudo-mean and the observed mean, weighted by the number of associated (pseudo-)observations.
- The precision was estimated from pseudo-observations (i.e. possibly a different number of pseudo-observations, to allow the variance of the mean and precision to be controlled separately) with sample mean and sample variance (i.e. with sum of squared deviations ).
- The posterior updates the number of pseudo-observations () simply by adding up the corresponding number of new observations ().
- The new sum of squared deviations is computed by adding the previous respective sums of squared deviations. However, a third "interaction term" is needed because the two sets of squared deviations were computed with respect to different means, and hence the sum of the two underestimates the actual total squared deviation.
As a consequence, if one has a prior mean of from samples and a prior precision of from samples, the prior distribution over and is
and after observing samples with mean and variance , the posterior probability is
Note that in some programming languages, such as Matlab, the gamma distribution is implemented with the inverse definition of , so the fourth argument of the Normal-Gamma distribution is .
Generating normal-gamma random variates
Generation of random variates is straightforward:
- Sample from a gamma distribution with parameters and
- Sample from a normal distribution with mean and variance
Related distributions
- The normal-inverse-gamma distribution is essentially the same distribution parameterized by variance rather than precision
- The normal-exponential-gamma distribution
Notes
43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.
References
- Bernardo, J.M.; Smith, A.F.M. (1993) Bayesian Theory, Wiley. ISBN 0-471-49464-X
- Dearden et al. "Bayesian Q-learning", Proceedings of the Fifteenth National Conference on Artificial Intelligence (AAAI-98), July 26–30, 1998, Madison, Wisconsin, USA.
55 yrs old Metal Polisher Records from Gypsumville, has interests which include owning an antique car, summoners war hack and spelunkering. Gets immense motivation from life by going to places such as Villa Adriana (Tivoli).
my web site - summoners war hack no survey ios
- ↑ Bernardo & Smith (1993, pages 136, 268, 434)