08744552 009588217 009284425 7 009518823 003923118

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: t$x2) ## the hat-matrix H = X %*% solve(t(X)%*%X) %*% t(X) ## leverages lev = rep(0,32) for (i in 1:32) {lev[i] = H[i,i]} lev [1] 0.08744552 0.09588217 0.09284425 [7] 0.09518823 0.03923118 0.03787039 [13] 0.04868743 0.08641369 0.12081830 [19] 0.08064726 0.12424110 0.08628552 [25] 0.24675474 0.08234239 0.28575991 [31] 0.44331534 0.08285421 0.03427419 0.08540308 0.04687672 0.04090267 0.11688367 0.09043677 0.15283620 0.33847490 0.23824594 0.16245986 0.17071980 0.07005945 0.09805986 0.07618809 0.14159717 OR > influence(clocks.fit2)$hat 1 2 3 0.08744552 0.09588217 0.09284425 8 9 10 0.03923118 0.03787039 0.08540308 15 16 17 0.12081830 0.04687672 0.33847490 22 23 24 0.04090267 0.23824594 0.07618809 29 30 31 0.16245986 0.14159717 0.44331534 > sum(lev) [1] 4 4 0.03427419 11 0.15283620 18 0.09805986 25 0.24675474 32 0.08285421 5 0.09043677 12 0.07005945 19 0.08064726 26 0.08234239 6 0.17071980 13 0.04868743 20 0.12424110 27 0.28575991 7 0.09518823 14 0.08641369 21 0.08628552 28 0.11688367...
View Full Document

Ask a homework question - tutors are online