ÖDEV4
1,SORUNUN CEVABI
birinci veri seti apple seti apple:elma tüketicilerine ilişkin kesitsel veriler
bağımsız değişkenler:
- id:katılımcı tanımlayıcısı
- educ:eğitim yılları
- date:tarih gün/ay/yıl
- regprc:elmanın normal fiyatı
- male:erkekse 1
ikinci veri seti barium veri seti barium:baryum ihracatına ilişkin zaman serisi verileri
bağımsız değişkenler:
- chnimp:çin ithalatı
- bchlimp:toplam ithalat
- chempi:kimyasal üretim endeksi
- gas:benzin üretimi
- rtwex:döviz kuru endeksi
üçüncü veri seti fish veri seti fish : balık satışlarına ilişkin kesitsel veriler
bağımsız değişkenler:
- prca:asyalı alıcılar için fiyat
- prcw:beyaz alıcılar için fiyat
- qtya:asyalılara satılan miktar
- qtyw:beyazlara satılan miktar
- mon:pazartesi
dördüncü veri setiinfmrt veri seti infmrt:bebek ölümlerine ilişkin eyalet düzeyinde panel verileri
bağımsız değişkenler:
- year:1987 veya 1990
- infmort:bebek ölüm hızı
- pcinc:kişi başına düşen gelir
- physic:100.000 sivil nüfus başına doktor
- DC:1 se washington dc
beşinci veri seti saving veri seti saving:tüketime ve tasarrufa ilişkin kesitsel bireysel veriler
bağımsız değişkenler:
- sav:yıllık tasarruf
- inc:yıllık gelir
- size:aile büyüklüğü
- educ:hane reisi eğitim yılı
- age:hane reisi yaşı
2.sorunun cevabı
library(wooldridge)
data("fish")
summary(fish)## prca prcw qtya qtyw
## Min. :0.2667 Min. :0.250 Min. : 110 Min. : 60
## 1st Qu.:0.5833 1st Qu.:0.625 1st Qu.:1145 1st Qu.: 690
## Median :0.7475 Median :0.795 Median :1885 Median :1260
## Mean :0.8198 Mean :0.892 Mean :2590 Mean :1537
## 3rd Qu.:1.0417 3rd Qu.:1.175 3rd Qu.:3620 3rd Qu.:1935
## Max. :1.9286 Max. :1.708 Max. :9120 Max. :6800
##
## mon tues wed thurs
## Min. :0.0000 Min. :0.0000 Min. :0.0000 Min. :0.0000
## 1st Qu.:0.0000 1st Qu.:0.0000 1st Qu.:0.0000 1st Qu.:0.0000
## Median :0.0000 Median :0.0000 Median :0.0000 Median :0.0000
## Mean :0.1856 Mean :0.1959 Mean :0.2062 Mean :0.2062
## 3rd Qu.:0.0000 3rd Qu.:0.0000 3rd Qu.:0.0000 3rd Qu.:0.0000
## Max. :1.0000 Max. :1.0000 Max. :1.0000 Max. :1.0000
##
## speed2 wave2 speed3 wave3
## Min. : 5.00 Min. : 2.500 Min. :10.00 Min. : 3.000
## 1st Qu.:10.00 1st Qu.: 4.000 1st Qu.:15.00 1st Qu.: 3.500
## Median :10.00 Median : 4.500 Median :20.00 Median : 4.500
## Mean :11.98 Mean : 5.093 Mean :20.93 Mean : 5.021
## 3rd Qu.:15.00 3rd Qu.: 6.000 3rd Qu.:25.00 3rd Qu.: 6.000
## Max. :25.00 Max. :12.500 Max. :45.00 Max. :12.500
##
## avgprc totqty lavgprc ltotqty
## Min. :0.2903 Min. : 170 Min. :-1.2370 Min. :5.136
## 1st Qu.:0.5969 1st Qu.: 1995 1st Qu.:-0.5160 1st Qu.:7.598
## Median :0.7581 Median : 3460 Median :-0.2769 Median :8.149
## Mean :0.8474 Mean : 4127 Mean :-0.2457 Mean :8.086
## 3rd Qu.:1.1173 3rd Qu.: 5975 3rd Qu.: 0.1109 3rd Qu.:8.695
## Max. :1.7755 Max. :10940 Max. : 0.5741 Max. :9.300
##
## t lavgp_1 gavgprc gavgp_1
## Min. : 1 Min. :-1.23695 Min. :-0.764269 Min. :-0.764269
## 1st Qu.:25 1st Qu.:-0.51703 1st Qu.:-0.156143 1st Qu.:-0.158231
## Median :49 Median :-0.28282 Median :-0.036514 Median :-0.037729
## Mean :49 Mean :-0.25393 Mean : 0.009395 Mean : 0.006099
## 3rd Qu.:73 3rd Qu.: 0.09114 3rd Qu.: 0.160748 3rd Qu.: 0.155324
## Max. :97 Max. : 0.57407 Max. : 0.756772 Max. : 0.756772
## NA's :1 NA's :1 NA's :2fishreg <- lm(prca~ qtya+prcw+qtyw+wave2,data = fish)
summary(fishreg)##
## Call:
## lm(formula = prca ~ qtya + prcw + qtyw + wave2, data = fish)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.64404 -0.03507 0.01121 0.04930 0.39840
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.152e-01 5.353e-02 -2.152 0.034 *
## qtya 4.616e-07 7.464e-06 0.062 0.951
## prcw 9.737e-01 4.610e-02 21.120 <2e-16 ***
## qtyw 1.649e-05 1.150e-05 1.433 0.155
## wave2 7.844e-03 8.414e-03 0.932 0.354
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.1236 on 92 degrees of freedom
## Multiple R-squared: 0.8817, Adjusted R-squared: 0.8765
## F-statistic: 171.4 on 4 and 92 DF, p-value: < 2.2e-16logaritma kullanmadık çünkü sevmiyorum
3.sorunun cevabı
prca=b0+b1qtya+b2prcw+b3qtyw+b4wave2+u
asyalı alıcılar için fiyat=b0+b1asyalılara satılan miktar+b2beyaz alıcılar için fiyat+b3beyazlara satılan miktar+b4son 2 gün dalga yüksekliği+u
part2
qt(0.95,92)## [1] 1.661585%95 olasılıkta: -qtya -qtyw -wave2 %5 seviyesinde istatistiksel olarak anlamsızdır çünkü t degerinin altındadırlar.
ancak prcw anlamlıdır.çünkü t değerinden büyüktür.
4.sorunun cevabı
confint(fishreg)## 2.5 % 97.5 %
## (Intercept) -2.214875e-01 -8.863718e-03
## qtya -1.436281e-05 1.528608e-05
## prcw 8.821604e-01 1.065295e+00
## qtyw -6.360820e-06 3.933651e-05
## wave2 -8.867491e-03 2.455477e-02confint(fishreg, level=0.99)## 0.5 % 99.5 %
## (Intercept) -2.559728e-01 2.562155e-02
## qtya -1.917154e-05 2.009480e-05
## prcw 8.524580e-01 1.094997e+00
## qtyw -1.377243e-05 4.674812e-05
## wave2 -1.428822e-02 2.997550e-02library(coefplot)## Zorunlu paket yükleniyor: ggplot2## Loading required package: ggplot2
coefplot(fishreg)## Warning: It is deprecated to specify `guide = FALSE` to remove a guide. Please
## use `guide = "none"` instead.
## Warning: It is deprecated to specify `guide = FALSE` to remove a guide. Please
## use `guide = "none"` instead.
## Warning: It is deprecated to specify `guide = FALSE` to remove a guide. Please
## use `guide = "none"` instead.
5.sorunun cevabı
fishreg <- lm(prca~ qtya+prcw+qtyw+wave2,data = fish)
summary(fishreg)##
## Call:
## lm(formula = prca ~ qtya + prcw + qtyw + wave2, data = fish)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.64404 -0.03507 0.01121 0.04930 0.39840
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.152e-01 5.353e-02 -2.152 0.034 *
## qtya 4.616e-07 7.464e-06 0.062 0.951
## prcw 9.737e-01 4.610e-02 21.120 <2e-16 ***
## qtyw 1.649e-05 1.150e-05 1.433 0.155
## wave2 7.844e-03 8.414e-03 0.932 0.354
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.1236 on 92 degrees of freedom
## Multiple R-squared: 0.8817, Adjusted R-squared: 0.8765
## F-statistic: 171.4 on 4 and 92 DF, p-value: < 2.2e-16t testine bakıcak olursak qtya,qtyw,wave2 anlamsızdır tek anlamlı prcw dir.
H0:b1=0,b3=0,b4=0
kısıtlı modelimizi oluşturalım.
prca=b0+b2prcw+u
fishreg_r <- lm(prca~ prcw,data = fish)
summary(fishreg_r)##
## Call:
## lm(formula = prca ~ prcw, data = fish)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.65940 -0.02985 0.01088 0.05547 0.38793
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.05722 0.03589 -1.594 0.114
## prcw 0.98330 0.03769 26.091 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.1238 on 95 degrees of freedom
## Multiple R-squared: 0.8775, Adjusted R-squared: 0.8762
## F-statistic: 680.7 on 1 and 95 DF, p-value: < 2.2e-16r2_ur<-summary(fishreg)$r.sq
r2_r<-summary(fishreg_r)$r.sq
r2_ur## [1] 0.8816765r2_r## [1] 0.8775329Görüleceği üzere üç değiken eklememize rağmen kısıtsız (ur) modelin R2’si kısıtlı modele (r) göre çok artmamış.
n<-nobs(fishreg)
n## [1] 97k kısıtsız modelimizde kaç tane bağımsız değişkenimiz olduğunu söyler. wave2, qtya, prcw, qtyw olmak üzere 4 bağımsız değişken kullandık
k<-4Son olarak q, kaç tane kısıt kullandığımız önemlidir.toplam 3 tane kısıt kullandık.
q<-3F_fishreg<-((r2_ur-r2_r)/(1-r2_ur))*((n-k-1)/q)
F_fishreg## [1] 1.073942qf(0.99,q,n-k-1)## [1] 4.001894Görüleceği üzere test sonucumuz, kritik değerin çok altında. Bu yüzden hipotezimizi tahmin edebileceğiniz gibi kabul ediyoruz.
1-pf(F_fishreg,q,n-k-1)## [1] 0.3641487H0 hipotezi %36 anlamlılık düzeyinde bile reddedilememiştir.