Cara menemukan residu dan plot mereka

Saya sudah diberi data

x = c(21,34,6,47,10,49,23,32,12,16,29,49,28,8,57,9,31,10,21,26,31,52,21,8,18,5,18,26,27,26,32,2,59,58,19,14,16,9,23,28,34,70,69,54,39,9,21,54,26)
y = c(47,76,33,78,62,78,33,64,83,67,61,85,46,53,55,71,59,41,82,56,39,89,31,43,29,55, 
     81,82,82,85,59,74,80,88,29,58,71,60,86,91,72,89,80,84,54,71,75,84,79)

Bagaimana saya bisa mendapatkan sisa dan plot versus ? Dan bagaimana saya bisa menguji apakah residu tampaknya mendekati normal? $x$

Saya tidak yakin apakah saya melakukan fit linear asli dengan benar karena saya mendapatkan persamaan tetapi catatan kuliah mengatakan bahwa garis regresi linier harus dalam bentuk . $y=6.9x-5.5$ $y_i=\beta_0+\beta_1x+\epsilon$

r regression tamu
sumber

Paket mana yang Anda gunakan? Misalnya fungsi 'regress' Matlab mengembalikan residu sebagai output dan Anda dapat membuat grafik menggunakan histogram

BGreene

Saya menggunakan Sagemath. Saya juga dapat menggunakan R melalui itu tetapi saya memiliki sedikit pengalaman tentang itu.

tamu

y = a + k x

$y = a +k x$

E [Y | X] = a + k X

$E[Y|X] = a+ k X$

Y = a + k X + ϵ

$Y = a + k X + \epsilon$

ϵ

$\epsilon$

\hat{β_{0}} = - 5.5

$\hat{\beta_0} = -5.5$

\hat{β_{1}} = 6.9

$\hat{\beta_1} = 6.9$

r_{i} = y_{y} - {\hat{y}}_{i} = y_{i} - (- 5.5 + 6.9 x_{i})

$r_i = y_y-\hat{y}_i = y_i - (-5.5 + 6.9 x_i)$

Jawaban:

EDIT: Anda memiliki Rtag, tetapi kemudian dalam komentar mengatakan Anda tidak tahu banyak tentang itu. Ini adalah Rkode. Saya tidak tahu apa-apa tentang Sage. Akhiri edit

Kamu bisa melakukan ini

x = c(21,34,6,47,10,49,23,32,12,16,29,49,28,8,57,9,31,10,21,
      26,31,52,21,8,18,5,18,26,27,26,32,2,59,58,19,14,16,9,23,
      28,34,70,69,54,39,9,21,54,26)
y = c(47,76,33,78,62,78,33,64,83,67,61,85,46,53,55,71,59,41,82,
      56,39,89,31,43,29,55, 81,82,82,85,59,74,80,88,29,58,71,60,
      86,91,72,89,80,84,54,71,75,84,79)

m1 <- lm(y~x)  #Create a linear model
resid(m1) #List of residuals
plot(density(resid(m1))) #A density plot
qqnorm(resid(m1)) # A quantile normal plot - good for checking normality
qqline(resid(m1))

Peter Flom - Pasang kembali Monica
sumber

+1 @guest, kode di atas adalah untuk R, yang tersedia secara bebas

BGreene

Baik. Jadi saya melihat gambar dengan caption density.default (x = resid (m1)). Haruskah kode ini menghasilkan dua grafik? Dan haruskah saya memeriksa dari grafik apakah residu tampaknya mendekati normal?

Tamu

Kode akan menampilkan dua grafik - satu adalah plot kerapatan (apakah terlihat seperti lonceng?) Yang lain adalah plot kuantil; jika residu normal, semua poin akan berada pada garis lurus.

Peter Flom - Reinstate Monica

Baik. Kode berfungsi jika Anda mengubah baris terakhir ke plot (qqnorm (resid (m1)))) dan plot (qqline (resid (m1)))). Jadi saya pikir residu tidak memenuhi distribusi normal karena ada titik lebih jauh di bawah garis daripada di atas garis. Apakah ada kriteria numerik untuk memeriksa normalitas?

Tamu