# Regresyon analizi
Regresyon analizi, iki ya da daha çok nicel değişken arasındaki ilişkiyi ölçmek için kullanılan analiz metodudur. Eğer tek bir değişken kullanılarak analiz yapılıyorsa buna tek değişkenli regresyon, birden çok değişken kullanılıyorsa çok değişkenli regresyon analizi olarak isimlendirilir. Regresyon analizi ile değişkenler arasındaki ilişkinin varlığı, eğer ilişki var ise bunun gücü hakkında bilgi edinilebilir. Regresyon terimi için öz Türkçe olarak bağlanım sözcüğü kullanılması teklif edilmiş ise de Türk ekonometriciler arasında bu kullanım yaygın değildir.
Örneğin, bir ziraat mühendisi buğday verimi ve gübre miktarı arasındaki ilişkiyi, bir mühendis, basınç ve sıcaklık, bir ekonomist gelir düzeyi ve tüketim harcamaları, bir eğitimci öğrencilerin devamsızlık gösterdiği gün sayıları ve başarı dereceleri arasındaki ilişkiyi bilmek isteyebilir. Regresyon, iki (ya da daha çok) değişken arasındaki doğrusal ilişkinin fonksiyonel şeklini, biri bağımlı diğeri bağımsız değişken olarak bir doğru denklemi olarak göstermekle kalmaz, değişkenlerden birinin değeri bilindiğinde diğeri hakkında kestirim yapılmasını sağlar. Genellikle bu iki (veya çok) değişkenlerin hepsinin niceliksel ölçekli olması zorunluluğu vardır.
Regresyonda, değişkenlerden biri bağımlı diğerleri bağımsız değişken olmalıdır. Buradaki mantık eşitliğin solunda yer alan değişkenin sağında yer alan değişkenlerden etkilenmesidir. Sağda yer alan değişkenlerse diğer değişkenlerden etkilenmemektedir. Burada etkilenmemek matematiksel anlamda bu değişkenleri bir doğrusal denkleme koyduğumuzda etki yapması anlamındadır. "Çoklu doğrusallık", "ardışık bağımlılık" sorunları kastedilmemektedir.


## Regresyon analizinin tarihi.
Regresyon yönteminin ilk şekli en küçük kareler prensibidir ve ilk olarak Adrien Marie Legendre tarafından 1805 yılında ortaya atılmıştır. Hemen takiben 1809 yılında C.F. Gauss aynı yöntemi açıklamıştır. "En küçük kareler" terimi Legendre tarafından "moindres carrés" olarak kullanılmış, ancak Gauss aynı yöntemi 1795 yılından beri kullandığını iddia etmiştir. Legendre ve Gauss bu yöntemi astronomik gözlemlerden uydularının güneş etrafındaki yörüngelerini tespit etmek için kullanırken ortaya çıkartmışlardır. 1748 yılında Eüler'in aynı problem üzerinde uğraştığı, fakat başarı sağlayamadığı bilinmektedir. En küçük kareler kuramında sonraki gelişme Gauss'in 1821 yılında yayınladığı bir makalede ortaya çıkartılmış ve bu yayında Gauss sonradan kendi adı verilen Gauss-Markov teoreminin bir şeklini açıklamıştır. 
"Regresyon" terimi 19. yüzyılda İngiliz istatistikçisi Francis Galton tarafından bir biyolojik inceleme için ortaya atılmıştır. Bu incelemenin ana konusu kalıtım olup, aile içinde baba ve annenin boyu ile çocukların boyu arasındaki bağlantıyı araştırmakta ve çocukların boylarının bir nesil içinde eski ata nesillerinin ortalamasına geri döndüklerini yani bir nesil içinde "ortalamaya geri dönüş" olduğu inceleme konusudur. Galton "geri dönüş" terimi için ilk yazısında İngilizce olarak "reversion" terimi kullanmışsa da sonradan aynı anlamda olan "regression" sözcüğü kullanmıştır. Bu çalışmalarında Galton istatistiksel 'regresyon' kavramını ve yöntemini de geliştirmiştir. Udny Yüle ve Karl Pearson bu yöntemi daha geniş genel istatistiksel alanlara uygulayıp geliştirmişlerdir... Bu yazılarda bağımlı ve bağımsız değişkenlerin normal dağılım gösterdiği varsayılmaktadır. Bu kısıtlayıcı varsayım R.A. Fisher 1922 ve 1925 yıllarındaki yayınları ile sadece bağımlı değişkenin koşullu dağılımının normal olduğu hallere uygulanmak üzere daha genişletilmiştir.).
Bu kavramları ve yöntemleri genel olarak, kalıtım konusu dışında "ortalamaya geri dönüş" ile hiçbir ilgisi olmayan konularda, kullanan istatistikçiler "regresyon" terimini kullanmakta devam etmişlerdir. Zamanımızda, bu terim, kavram ve yöntemin Galton'un konusu ile bütün ilişkisi yok olmuştur ve artık "regresyon" terimi doğrusal bağlantı bulunması ve eğri uydurma ile eş anlamlar vermektedir.


