Lewati ke konten utama
Logo Memahami Aljabar Linear

Subbab 6.5 Kuadrat terkecil ortogonal

Pembahasan kita sejauh ini telah menghasilkan seperangkat alat yang kaya untuk memanfaatkan keortogonalan dalam menyederhanakan berbagai persoalan. Seperti yang akan kita lihat dalam bagian ini, alat-alat tersebut secara alami mengarah pada beberapa teknik penting untuk mengolah data, seperti regresi linear.

Subbagian 6.5.1 Pendahuluan

Misalkan kita mengumpulkan sejumlah data ketika melakukan percobaan, lalu memplotnya seperti yang tampak di sebelah kiri Gambar 6.5.1. Perhatikan bahwa tidak ada satu garis yang memuat semua titik; bahkan, keberadaan garis semacam itu akan mengejutkan karena kita dapat mengharapkan adanya ketidakpastian dalam pengukuran yang dicatat. Namun, tampaknya ada sebuah garis, seperti yang diperlihatkan di sebelah kanan, yang hampir memuat titik-titik tersebut.
Plot pencar titik-titik data dua dimensi yang kira-kira terletak pada sebuah garis.
Plot pencar titik-titik data dua dimensi dan sebuah garis yang kira-kira memuat titik-titik tersebut.
Kontrol papan ketik untuk menjelajahi diagram
Tombol Tindakan
Enter, A Aktifkan penjelajahan dengan papan ketik
B Aktifkan penjelajahan dengan menu
Escape Keluar dari mode penjelajahan
Panah bawah Jelajahi tingkat berikutnya di bawah
Panah atas Jelajahi tingkat berikutnya di atas
Panah kanan Jelajahi elemen berikutnya pada tingkat ini
Panah kiri Jelajahi elemen sebelumnya pada tingkat ini
X Aktifkan atau nonaktifkan mode ahli
W Tampilkan rincian tambahan jika tersedia
Spasi Ulangi keluaran suara
M Aktifkan pembesaran bertahap
Koma Aktifkan pembesaran langsung
N Nonaktifkan pembesaran
Z Tampilkan atau sembunyikan subtitel
C Gilir pengaturan kontras
T Gunakan warna monokrom
L Ganti bahasa jika tersedia
K Hentikan suara saat ini
Y Hentikan keluaran suara
O Mulai atau hentikan sonifikasi
P Ulangi keluaran sonifikasi
Gambar 6.5.1. Sekumpulan titik dan sebuah garis yang menghampiri hubungan linear yang ditunjukkan oleh titik-titik tersebut.
Dalam bagian ini, kita akan menelaah bagaimana teknik-teknik yang dikembangkan dalam bab ini memungkinkan kita menentukan garis yang paling baik menghampiri data. Secara lebih khusus, kita akan melihat bahwa pencarian garis yang melalui titik-titik data menghasilkan sistem tak konsisten \(A\xvec=\bvec\text{.}\) Karena kita tidak dapat menemukan solusi, sebagai gantinya kita mencari vektor \(\xvec\) sedemikian sehingga \(A\xvec\) sedekat mungkin dengan \(\bvec\text{.}\) Proyeksi ortogonal menyediakan alat yang tepat untuk melakukannya.

Aktivitas Pratinjau 6.5.1.

(a)
Apakah persamaan \(A\xvec=\bvec\) memiliki solusi jika \(A\) dan \(\bvec\) diberikan oleh
\begin{equation*} \begin{bmatrix} 1 \amp 2 \\ 2 \amp 5 \\ -1 \amp 0 \\ \end{bmatrix} \xvec = \threevec5{-3}{-1}\text{.} \end{equation*}
(b)
Kita mengetahui bahwa \(\threevec12{-1}\) dan \(\threevec250\) membentuk basis bagi \(\col(A)\text{.}\) Tentukan sebuah basis ortogonal bagi \(\col(A)\text{.}\)
(c)
Tentukan proyeksi ortogonal \(\widehat\bvec\) dari \(\bvec\) pada \(\col(A)\text{.}\)
(d)
Jelaskan mengapa persamaan \(A\xvec=\widehat\bvec\) pasti konsisten, lalu tentukan solusinya.

Subbagian 6.5.2 Contoh pertama

Ketika sebelumnya menjumpai sistem tak konsisten, kita hanya menyatakan bahwa tidak ada solusi, lalu beralih ke persoalan lain. Namun, aktivitas pratinjau memperlihatkan cara mencari solusi hampiran bagi sistem tak konsisten: jika \(A\xvec = \bvec\) tidak memiliki solusi, sebagai gantinya kita menyelesaikan sistem konsisten \(A\xvec = \bhat\text{,}\) dengan ruas kanannya merupakan proyeksi ortogonal \(\bvec\) pada \(\col(A)\text{.}\) Seperti yang akan kita lihat, dalam arti tertentu solusi ini adalah solusi terbaik yang mungkin.

Kegiatan 6.5.2.

