Thin plate spline: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
 
No edit summary
Line 1: Line 1:
Nice to meet you, my title is Ling and I completely dig that name. Playing crochet is something that I've done for many years. Bookkeeping is what I do for a residing. Kansas is where her home is but she needs to move because of her family members.<br><br>Also visit my blog post - [http://Pcresourcesllc.com/UserProfile/tabid/61/userId/94025/Default.aspx Pcresourcesllc.com]
{{Probability distribution |
  name      =Inverse-Wishart|
  type      =density|
  pdf_image  =|
  cdf_image  =|
  notation  =<math> \mathcal{W}^{-1}({\mathbf\Psi},\nu)</math>|
  parameters =<math> \nu > p-1  </math> [[degrees of freedom (statistics)|degrees of freedom]] ([[real numbers|real]])<br /><math>\mathbf{\Psi} > 0</math> [[scale matrix]] ([[positive-definite matrix|pos. def]])|
  support    =<math>\mathbf{X}</math> is [[positive-definite matrix|positive definite]]|
  pdf        =<math>\frac{\left|{\mathbf\Psi}\right|^{\frac{\nu}{2}}}{2^{\frac{\nu p}{2}}\Gamma_p(\frac{\nu}{2})} \left|\mathbf{X}\right|^{-\frac{\nu+p+1}{2}}e^{-\frac{1}{2}\operatorname{tr}({\mathbf\Psi}\mathbf{X}^{-1})}</math>
*<math>\Gamma_p</math> is the [[multivariate gamma function]]
*<math>\mathrm{tr}</math> is the [[trace (linear algebra)|trace]] function
|
  cdf        =|
  mean      = <math>\frac{\mathbf{\Psi}}{\nu - p - 1}</math>|
  median    =|
  mode      = <math>\frac{\mathbf{\Psi}}{\nu + p + 1}</math><ref name="Ohagan">{{Cite book
| author = A. O'Hagan, and J. J. Forster
| title = Kendall's Advanced Theory of Statistics: Bayesian Inference
| volume = 2B
| edition = 2
| publisher = Arnold
| year = 2004
| isbn = 0-340-80752-0
  }}</ref>{{rp|406}}|
  variance  =see below|
  skewness  =|
  kurtosis  =|
  entropy    =|
  mgf        =|
  char      =|
}}
 
In [[statistics]], the '''inverse Wishart distribution''', also called the '''inverted Wishart distribution''', is a [[probability distribution]] defined on real-valued [[positive-definite matrix|positive-definite]] [[matrix (mathematics)|matrices]]. In [[Bayesian statistics]] it is used as the [[conjugate prior]] for the covariance matrix of a
[[multivariate normal]] distribution.
 
We say <math>\mathbf{X}</math> follows an inverse Wishart distribution, denoted as <math> \mathbf{X}\sim \mathcal{W}^{-1}({\mathbf\Psi},\nu)</math>, if its [[matrix inverse|inverse]] <math> \mathbf{X}^{-1}</math> has a [[Wishart distribution]] <math> \mathcal{W}({\mathbf \Psi}^{-1}, \nu) </math>.
 
==Density==
The [[probability density function]] of the inverse Wishart is:
 
:<math>
\frac{\left|{\mathbf\Psi}\right|^{\frac{\nu}{2}}}{2^{\frac{\nu p}{2}}\Gamma_p(\frac{\nu}{2})} \left|\mathbf{X}\right|^{-\frac{\nu+p+1}{2}}e^{-\frac{1}{2}\operatorname{tr}({\mathbf\Psi}\mathbf{X}^{-1})}
</math>
 
where <math>\mathbf{X}</math> and <math>{\mathbf\Psi}</math> are <math>p\times p</math> [[positive-definite matrix|positive definite]] matrices, and &Gamma;<sub>''p''</sub>(&middot;) is the [[multivariate gamma function]].
 
==Theorems==
===Distribution of the inverse of a Wishart-distributed matrix===
 