## Temeldeki varsayımlar.
Doğrusal regresyon yöntemini kullanmak için temelde şu varsayımların bulunduğu kabul edilmektedir:
  - Çıkarımsal yöntem olduğu için kullanılan iki değişkenli örneklemin ya istatistiksel rastgele örneklem olduğu ya da anakütleyi çok iyi temsil ettiği bilinmektedir.
  - Bağımlı değişken içinde hata bulunmaktadır. Bu hatanın bir rassal değişken olduğu ve ortalama hatanın sıfır olduğudur. Sistematik hata da bulunması mümkündür ama bu hatanın incelemeye alınması regresyon analizi kapsamı dışındadır.
  - Bağımsız değişken hatasızdır. Eğer bağımsız değişkende hata bulunduğu varsayılırsa özel bir yöntem şekli olan değişkenler-içinde-hata modeli teknikler kullanılarak model kurulmalıdır.
  - Hatalar zaman içinde ve kendi aralarında birbirine bağımlı değildir. Buna otokorelasyon veya serisel korelasyon bulunmaması varsayımı adı verilir.
  - Hata varyansı sabittir ve veriler arasında hiç değişmediği varsayılır. Bu eşvaryanslılık veya "homoskedastisite" varsayımı adı ile anılır. Eger bu varsayim uygun degilse ağırlıklı en küçük kareler yöntemi uygulanabilir.
  - Hataların varyans-kovaryans matrisinin çapraz elamanları sabit hata varyansı olur ve matrisin diğer çapraz dışı elemanları 0 olur.
  - Eğer çoklu regresyon analizi yapılıyor ve üç veya daha çok parametre için kestirim isteniyorsa, bağımsız değişkenlerin birbirleri ile bağlantısının olmaması gereklidir. Buna çoklu doğrusallık (multicollinearity) olmaması varsayımı adı verilir. 
  - Hatalar bir normal dağılım gösterir. Eğer bu hataların normalliği varsayımı uygun değilse genelleştirilmiş doğrusal model uygulanabilir.


