ฉันชอบที่จะควบคุมวัตถุที่ฉันสร้างแม้ว่าพวกเขาอาจจะไม่ได้ตั้งใจก็ตาม
พิจารณาแล้วว่าเป็นไปได้ทั้งหมดแปรปรวนเมทริกซ์Σสามารถแสดงในรูปแบบn×nΣ
Σ=P′ Diagonal(σ1,σ2,…,σn) P
ที่เป็นมุมฉากเมทริกซ์และσ 1 ≥ σ 2 ≥ ⋯ ≥ σ n ≥ 0Pσ1≥σ2≥⋯≥σn≥0
เรขาคณิตนี้อธิบายถึงโครงสร้างความแปรปรวนที่มีช่วงขององค์ประกอบหลักที่มีขนาดฉัน ส่วนประกอบเหล่านี้ชี้ไปในทิศทางของแถวของP ดูตัวเลขที่ทำให้ความรู้สึกของการวิเคราะห์องค์ประกอบหลัก, eigenvectors และค่าลักษณะเฉพาะตัวอย่างกับn = 3 การตั้งค่าσ iจะกำหนดขนาดของโควาเรียสและขนาดสัมพัทธ์ แถวของPปรับแกนของรูปร่างตามที่คุณต้องการσiPn=3σiP
ประโยชน์เชิงพีชคณิตและการคำนวณอย่างหนึ่งของวิธีการนี้คือเมื่อ , Σกลับด้านได้อย่างง่ายดาย (ซึ่งเป็นการดำเนินการทั่วไปในเมทริกซ์ความแปรปรวนร่วม):σn>0Σ
Σ−1=P′ Diagonal(1/σ1,1/σ2,…,1/σn) P.
ไม่สนใจเกี่ยวกับทิศทาง แต่เพียงช่วงของขนาดที่ ? ไม่เป็นไร: คุณสามารถสร้างเมทริกซ์แบบ orthogonal แบบสุ่มได้อย่างง่ายดาย เพียงแค่ห่อn 2 iid มาตรฐานค่าปกติลงในเมทริกซ์จตุรัสแล้วปรับมุมมัน มันเกือบจะแน่นอนจะทำงาน (ให้nไม่ได้เป็นขนาดใหญ่) การย่อยสลาย QR จะทำเช่นนั้นในรหัสนี้σin2n
n <- 5
p <- qr.Q(qr(matrix(rnorm(n^2), n)))
สิ่งนี้ได้ผลเพราะการกระจายแบบพหุคูณ -variate ที่สร้างขึ้นจึงเป็น "รูปไข่": มันคงที่ภายใต้การหมุนและการสะท้อนทั้งหมด (ผ่านจุดกำเนิด) ดังนั้นเมทริกซ์ฉากมุมฉากทั้งหมดจะถูกสร้างขึ้นอย่างสม่ำเสมอตามที่ถกเถียงกันอยู่ว่าจะสร้างจุดกระจายอย่างสม่ำเสมอบนพื้นผิวของทรงกลม 3 มิติได้อย่างไร .n
วิธีที่รวดเร็วในการขอรับจากPและσ ฉันเมื่อคุณได้ระบุหรือสร้างพวกใช้และการหาประโยชน์'s กลับมาใช้ของอาร์เรย์ในการดำเนินการทางคณิตศาสตร์เช่นในตัวอย่างนี้กับσ = ( σ 1 , ... , σ 5 ) = ( 5 , 4 , 3 , 2 , 1 ) :ΣPσicrossprod
R
σ=(σ1,…,σ5)=(5,4,3,2,1)
Sigma <- crossprod(p, p*(5:1))
ขณะที่การตรวจสอบ, การสลายตัวมูลค่าเอกพจน์ควรกลับทั้งและP ' คุณสามารถตรวจสอบได้ด้วยคำสั่งσP′
svd(Sigma)
The inverse of Sigma
of course is obtained merely by changing the multiplication by σ into a division:
Tau <- crossprod(p, p/(5:1))
You may verify this by viewing zapsmall(Sigma %*% Tau)
, which should be the n×n identity matrix. A generalized inverse (essential for regression calculations) is obtained by replacing any σi≠0 by 1/σi, exactly as above, but keeping any zeros among the σi as they were.