Misalkan kita mempunyai tiga titik data \((1,1)\text{,}\) \((2,1)\text{,}\) dan \((3,3)\text{,}\) serta ingin menentukan sebuah garis yang melalui ketiganya.
  1. Plot ketiga titik ini pada Gambar 6.5.2. Dapatkah Anda menggambar sebuah garis yang melalui ketiga titik tersebut?
    Kontrol papan ketik untuk menjelajahi diagram
    Tombol Tindakan
    Enter, A Aktifkan penjelajahan dengan papan ketik
    B Aktifkan penjelajahan dengan menu
    Escape Keluar dari mode penjelajahan
    Panah bawah Jelajahi tingkat berikutnya di bawah
    Panah atas Jelajahi tingkat berikutnya di atas
    Panah kanan Jelajahi elemen berikutnya pada tingkat ini
    Panah kiri Jelajahi elemen sebelumnya pada tingkat ini
    X Aktifkan atau nonaktifkan mode ahli
    W Tampilkan rincian tambahan jika tersedia
    Spasi Ulangi keluaran suara
    M Aktifkan pembesaran bertahap
    Koma Aktifkan pembesaran langsung
    N Nonaktifkan pembesaran
    Z Tampilkan atau sembunyikan subtitel
    C Gilir pengaturan kontras
    T Gunakan warna monokrom
    L Ganti bahasa jika tersedia
    K Hentikan suara saat ini
    Y Hentikan keluaran suara
    O Mulai atau hentikan sonifikasi
    P Ulangi keluaran sonifikasi
    Kisi koordinat baku \(1\times1\) beserta sumbu-sumbunya. Rentang horizontal dan vertikal membentang dari \(-1\) hingga \(4\text{.}\)
    Gambar 6.5.2. Plot ketiga titik data di sini.
  2. Ingat bahwa persamaan garis dapat ditulis sebagai \(b + mx=y\text{,}\) dengan \(m\) sebagai gradien dan \(b\) sebagai titik potong sumbu \(y\text{.}\) Kita akan mencari \(b\) dan \(m\) sedemikian sehingga ketiga titik terletak pada garis itu.
    Titik data pertama \((1,1)\) memberikan sebuah persamaan bagi \(b\) dan \(m\text{.}\) Secara khusus, kita mengetahui bahwa ketika \(x=1\text{,}\) maka \(y=1\text{,}\) sehingga kita memperoleh \(b + m(1) = 1\) atau \(b + m = 1\text{.}\) Gunakan dua titik data lainnya untuk membentuk sistem linear yang memuat \(m\) dan \(b\text{.}\)
  3. Kita telah memperoleh sistem linear dengan tiga persamaan, masing-masing berasal dari satu titik data, untuk dua peubah tak diketahui \(b\) dan \(m\text{.}\) Tentukan matriks \(A\) dan vektor \(\bvec\) sehingga sistem tersebut berbentuk \(A\xvec=\bvec\text{,}\) dengan \(\xvec=\ctwovec bm\text{.}\)
    Perhatikan bahwa vektor tak diketahui \(\xvec=\ctwovec bm\) menyatakan garis yang kita cari.
  4. Apakah sistem linear ini memiliki solusi? Bagaimana pertanyaan ini berkaitan dengan upaya Anda menggambar sebuah garis yang melalui ketiga titik di atas?
  5. Karena sistem ini tak konsisten, kita mengetahui bahwa \(\bvec\) tidak berada dalam ruang kolom \(\col(A)\text{.}\) Tentukan basis ortogonal bagi \(\col(A)\text{,}\) lalu gunakan basis itu untuk menentukan proyeksi ortogonal \(\widehat\bvec\) dari \(\bvec\) pada \(\col(A)\text{.}\)
  6. Karena \(\widehat\bvec\) berada dalam \(\col(A)\text{,}\) persamaan \(A\xvec = \widehat\bvec\) konsisten. Tentukan solusinya, \(\xvec = \ctwovec{b}{m}\text{,}\) lalu sketsakan garis \(y=b + mx\) pada Gambar 6.5.2. Garis ini disebut garis kecocokan terbaik.
Aktivitas ini menggambarkan gagasan yang melandasi teknik bernama kuadrat terkecil ortogonal, yang telah kita persiapkan sepanjang bab ini. Jika titik-titik data dinyatakan sebagai \((x_i, y_i)\text{,}\) kita membentuk matriks \(A\) dan vektor \(\bvec\) sebagai
\begin{equation*} A = \begin{bmatrix} 1 \amp x_1 \\ 1 \amp x_2 \\ 1 \amp x_3 \\ \end{bmatrix},\hspace{24pt} \bvec = \threevec{y_1}{y_2}{y_3}\text{.} \end{equation*}
Dengan vektor \(\xvec=\ctwovec bm\) menyatakan garis \(b+mx = y\text{,}\) kita melihat bahwa persamaan \(A\xvec=\bvec\) menyatakan garis yang melalui semua titik data. Dalam aktivitas kita, secara visual tampak jelas bahwa tidak ada garis semacam itu, selaras dengan fakta bahwa persamaan \(A\xvec=\bvec\) tak konsisten.
Ingat bahwa \(\bhat\text{,}\) yaitu proyeksi ortogonal \(\bvec\) pada \(\col(A)\text{,}\) merupakan vektor dalam \(\col(A)\) yang paling dekat dengan \(\bvec\text{.}\) Oleh karena itu, ketika kita menyelesaikan persamaan \(A\xvec=\bhat\text{,}\) kita mencari vektor \(\xvec\) sedemikian sehingga \(A\xvec = \threevec{b+mx_1}{b+mx_2}{b+mx_3}\) sedekat mungkin dengan \(\bvec=\threevec{y_1}{y_2}{y_3}\text{.}\) Mari kita pikirkan makna hal ini dalam konteks persoalan tersebut.
Selisihnya adalah \(\bvec-A\xvec = \threevec{y_1-(b+mx_1)}{y_2-(b+mx_2)}{y_3-(b+mx_3)}\text{,}\) sehingga kuadrat jarak antara \(A\xvec\) dan \(\bvec\) adalah
\begin{align*} \len{\bvec - A\xvec}^2 \amp =\\ \amp \left(y_1-(b+mx_1)\right)^2 + \left(y_2-(b+mx_2)\right)^2 + \left(y_3-(b+mx_3)\right)^2\text{.} \end{align*}
Pendekatan kita menentukan nilai \(b\) dan \(m\) yang membuat jumlah kuadrat ini sekecil mungkin; karena itulah persoalan ini disebut persoalan kuadrat terkecil.
Ketika garis yang ditentukan oleh vektor \(\xvec=\ctwovec bm\) digambar, besaran \(y_i - (b + mx_i)\) menyatakan jarak vertikal antara garis tersebut dan titik data \((x_i, y_i)\text{,}\) seperti diperlihatkan pada Gambar 6.5.5. Dengan sudut pandang ini, kuadrat jarak \(\len{\bvec-A\xvec}^2\) mengukur seberapa jauh garis yang ditentukan oleh vektor \(\xvec\) meleset dari titik-titik data. Solusi persoalan kuadrat terkecil adalah garis yang meleset dari titik-titik data sekecil mungkin.
Kontrol papan ketik untuk menjelajahi diagram
Tombol Tindakan
Enter, A Aktifkan penjelajahan dengan papan ketik
B Aktifkan penjelajahan dengan menu
Escape Keluar dari mode penjelajahan
Panah bawah Jelajahi tingkat berikutnya di bawah
Panah atas Jelajahi tingkat berikutnya di atas
Panah kanan Jelajahi elemen berikutnya pada tingkat ini
Panah kiri Jelajahi elemen sebelumnya pada tingkat ini
X Aktifkan atau nonaktifkan mode ahli
W Tampilkan rincian tambahan jika tersedia
Spasi Ulangi keluaran suara
M Aktifkan pembesaran bertahap
Koma Aktifkan pembesaran langsung
N Nonaktifkan pembesaran
Z Tampilkan atau sembunyikan subtitel
C Gilir pengaturan kontras
T Gunakan warna monokrom
L Ganti bahasa jika tersedia
K Hentikan suara saat ini
Y Hentikan keluaran suara
O Mulai atau hentikan sonifikasi
P Ulangi keluaran sonifikasi
Tiga titik data dari aktivitas sebelumnya dan garis yang diberikan sebagai solusi persoalan kuadrat terkecil. Terdapat ruas-ruas garis vertikal dari setiap titik data menuju garis kuadrat terkecil yang menunjukkan seberapa jauh garis tersebut meleset dari titik-titik itu.
Gambar 6.5.5. Solusi persoalan kuadrat terkecil dan jarak vertikal antara garis dengan titik-titik data.

Subbagian 6.5.3 Menyelesaikan persoalan kuadrat terkecil

