Kalıntılar nasıl bulunur ve çizilir

Bana veri verildi

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)

Kalıntıları nasıl elde edebilirim ve nasıl çizebilirim ? Ve artıkların yaklaşık normal olup olmadığını nasıl test edebilirim? $x$

denklemini aldığım için orijinal doğrusal uyumu doğru şekilde yapıp yapmadığımdan emin değilim ama ders notları doğrusal regresyon çizgisinin biçiminde olması gerektiğini söylüyor . $y=6.9x-5.5$ $y_i=\beta_0+\beta_1x+\epsilon$

r regression

— misafir
kaynak

Hangi paketi kullanıyorsunuz? Örneğin Matlab'ın 'regress' fonksiyonu artıkları çıktı olarak döndürür ve histogram kullanarak grafik

— çizebilirsiniz

Sagemath kullanıyorum. Ayrıca R'yi de kullanabilirim ama çok az tecrübem var.

— misafir

Burada bulunan 2 denklemle ilgili olarak. (Doğrusal fonksiyonu olarak) regresyon çizgisi biçimi ise sonra doğrusal model ve bu da bir hata terimlerini kullanarak burada sıfır beklentili bir hata terimidir. Bu iki denklemin birbirine uyduğu duygusudur.

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$

— Ric

Elindeki denklem olduğu ile notlarınıza belirtilen formun ve . sadece

\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)$

— Glen_b-Monica

DÜZENLEME: Bir Retiketiniz var, ancak bir yorumda bunun hakkında fazla bir şey bilmediğinizi söylüyor. Bu Rkod. Sage hakkında hiçbir şey bilmiyorum. Düzenlemeyi bitir

Bunu yapabilirsiniz

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 - Monica'yı eski durumuna döndürün
kaynak

+1 @guest, yukarıdaki kod serbestçe kullanılabilen R içindir

— BGreene

Tamam. Bu yüzden resmi altyazı yoğunluğu ile gördüm. Varsayılan (x = kalıntı (m1)). Bu kod iki grafik çıkarmalı mı? Grafikten baktığımda artıkların normal görünüp görünmediğini kontrol etmeli miyim?

— misafir

Kod iki grafik çıkarır - biri yoğunluk grafiği (çan şeklinde mi görünüyor?) Diğeri ise kantil grafiktir; artıklar tamamen normal olsaydı, noktaların hepsi düz çizgide uzanırdı.

— Peter Flom - Monica'yı eski durumuna döndürün

Sağ. Son satırları çizim (qqnorm (artık (m1))) ve çizim (qqline (artık (m1))) olarak değiştirirseniz kod çalışır. Sanırım artıklar çizginin altında, çizginin üzerinde olduğundan daha fazla nokta olduğu için normal dağılımı tatmin etmiyor. Normalliği kontrol etmek için herhangi bir sayısal kriter var mı?

— misafir