Etkileşimli tek bir model

Python ile statsmodels'da Orta Düzey Regresyon

Maarten Van den Broeck

Content Developer at DataCamp

Etkileşim nedir?

Balık veri setinde

Farklı türlerin kütle/uzunluk oranları farklıdır.
Uzunluğun beklenen kütle üzerindeki etkisi türe göre değişir.

Daha genel olarak

Bir açıklayıcı değişkenin beklenen yanıta etkisi, başka bir açıklayıcı değişkenin değerine bağlı olarak değişir.

Etkileşimleri belirtme

Etkileşim yok

response ~ explntry1 + explntry2

Etkileşimli (örtük)

response_var ~ explntry1 * explntry2

Etkileşimli (açık)

response ~ explntry1 + explntry2 + explntry1:explntry2

Etkileşim yok

mass_g ~ length_cm + species

Etkileşimli (örtük)

mass_g ~ length_cm * species

Etkileşimli (açık)

mass_g ~ length_cm + species + length_cm:species

Modeli çalıştırma

mdl_mass_vs_both = ols("mass_g ~ length_cm * species", data=fish).fit()

print(mdl_mass_vs_both.params)

Intercept                    -1035.3476
species[T.Perch]               416.1725
species[T.Pike]               -505.4767
species[T.Roach]               705.9714
length_cm                       54.5500
length_cm:species[T.Perch]     -15.6385
length_cm:species[T.Pike]       -1.3551
length_cm:species[T.Roach]     -31.2307

Katsayıları okumak daha kolay

mdl_mass_vs_both_inter = ols("mass_g ~ species + species:length_cm + 0", data=fish).fit()

print(mdl_mass_vs_both_inter.params)

species[Bream]             -1035.3476
species[Perch]              -619.1751
species[Pike]              -1540.8243
species[Roach]              -329.3762
species[Bream]:length_cm      54.5500
species[Perch]:length_cm      38.9115
species[Pike]:length_cm       53.1949
species[Roach]:length_cm      23.3193

Tanıdık sayılar

print(mdl_mass_vs_both_inter.params)

species[Bream]             -1035.3476
species[Perch]              -619.1751
species[Pike]              -1540.8243
species[Roach]              -329.3762
species[Bream]:length_cm      54.5500
species[Perch]:length_cm      38.9115
species[Pike]:length_cm       53.1949
species[Roach]:length_cm      23.3193

print(mdl_bream.params)

Intercept   -1035.3476
length_cm      54.5500

Haydi pratik yapalım!

Python ile statsmodels'da Orta Düzey Regresyon