Kuis 2 1. Panggil Data Wages Penyelesaian: Script
Kuis 2 1. Panggil Data Wages Penyelesaian: Script
Kuis 2 1. Panggil Data Wages Penyelesaian: Script
Npm : 1817031056
Ekonomatrika B
KUIS 2
Script :
> library(plm)
> data("Wages",package="plm")
> data1<-Wages
> data1
Output :
exp wks bluecol ind south smsa married sex union ed black lwage
1 3 32 no 0 yes no yes male no 9 no 5.56068
2 4 43 no 0 yes no yes male no 9 no 5.72031
3 5 40 no 0 yes no yes male no 9 no 5.99645
4 6 39 no 0 yes no yes male no 9 no 5.99645
5 7 42 no 1 yes no yes male no 9 no 6.06146
6 8 35 no 1 yes no yes male no 9 no 6.17379
7 9 32 no 1 yes no yes male no 9 no 6.24417
8 30 34 yes 0 no no yes male no 11 no 6.16331
9 31 27 yes 0 no no yes male no 11 no 6.21461
10 32 33 yes 1 no no yes male yes 11 no 6.26340
11 33 30 yes 1 no no yes male no 11 no 6.54391
12 34 30 yes 1 no no yes male no 11 no 6.69703
13 35 37 yes 1 no no yes male no 11 no 6.79122
14 36 30 yes 1 no no yes male no 11 no 6.81564
15 6 50 yes 1 no no yes male yes 12 no 5.65249
16 7 51 yes 1 no no yes male yes 12 no 6.43615
17 8 50 yes 1 no no yes male yes 12 no 6.54822
18 9 52 yes 1 no no yes male yes 12 no 6.60259
19 10 52 yes 1 no no yes male yes 12 no 6.69580
20 11 52 yes 1 no no no male yes 12 no 6.77878
21 12 46 yes 1 no no no male yes 12 no 6.86066
22 31 52 yes 0 no yes no female no 10 yes 6.15698
23 32 46 yes 0 no yes no female no 10 yes 6.23832
24 33 46 yes 0 no yes no female no 10 yes 6.30079
25 34 49 yes 0 no yes no female no 10 yes 6.35957
26 35 44 yes 0 no yes no female no 10 yes 6.46925
27 36 52 yes 0 no yes no female no 10 yes 6.56244
28 37 46 yes 0 no yes no female no 10 yes 6.62141
29 10 50 yes 0 no no yes male yes 16 no 6.43775
30 11 46 yes 0 no no yes male yes 16 no 6.62007
31 12 40 yes 0 no no yes male yes 16 no 6.63332
32 13 50 no 0 no no yes male no 16 no 6.98286
33 14 47 yes 0 no yes yes male no 16 no 7.04752
34 15 47 no 0 no no yes male no 16 no 7.31322
35 16 49 no 0 no no yes male no 16 no 7.29574
36 26 44 yes 1 no yes yes male no 12 no 6.90575
37 27 47 yes 1 no yes yes male no 12 no 6.90575
38 28 47 yes 1 no yes yes male no 12 no 6.90776
39 29 47 yes 1 no yes yes male no 12 no 7.00307
40 30 44 yes 1 no yes yes male no 12 no 7.06902
41 31 45 yes 1 no yes yes male no 12 no 7.52023
42 32 47 yes 1 no yes yes male no 12 no 7.33889
43 15 46 yes 0 no no yes male yes 12 no 6.13340
44 16 48 yes 0 no no yes male yes 12 no 6.17379
45 17 49 yes 0 no no yes male yes 12 no 6.21261
46 18 46 yes 0 no no yes male yes 12 no 6.31355
47 19 47 yes 0 no no yes male yes 12 no 6.37502
48 20 47 yes 0 no no yes male yes 12 no 6.44572
49 21 48 yes 0 no no yes male yes 12 no 6.52209
50 23 51 yes 1 yes no yes male yes 10 no 6.33150
51 24 50 yes 1 yes no yes male yes 10 no 6.40357
52 25 50 yes 1 yes no yes male yes 10 no 6.54391
53 26 50 yes 1 yes no yes male yes 10 no 6.56244
54 27 44 yes 1 no yes yes male yes 10 no 6.59167
55 28 51 yes 1 no yes yes male yes 10 no 6.81783
56 29 51 yes 1 no yes yes male yes 10 no 6.89163
57 3 50 no 0 yes yes yes male no 16 no 6.55108
58 4 48 no 0 yes yes yes male no 16 no 6.55108
59 5 50 no 0 yes yes yes male no 16 no 6.80239
60 6 48 no 0 yes yes yes male no 16 no 6.90776
61 7 48 no 0 yes yes yes male no 16 no 7.09008
62 8 44 no 0 yes yes yes male no 16 no 7.17012
63 9 48 no 0 yes yes yes male no 16 no 7.20786
64 3 49 no 0 yes yes yes male no 16 no 6.39693
65 4 47 no 0 yes yes yes male no 16 no 6.43775
66 5 46 no 0 yes yes yes male no 16 no 6.43775
67 6 44 no 0 yes yes yes male no 16 no 6.43775
68 7 43 no 0 yes yes yes male no 16 no 6.52209
69 8 34 no 0 yes yes yes male no 16 no 6.61338
70 9 40 no 0 yes yes yes male no 16 no 6.73934
71 24 47 yes 0 no yes no male yes 12 no 6.65801
72 25 48 yes 0 no yes no male yes 12 no 6.72623
73 26 45 yes 0 no yes no male yes 12 no 6.80239
74 27 45 yes 0 no yes yes male yes 12 no 6.90776
75 28 47 yes 0 no yes yes male yes 12 no 7.03966
76 29 17 yes 0 no yes yes male yes 12 no 7.12769
77 30 47 yes 0 no yes yes male yes 12 no 7.24423
78 21 47 yes 0 yes yes yes male yes 12 no 6.55108
79 22 46 yes 0 yes yes yes male yes 12 no 6.62936
80 23 47 yes 0 yes yes yes male yes 12 no 6.72263
81 24 47 yes 0 yes yes yes male yes 12 no 6.73340
82 25 47 yes 0 yes yes yes male yes 12 no 6.72263
83 26 47 yes 0 yes yes yes male yes 12 no 6.95177
[ reached 'max' / getOption("max.print") -- omitted 4082 rows ]
4. Struktur data
Penyelesaian :
> str(data1)
$ ed : int 9 9 9 9 9 9 9 11 11 11 ...
Interpretasi :
Data tersebut merupakan jenis data frame yang terdisi dari 4165 data dan memiliki 14
variable yaitu i, time, exp, wks, bluecol,ind,south,smsa,married,sex,union,ed,black dan
iwage. Dengan lwage adalah merupakan variabel y atau yang merupakan variabel numerik
dan yang lainnya yaitu exp, wks, buecol, ind, south, smsa, married, sex, union, ed, black
merupakan variabel kategori
> summary(data1)
(Other):4123 7:595
Median :12.00
Mean :12.85
3rd Qu.:16.00
Max. :17.00
black lwage
Median :6.685
Mean :6.676
3rd Qu.:6.953
Max. :8.537
> names(data1)
Penyelesaian :
Script :
plm(lwage~exp+wks+bluecol+ind+south+smsa+married+sex+union+ed+black, data
summary(data1)
Output :
> data1<-
+ plm(lwage~exp+wks+bluecol+ind+south+smsa+married+sex+union+ed+black, data
> summary(data1)
Oneway (individual) effect Within Model
Call:
Residuals:
Coefficients:
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Penyelesaian :
Script :
data1<-
plm(lwage~exp+wks+bluecol+ind+south+smsa+married+sex+union+ed+black, data
summary(data1)
Output :
> data1<-
+ plm(lwage~exp+wks+bluecol+ind+south+smsa+married+sex+union+ed+black, data
> summary(data1)
Call:
Residuals:
Min. 1st Qu. Median 3rd Qu. Max.
Coefficients:
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
R-Squared: 0.65251
Penyelesaian :
> #Nomor 5c
> summary(data1)
Call:
Residuals:
Coefficients:
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
R-Squared: 0.0049904
Penyelesaiaan :
> #Nomor 5d
> summary(Data1.bt)
Call:
Residuals:
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
R-Squared: 0.5215
Penyelesaian :
> #Nomor 5e
> summary(DataW.ra)
(Swamy-Arora's transformation)
Call:
Effects:
theta: 0.7896
Residuals:
Coefficients:
Estimate Std. Error z-value Pr(>|z|)
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
R-Squared: 0.36909
> #Nomor 6a
> #Uji chow
> pFtest(1ata1.fe , data1.ce)
data: lwage ~ exp + wks + bluecol + ind + south + smsa + married + ...
F = 32.269, df1 = 591, df2 = 3562, p-value < 2.2e-16
alternative hypothesis: significant effects
b. Uji Hausman
Penyelesaian :
> #Nomor 6b
> #uji hausman
> phtest(data1.ra,data1.fe)
Hausman Test
data: lwage ~ exp + wks + bluecol + ind + south + smsa + married + ...
chisq = 6529.5, df = 8, p-value < 2.2e-16
alternative hypothesis: one model is inconsistent