## Doğrusal regresyon.
### Anakütle doğrusal regresyon modeli.
Doğrusal regresyonda, anakütle model belirlenmesine göre bağımlı değişken y_i "parametreler"in bir doğrusal birleşiği olur. Dikkat edilirse "parametrelerden" bahis edilmektedir, çünkü bağımsız değişkenlerin bir doğrusal bileşiği olması gerekli değildir. Örneğin, tek bir bağımsız değişkenli ( x_i ) ve iki parametreli (β_0 ve β_1): 
doğru: y_i=β_0 +β_1 x_i +ϵ_i, i=1,n
Burada ϵ_i bir "hata terimi"dir ve i alt-indeksleri bir belirlenmiş mümkün gözlemi ifade eder. Ayrıca ϵ_i, i'nci gözlemin regresyon doğrusuna olan uzaklığını ifade etmekte olup ortalaması 0 ve varyansı σ^2 olan bir normal dağılış gösterir.
Çoklu doğrusal regresyonda iki veya daha çok sayıda bağımsız değişken veya bağımsız değişken fonksiyonu bulunur. Örneğin, önce verilmiş olan regresyon modeli yeni bir terim "xi"2 eklenerek değiştirilirse; şu anakütle çoklu doğrusal regresyon modeli ortaya çıkar:
parabol: y_i=β_0 +β_1 x_i +β_2 x_i^2+ϵ_i, i=1,m
Denklemin sağ tarafındaki bağımsız değişken için bir ikinci derece (kuadratik) ifade bulunmasına rağmen bu model hala doğrusal regresyon modelidir; çünkü üç tane parametre, yani β_0, β_1 ve β_2 ile çoklu değişkenli doğrusaldır.
Daha genel çoklu doğrusal regresyon modelinde "p" tane bağımsız değişken olduğu varsayılır ve anakütle modeli şöyle ifade edilir:
 y_i = β_0 + β_1 x_1i + ⋯+ β_p x_pi + ε_i ,
veya biraz daha kısa olarak
y_i = β_0 + ∑_j=1^p X_ij β_j + ε_i


### İki Değişkenli regresyon katsayı kestirimleri.
Anakütleden bir rastgele örneklem elde edilirse, bu β_i anakütle parametreleri için örneklem tahminleri β_i bulunur ve şu örneklem doğrusal regresyon denklemi elde edilir: 
 y_i = β_0 + β_1 X_i + e_i 
Burada e_i terimi örneklemden elde edilen "artık" olur; ve 
 e_i = y_i - y_i 
olur.
Bu örneklem regresyon denkleminin kestirimini elde etmenin yöntemi adi en küçük kareler (Ordinary least squares) olarak adlandırılır. Bu yönteme göre parametre kestirimleri artıkların kare toplamının minimum (en küçük) değerini bulmakla elde edilir. Artıklar minimum toplamı AKT şudur: 
AKT=∑_i=1^Ne_i^2
Bu fonksiyonun minimum değerini bulmak her parametre ile birinci kısımsal türevi alınarak sıfıra eşitlenir. Böylece her kısimsal türeve göre bir denklem elde edilir. Örneğin iki parametreli doğrusal regresyon için iki değişkenli iki denklem elde edilir. Bu simultane denklem sitemine normal denklemler adı verilir ve bu simultane denklemler sistemi birlikte çözülüp her anakütle parametresi için bir kestrim değeri bulunur.
İki parametreli basit doğrusal regresyon için iki anakütle parametre kestirimleri olan β_0, β_1 şu normal denklemlerin birlikte çözülmesi ile elde edilir:
[ m α+ ∑x_i β=∑y_i; ∑x_i α+ ∑x_i^2 β=∑x_iy_i ]
Burada bütün toplamlar "i"=1 den "i"="n" kadardır ve örneklem toplamları, toplam kareleri ve toplam karşı çarpımı ile elde edilirler.
Bu iki değiskenli iki doğrusal denklem için çeşitli çözüm yolları vardır. Burada Cramer'in kuralı kullanılıp şu ifade elde edilir:
β̂= m ∑x_iy_i - ∑x_i ∑y_i/Δ
=∑(x_i-x̅)(y_i-y̅)/∑(x_i-x̅)^2
α̂= ∑x_i^2 ∑y_i - ∑x_i ∑x_iy_i/Δ= y̅-x̅ β̂
burada
Δ= m∑x_i^2 - (∑x_i)^2
Böylece iki parametreli doğrusal regresyon icin en küçük kareler parametre tahminleri için formüller şöyledir:
β_1=∑(x_i-x̅)(y_i-y̅)/∑(x_i-x̅)^2 a
ve
β̂_̂0̂=y̅-β_1x̅
Burada x̅ x gözlem değerlerinin ortalaması ve y̅ ise y gözlem değerlerinin ortalamasıdır.