Setelah melihat sebuah contoh mengenai tujuan yang hendak kita capai, mari kita tempatkan teknik ini dalam kerangka yang lebih umum.
Untuk sistem tak konsisten \(A\xvec = \bvec\text{,}\) kita mencari vektor \(\xvec\) yang meminimumkan jarak dari \(A\xvec\) ke \(\bvec\text{.}\) Dengan kata lain, \(\xvec\) memenuhi \(A\xvec = \widehat\bvec\text{,}\) dengan \(\bhat\) merupakan proyeksi ortogonal \(\bvec\) pada ruang kolom \(\col(A)\text{.}\) Kita mengetahui bahwa persamaan \(A\xvec=\bhat\) konsisten karena \(\bhat\) berada dalam \(\col(A)\text{,}\) dan solusinya tunggal jika kita mengasumsikan kolom-kolom \(A\) bebas linear.
Biasanya kita menyatakan solusi \(A\xvec = \bhat\) dengan \(\xhat\) dan menyebut vektor ini solusi hampiran kuadrat terkecil dari \(A\xvec=\bvec\text{,}\) untuk membedakannya dari solusi \(A\xvec=\bvec\) yang mungkin saja tidak ada.
Ada metode alternatif untuk menentukan \(\xhat\) tanpa terlebih dahulu menentukan proyeksi ortogonal \(\bhat\text{.}\) Ingat bahwa \(\bhat\) didefinisikan oleh fakta bahwa \(\widehat\bvec - \bvec\) ortogonal terhadap \(\col(A)\text{.}\) Dengan kata lain, \(\bhat-\bvec\) berada dalam komplemen ortogonal \(\col(A)^\perp\text{,}\) yang menurut Proposisi 6.2.10 sama dengan \(\nul(A^T)\text{.}\) Karena \(\bhat-\bvec\) berada dalam \(\nul(A^T)\text{,}\) diperoleh
\begin{equation*} A^T(\widehat\bvec-\bvec) = \zerovec\text{.} \end{equation*}
Karena solusi hampiran kuadrat terkecil adalah vektor \(\xhat\) sedemikian sehingga \(A\xhat = \bhat\text{,}\) kita dapat menata ulang persamaan ini untuk memperoleh
\begin{align*} A^T(A\xhat - \bvec) \amp = \zerovec\\ A^TA\xhat - A^T\bvec \amp = \zerovec\\ A^TA\xhat \amp = A^T\bvec\text{.} \end{align*}
Persamaan ini disebut persamaan normal, dan kita mempunyai proposisi berikut.

Contoh 6.5.7.

Perhatikan persamaan
\begin{equation*} \begin{bmatrix} 2 \amp 1 \\ 2 \amp 0 \\ -1 \amp 3 \\ \end{bmatrix} \xvec = \threevec{16}{-1}7 \end{equation*}
dengan matriks \(A\) dan vektor \(\bvec\text{.}\) Karena persamaan ini tak konsisten, kita akan menentukan solusi hampiran kuadrat terkecil \(\xhat\) dengan menyelesaikan persamaan normal \(A^TA\xhat = A^T\bvec\text{,}\) yang berbentuk
\begin{equation*} A^TA\xhat = \begin{bmatrix} 9 \amp -1 \\ -1 \amp 10 \\ \end{bmatrix} \xhat = \twovec{23}{37} = A^T\bvec \end{equation*}
dan memiliki solusi \(\xhat=\twovec34\text{.}\)

Kegiatan 6.5.3.

Laju bunyi jangkrik berkaitan dengan suhu luar ruangan, sebagaimana tercermin dalam sejumlah data percobaan yang akan kita kaji dalam aktivitas ini. Laju bunyi \(C\) dinyatakan dalam bunyi per detik, sedangkan suhu \(T\) dinyatakan dalam derajat Fahrenheit. Evaluasi sel berikut untuk memuat data:
Evaluasi sel ini juga menyediakan:
  • vektor chirps dan temps yang dibentuk dari kolom-kolom himpunan data.
  • perintah onesvec(n), yang membentuk vektor berdimensi \(n\) dengan semua entrinya bernilai satu.
  • Ingat bahwa Anda dapat membentuk matriks yang kolom-kolomnya adalah vektor v1 dan v2 dengan matrix([v1, v2]).T.
Kita ingin menyatakan hubungan ini dengan fungsi linear
\begin{equation*} \beta_0 + \beta_1 C = T\text{.} \end{equation*}
  1. Gunakan titik data pertama \((C_1,T_1)=(20.0,88.6)\) untuk menulis persamaan yang memuat \(\beta_0\) dan \(\beta_1\text{.}\)
  2. Misalkan peubah-peubah tak diketahui dinyatakan dengan vektor \(\xvec = \twovec{\beta_0}{\beta_1}\text{.}\) Gunakan 15 titik data untuk membentuk matriks \(A\) dan vektor \(\bvec\) sedemikian sehingga sistem linear \(A\xvec= \bvec\) menyatakan vektor tak diketahui \(\xvec\text{.}\)
  3. Tuliskan persamaan normal \(A^TA\xhat = A^T\bvec\text{;}\) dengan kata lain, tentukan matriks \(A^TA\) dan vektor \(A^T\bvec\text{.}\)
  4. Selesaikan persamaan normal untuk menentukan \(\xhat\text{,}\) yaitu solusi hampiran kuadrat terkecil bagi persamaan \(A\xvec=\bvec\text{.}\) Namai solusi Anda xhat karena x memiliki arti lain dalam Sage.
    Berapakah nilai \(\beta_0\) dan \(\beta_1\) yang Anda peroleh?
  5. Jika laju bunyi jangkrik adalah 22 bunyi per detik, berapakah prediksi Anda untuk suhunya?
    Dengan asumsi solusi Anda dinamai xhat, Anda dapat memplot data dan garis tersebut menggunakan sel di bawah ini.
Contoh ini memperlihatkan pendekatan yang disebut regresi linear, yaitu sekumpulan data dimodelkan menggunakan fungsi linear yang diperoleh dengan menyelesaikan persoalan kuadrat terkecil. Setelah memperoleh fungsi linear yang paling cocok dengan data, kita dapat membuat prediksi untuk keadaan yang belum muncul dalam data.
Jika fungsi tersebut akan kita gunakan untuk membuat prediksi, wajar untuk menanyakan seberapa besar keyakinan kita terhadap prediksi itu. Pertanyaan statistis ini mengarah pada teori yang kaya dan telah berkembang dengan baik
 1 