If <math>{\mathbf A}\sim \mathcal{W}({\mathbf\Sigma},\nu)</math> and <math>{\mathbf\Sigma}</math> is of size <math>p \times p</math>, then <math>\mathbf{X}={\mathbf A}^{-1}</math> has an inverse Wishart distribution <math>\mathbf{X}\sim \mathcal{W}^{-1}({\mathbf\Sigma}^{-1},\nu)</math> .<ref name="MardiaK1979Multivariate">{{Cite book
| author = [[Kanti V. Mardia]], J. T. Kent and J. M. Bibby
| title = Multivariate Analysis
| publisher = [[Academic Press]]
| year = 1979
| isbn = 0-12-471250-9
}}</ref>
 
===Marginal and conditional distributions from an inverse Wishart-distributed matrix===
 
Suppose <math>{\mathbf A}\sim \mathcal{W}^{-1}({\mathbf\Psi},\nu)</math> has an inverse Wishart distribution. Partition the matrices <math> {\mathbf A} </math> and <math> {\mathbf\Psi} </math> [[Conformable matrix|conformably]] with each other
:<math>
    {\mathbf{A}} = \begin{bmatrix} \mathbf{A}_{11} & \mathbf{A}_{12} \\ \mathbf{A}_{21} & \mathbf{A}_{22} \end{bmatrix}, \;
    {\mathbf{\Psi}} = \begin{bmatrix} \mathbf{\Psi}_{11} & \mathbf{\Psi}_{12} \\ \mathbf{\Psi}_{21} & \mathbf{\Psi}_{22} \end{bmatrix}
</math>
where <math>{\mathbf A_{ij}}</math> and <math>{\mathbf \Psi_{ij}} </math> are <math> p_{i}\times p_{j}</math> matrices, then we have
 
i) <math> {\mathbf A_{11} } </math> is independent of <math> {\mathbf A}_{11}^{-1}{\mathbf A}_{12} </math> and <math> {\mathbf A}_{22\cdot 1} </math>, where <math>{\mathbf A_{22\cdot 1}} = {\mathbf A}_{22} - {\mathbf A}_{21}{\mathbf A}_{11}^{-1}{\mathbf A}_{12}</math> is the [[Schur complement]] of <math> {\mathbf A_{11} } </math> in <math> {\mathbf A} </math>;
 
ii) <math> {\mathbf A_{11} } \sim \mathcal{W}^{-1}({\mathbf \Psi_{11} }, \nu-p_{2}) </math>;
 
iii) <math> {\mathbf A}_{11}^{-1} {\mathbf A}_{12}| {\mathbf A}_{22\cdot 1} \sim MN_{p_{1}\times p_{2}}
( {\mathbf \Psi}_{11}^{-1} {\mathbf \Psi}_{12},  {\mathbf A}_{22\cdot 1} \otimes  {\mathbf \Psi}_{11}^{-1}) </math>, where <math> MN_{p\times q}(\cdot,\cdot) </math> is a [[matrix normal distribution]];
 
iv) <math> {\mathbf A}_{22\cdot 1} \sim  \mathcal{W}^{-1}({\mathbf \Psi}_{22\cdot 1}, \nu) </math>, where <math>{\mathbf \Psi_{22\cdot 1}} = {\mathbf \Psi}_{22} - {\mathbf \Psi}_{21}{\mathbf \Psi}_{11}^{-1}{\mathbf \Psi}_{12}</math>;
 
===Conjugate distribution===
 
Suppose we wish to make inference about a covariance matrix <math>{\mathbf{\Sigma}}</math> whose [[prior probability|prior]] <math>{p(\mathbf{\Sigma})}</math> has a <math>\mathcal{W}^{-1}({\mathbf\Psi},\nu)</math> distribution.  If the observations <math>\mathbf{X}=[\mathbf{x}_1,\ldots,\mathbf{x}_n]</math> are independent p-variate Gaussian variables drawn from a <math>N(\mathbf{0},{\mathbf \Sigma})</math> distribution, then the conditional distribution  <math>{p(\mathbf{\Sigma}|\mathbf{X})}</math> has a <math>\mathcal{W}^{-1}({\mathbf A}+{\mathbf\Psi},n+\nu)</math> distribution, where <math>{\mathbf{A}}=\mathbf{X}\mathbf{X}^T</math> is <math>n-1</math> times the sample covariance matrix.
 