### Çok değişkenli regresyon katsayı kestirimleri.
Çok değişkenli regresyon için regresyon katsayısı kestrimi de iki değişkenli regresyon hâli gibi "en küçük kareler" yöntemi ile yapılır. Bu demektir ki kestirmde ortaya çıkacak artıkların karelerinin toplamı minimize edilecektir. Artıklar "ri" olur ve minimize edilecek objektif fonksiyon soyle ifade edilir: 
S=∑_i=1^mr_i^2,
Her bir artık, gözlemi yapılan bir değer ile kestirim modelini kullanarak elde edilen bir kestirim değeri arasındaki farktır:
r_i = y_i-∑_j=1^n X_ij β̂_j 
En küçük kareler yöntemine göre "S", toplam artıklar karesi, minimize edilecektir. Belirli koşullar geçerli ise, elde edilecek parametre kestrimleri (Gauss-Markov teoremine göre) en küçük varyans gösteren kestirim ve hatta maksimum değişirlilik özelliği gösteren dogrusal kestirim olabilirler.
Anakütle regresyon katsayılarının "kestirim regresyon katsayıları" için bir "^ (şapka)" notasyonu kullanılır: yani β̂_i. Genellikle çok değişkenli regresyon için normal denklemler şöyle ifade edilir
∑_i=1^N∑_p=1^p X_ijX_ikβ̂_k=∑_i=1^N X_ijy_i, j=1,p 
Matris notasyonu kullanılırsa, normal denklemler şöyle yazılabilir:
(𝐗^𝐓𝐗)β=𝐗^𝐓𝐲
Eğer X^TX matrisi singuler değilse 
 β̂ = (𝐗^𝐓𝐗)^-1 𝐗^𝐓𝐲,
Bu doğru uydurma tahminleridir.


### Hatalar varyansı ve toplam kareler.
Anakütle hatalarının normal dağılım gösterdiğine dair bir diğer varsayımı da kullanarak incelemeci önce "hatalar varyansı" ve "toplam kareler" değerlerini bulur ve bunları kullanarak tahmin edilen denklem ve parametreler üzerinde çıkarımsal istatistik sonuçlara varabilir.
Anakütle hata teriminin sabit bir varyansı bulunduğu varsayımına göre, "hatalar varyansı" kestirimi şöyle bulunur:
 σ̂_̂ε̂ = √(SSE/N-2) 
Bu ifadeye regresyon için kare kök ortalama hata karesi adı verilir. Parametre kestirimleri için standart hata şöyle bulunur:
σ̂_β_0=σ̂_ε √(1/N + x̅^2/∑(x_i-x̅)^2)
σ̂_β_1=σ̂_ε √(1/∑(x_i-x̅)^2)
Örneklem veri serisinin "değişebilirliği" değişik toplam kareler suretiyle ifade edilebilirler.
  - Tüm toplam kareler (örneklem varyansına orantılı olur):
SS_tot=∑_i (y_i-y̅)^2.
  - Regresyon toplam kareler: Bazen "açıklanan toplam kareler" diye anılır.
SS_reg=∑_i (f_i-f̅)^2
  - Toplam hatalar karesi: "Artıklar toplam karesi" olarak da isimlendirilir.
SS_err=∑_i (y_i - f_i)^2 
Bir regresyon modeli için parametreler, en küçük kareler yöntemi kullanılarak, tahmin edildikten sonra regresyon kontrol etme yöntemleri (yani çıkarımsal kontrol) uygulamak gerekir. Bunlar bulunan tahmin edilmiş genel doğrusal denklemin örnekleme uyum iyiliği sınaması ve kestirimci regresyon katsayılarının tek tek istatistiksel anlamlılığını sınamak veya güvenlik aralığı sağlamak olarak özetlenir.