Sebagai contoh, lihat Gareth James, Daniela Witten, Trevor Hastie, Robert Tibshirani. An Introduction to Statistical Learning: with Applications in R. Springer, 2013.
, yang tidak akan kita bahas secara terperinci di sini. Namun, ada satu ukuran sederhana mengenai seberapa baik fungsi linear kita cocok dengan data, yaitu koefisien determinasi yang dinyatakan dengan \(R^2\text{.}\)
Kita telah melihat bahwa kuadrat jarak \(\len{\bvec-A\xvec}^2\) mengukur besarnya kegagalan garis melewati titik-titik data. Ketika garis dekat dengan titik-titik data, kita mengharapkan bilangan ini kecil. Namun, besar ukuran ini bergantung pada skala data. Sebagai contoh, dua garis pada Gambar 6.5.8 tampak sama-sama cocok dengan data, tetapi \(|\bvec-A\xhat|^2\) di sebelah kanan 100 kali lebih besar.
Tiga titik data dari aktivitas, garis kuadrat terkecil, dan ruas-ruas garis vertikal dari titik-titik data menuju garis yang menunjukkan seberapa jauh garis tersebut meleset dari titik-titik itu.
Diagram ini serupa dengan diagram sebelumnya, tetapi koordinat titik-titik data telah diskalakan dengan faktor \(10\text{.}\) Sebagai contoh, titik \((1,1)\) dalam himpunan data semula diskalakan menjadi \((10,10)\) dalam himpunan ini. Garis kuadrat terkecil dan ruas-ruas garis vertikal dari titik menuju garis juga ditampilkan. Pengamatan utamanya ialah kedua diagram tampak persis sama, kecuali skala pada sumbu.
Kontrol papan ketik untuk menjelajahi diagram
Tombol Tindakan
Enter, A Aktifkan penjelajahan dengan papan ketik
B Aktifkan penjelajahan dengan menu
Escape Keluar dari mode penjelajahan
Panah bawah Jelajahi tingkat berikutnya di bawah
Panah atas Jelajahi tingkat berikutnya di atas
Panah kanan Jelajahi elemen berikutnya pada tingkat ini
Panah kiri Jelajahi elemen sebelumnya pada tingkat ini
X Aktifkan atau nonaktifkan mode ahli
W Tampilkan rincian tambahan jika tersedia
Spasi Ulangi keluaran suara
M Aktifkan pembesaran bertahap
Koma Aktifkan pembesaran langsung
N Nonaktifkan pembesaran
Z Tampilkan atau sembunyikan subtitel
C Gilir pengaturan kontras
T Gunakan warna monokrom
L Ganti bahasa jika tersedia
K Hentikan suara saat ini
Y Hentikan keluaran suara
O Mulai atau hentikan sonifikasi
P Ulangi keluaran sonifikasi
Gambar 6.5.8. Garis-garis tampak sama-sama cocok meskipun \(\len{\bvec-A\xhat}^2\) berbeda dengan faktor 100.
Koefisien determinasi \(R^2\) didefinisikan dengan menormalkan \(|\bvec-A\xhat|^2\) agar tidak bergantung pada skala. Ingat bahwa pada Subbab 6.1 kita telah menjelaskan cara memusatkan sebuah vektor: untuk vektor \(\vvec\text{,}\) kita memperoleh \(\widetilde{\vvec}\) dengan mengurangkan rata-rata komponennya dari setiap komponen.

Definisi 6.5.9. Koefisien determinasi.

Koefisien determinasi adalah
\begin{equation*} R^2 = 1 - \frac{|\bvec - A\xhat|^2} {|\widetilde{\bvec}|^2}, \end{equation*}
dengan \(\widetilde{\bvec}\) merupakan vektor yang diperoleh dengan memusatkan \(\bvec\text{.}\)
Penjelasan yang lebih lengkap tentang definisi ini bergantung pada konsep varians, yang kita telaah dalam Latihan 6.5.7.12 dan bab berikutnya. Untuk saat ini, cukup diketahui bahwa \(0\leq R^2 \leq 1\text{,}\) dan semakin dekat \(R^2\) ke 1, semakin baik garis tersebut cocok dengan data. Dalam contoh awal yang digambarkan pada Gambar 6.5.8, kita memperoleh \(R^2 = 0.75\text{,}\) sedangkan dalam kajian laju bunyi jangkrik kita memperoleh \(R^2=0.69\text{.}\) Namun, menilai keyakinan terhadap prediksi yang dibuat dengan menyelesaikan persoalan kuadrat terkecil dapat memerlukan pertimbangan yang matang; mengandalkan nilai \(R^2\) saja merupakan tindakan yang naif.

Subbagian 6.5.4 Menggunakan faktorisasi \(QR\)

Seperti yang telah kita lihat, solusi hampiran kuadrat terkecil \(\xhat\) bagi \(A\xvec=\bvec\) dapat ditentukan dengan menyelesaikan persamaan normal \(A^TA\xhat = A^T\bvec\text{.}\) Strategi ini praktis untuk beberapa persoalan. Namun, pendekatan tersebut dapat bermasalah karena galat pembulatan kecil dapat menumpuk dan menghasilkan hasil akhir yang tidak akurat.
Seperti ditunjukkan oleh aktivitas berikutnya, terdapat metode lain untuk menentukan solusi hampiran kuadrat terkecil \(\xhat\) menggunakan faktorisasi \(QR\) dari matriks \(A\text{.}\) Metode ini lebih disukai karena lebih andal secara numerik.