Because the prior and posterior distributions are the same family, we say the inverse Wishart distribution is [[conjugate prior|conjugate]] to the multivariate Gaussian.
 
Due to its conjugacy to the multivariate Gaussian, it is possible to [[marginalize out]] (integrate out) the Gaussian's parameter <math>\mathbf{\Sigma}</math>.
 
<math>P(\mathbf{X}|\mathbf{\Psi},\nu) = \int P(\mathbf{X}|\mathbf{\Sigma})P(\mathbf{\Sigma}|\mathbf{\Psi},\nu) d\mathbf{\Sigma} = \frac{|\mathbf{\Psi}|^{\frac{\nu}{2}}\Gamma_p\left(\frac{\nu+n}{2}\right)}{\pi^{\frac{np}{2}}|\mathbf{\Psi}+\mathbf{A}|^{\frac{\nu+n}{2}}\Gamma_p(\frac{\nu}{2})}</math>
 
(this is useful because the variance matrix <math>\mathbf{\Sigma}</math> is not known in practice, but because <math>{\mathbf\Psi}</math> is known ''a priori'', and <math>{\mathbf A}</math> can be obtained from the data, the right hand side can be evaluated directly).
 
===Moments===
The following is based on Press, S. J. (1982) "Applied Multivariate Analysis", 2nd ed. (Dover Publications, New York), after reparameterizing the degree of freedom to be consistent with the p.d.f. definition above.
 
The mean:<ref name="MardiaK1979Multivariate"/>{{rp|85}}
:<math>
E(\mathbf X) = \frac{\mathbf\Psi}{\nu-p-1}.</math>
 
The variance of each element of <math>\mathbf{X}</math>:
:<math>
\operatorname{Var}(x_{ij}) = \frac{(\nu-p+1)\psi_{ij}^2 + (\nu-p-1)\psi_{ii}\psi_{jj}}
{(\nu-p)(\nu-p-1)^2(\nu-p-3)}</math>
The variance of the diagonal uses the same formula as above with <math>i=j</math>, which simplifies to:
:<math>
\operatorname{Var}(x_{ii}) = \frac{2\psi_{ii}^2}{(\nu-p-1)^2(\nu-p-3)}.</math>
The covariance of elements of <math>\mathbf{X}</math> are given by:
:<math>
\operatorname{Cov}(x_{ij},x_{kl}) = \frac{2\psi_{ij}\psi_{kl} + (\nu-p-1) (\psi_{ik}\psi_{jl} + \psi_{il}\psi_{kj})}{(\nu-p)(\nu-p-1)^2(\nu-p-3)}</math>
 
== Related distributions ==
A [[univariate]] specialization of the inverse-Wishart distribution is the [[inverse-gamma distribution]]. With <math>p=1</math> (i.e. univariate) and <math>\alpha = \nu/2</math>, <math>\beta = \mathbf{\Psi}/2</math> and <math>x=\mathbf{X}</math> the [[probability density function]] of the inverse-Wishart distribution becomes
 
: <math>p(x|\alpha, \beta) = \frac{\beta^\alpha\, x^{-\alpha-1} \exp(-\beta/x)}{\Gamma_1(\alpha)}.</math>
 
i.e., the inverse-gamma distribution, where <math>\Gamma_1(\cdot)</math> is the ordinary [[Gamma function]].
 
A generalization is the [[inverse multivariate gamma distribution]].
 
Another generalization has been termed the generalized inverse Wishart distribution, <math>\mathcal{GW}^{-1}</math>. A <math> p \times p</math> positive definite matrix <math>\mathbf{X}</math> is said to be distributed as <math>\mathcal{GW}^{-1}(\mathbf{\Psi},\nu,\mathbf{S})</math> if <math>\mathbf{Y} = \mathbf{X}^{1/2}\mathbf{S}^{-1}\mathbf{X}^{1/2}</math> is distributed as <math>\mathcal{W}^{-1}(\mathbf{\Psi},\nu)</math>. Here <math>\mathbf{X}^{1/2}</math> denotes the symmetric matrix square root of <math>\mathbf{X}</math>, the parameters <math>\mathbf{\Psi},\mathbf{S}</math> are <math> p \times p</math> positive definite matrices, and the parameter <math>\nu</math> is a positive scalar larger than <math>2p</math>. Note that when <math>\mathbf{S}</math> is equal to an identity matrix, <math>\mathcal{GW}^{-1}(\mathbf{\Psi},\nu,\mathbf{S}) = \mathcal{W}^{-1}(\mathbf{\Psi},\nu)</math>. This generalized inverse Wishart distribution has been applied to estimating the distributions of multivariate autoregressive processes.<ref name="Triantafyllopoulos2011">{{Cite doi|10.1111/j.1467-9892.2010.00686.x}}</ref>
 