### Kestirim denklemin genel uyum iyiliğinin çıkarımsal kontrolü.
Tahmin edilen doğrusal regresyon denkleminin genel olarak örnekleme uygun olup olmadığı için uygulanan çıkarımsal istatistik işlemleri iki türde olabilir:
  - belirleme katsayısı yani R-kare değerinin bulunması  
  - regresyon denklemine F-sınaması uygulaması.
Bu iki çıkarımsal kontrol aynı konunun çıkarımsal kontrolü için kullanılır: tahmin regresyon denkleminin tümünün uygun olup olmadığı. Bu nedenle iki yöntemden tek birini kullanmak yeterlidir. Buna rağmen birçok pratik analizde her iki çıkarımsal analiz de kullanılmaktadır; çünkü ikisinde hesaplanması için nerede ise aynı yöntem kullanılır ve eğer istatistik veya ekonometrik kompüter paketi kullanılırsa her iki yöntem için gerekli sonuçlar birlikte verilir.


#### Belirleme katsayısı yani R-kare (R2) değeri.
Belirleme katsayısı yani (R2) ) R-kare) için en genel tanımlama formülü sudur:R^2 ≡1-SS_errSS_tot.
yani (1-"toplam hatalar karesi") ile ("tüm toplam kareler") oranı;
veya (1-"Artıklar toplam karesi") ile ("tüm toplam kareler") oranı.
Çıkarımsal analizde R-kare değeri bulunur ve bulunan değer doğru hesaplanmışsa 0 ile 1 arasında olmalıdır. Yapılan bu analiz çeşitli hallerde açıklanabilir: 
  - Eğer sadece iki değişkenli (Y ve X0) regresyon analizi yapılıyorsa, geometrik olarak X-Y düzeyinde doğrusal regresyon tahmini bu düzeyde olan noktalara en-küçük-kareler prensibine göre en iyi uyan bir doğru uydurmaktır. Eğer bütün noktalar bu tahmin edilen doğru üzerinde ise tahmin edilen doğru uyumu bu veri noktaları için hiç hatasızdır ve veri noktaları doğrudan hiç ayrılık göstermez. Hat varyasyonu bu halde sıfır olur. Bu veri noktaları tahmin edilen doğrudan uzaklaştıkça hatalar varyasyonu çoğalmaktadır.
  - Üç değişkenli (Y ve X1 ve X2 değişkenli) çoklu regresyonda ise tahmin edilen bir üç boyutlu düzeydir. EĞer bu düzey veri noktalarına tam olarak uyarsa bütün veri noktaları tahmin edilen düzey üstüne düşerler. Veri noktaları tahmin edilen düzey etrafında dağılmaya başlarlarsa, hatalar varyasyonu artmaya başlar.
  - Üç değişkenden daha çoklu değişkenli regresyon analizi ise grafik olarak verilemez  çok boyutlu uzayı gösterir ve tahmin edilen regresyon katsayıları ise bu Çok boyutlu uzayda bir "hiper-düzey" ortaya çıkartır ve bu "hiper-düzey"den ayrılmalar hata varyasyonudur.
Genel olarak:
   - Eğer R2 değeri sıfıra "yakınsa", uyum iyiliği uygun olmadığı kabul edilir. Bu sonuç ortaya çıkarsa toplanan verilere kullanılan modelin uygun olmadığı sonucu çıkarılır ve bu uygunsuzluk modelinin değiştirilmesini gerektirir. Bu demektir ki model ile açıklanan varyasyon tüm varyansyonu %0ini açıklamakta ve geometrik olarak örneklem verileri regresyon ile elde edilen hiperdüzeyin etrafına çok dağılmış olarak bulunmaktadırlar. Bu çıkarıma varılırsa bu basamağa kadar yapılmış olan analizin bir kenara bırakılması ve diğer bir modelin bulunup kullanılması gerekir.
   - Eğer R2 değeri bire "yakınsa", uyum iyiliği uygun olduğu sonucu çıkarılır. Geometrik olarak çok değişken boyutlu uzayda, örneklem veri noktalarının hepsi uzayda dağılma göstermeyip hemen hepsi regresyon ile elde edilen hiperdüzey üzerinde bulunmaktadır. Bu halde çıkarımsal kontrol analizinin diğer basamağına geçilir.