Kegiatan 6.5.4.

  1. Misalkan kita ingin menentukan solusi hampiran kuadrat terkecil bagi persamaan \(A\xvec = \bvec\) dan mempunyai faktorisasi \(QR\text{,}\) yaitu \(A=QR\text{.}\) Jelaskan mengapa solusi hampiran kuadrat terkecil diberikan dengan menyelesaikan
    \begin{align*} A\xhat \amp = QQ^T\bvec \\\\ QR\xhat \amp = QQ^T\bvec \\ \end{align*}
  2. Kalikan kedua ruas persamaan kedua dengan \(Q^T\text{,}\) lalu jelaskan mengapa
    \begin{equation*} R\xhat = Q^T\bvec. \end{equation*}
    Karena \(R\) berbentuk segitiga atas, persamaan ini relatif mudah diselesaikan dengan substitusi mundur, seperti yang kita lihat dalam Subbab 5.1. Oleh karena itu, kita menulis solusi hampiran kuadrat terkecil sebagai
    \begin{equation*} \xhat = R^{-1}Q^T\bvec, \end{equation*}
    dan menerapkannya dalam konteks berikut.
  3. Rumus Brozak, yang digunakan untuk menghitung indeks lemak tubuh seseorang \(BFI\text{,}\) adalah
    \begin{equation*} BFI = 100 \left(\frac{4.57}{\rho} - 4.142\right) \end{equation*}
    dengan \(\rho\) menyatakan kerapatan tubuh seseorang dalam gram per sentimeter kubik. Namun, memperoleh pengukuran \(\rho\) yang akurat sulit karena orang tersebut harus dibenamkan ke dalam air dan volume air yang dipindahkan harus diukur. Sebagai gantinya, kita akan mengumpulkan beberapa ukuran tubuh lain yang lebih mudah diperoleh dan menggunakannya untuk memprediksi \(BFI\text{.}\)
    Sebagai contoh, misalkan kita mengambil 10 pasien dan mengukur berat badan \(w\) dalam pon, tinggi badan \(h\) dalam inci, lingkar perut \(a\) dalam sentimeter, lingkar pergelangan tangan \(r\) dalam sentimeter, lingkar leher \(n\) dalam sentimeter, dan \(BFI\text{.}\) Evaluasi sel berikut akan memuat dan menampilkan data.
    Selain itu, sel tersebut menyediakan:
    1. vektor weight, height, abdomen, wrist, neck, dan BFI yang dibentuk dari kolom-kolom himpunan data.
    2. perintah onesvec(n), yang mengembalikan vektor berdimensi \(n\) dengan semua entrinya bernilai satu.
    3. perintah QR(A), yang mengembalikan faktorisasi \(QR\) dari \(A\) sebagai Q, R = QR(A).
    4. perintah demean(v), yang mengembalikan vektor terpusat \(\widetilde{\vvec}\text{.}\)
    Kita ingin menentukan fungsi linear
    \begin{equation*} \beta_0 + \beta_1w + \beta_2h + \beta_3a + \beta_4r + \beta_5n = BFI \end{equation*}
    yang paling cocok dengan data.
    Gunakan titik data pertama untuk menulis persamaan bagi parameter \(\beta_0,\beta_1,\ldots,\beta_5\text{.}\)
  4. Jelaskan sistem linear \(A\xvec = \bvec\) untuk parameter tersebut. Secara lebih khusus, jelaskan bagaimana matriks \(A\) dan vektor \(\bvec\) dibentuk.
  5. Bentuk matriks \(A\) dan tentukan faktorisasi \(QR\)-nya dalam sel di bawah ini.
  6. Tentukan solusi hampiran kuadrat terkecil \(\xhat\) dengan menyelesaikan persamaan \(R\xhat = Q^T\bvec\text{.}\) Anda dapat menggunakan N(xhat) untuk menampilkan hampiran desimal vektor tersebut. Berapakah parameter \(\beta_0,\beta_1,\ldots,\beta_5\) yang paling cocok dengan data?
  7. Tentukan koefisien determinasi \(R^2\) untuk parameter Anda. Apa maknanya bagi kualitas kecocokan tersebut?
  8. Misalkan ukuran seseorang adalah: berat badan 190, tinggi badan 70, lingkar perut 90, lingkar pergelangan tangan 18, dan lingkar leher 35. Perkirakan \(BFI\) orang tersebut.
Sebagai rangkuman, kita telah melihat bahwa

Subbagian 6.5.5 Regresi polinomial

Dalam contoh-contoh sejauh ini, kita telah mencocokkan fungsi linear dengan suatu himpunan data. Namun, terkadang polinomial seperti fungsi kuadrat lebih sesuai. Ternyata teknik-teknik yang kita kembangkan dalam bagian ini tetap berguna, sebagaimana diperlihatkan oleh aktivitas berikut.

Kegiatan 6.5.5.

  1. Misalkan kita mempunyai himpunan data kecil yang memuat titik-titik \((0,2)\text{,}\) \((1,1)\text{,}\) \((2,3)\text{,}\) dan \((3,3)\text{,}\) seperti yang tampil ketika sel berikut dievaluasi.
    Selain memuat dan memplot data, evaluasi sel tersebut menyediakan perintah-perintah berikut:
    • Q, R = QR(A) mengembalikan faktorisasi \(QR\) dari \(A\text{.}\)
    • demean(v) mengembalikan vektor terpusat \(\widetilde{\vvec}\text{.}\)
    Mari kita cocokkan fungsi kuadrat berbentuk
    \begin{equation*} \beta_0 + \beta_1 x + \beta_2 x^2 = y \end{equation*}
    dengan himpunan data ini.
    Tuliskan empat persamaan, satu untuk setiap titik data, yang menyatakan koefisien \(\beta_0\text{,}\) \(\beta_1\text{,}\) dan \(\beta_2\text{.}\)
  2. Nyatakan keempat persamaan ini sebagai sistem linear \(A\xvec = \bvec\text{,}\) dengan \(\xvec = \threevec{\beta_0}{\beta_1}{\beta_2}\text{.}\)
    Tentukan faktorisasi \(QR\) dari \(A\text{,}\) lalu gunakan faktorisasi itu untuk menentukan solusi hampiran kuadrat terkecil \(\xhat\text{.}\)
  3. Gunakan parameter \(\beta_0\text{,}\) \(\beta_1\text{,}\) dan \(\beta_2\) yang Anda peroleh untuk menulis fungsi kuadrat yang cocok dengan data. Anda dapat memplot fungsi ini bersama data dengan memasukkan fungsi Anda pada tempat yang ditunjukkan di bawah ini.
  4. Berapakah nilai \(y\) yang Anda prediksi ketika \(x=1.5\text{?}\)
  5. Tentukan koefisien determinasi \(R^2\) untuk fungsi kuadrat. Apa yang ditunjukkannya tentang kualitas kecocokan?
  6. Sekarang cocokkan polinomial kubik berbentuk
    \begin{equation*} \beta_0 + \beta_1x + \beta_2 x^2 + \beta_3x^3 = y \end{equation*}
    dengan himpunan data ini.
  7. Tentukan koefisien determinasi \(R^2\) untuk fungsi kubik. Apa yang ditunjukkannya tentang kualitas kecocokan?
  8. Apa yang Anda perhatikan ketika memplot fungsi kubik bersama data? Bagaimana hal itu mencerminkan nilai \(R^2\) yang Anda peroleh?
Matriks \(A\) yang Anda bentuk dalam aktivitas sebelumnya ketika mencocokkan fungsi kuadrat dan kubik dengan himpunan data mempunyai bentuk khusus. Secara khusus, jika titik-titik data diberi label \((x_i, y_i)\) dan kita mencari polinomial berderajat \(k\text{,}\) maka
\begin{equation*} A = \begin{bmatrix} 1 \amp x_1 \amp x_1^2 \amp \ldots \amp x_1^k \\ 1 \amp x_2 \amp x_2^2 \amp \ldots \amp x_2^k \\ \vdots \amp \vdots \amp \vdots \amp \ddots \amp \vdots \\ 1 \amp x_m \amp x_m^2 \amp \ldots \amp x_m^k \\ \end{bmatrix}. \end{equation*}
Matriks ini disebut matriks Vandermonde berderajat \(k\text{.}\)

Kegiatan 6.5.6.