A different type of generalization is the [[normal-inverse-Wishart distribution]], essentially the product of a [[multivariate normal distribution]] with an inverse Wishart distribution.
 
==See also==
*[[Inverse multivariate gamma distribution]]
*[[Matrix normal distribution]]
*[[Wishart distribution]]
 
== References ==
<div class="references-small"><references/></div>
 
{{ProbDistributions|multivariate}}
 
[[Category:Continuous distributions]]
[[Category:Multivariate continuous distributions]]
[[Category:Conjugate prior distributions]]
[[Category:Probability distributions]]
[[Category:Exponential family distributions]]

Revision as of 17:53, 4 November 2013

Template:Probability distribution

In statistics, the inverse Wishart distribution, also called the inverted Wishart distribution, is a probability distribution defined on real-valued positive-definite matrices. In Bayesian statistics it is used as the conjugate prior for the covariance matrix of a multivariate normal distribution.

We say follows an inverse Wishart distribution, denoted as , if its inverse has a Wishart distribution .

Density

The probability density function of the inverse Wishart is:

where and are positive definite matrices, and Γp(·) is the multivariate gamma function.

Theorems

Distribution of the inverse of a Wishart-distributed matrix

If and is of size , then has an inverse Wishart distribution .[1]

Marginal and conditional distributions from an inverse Wishart-distributed matrix

Suppose has an inverse Wishart distribution. Partition the matrices and conformably with each other

where and are matrices, then we have

i) is independent of and , where is the Schur complement of in ;

ii) ;

iii) , where is a matrix normal distribution;

iv) , where ;

Conjugate distribution

Suppose we wish to make inference about a covariance matrix whose prior has a distribution. If the observations are independent p-variate Gaussian variables drawn from a distribution, then the conditional distribution has a distribution, where is times the sample covariance matrix.

Because the prior and posterior distributions are the same family, we say the inverse Wishart distribution is conjugate to the multivariate Gaussian.

Due to its conjugacy to the multivariate Gaussian, it is possible to marginalize out (integrate out) the Gaussian's parameter .

(this is useful because the variance matrix is not known in practice, but because is known a priori, and can be obtained from the data, the right hand side can be evaluated directly).

Moments

The following is based on Press, S. J. (1982) "Applied Multivariate Analysis", 2nd ed. (Dover Publications, New York), after reparameterizing the degree of freedom to be consistent with the p.d.f. definition above.

The mean:[1]Primarily based on the most recent URA personal property value index (PPPI) flash estimates, we know that the PPPI, which represents the overall real property price development, has dipped in 2013Q4. That is the first dip the market has seen within the final 2 years.

To give you some perspective, the entire number of personal properties in Singapore (together with govt condominiums) is 297,689 in 2013Q3. Primarily based on the projection that there will be 19,302 units accomplished in 2014, the rise in residential models works out to be more than 6%. With a lot New Ec Launch Singapore provide, buyers might be spoilt for alternative and this in flip will lead to their reluctance to pay a premium for potential models. The complete textual content of the Copyright Act (Cap sixty three) and different statutes regarding IPR might be found on the Singapore Statutes Online Website online The Group's income jumped forty.1 p.c to $324.5 million from $231.6 million in FY 2013, lifted by increased development income and sales of growth properties in Singapore and China. Actual Estate Shopping for