#### F-testi.
Regresyon denklem uyum iyiliği için F-testi için sıfır hipotez (H0) anakütle model katsayılarının hepsinin değerinin 0 (sıfır)a eşit olduğudur. Yani sıfır hipotez gerçekse hesap ile elde edilen katsayı kestirimlerinin hepsi anakütle için 0dir yani hiçbir etki ve bağlantı anlamı vermez. Bu basamağa kadar yapılmış olan analizin bir kenara bırakılması ve diğer bir modelin bulunup kullanılması gerekir. Eğer sıfır hipotez reddedilirse regresyon kestirimleri genellikle uygundur ve cikarimsal kontrol analizine devam edilir.


#### Kestirimi yapılan her tek regresyon parametresinin değerinin çıkarımsal kontrolü.
Yani bu tahmin edilmiş standart hataları kullanarak her tek anakütle regresyon parametresi hakkında "kestirim güvenlik aralıkları" kurabilir ve "hipotez sınamaları" yapılabilir.


### İnterpolasyon ve ekstrapolasyon.
Regresyon modelleri kullanarak kestirim yapılmak istenirse, bağımsız değişken olan x değişken veri değerleri verilirse bağımlı değişken y için kestrim değerleri (f_i) tahmin etmek için kullanabilirler.
Eğer bu kestirim, modeli kurmak için kullanılan x değişken değerleri için yapılıyorsa, bu işlem interpolasyon olarak adlandırılır. Eğer kestirim modeli kurmak için kullanılan bağımsız değişken değerlerinin dışındaki değerler ile yapılırsa, ekstrapolasyon olarak adlandırılır ve ekstrapolasyon çok daha yanlı olabilir.


## Diğer yaklaşımlar.
### Ağırlıklı en küçük kareler yöntemi.
En küçük kareler yöntemi kullanılırken yapılan temel varsayımlarından biri hata terimi varyanslarının birbirine eşit olduğudur. Eğer gözlem hataları birbirine eşit olmayan hata terimi varyansları gösteriyorlarsa, en küçük kareler yönteminin bu önemli varsayımı ihlal edilmiş olur ve "en küçük kareler" yöntemi ile elde edilen regresyon kestirimleri anlamlı olmayabilir. Bu sorunu çözümlemek için her gözleme ayrı ağırlık vererek en küçük kareler yöntemi uygulamak imkânı vardır ve bu genelleştirmeye ağırlıklı en küçük kareler adı verilir. Ağırlıklı En Küçük Kareler Yöntemi,değişkenlere ağırlık vererek veya değişkenlerin önem derecesini değiştirerek uygulanan en küçük kareler yöntemidir.


### Değişkenlerde-hatalar modeli.
En küçük kareler yöntemi kullanılırken yapılan temel varsayımlarından biri de gözlem hatalarının yalnızca bağımlı değişkenlerde yapıldığı ve bağımsız değişkende gözlem hatası bulunmadığıdır. Eğer bağımlı değişkende hata bulunduğu ve bunun elde edilen regresyon kestirim değerlerine çok etki yaptığı bilinirse, değişkenlerde-hatalar-modeli veya diğer ismi ile total en küçük kareler modeli ve bu modelle ilişkili kestirim yöntemleri kullanılabilir.