Aktivitas ini menelaah himpunan data mengenai es laut Arktik yang berasal dari Sustainability Math.
Evaluasi sel di bawah ini akan memplot luas es laut Arktik, dalam jutaan kilometer persegi, selama dua belas bulan pada tahun 2012.
Selain itu, Anda mempunyai akses ke beberapa variabel dan perintah khusus:
  • month adalah vektor nilai bulan dan ice adalah vektor nilai es laut dari tabel di atas.
  • vandermonde(x, k) membentuk matriks Vandermonde berderajat \(k\) menggunakan titik-titik dalam vektor x.
  • Q, R = QR(A) memberikan faktorisasi \(QR\) dari \(A\text{.}\)
  • demean(v) mengembalikan vektor terpusat \(\widetilde{\vvec}\text{.}\)
  1. Tentukan vektor \(\xhat\text{,}\) yaitu solusi hampiran kuadrat terkecil bagi sistem linear yang dihasilkan ketika polinomial berderajat 5 dicocokkan dengan data.
  2. Jika hasil Anda disimpan dalam variabel xhat, Anda dapat memplot polinomial dan data bersama-sama menggunakan sel berikut.
  3. Tentukan koefisien determinasi \(R^2\) untuk kecocokan polinomial ini.
  4. Ulangi langkah-langkah ini untuk mencocokkan polinomial berderajat 8 dengan data, memplot polinomial bersama data, dan menentukan \(R^2\text{.}\)
  5. Ulangi sekali lagi dengan mencocokkan polinomial berderajat 11 dengan data, membuat plot, dan menentukan \(R^2\text{.}\)
    Memang benar bahwa polinomial berderajat lebih tinggi lebih cocok dengan data, sebagaimana terlihat dari meningkatnya nilai \(R^2\text{,}\) tetapi hal itu tidak selalu baik. Sebagai contoh, ketika \(k=11\text{,}\) Anda mungkin memperhatikan bahwa grafik polinomial berombak sedikit lebih banyak daripada yang kita harapkan. Dalam hal ini, polinomial tersebut berusaha terlalu keras mencocokkan data, yang biasanya mengandung ketidakpastian, terutama jika diperoleh dari pengukuran. Galat yang terkandung dalam data disebut derau, dan kehadirannya berarti kita tidak patut mengharapkan polinomial cocok sempurna dengan data. Ketika kita memilih polinomial yang derajatnya terlalu tinggi, derau diberi bobot terlalu besar dalam model, sehingga menimbulkan perilaku yang tidak diinginkan, seperti gelombang pada grafik.
    Mencocokkan data dengan polinomial yang derajatnya terlalu tinggi disebut pencocokan berlebih, suatu fenomena yang dapat muncul dalam banyak aplikasi pembelajaran mesin. Secara umum, kita ingin memilih \(k\) yang cukup besar untuk menangkap ciri penting data, tetapi tidak begitu besar sehingga terjadi pencocokan berlebih dan derau terbangun ke dalam model. Ada berbagai cara menentukan nilai optimal \(k\text{,}\) tetapi hal itu tidak akan kita bahas di sini.
  6. Dengan memilih nilai \(k\) yang wajar, perkirakan luas es laut Arktik pada bulan ke-6,5, kira-kira saat titik balik matahari musim panas.

Subbagian 6.5.6 Rangkuman

Bagian ini memperkenalkan beberapa jenis persoalan kuadrat terkecil dan suatu kerangka untuk menanganinya.
  • Untuk sistem tak konsisten \(A\xvec=\bvec\text{,}\) kita menentukan solusi hampiran kuadrat terkecil \(\xhat\) dengan mensyaratkan agar \(A\xhat\) sedekat mungkin dengan \(\bvec\text{.}\) Dengan kata lain, \(A\xhat = \bhat\text{,}\) dengan \(\bhat\) merupakan proyeksi ortogonal \(\bvec\) pada \(\col(A)\text{.}\)
  • Salah satu cara menentukan \(\xhat\) adalah menyelesaikan persamaan normal \(A^TA\xhat = A^T\bvec.\) Metode ini tidak kita utamakan karena dapat menimbulkan masalah numerik.
  • Cara kedua menentukan \(\xhat\) menggunakan faktorisasi \(QR\) dari \(A\text{.}\) Jika \(A=QR\text{,}\) maka \(\xhat = R^{-1}Q^T\bvec\text{,}\) dan menentukan \(R^{-1}\) layak secara komputasional karena \(R\) berbentuk segitiga atas.
  • Teknik ini dapat diterapkan secara luas dan berguna untuk memodelkan data. Dalam bagian ini kita melihat contoh ketika fungsi linear dari beberapa variabel masukan dan polinomial memberikan model yang efektif bagi himpunan data yang berbeda.
  • Salah satu ukuran sederhana bagi kualitas kecocokan adalah koefisien determinasi \(R^2\text{,}\) meskipun penerapan nyata memerlukan pertimbangan tambahan.

Latihan 6.5.7 Latihan

Evaluasi sel berikut memuat beberapa perintah yang akan berguna dalam latihan-latihan berikut. Secara khusus, terdapat perintah:
  • QR(A), yang mengembalikan faktorisasi \(QR\) dari A sebagai Q, R = QR(A),
  • onesvec(n), yang mengembalikan vektor berdimensi \(n\) dengan semua entrinya bernilai 1,
  • demean(v), yang memusatkan vektor v,
  • vandermonde(x, k), yang mengembalikan matriks Vandermonde berderajat \(k\) yang dibentuk dari komponen-komponen vektor x, dan
  • plot_model(xhat, data), yang memplot data dan model xhat.

1.

Misalkan kita menulis sistem linear
\begin{equation*} \begin{bmatrix} 1 \amp -1 \\ 2 \amp -1 \\ -1 \amp 3 \end{bmatrix} \xvec = \threevec{-8}5{-10} \end{equation*}
sebagai \(A\xvec=\bvec\text{.}\)
  1. Tentukan basis ortogonal bagi \(\col(A)\text{.}\)
  2. Tentukan \(\bhat\text{,}\) yaitu proyeksi ortogonal \(\bvec\) pada \(\col(A)\text{.}\)
  3. Tentukan solusi sistem linear \(A\xvec = \bhat\text{.}\)

2.

Perhatikan data dalam Tabel 6.5.11.
Tabel 6.5.11. Himpunan data dengan empat titik.
\(x\) \(y\)
1 1
2 1
3 1
4 2
  1. Susun sistem linear \(A\xvec=\bvec\) yang menyatakan garis \(b + mx = y\) yang melalui titik-titik tersebut.
  2. Tuliskan persamaan normal yang menyatakan solusi hampiran kuadrat terkecil bagi \(A\xvec=\bvec\text{.}\)
  3. Tentukan solusi hampiran kuadrat terkecil \(\xhat\text{,}\) lalu plot data dan garis yang dihasilkan.
  4. Berapakah nilai \(y\) yang Anda prediksi ketika \(x=3.5\text{?}\)
  5. Tentukan koefisien determinasi \(R^2\text{.}\)

3.

