<?xml version="1.0" encoding="UTF-8"?><rss version="2.0"
	xmlns:content="http://purl.org/rss/1.0/modules/content/"
	xmlns:dc="http://purl.org/dc/elements/1.1/"
	xmlns:atom="http://www.w3.org/2005/Atom"
	xmlns:sy="http://purl.org/rss/1.0/modules/syndication/"
	
	xmlns:georss="http://www.georss.org/georss"
	xmlns:geo="http://www.w3.org/2003/01/geo/wgs84_pos#"
	
	>
<channel>
	<title>
	Comments on: Principal Component Regression: Part 1 &#8211; The Magic of the SVD	</title>
	<atom:link href="https://www.centerspace.net/theoretical-motivation-behind-pcr/feed" rel="self" type="application/rss+xml" />
	<link>https://www.centerspace.net/theoretical-motivation-behind-pcr</link>
	<description>.NET numerical class libraries</description>
	<lastBuildDate>Sun, 03 May 2020 15:30:07 +0000</lastBuildDate>
	<sy:updatePeriod>
	hourly	</sy:updatePeriod>
	<sy:updateFrequency>
	1	</sy:updateFrequency>
	<generator>https://wordpress.org/?v=6.1.1</generator>
	<item>
		<title>
		By: Principal Components Regression: Part 3 – The NIPALS Algorithm - CenterSpace		</title>
		<link>https://www.centerspace.net/theoretical-motivation-behind-pcr#comment-111331</link>

		<dc:creator><![CDATA[Principal Components Regression: Part 3 – The NIPALS Algorithm - CenterSpace]]></dc:creator>
		<pubDate>Tue, 29 Nov 2016 22:05:10 +0000</pubDate>
		<guid isPermaLink="false">http://www.centerspace.net/blog/?p=1307#comment-111331</guid>

					<description><![CDATA[[&#8230;] in this series, &#8220;Principal Component Regression: Part 1 – The Magic of the SVD&#8221; is here. And the second, &#8220;Principal Components Regression: Part 2 – The Problem With Linear [&#8230;]]]></description>
			<content:encoded><![CDATA[<p>[&#8230;] in this series, &#8220;Principal Component Regression: Part 1 – The Magic of the SVD&#8221; is here. And the second, &#8220;Principal Components Regression: Part 2 – The Problem With Linear [&#8230;]</p>
]]></content:encoded>
		
			</item>
		<item>
		<title>
		By: Willem		</title>
		<link>https://www.centerspace.net/theoretical-motivation-behind-pcr#comment-4505</link>

		<dc:creator><![CDATA[Willem]]></dc:creator>
		<pubDate>Mon, 02 Jul 2012 11:14:45 +0000</pubDate>
		<guid isPermaLink="false">http://www.centerspace.net/blog/?p=1307#comment-4505</guid>

					<description><![CDATA[Hi Steve,

thank you very much for this article. You mentioned a third part in which you discuss the implementation of PCR. Is it correct that the third part is missing?

Greetings,
Willem]]></description>
			<content:encoded><![CDATA[<p>Hi Steve,</p>
<p>thank you very much for this article. You mentioned a third part in which you discuss the implementation of PCR. Is it correct that the third part is missing?</p>
<p>Greetings,<br />
Willem</p>
]]></content:encoded>
		
			</item>
		<item>
		<title>
		By: Trevor Misfeldt		</title>
		<link>https://www.centerspace.net/theoretical-motivation-behind-pcr#comment-3943</link>

		<dc:creator><![CDATA[Trevor Misfeldt]]></dc:creator>
		<pubDate>Fri, 30 Sep 2011 23:21:06 +0000</pubDate>
		<guid isPermaLink="false">http://www.centerspace.net/blog/?p=1307#comment-3943</guid>

					<description><![CDATA[In reply to &lt;a href=&quot;https://www.centerspace.net/theoretical-motivation-behind-pcr#comment-3845&quot;&gt;Eran Zadok&lt;/a&gt;.

1. Look at the two dimension case. Consider the image of the unit sphere S={v:‖v‖≤1} under X. Since S is compact and X is continuous there exists a vector v_1 such that ‖Xv_1 ‖ is maximal. That is 

‖Xv_1 ‖≥‖Xv‖  ∀ v∈S

Set 〖σu〗_1=Xv_1, where u_1 is a unit vector so that ‖Xv_1 ‖=σ. We will show that for any vector v_2 orthogonal to v_1its image u_2=Xv_2is orthogonal to u_1.

Choose a unit vector u_2  orthogonal to  u_1 (ala Gram-Schmidt, for example), and write

 Xv_2=αu_1+βu_2

We will show that α=0. Set

v=(σv_1+αv_2)/(σ^2+α^2 )^(1⁄2) 

Then

Xv=(σ^2 u_1+α^2 u_1+αβu_2)/(σ^2+α^2 )^(1⁄2) =((σ^2+α^2 ) u_1+αβu_2)/(σ^2+α^2 )^(1⁄2) 


Taking the norm of both sides and using the fact that u_1 and u_2 are orthonormal vectors* yields

‖Xv‖=(‖(σ^2+α^2 ) u_1 ‖^2+‖αβu_2 ‖^2 )^(1⁄2)/(σ^2+α^2 )^(1⁄2) =((σ^2+α^2 )^2+(αβ)^2 )^(1⁄2)/(σ^2+α^2 )^(1⁄2) ≥(σ^2+α^2)/(σ^2+α^2 )^(1⁄2) 

So that
〖‖Xv‖≥(σ^2+α^2 )〗^(1⁄2)≥σ

Since ‖Xv‖ cannot be strictly greater than σ, equality must hold which means α=0.

*For orthogonal vectors u and v we have:

‖u+v‖^2=(u+v)∙(u+v)=u∙u+u∙v+v∙u+v∙v=‖u‖^2+‖v‖^2

2. The terms “left” and “right” do seem a bit awkward. The come from the positions of the factors U and V in the decomposition X=UΣV^T.]]></description>
			<content:encoded><![CDATA[<p>In reply to <a href="https://www.centerspace.net/theoretical-motivation-behind-pcr#comment-3845">Eran Zadok</a>.</p>
<p>1. Look at the two dimension case. Consider the image of the unit sphere S={v:‖v‖≤1} under X. Since S is compact and X is continuous there exists a vector v_1 such that ‖Xv_1 ‖ is maximal. That is </p>
<p>‖Xv_1 ‖≥‖Xv‖  ∀ v∈S</p>
<p>Set 〖σu〗_1=Xv_1, where u_1 is a unit vector so that ‖Xv_1 ‖=σ. We will show that for any vector v_2 orthogonal to v_1its image u_2=Xv_2is orthogonal to u_1.</p>
<p>Choose a unit vector u_2  orthogonal to  u_1 (ala Gram-Schmidt, for example), and write</p>
<p> Xv_2=αu_1+βu_2</p>
<p>We will show that α=0. Set</p>
<p>v=(σv_1+αv_2)/(σ^2+α^2 )^(1⁄2) </p>
<p>Then</p>
<p>Xv=(σ^2 u_1+α^2 u_1+αβu_2)/(σ^2+α^2 )^(1⁄2) =((σ^2+α^2 ) u_1+αβu_2)/(σ^2+α^2 )^(1⁄2) </p>
<p>Taking the norm of both sides and using the fact that u_1 and u_2 are orthonormal vectors* yields</p>
<p>‖Xv‖=(‖(σ^2+α^2 ) u_1 ‖^2+‖αβu_2 ‖^2 )^(1⁄2)/(σ^2+α^2 )^(1⁄2) =((σ^2+α^2 )^2+(αβ)^2 )^(1⁄2)/(σ^2+α^2 )^(1⁄2) ≥(σ^2+α^2)/(σ^2+α^2 )^(1⁄2) </p>
<p>So that<br />
〖‖Xv‖≥(σ^2+α^2 )〗^(1⁄2)≥σ</p>
<p>Since ‖Xv‖ cannot be strictly greater than σ, equality must hold which means α=0.</p>
<p>*For orthogonal vectors u and v we have:</p>
<p>‖u+v‖^2=(u+v)∙(u+v)=u∙u+u∙v+v∙u+v∙v=‖u‖^2+‖v‖^2</p>
<p>2. The terms “left” and “right” do seem a bit awkward. The come from the positions of the factors U and V in the decomposition X=UΣV^T.</p>
]]></content:encoded>
		
			</item>
		<item>
		<title>
		By: Eran Zadok		</title>
		<link>https://www.centerspace.net/theoretical-motivation-behind-pcr#comment-3845</link>

		<dc:creator><![CDATA[Eran Zadok]]></dc:creator>
		<pubDate>Tue, 26 Jul 2011 23:32:40 +0000</pubDate>
		<guid isPermaLink="false">http://www.centerspace.net/blog/?p=1307#comment-3845</guid>

					<description><![CDATA[Hi Ken,

two questions:

1.In the section &#039;a geometric fact&#039; you say &#039;...the columns of  V are orthonormal by construction&#039;.
 It may be so, but it is not obvious, since V columns are defined &#039;the pre-images of the principal semiaxes of XS&#039;.

2. In the same section you define &#039;u&#039; as left singular vectors and &#039;v&#039; right singular vectors. Shouldn&#039;t it be the other way around?

Regards,
Eran.]]></description>
			<content:encoded><![CDATA[<p>Hi Ken,</p>
<p>two questions:</p>
<p>1.In the section &#8216;a geometric fact&#8217; you say &#8216;&#8230;the columns of  V are orthonormal by construction&#8217;.<br />
 It may be so, but it is not obvious, since V columns are defined &#8216;the pre-images of the principal semiaxes of XS&#8217;.</p>
<p>2. In the same section you define &#8216;u&#8217; as left singular vectors and &#8216;v&#8217; right singular vectors. Shouldn&#8217;t it be the other way around?</p>
<p>Regards,<br />
Eran.</p>
]]></content:encoded>
		
			</item>
		<item>
		<title>
		By: Ken Baldwin		</title>
		<link>https://www.centerspace.net/theoretical-motivation-behind-pcr#comment-3524</link>

		<dc:creator><![CDATA[Ken Baldwin]]></dc:creator>
		<pubDate>Mon, 28 Mar 2011 17:33:00 +0000</pubDate>
		<guid isPermaLink="false">http://www.centerspace.net/blog/?p=1307#comment-3524</guid>

					<description><![CDATA[Fixed, thanks]]></description>
			<content:encoded><![CDATA[<p>Fixed, thanks</p>
]]></content:encoded>
		
			</item>
		<item>
		<title>
		By: Lei Fang		</title>
		<link>https://www.centerspace.net/theoretical-motivation-behind-pcr#comment-3394</link>

		<dc:creator><![CDATA[Lei Fang]]></dc:creator>
		<pubDate>Wed, 16 Mar 2011 20:57:48 +0000</pubDate>
		<guid isPermaLink="false">http://www.centerspace.net/blog/?p=1307#comment-3394</guid>

					<description><![CDATA[In matrix form we have
AV = U \Sigma

should be
XV = U \Sigma]]></description>
			<content:encoded><![CDATA[<p>In matrix form we have<br />
AV = U \Sigma</p>
<p>should be<br />
XV = U \Sigma</p>
]]></content:encoded>
		
			</item>
		<item>
		<title>
		By: CenterSpace Blog &#187; Blog Archive &#187; Principal Components Regression: Part 2 The Problem With Linear Regression		</title>
		<link>https://www.centerspace.net/theoretical-motivation-behind-pcr#comment-1772</link>

		<dc:creator><![CDATA[CenterSpace Blog &#187; Blog Archive &#187; Principal Components Regression: Part 2 The Problem With Linear Regression]]></dc:creator>
		<pubDate>Wed, 10 Mar 2010 17:23:20 +0000</pubDate>
		<guid isPermaLink="false">http://www.centerspace.net/blog/?p=1307#comment-1772</guid>

					<description><![CDATA[[...] This is the second part in a three part series on PCR, the first article on the topic can be found here. [...]]]></description>
			<content:encoded><![CDATA[<p>[&#8230;] This is the second part in a three part series on PCR, the first article on the topic can be found here. [&#8230;]</p>
]]></content:encoded>
		
			</item>
		<item>
		<title>
		By: Brand Hunt		</title>
		<link>https://www.centerspace.net/theoretical-motivation-behind-pcr#comment-1649</link>

		<dc:creator><![CDATA[Brand Hunt]]></dc:creator>
		<pubDate>Wed, 10 Feb 2010 22:47:07 +0000</pubDate>
		<guid isPermaLink="false">http://www.centerspace.net/blog/?p=1307#comment-1649</guid>

					<description><![CDATA[Great article; I also dig that book.  Also, it looks as if your code snippets aren&#039;t using character entity references correctly (or your doubling up?) -- the loop has a the &quot;lt&quot; escape instead of the less than character.]]></description>
			<content:encoded><![CDATA[<p>Great article; I also dig that book.  Also, it looks as if your code snippets aren&#8217;t using character entity references correctly (or your doubling up?) &#8212; the loop has a the &#8220;lt&#8221; escape instead of the less than character.</p>
]]></content:encoded>
		
			</item>
	</channel>
</rss>