### Genelleştirilmiş doğrusal model.
Eğer anakütle regresyon modeli içindeki hatalar için olasılık dağılım fonksiyonu bir normal dağılım göstermiyorsa genelleştirilmiş doğrusal model kullanılabilir. Bu halde hataların olasılık dağılım fonksiyonu ussel dağılım, gamma dağılımı, ters Gauss tipi dağılım, Poisson dağılımı, binom dağılım, multinom dağılım vb. olabilir.


### Doğrusal olmayan regresyon.
Eğer değişkenler hakkında kurulan teori dolayasıyla anakütle modeli parametreleri ile doğrusal değilse, kestirim yapılırken toplam kareleri tekrarlama usulü kullanarak minimize edilmesi gerekir. Bu kullanılan tekrarlama yöntemi birçok karışık sorunlar ortaya çıkarır. Bunların özet olarak incelenmesi için doğrusal olmayan regresyon#Dogrusal olmayan ve dogrusal en kucuk kareler arasindaki farklar maddesine bakiniz.


## İçsel kaynaklar.
  - Parçalar için regresyon
  - Güvenlik aralığı
  - Güvenlik yöresi
  - Ekstrapolasyon
  - Krigleme
  - Tahmin etme
  - Tahmin aralığı
  - İstatistik
  - Trend kestirimi
  - Güçlü regresyon
  - Çokdeğişirli normal dağılım


## Bibliyografya.
  - Gujarati, Damodar (çev. Ümit Şenesen, Gülay Günlük Şenesen) (2008) "Temel Ekonometri", Literatür Yayınları ISBN 975-7860-99-9.
  - Audi, R., Ed. (1996). "Curve fitting problem," "The Cambridge Dictionary of Philosophy". Cambridge, Cambridge University Press. pp. 172–173.  YÖNLENDİRME 
  - William H. Kruskal and Judith M. Tanur, ed. (1978), "Linear Hypotheses," "International Encyclopedia of Statistics". Free Press, cilt 1  YÖNLENDİRME 
Evan J. Williams, "I. Regression," say. 523-41. 
Julian C. Stanley, "II. Analysis of Variance," pp. 541-554.
  - Lindley, D.V. (1987). "Regression and correlation analysis," , Cilt. 4, say. 120-23.  YÖNLENDİRME 
  - Birkes, David and Yadolah Dodge, "Alternative Methods of Regression". ISBN 0-471-56881-3  YÖNLENDİRME 
  - Chatfield, C. (1993) "Calculating Interval Forecasts," "Journal of Business and Economic Statistics," 11. pp. 121–135.  YÖNLENDİRME 
  - Dinçer, F. "Regresyon Kelimesinin Tarihi", 
  - Draper, N.R. ve Smith, H. (1998)."Applied Regression Analysis" Wiley Series in Probability and Statistics  YÖNLENDİRME 
  - Fox, J. (1997). "Applied Regression Analysis, Linear Models and Related Methods." Sage  YÖNLENDİRME 
  - Hardle, W., "Applied Nonparametric Regression" (1990), ISBN 0-521-42950-1  YÖNLENDİRME 
  - Meade, N. ve T. Islam (1995) "Prediction Intervals for Growth Curve Forecasts," "Journal of Forecasting," 14, say. 413-430.  YÖNLENDİRME 
  - Munro, Barbara Hazard (2005) "Statistical Methods for Health Care Research 5th ed." Lippincott Williams & Wilkins,  YÖNLENDİRME 
  - Sykes, A.O. "An Introduction to Regression Analysis" (Innaugural Coase Lecture)
  - Kotsiantis, S., Kanellopoulos, D. ve Pintelas, P. (2006) "Local Additive Regression of Decision Stumps", "Lecture Notes in Artificial Intelligence", Springer-Verlag, Vol. 3955, SETN 2006, say.. 148 – 157  YÖNLENDİRME 
  - Kotsiantis, S. ve Pintelas, P. (2005) "Selective Averaging of Regression Models", "Annals of Mathematics, Computing & TeleInformatics," Cilt 1, No 3, say. 66-75  YÖNLENDİRME