Perhatikan empat titik dalam Tabel 6.5.11.
  1. Susun sistem linear \(A\xvec = \bvec\) yang menyatakan fungsi kuadrat
    \begin{equation*} \beta_0+\beta_1x+\beta_2x^2 = y \end{equation*}
    yang melalui titik-titik tersebut.
  2. Gunakan faktorisasi \(QR\) untuk menentukan solusi hampiran kuadrat terkecil \(\xhat\text{,}\) lalu plot data dan grafik fungsi kuadrat yang dihasilkan.
  3. Berapakah nilai \(y\) yang Anda prediksi ketika \(x=3.5\text{?}\)
  4. Tentukan koefisien determinasi \(R^2\text{.}\)

4.

Perhatikan data dalam Tabel 6.5.12.
Tabel 6.5.12. Himpunan data sederhana
\(x_1\) \(x_2\) \(y\)
1 1 4.2
1 2 3.3
2 1 5.9
2 2 5.1
3 2 7.5
3 3 6.3
  1. Susun sistem linear \(A\xvec = \bvec\) yang menyatakan hubungan
    \begin{equation*} \beta_0 + \beta_1 x_1 + \beta_2 x_2 = y. \end{equation*}
  2. Tentukan solusi hampiran kuadrat terkecil \(\xhat\text{.}\)
  3. Berapakah nilai \(y\) yang Anda prediksi ketika \(x_1 = 2.4\) dan \(x_2=2.9\text{?}\)
  4. Tentukan koefisien determinasi \(R^2\text{.}\)

5.

Tentukan apakah pernyataan-pernyataan berikut benar atau salah, serta jelaskan penalaran Anda.
  1. Jika \(A\xvec=\bvec\) konsisten, maka \(\xhat\) merupakan solusi \(A\xvec=\bvec\text{.}\)
  2. Jika \(R^2=1\text{,}\) maka solusi hampiran kuadrat terkecil \(\xhat\) juga merupakan solusi persamaan semula \(A\xvec=\bvec\text{.}\)
  3. Dengan faktorisasi \(QR\text{,}\) yaitu \(A=QR\text{,}\) kita mempunyai \(A\xhat=Q^TQ\bvec\text{.}\)
  4. Faktorisasi \(QR\) menyediakan metode untuk menentukan solusi hampiran kuadrat terkecil bagi \(A\xvec=\bvec\) yang lebih andal daripada menyelesaikan persamaan normal.
  5. Solusi \(AA^T\xvec = A\bvec\) merupakan solusi hampiran kuadrat terkecil bagi \(A\xvec = \bvec\text{.}\)

6.

Jelaskan jawaban Anda atas pertanyaan-pertanyaan berikut.
  1. Jika \(\xhat=\zerovec\text{,}\) apa artinya bagi vektor \(\bvec\text{?}\)
  2. Jika kolom-kolom \(A\) ortonormal, bagaimana Anda dapat dengan mudah menentukan solusi hampiran kuadrat terkecil bagi \(A\xvec=\bvec\text{?}\)

7.

Sel berikut memuat data yang menunjukkan jumlah penduduk Bangladesh yang hidup tanpa listrik selama 27 tahun. Sel ini juga mendefinisikan vektor year, yang mencatat tahun-tahun dalam himpunan data, dan people, yang mencatat jumlah penduduk.
  1. Misalkan kita ingin menulis
    \begin{equation*} N = \beta_0 + \beta_1 t \end{equation*}
    dengan \(t\) menyatakan tahun dan \(N\) menyatakan jumlah penduduk. Bentuk matriks \(A\) dan vektor \(\bvec\) sedemikian sehingga sistem linear \(A\xvec=\bvec\) menyatakan vektor \(\xvec=\twovec{\beta_0}{\beta_1}\text{.}\)
  2. Dengan menggunakan faktorisasi \(QR\) dari \(A\text{,}\) tentukan nilai \(\beta_0\) dan \(\beta_1\) dalam solusi hampiran kuadrat terkecil \(\xhat\text{.}\)
  3. Berapakah koefisien determinasi \(R^2\text{,}\) dan apa yang dinyatakannya tentang kualitas hampiran tersebut?
  4. Berapakah prediksi Anda untuk jumlah penduduk yang hidup tanpa listrik pada tahun 1985?
  5. Perkirakan tahun ketika tidak ada lagi penduduk yang hidup tanpa listrik.

8.

Persoalan ini berkenaan dengan himpunan data yang menggambarkan planet-planet dalam Tata Surya kita. Untuk setiap planet, kita mempunyai panjang sumbu semi-mayor \(L\text{,}\) yang pada dasarnya adalah jarak planet ke Matahari dalam SA (satuan astronomi), dan periode \(P\text{,}\) yaitu waktu dalam tahun yang diperlukan untuk menyelesaikan satu orbit mengelilingi Matahari.
Kita ingin memodelkan data ini menggunakan fungsi \(P = CL^r\text{,}\) dengan \(C\) dan \(r\) sebagai parameter yang perlu ditentukan. Karena fungsi ini tidak linear, kita mengubah hubungan tersebut dengan mengambil logaritma natural kedua ruas untuk memperoleh
\begin{equation*} \ln(P) = \ln(C) + r\ln(L). \end{equation*}
Evaluasi sel berikut memuat himpunan data dan mendefinisikan dua vektor: logaxis, yang komponennya adalah \(\ln(L)\text{,}\) dan logperiod, yang komponennya adalah \(\ln(P)\text{.}\)
  1. Bentuk matriks \(A\) dan vektor \(\bvec\) sedemikian sehingga solusi \(A\xvec=\bvec\) adalah vektor \(\xvec=\ctwovec{\ln(C)}r\text{.}\)
  2. Tentukan solusi hampiran kuadrat terkecil \(\xhat\text{.}\) Nilai berapakah yang diperoleh untuk \(C\) dan \(r\text{?}\)
  3. Tentukan koefisien determinasi \(R^2\text{.}\) Apa yang dinyatakannya tentang kualitas hampiran tersebut?
  4. Misalkan orbit sebuah asteroid mempunyai sumbu semi-mayor dengan panjang \(L=4.0\) SA. Perkirakan periode \(P\) orbit asteroid itu.
  5. Komet Halley mempunyai periode \(P=75\) tahun. Perkirakan panjang sumbu semi-mayornya.

9.