One factor we've on this nation is a big group of "economists," and "market analysts." What's interesting about this group of real property market-watchers is that there are two very other ways wherein they predict Boomers will affect housing markets over the subsequent decade. Let's check out those two opposites and see how every can change the best way real property investors strategy their markets. The good news is that actual property buyers are prepared for either state of affairs, and there's profit in being ready. I'm excited and searching ahead to the alternatives both or each of these conditions will supply; thank you Boomers! Mapletree to further broaden past Asia Why fortune will favour the brave in Asia's closing real property frontier

The story of the 23.2 home begins with a stack of Douglas fir beams salvaged from varied demolished warehouses owned by the consumer's household for a number of generations. Design and structure innovator Omer Arbel, configured them to type a triangulated roof, which makes up one of the placing features of the home. The transient from the entrepreneur-proprietor was not solely to design a house that integrates antique wood beams, however one which erases the excellence between inside and exterior. Built on a gentle slope on a large rural acreage surrounded by two masses of previous-development forests, the indoors movement seamlessly to the outdoors and, from the within looking, one enjoys unobstructed views of the existing natural panorama which is preserved

First, there are typically extra rental transactions than gross sales transactions, to permit AV to be decided for each property primarily based on comparable properties. Second, movements in sale costs are more unstable than leases. Hence, utilizing rental transactions to derive the AV helps to maintain property tax more steady for property homeowners. If you're buying or trying to lease a property. It's tiring to call up individual property agent, organize appointments, coordinate timing and to go for particular person property viewing. What most individuals do is to have a property agent representing them who will arrange and coordinate the viewings for all the properties out there based mostly on your necessities & most well-liked timing. Rent Property District 12 Rent Property District thirteen

The Annual Worth of a property is mostly derived based mostly on the estimated annual hire that it may well fetch if it have been rented out. In determining the Annual Worth of a property, IRAS will think about the leases of similar properties within the vicinity, dimension and condition of the property, and different relevant components. The Annual Worth of a property is determined in the identical method regardless of whether the property is let-out, proprietor-occupied or vacant. The Annual Worth of land is determined at 5% of the market price of the land. When a constructing is demolished, the Annual Worth of the land is assessed by this method. Property Tax on Residential Properties Buyer Stamp Responsibility on Buy of Properties – Business and residential properties Rent House District 01

Within the event the Bank's valuation is decrease than the acquisition price, the purchaser has to pay the distinction between the purchase value and the Bank's valuation utilizing money. As such, the money required up-front might be increased so it's at all times essential to know the valuation of the property before making any offer. Appoint Lawyer The Bank will prepare for a proper valuation of the property by way of physical inspection The completion statement will present you the balance of the acquisition price that you must pay after deducting any deposit, pro-rated property tax and utility costs, upkeep prices, and different relevant expenses in addition to any fees payable to the agent and the lawyer. Stamp Responsibility Primarily based on the Purchase Price or Market Value, whichever is larger

The variance of each element of :

The variance of the diagonal uses the same formula as above with , which simplifies to:

The covariance of elements of are given by:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \operatorname {Cov} (x_{ij},x_{kl})={\frac {2\psi _{ij}\psi _{kl}+(\nu -p-1)(\psi _{ik}\psi _{jl}+\psi _{il}\psi _{kj})}{(\nu -p)(\nu -p-1)^{2}(\nu -p-3)}}}

Related distributions

A univariate specialization of the inverse-Wishart distribution is the inverse-gamma distribution. With (i.e. univariate) and , and the probability density function of the inverse-Wishart distribution becomes

i.e., the inverse-gamma distribution, where is the ordinary Gamma function.

A generalization is the inverse multivariate gamma distribution.

Another generalization has been termed the generalized inverse Wishart distribution, . A positive definite matrix is said to be distributed as if is distributed as . Here denotes the symmetric matrix square root of , the parameters are positive definite matrices, and the parameter is a positive scalar larger than . Note that when is equal to an identity matrix, . This generalized inverse Wishart distribution has been applied to estimating the distributions of multivariate autoregressive processes.[2]

A different type of generalization is the normal-inverse-Wishart distribution, essentially the product of a multivariate normal distribution with an inverse Wishart distribution.

See also

References

  1. 1.0 1.1 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.

    My blog: http://www.primaboinca.com/view_profile.php?userid=5889534
  2. Template:Cite doi

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