Evaluasi sel berikut memuat himpunan data yang menggambarkan suhu atmosfer Bumi pada berbagai ketinggian. Terdapat pula dua vektor: altitude, yang dinyatakan dalam kilometer, dan temperature, dalam derajat Celsius.
  1. Jelaskan cara membentuk matriks \(A\) dan vektor \(\bvec\) sedemikian sehingga sistem linear \(A\xvec=\bvec\) menyatakan polinomial berderajat \(k\) yang cocok dengan data.
  2. Setelah memilih nilai \(k\text{,}\) bentuk matriks \(A\) dan vektor \(\bvec\text{,}\) lalu tentukan solusi hampiran kuadrat terkecil \(\xhat\text{.}\)
  3. Plot polinomial dan data menggunakan plot_model(xhat, data).
  4. Sekarang amati apa yang terjadi ketika Anda mengubah derajat polinomial \(k\text{.}\) Pilih nilai \(k\) yang sesuai, yang tampaknya menangkap ciri terpenting data sekaligus menghindari pencocokan berlebih, lalu jelaskan pilihan Anda.
  5. Gunakan nilai \(k\) pilihan Anda untuk memperkirakan suhu pada ketinggian 55 kilometer.

10.

Sel berikut memuat data mengenai 1057 rumah dalam suatu pasar properti. Untuk setiap rumah, dicatat luas hunian dalam kaki persegi, luas lahan dalam acre, usia dalam tahun, dan harga dalam dolar. Sel tersebut juga mendefinisikan variabel area, size, age, dan price.
Kita akan menggunakan regresi linear untuk memprediksi harga rumah berdasarkan luas hunian, luas lahan, dan usianya:
\begin{equation*} \beta_0 + \beta_1~\text{Luas Hunian} + \beta_2~\text{Luas Lahan} + \beta_3~\text{Usia} = \text{Harga}. \end{equation*}
  1. Gunakan faktorisasi \(QR\) untuk menentukan solusi hampiran kuadrat terkecil \(\xhat\text{.}\)
  2. Bahas arti tanda \(\beta_1\text{,}\) \(\beta_2\text{,}\) dan \(\beta_3\text{.}\)
  3. Jika dua rumah identik kecuali usianya berbeda satu tahun, bagaimana Anda memprediksi perbandingan harga keduanya?
  4. Tentukan koefisien determinasi \(R^2\text{.}\) Apa yang dinyatakannya tentang kualitas kecocokan?
  5. Prediksi harga rumah yang luas huniannya 2000 kaki persegi, luas lahannya 1,5 acre, dan usianya 50 tahun.

11.

Kita telah mengamati bahwa jika kolom-kolom \(A\) bebas linear, maka persamaan \(A\xvec=\bvec\) mempunyai solusi hampiran kuadrat terkecil tunggal karena persamaan \(A\xhat=\bhat\) mempunyai solusi tunggal. Kita juga menyatakan bahwa \(\xhat\) merupakan solusi tunggal persamaan normal \(A^TA\xhat = A^T\bvec\) tanpa menjelaskan mengapa persamaan ini mempunyai solusi tunggal. Latihan ini memberikan penjelasannya.
Dengan mengasumsikan bahwa kolom-kolom \(A\) bebas linear, kita ingin menyimpulkan bahwa persamaan \(A^TA\xhat=A^T\bvec\) mempunyai solusi tunggal.
  1. Misalkan \(\xvec\) adalah vektor yang memenuhi \(A^TA\xvec = \zerovec\text{.}\) Jelaskan mengapa argumen berikut sah dan memungkinkan kita menyimpulkan bahwa \(A\xvec = \zerovec\text{.}\)
    \begin{equation*} \begin{aligned} A^TA\xvec \amp = \zerovec \\ \xvec\cdot A^TA\xvec \amp = \xvec\cdot\zerovec = 0 \\ (A\xvec)\cdot(A\xvec) \amp = 0 \\ \len{A\xvec}^2 \amp = 0. \\ \end{aligned} \end{equation*}
    Dengan kata lain, jika \(A^TA\xvec = \zerovec\text{,}\) kita mengetahui bahwa \(A\xvec = \zerovec\text{.}\)
  2. Jika kolom-kolom \(A\) bebas linear dan \(A\xvec = \zerovec\text{,}\) apa yang kita ketahui tentang vektor \(\xvec\text{?}\)
  3. Jelaskan mengapa \(A^TA\xvec = \zerovec\) hanya mungkin terjadi ketika \(\xvec = \zerovec\text{.}\)
  4. Dengan mengasumsikan bahwa kolom-kolom \(A\) bebas linear, jelaskan mengapa \(A^TA\xhat=A^T\bvec\) mempunyai solusi tunggal.

12.

Persoalan ini membahas makna koefisien determinasi \(R^2\) dan hubungannya dengan varians, topik yang muncul dalam bagian berikutnya. Sepanjang persoalan ini, kita meninjau sistem linear \(A\xvec=\bvec\) dan solusi hampiran kuadrat terkecil \(\xhat\text{,}\) dengan \(A\xhat=\bhat\text{.}\) Misalkan \(A\) adalah matriks \(m\times n\text{,}\) dan vektor berdimensi \(m\) dinyatakan dengan \(\onevec = \fourvec11{\vdots}1\text{.}\)
  1. Jelaskan mengapa \(\bbar\text{,}\) yaitu rata-rata komponen \(\bvec\text{,}\) dapat ditentukan sebagai hasil kali titik
    \begin{equation*} \bbar = \frac 1m \bvec\cdot\onevec. \end{equation*}
  2. Dalam contoh-contoh yang telah kita lihat dalam bagian ini, jelaskan mengapa \(\onevec\) berada dalam \(\col(A)\text{.}\)
  3. Jika kita menulis \(\bvec = \bhat + \bvec^\perp\text{,}\) jelaskan mengapa
    \begin{equation*} \bvec^\perp\cdot\onevec = 0 \end{equation*}
    dan karena itu mengapa rata-rata komponen \(\bvec^\perp\) adalah nol.
  4. Varians vektor berdimensi \(m\text{,}\) yaitu \(\vvec\text{,}\) adalah \(\var(\vvec) = \frac1m \len{\widetilde{\vvec}}^2\text{,}\) dengan \(\widetilde{\vvec}\) adalah vektor yang diperoleh dengan memusatkan \(\vvec\text{.}\)
    Jelaskan mengapa
    \begin{equation*} \var(\bvec) = \var(\bhat) + \var(\bvec^\perp). \end{equation*}
  5. Jelaskan mengapa
    \begin{equation*} \frac{\len{\bvec - A\xhat}^2}{\len{\widetilde{\bvec}}^2} = \frac{\var(\bvec^\perp)}{\var(\bvec)} \end{equation*}
    dan dengan demikian
    \begin{equation*} R^2 = \frac{\var(\bhat)}{\var(\bvec)} = \frac{\var(A\xhat)}{\var(\bvec)}. \end{equation*}
    Pernyataan-pernyataan ini menunjukkan mengapa \(R^2\) terkadang dikatakan mengukur “fraksi varians yang dijelaskan” oleh fungsi yang kita gunakan untuk mencocokkan data. Seperti terlihat dalam latihan sebelumnya, mungkin terdapat ciri lain yang tidak dicatat dalam himpunan data tetapi memengaruhi besaran yang ingin kita prediksi.
  6. Jelaskan mengapa \(0\leq R^2 \leq 1\text{.}\)