Dane pełne 16 zmiennych
## Warning: pakiet 'readxl' został zbudowany w wersji R 4.3.3
dane2011<- read_excel("C:/Users/HP/Downloads/WAD projekt dane.xlsx",
sheet = "2011")
head(dane2011)
## # A tibble: 6 × 17
## `#` Eksport towarów i usług (…¹ import towarów i usł…² `Bezrobocie gółem`
## <chr> <dbl> <dbl> <dbl>
## 1 Austria 53.9 51.2 4.6
## 2 Belgium 80.7 80.8 7.1
## 3 Bulgaria 58.8 58.6 11.3
## 4 Croatia 38.1 39.8 13.7
## 5 Cyprus 53.5 57.2 7.9
## 6 Czechia 70.8 67 6.7
## # ℹ abbreviated names: ¹`Eksport towarów i usług (% PKB)`,
## # ²`import towarów i usług (% PKB)`
## # ℹ 13 more variables: `Bezrobocie, ogółem wśród młodzieży` <dbl>,
## # `Bezpośrednie inwestycje zagraniczne` <dbl>,
## # `Dług rządu centralnego` <dbl>, `rządowe wydatki konsumpcyjne` <dbl>,
## # `wydatki konsumpcyjne ogółem (%PKB)` <dbl>, `PKB per capita` <dbl>,
## # Inflacja <dbl>, `Oczekiwana długość życia` <dbl>, …
rownames(dane2011) <- dane2011[[1]]
## Warning: Setting row names on a tibble is deprecated.
# Remove the first column (now it's in row names)
dane2011 <- dane2011[, -1]
analiza składowych głównych
wynik = prcomp(dane2011, scale=T)
wynik
## Standard deviations (1, .., p=16):
## [1] 2.20619678 1.88443361 1.56627696 1.20611141 0.98455059 0.86970814
## [7] 0.76117664 0.65140445 0.57702785 0.49144475 0.40872761 0.28602840
## [13] 0.25770415 0.18755162 0.13132518 0.04521953
##
## Rotation (n x k) = (16 x 16):
## PC1 PC2
## Eksport towarów i usług (% PKB) 0.24901799 0.330226024
## import towarów i usług (% PKB) 0.20341504 0.360343085
## Bezrobocie gółem -0.26154124 0.007548019
## Bezrobocie, ogółem wśród młodzieży -0.27612385 -0.005380294
## Bezpośrednie inwestycje zagraniczne 0.02590095 -0.124860117
## Dług rządu centralnego -0.06266567 -0.373699780
## rządowe wydatki konsumpcyjne 0.12466293 -0.378969212
## wydatki konsumpcyjne ogółem (%PKB) -0.34183644 -0.279171625
## PKB per capita 0.36600407 -0.059140262
## Inflacja -0.26139967 0.259335067
## Oczekiwana długość życia 0.25063831 -0.369794413
## Eksport żywności -0.14833096 -0.255079290
## Eksport wyrobów przemysłowych 0.19617492 0.206234406
## Eksport usług 0.24191670 -0.217244010
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana -0.35974037 0.071554613
## Zatrudnienie w rolnictwie -0.31849882 0.129647320
## PC3 PC4
## Eksport towarów i usług (% PKB) 0.273747455 -0.246677762
## import towarów i usług (% PKB) 0.278769491 -0.252875331
## Bezrobocie gółem -0.155590658 -0.617751988
## Bezrobocie, ogółem wśród młodzieży -0.179166190 -0.573216328
## Bezpośrednie inwestycje zagraniczne 0.383282009 -0.047559906
## Dług rządu centralnego 0.155827274 -0.127716371
## rządowe wydatki konsumpcyjne 0.005427735 -0.004955756
## wydatki konsumpcyjne ogółem (%PKB) -0.105930127 0.020891432
## PKB per capita 0.177254310 -0.052383921
## Inflacja 0.202940611 0.098449182
## Oczekiwana długość życia -0.019143178 -0.060625314
## Eksport żywności 0.489587657 -0.056314233
## Eksport wyrobów przemysłowych -0.442682390 0.052290369
## Eksport usług -0.239700038 0.012566443
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana 0.188930079 0.229935843
## Zatrudnienie w rolnictwie -0.082572009 0.266325958
## PC5 PC6
## Eksport towarów i usług (% PKB) -0.06128878 0.02207741
## import towarów i usług (% PKB) 0.01899040 -0.04339067
## Bezrobocie gółem -0.08872403 -0.01064980
## Bezrobocie, ogółem wśród młodzieży -0.02337106 0.13661373
## Bezpośrednie inwestycje zagraniczne 0.73888207 0.00715146
## Dług rządu centralnego -0.02832343 0.37888044
## rządowe wydatki konsumpcyjne -0.32372396 -0.49877530
## wydatki konsumpcyjne ogółem (%PKB) 0.22791455 -0.10960611
## PKB per capita -0.29421767 0.32368879
## Inflacja -0.14993962 0.24511471
## Oczekiwana długość życia -0.12254799 0.32304791
## Eksport żywności 0.00312608 0.06366561
## Eksport wyrobów przemysłowych 0.24604914 0.15673129
## Eksport usług 0.06091649 0.30393060
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana -0.30893088 0.01970972
## Zatrudnienie w rolnictwie -0.01451599 0.42879526
## PC7 PC8
## Eksport towarów i usług (% PKB) -0.136122610 0.16425280
## import towarów i usług (% PKB) -0.163099807 0.17192316
## Bezrobocie gółem 0.144971255 0.06574988
## Bezrobocie, ogółem wśród młodzieży 0.006257848 -0.22307968
## Bezpośrednie inwestycje zagraniczne 0.057319782 -0.19139196
## Dług rządu centralnego -0.555188878 0.42431767
## rządowe wydatki konsumpcyjne -0.177520789 -0.09819757
## wydatki konsumpcyjne ogółem (%PKB) -0.061601745 -0.06355955
## PKB per capita 0.076635349 -0.36312228
## Inflacja 0.243404100 -0.23493008
## Oczekiwana długość życia -0.056906865 -0.33800731
## Eksport żywności 0.179764851 0.12565494
## Eksport wyrobów przemysłowych -0.291710965 0.06345987
## Eksport usług 0.588160620 0.48319033
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana 0.031303179 0.28362485
## Zatrudnienie w rolnictwie -0.233032600 -0.15487338
## PC9 PC10
## Eksport towarów i usług (% PKB) -0.11279341 -0.11601521
## import towarów i usług (% PKB) -0.04390247 -0.18317643
## Bezrobocie gółem -0.05619780 -0.01851305
## Bezrobocie, ogółem wśród młodzieży -0.07862362 -0.06176949
## Bezpośrednie inwestycje zagraniczne -0.03657059 -0.19799968
## Dług rządu centralnego 0.31528753 -0.03549657
## rządowe wydatki konsumpcyjne 0.01213635 -0.61758122
## wydatki konsumpcyjne ogółem (%PKB) 0.08003094 -0.06370672
## PKB per capita -0.10721121 -0.04879448
## Inflacja 0.69771833 -0.31896644
## Oczekiwana długość życia 0.01816851 0.26915687
## Eksport żywności -0.18203604 0.10546843
## Eksport wyrobów przemysłowych 0.16388473 -0.05599493
## Eksport usług -0.05325955 -0.35332729
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana -0.12850315 0.15304596
## Zatrudnienie w rolnictwie -0.53919783 -0.42873245
## PC11 PC12
## Eksport towarów i usług (% PKB) 0.06085436 0.152378993
## import towarów i usług (% PKB) 0.18699434 0.223048443
## Bezrobocie gółem -0.10423538 0.311824645
## Bezrobocie, ogółem wśród młodzieży -0.14657645 -0.395801486
## Bezpośrednie inwestycje zagraniczne -0.32026931 -0.001954585
## Dług rządu centralnego 0.07255776 -0.126777862
## rządowe wydatki konsumpcyjne -0.24186895 0.022792763
## wydatki konsumpcyjne ogółem (%PKB) 0.54973853 0.091420766
## PKB per capita 0.17090928 -0.444811343
## Inflacja -0.04118267 0.113051428
## Oczekiwana długość życia -0.19744413 0.606883998
## Eksport żywności -0.09470535 -0.135452506
## Eksport wyrobów przemysłowych -0.40655189 -0.140168913
## Eksport usług 0.03084750 0.032532969
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana -0.45939602 -0.012342054
## Zatrudnienie w rolnictwie 0.06549889 0.173296065
## PC13 PC14
## Eksport towarów i usług (% PKB) -0.121034125 0.07523805
## import towarów i usług (% PKB) -0.214827651 0.16587221
## Bezrobocie gółem 0.474043957 -0.34138181
## Bezrobocie, ogółem wśród młodzieży -0.461860409 0.25062949
## Bezpośrednie inwestycje zagraniczne -0.092365185 -0.30876380
## Dług rządu centralnego 0.034986020 -0.18662938
## rządowe wydatki konsumpcyjne 0.004909706 0.06201521
## wydatki konsumpcyjne ogółem (%PKB) -0.197778127 0.11773323
## PKB per capita 0.111286432 -0.28678716
## Inflacja 0.047287828 0.09623750
## Oczekiwana długość życia -0.221541543 0.16772551
## Eksport żywności 0.397031232 0.61706918
## Eksport wyrobów przemysłowych 0.218306497 0.30559779
## Eksport usług -0.175152916 0.02581040
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana -0.369267150 -0.21147559
## Zatrudnienie w rolnictwie 0.127442957 0.02008734
## PC15 PC16
## Eksport towarów i usług (% PKB) 0.05224459 -0.7531647511
## import towarów i usług (% PKB) 0.18624627 0.6325848656
## Bezrobocie gółem 0.19697929 0.0269731207
## Bezrobocie, ogółem wśród młodzieży -0.17251697 0.0017519760
## Bezpośrednie inwestycje zagraniczne 0.04034117 -0.0068110648
## Dług rządu centralnego -0.15679398 0.0235098040
## rządowe wydatki konsumpcyjne 0.03816247 -0.0125202567
## wydatki konsumpcyjne ogółem (%PKB) 0.56711988 -0.1616301646
## PKB per capita 0.40204475 0.0457958185
## Inflacja -0.01591706 -0.0419446619
## Oczekiwana długość życia 0.04671937 0.0078900311
## Eksport żywności 0.08814406 0.0119469495
## Eksport wyrobów przemysłowych 0.43537040 -0.0230855492
## Eksport usług 0.01205254 -0.0004394145
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana 0.40906181 0.0103517784
## Zatrudnienie w rolnictwie -0.12592521 0.0163269604
# pomiar rzetelności skali miarą alfa-Cronbacha oraz pomiar zasadności stosowania PCA miarą KMO
#?install.packages("psych")
library(psych)
## Warning: pakiet 'psych' został zbudowany w wersji R 4.3.3
#install.packages("ltm")
library(ltm)
## Warning: pakiet 'ltm' został zbudowany w wersji R 4.3.3
## Ładowanie wymaganego pakietu: MASS
## Warning: pakiet 'MASS' został zbudowany w wersji R 4.3.2
## Ładowanie wymaganego pakietu: msm
## Warning: pakiet 'msm' został zbudowany w wersji R 4.3.3
## Ładowanie wymaganego pakietu: polycor
## Warning: pakiet 'polycor' został zbudowany w wersji R 4.3.3
##
## Dołączanie pakietu: 'polycor'
## Następujący obiekt został zakryty z 'package:psych':
##
## polyserial
##
## Dołączanie pakietu: 'ltm'
## Następujący obiekt został zakryty z 'package:psych':
##
## factor.scores
miara KMO
KMO(dane2011) #>0.5 wówczas uzasadnione jest stosowanie analizy głównych składowych
## Kaiser-Meyer-Olkin factor adequacy
## Call: KMO(r = dane2011)
## Overall MSA = 0.59
## MSA for each item =
## Eksport towarów i usług (% PKB)
## 0.57
## import towarów i usług (% PKB)
## 0.57
## Bezrobocie gółem
## 0.50
## Bezrobocie, ogółem wśród młodzieży
## 0.59
## Bezpośrednie inwestycje zagraniczne
## 0.35
## Dług rządu centralnego
## 0.45
## rządowe wydatki konsumpcyjne
## 0.74
## wydatki konsumpcyjne ogółem (%PKB)
## 0.56
## PKB per capita
## 0.58
## Inflacja
## 0.65
## Oczekiwana długość życia
## 0.84
## Eksport żywności
## 0.59
## Eksport wyrobów przemysłowych
## 0.46
## Eksport usług
## 0.88
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana
## 0.57
## Zatrudnienie w rolnictwie
## 0.67
Cronbach’s alpha
cronbach.alpha(dane2011) #>0.8 skale rzetelne, mierzą to samo #>0.8 skale rzetelne, mierzą to samo
##
## Cronbach's alpha for the 'dane2011' data-set
##
## Items: 16
## Sample units: 27
## alpha: 0
## Number of categories should be increased in order to count frequencies.
## Warning in alpha(dane2011): Some items were negatively correlated with the first principal component and probably
## should be reversed.
## To do this, run the function again with the 'check.keys=TRUE' option
## Some items ( Eksport towarów i usług (% PKB) import towarów i usług (% PKB) Bezpośrednie inwestycje zagraniczne rządowe wydatki konsumpcyjne PKB per capita Oczekiwana długość życia Eksport wyrobów przemysłowych Eksport usług ) were negatively correlated with the first principal component and
## probably should be reversed.
## To do this, run the function again with the 'check.keys=TRUE' option
##
## Reliability analysis
## Call: alpha(x = dane2011)
##
## raw_alpha std.alpha G6(smc) average_r S/N ase mean sd median_r
## 2.5e-07 -0.51 0.78 -0.022 -0.34 1.2e-07 3.8e+09 4.4e+09 -0.042
##
## 95% confidence boundaries
## lower alpha upper
## Feldt -0.65 0 0.47
## Duhachek 0.00 0 0.00
##
## Reliability if an item is dropped:
## raw_alpha std.alpha G6(smc)
## Eksport towarów i usług (% PKB) 2.5e-07 -0.2849 0.77
## import towarów i usług (% PKB) 2.5e-07 -0.2770 0.78
## Bezrobocie gółem 2.5e-07 -0.6812 0.69
## Bezrobocie, ogółem wśród młodzieży 2.5e-07 -0.7379 0.70
## Bezpośrednie inwestycje zagraniczne 2.5e-07 -0.6658 0.75
## Dług rządu centralnego 2.5e-07 -1.1429 0.67
## rządowe wydatki konsumpcyjne 2.5e-07 -0.4106 0.80
## wydatki konsumpcyjne ogółem (%PKB) 2.5e-07 -0.5187 0.73
## PKB per capita -2.8e-10 -0.3918 0.76
## Inflacja 2.5e-07 -0.4251 0.79
## Oczekiwana długość życia 2.5e-07 -0.5389 0.76
## Eksport żywności 2.5e-07 -1.0714 0.64
## Eksport wyrobów przemysłowych 2.5e-07 0.0069 0.82
## Eksport usług 1.7e-03 -0.2904 0.83
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana 2.5e-07 -0.4270 0.75
## Zatrudnienie w rolnictwie 2.5e-07 -0.3766 0.80
## average_r S/N alpha se
## Eksport towarów i usług (% PKB) -0.01500 -0.222 1.2e-07
## import towarów i usług (% PKB) -0.01467 -0.217 1.2e-07
## Bezrobocie gółem -0.02776 -0.405 1.2e-07
## Bezrobocie, ogółem wśród młodzieży -0.02913 -0.425 1.2e-07
## Bezpośrednie inwestycje zagraniczne -0.02737 -0.400 1.2e-07
## Dług rządu centralnego -0.03687 -0.533 1.2e-07
## rządowe wydatki konsumpcyjne -0.01979 -0.291 1.2e-07
## wydatki konsumpcyjne ogółem (%PKB) -0.02330 -0.342 1.2e-07
## PKB per capita -0.01912 -0.281 5.9e-09
## Inflacja -0.02029 -0.298 1.2e-07
## Oczekiwana długość życia -0.02390 -0.350 1.2e-07
## Eksport żywności -0.03571 -0.517 1.2e-07
## Eksport wyrobów przemysłowych 0.00046 0.007 1.2e-07
## Eksport usług -0.01523 -0.225 1.2e-03
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana -0.02035 -0.299 1.2e-07
## Zatrudnienie w rolnictwie -0.01858 -0.274 1.2e-07
## var.r med.r
## Eksport towarów i usług (% PKB) 0.12 -0.0030
## import towarów i usług (% PKB) 0.13 -0.0030
## Bezrobocie gółem 0.13 -0.0390
## Bezrobocie, ogółem wśród młodzieży 0.13 -0.0390
## Bezpośrednie inwestycje zagraniczne 0.14 -0.1032
## Dług rządu centralnego 0.13 -0.1073
## rządowe wydatki konsumpcyjne 0.13 -0.0390
## wydatki konsumpcyjne ogółem (%PKB) 0.12 -0.0572
## PKB per capita 0.12 -0.0572
## Inflacja 0.13 -0.0572
## Oczekiwana długość życia 0.12 -0.0390
## Eksport żywności 0.13 -0.1073
## Eksport wyrobów przemysłowych 0.13 -0.0016
## Eksport usług 0.13 -0.0030
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana 0.12 -0.0390
## Zatrudnienie w rolnictwie 0.13 -0.0352
##
## Item statistics
## n raw.r std.r r.cor
## Eksport towarów i usług (% PKB) 27 -0.124 -7.4e-03 -0.0096
## import towarów i usług (% PKB) 27 -0.216 -1.8e-02 -0.0221
## Bezrobocie gółem 27 -0.181 4.0e-01 0.4225
## Bezrobocie, ogółem wśród młodzieży 27 -0.229 4.4e-01 0.4662
## Bezpośrednie inwestycje zagraniczne 27 -0.057 3.9e-01 0.3626
## Dług rządu centralnego 27 0.101 6.9e-01 0.7065
## rządowe wydatki konsumpcyjne 27 0.273 1.5e-01 0.0402
## wydatki konsumpcyjne ogółem (%PKB) 27 -0.161 2.6e-01 0.2823
## PKB per capita 27 0.384 1.2e-01 0.1196
## Inflacja 27 -0.529 1.6e-01 0.0877
## Oczekiwana długość życia 27 0.549 2.8e-01 0.2664
## Eksport żywności 27 -0.176 6.5e-01 0.7183
## Eksport wyrobów przemysłowych 27 0.295 -5.0e-01 -0.5904
## Eksport usług 27 1.000 -7.2e-05 -0.1467
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana 27 -0.540 1.6e-01 0.1623
## Zatrudnienie w rolnictwie 27 -0.390 1.1e-01 0.0304
## r.drop mean sd
## Eksport towarów i usług (% PKB) -0.124 6.0e+01 3.0e+01
## import towarów i usług (% PKB) -0.216 5.8e+01 2.5e+01
## Bezrobocie gółem -0.181 1.0e+01 4.3e+00
## Bezrobocie, ogółem wśród młodzieży -0.229 2.4e+01 9.6e+00
## Bezpośrednie inwestycje zagraniczne -0.057 1.1e+01 2.8e+01
## Dług rządu centralnego 0.101 6.8e+01 3.5e+01
## rządowe wydatki konsumpcyjne 0.273 2.1e+01 3.2e+00
## wydatki konsumpcyjne ogółem (%PKB) -0.161 7.6e+01 7.6e+00
## PKB per capita 0.384 3.0e+04 2.1e+04
## Inflacja -0.529 3.3e+00 9.7e-01
## Oczekiwana długość życia 0.549 7.9e+01 3.0e+00
## Eksport żywności -0.176 1.2e+01 9.1e+00
## Eksport wyrobów przemysłowych 0.295 6.8e+01 1.7e+01
## Eksport usług 0.384 6.0e+10 7.0e+10
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana -0.540 2.5e+00 1.5e+00
## Zatrudnienie w rolnictwie -0.390 6.4e+00 5.8e+00
Wyniki
# zamiast obrazka możemy tę informację mieć przedstawioną jako ramkę danych
summary(wynik)
## Importance of components:
## PC1 PC2 PC3 PC4 PC5 PC6 PC7
## Standard deviation 2.2062 1.8844 1.5663 1.20611 0.98455 0.86971 0.76118
## Proportion of Variance 0.3042 0.2219 0.1533 0.09092 0.06058 0.04727 0.03621
## Cumulative Proportion 0.3042 0.5262 0.6795 0.77040 0.83098 0.87825 0.91447
## PC8 PC9 PC10 PC11 PC12 PC13 PC14
## Standard deviation 0.65140 0.57703 0.49144 0.40873 0.28603 0.25770 0.1876
## Proportion of Variance 0.02652 0.02081 0.01509 0.01044 0.00511 0.00415 0.0022
## Cumulative Proportion 0.94099 0.96180 0.97689 0.98733 0.99245 0.99660 0.9988
## PC15 PC16
## Standard deviation 0.13133 0.04522
## Proportion of Variance 0.00108 0.00013
## Cumulative Proportion 0.99987 1.00000
par(mfcol=c(1,2))
# ten wykres przedstawia ile wariancji jest wyjaśnione przez kolejne zmienne
plot(wynik, main= "Ile wariancji jest wyjaśnione przez kolejne zmienne")
# narysujmy biplot dla tych wyników
biplot(wynik)
#narysuj wykres osypiska
screeplot(wynik, type="line", main="Wykres osypiska")
# macierz przekształcenia
wynik$rotation[,1:4]
## PC1 PC2
## Eksport towarów i usług (% PKB) 0.24901799 0.330226024
## import towarów i usług (% PKB) 0.20341504 0.360343085
## Bezrobocie gółem -0.26154124 0.007548019
## Bezrobocie, ogółem wśród młodzieży -0.27612385 -0.005380294
## Bezpośrednie inwestycje zagraniczne 0.02590095 -0.124860117
## Dług rządu centralnego -0.06266567 -0.373699780
## rządowe wydatki konsumpcyjne 0.12466293 -0.378969212
## wydatki konsumpcyjne ogółem (%PKB) -0.34183644 -0.279171625
## PKB per capita 0.36600407 -0.059140262
## Inflacja -0.26139967 0.259335067
## Oczekiwana długość życia 0.25063831 -0.369794413
## Eksport żywności -0.14833096 -0.255079290
## Eksport wyrobów przemysłowych 0.19617492 0.206234406
## Eksport usług 0.24191670 -0.217244010
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana -0.35974037 0.071554613
## Zatrudnienie w rolnictwie -0.31849882 0.129647320
## PC3 PC4
## Eksport towarów i usług (% PKB) 0.273747455 -0.246677762
## import towarów i usług (% PKB) 0.278769491 -0.252875331
## Bezrobocie gółem -0.155590658 -0.617751988
## Bezrobocie, ogółem wśród młodzieży -0.179166190 -0.573216328
## Bezpośrednie inwestycje zagraniczne 0.383282009 -0.047559906
## Dług rządu centralnego 0.155827274 -0.127716371
## rządowe wydatki konsumpcyjne 0.005427735 -0.004955756
## wydatki konsumpcyjne ogółem (%PKB) -0.105930127 0.020891432
## PKB per capita 0.177254310 -0.052383921
## Inflacja 0.202940611 0.098449182
## Oczekiwana długość życia -0.019143178 -0.060625314
## Eksport żywności 0.489587657 -0.056314233
## Eksport wyrobów przemysłowych -0.442682390 0.052290369
## Eksport usług -0.239700038 0.012566443
## Rolnictwo, leśnictwo i rybołówstwo, wartość dodana 0.188930079 0.229935843
## Zatrudnienie w rolnictwie -0.082572009 0.266325958
# ładunki, przystępny sposób wyboru ładunków do tworzenia głównej składowej
loadings(wynik)
## NULL
DANE PO REDUKCJI - 5 ZMIENNYCH
library(readxl)
redukcja2011 <- read_excel("C:/Users/HP/Downloads/WAD projekt dane.xlsx",
sheet = "2011 redukcja")
rownames(redukcja2011) <- redukcja2011[[1]]
## Warning: Setting row names on a tibble is deprecated.
redukcja2011$kraje<- as.factor(redukcja2011$kraje)
# Remove the first column (now it's in row names)
redukcja2011 <- redukcja2011[, -1]
library(readxl)
redukcja11 <- read_excel("C:/Users/HP/Downloads/Zeszyt1_with_index.xlsx")
wynik = prcomp(redukcja2011, scale=T)
miara KMO i Cronbach’s alpha
pomiar rzetelności skali miarą alfa-Cronbacha oraz pomiar zasadności
stosowania PCA miarą KMO
#install.packages("psych")
library(psych)
KMO(redukcja2011) #>0.5 wówczas uzasadnione jest stosowanie analizy głównych składowych
## Kaiser-Meyer-Olkin factor adequacy
## Call: KMO(r = redukcja2011)
## Overall MSA = 0.47
## MSA for each item =
## import towarów i usług (% PKB) Bezrobocie gółem
## 0.57 0.55
## rządowe wydatki konsumpcyjne Eksport żywności
## 0.48 0.31
## Inflacja
## 0.47
cronbach.alpha(redukcja2011) #>0.8 skale rzetelne, mierzą to samo #>0.8 skale rzetelne, mierzą to samo
##
## Cronbach's alpha for the 'redukcja2011' data-set
##
## Items: 5
## Sample units: 27
## alpha: -0.252
alpha(redukcja2011, check.keys=TRUE)
## Number of categories should be increased in order to count frequencies.
## Warning in alpha(redukcja2011, check.keys = TRUE): Some items were negatively correlated with the first principal component and were automatically reversed.
## This is indicated by a negative sign for the variable name.
##
## Reliability analysis
## Call: alpha(x = redukcja2011, check.keys = TRUE)
##
## raw_alpha std.alpha G6(smc) average_r S/N ase mean sd median_r
## 0.15 0.42 0.48 0.13 0.72 0.15 64 5.7 0.14
##
## 95% confidence boundaries
## lower alpha upper
## Feldt -0.49 0.15 0.57
## Duhachek -0.14 0.15 0.44
##
## Reliability if an item is dropped:
## raw_alpha std.alpha G6(smc) average_r S/N
## import towarów i usług (% PKB) 0.128 0.39 0.45 0.136 0.628
## Bezrobocie gółem 0.208 0.48 0.52 0.188 0.928
## rządowe wydatki konsumpcyjne- 0.055 0.05 0.10 0.013 0.053
## Eksport żywności- 0.068 0.49 0.50 0.192 0.951
## Inflacja 0.150 0.30 0.31 0.098 0.432
## alpha se var.r med.r
## import towarów i usług (% PKB) 0.22 0.057 0.17
## Bezrobocie gółem 0.14 0.051 0.17
## rządowe wydatki konsumpcyjne- 0.16 0.025 0.03
## Eksport żywności- 0.11 0.047 0.17
## Inflacja 0.16 0.030 0.14
##
## Item statistics
## n raw.r std.r r.cor r.drop mean sd
## import towarów i usług (% PKB) 27 0.923 0.53 0.30 0.16 57.7 24.61
## Bezrobocie gółem 27 0.029 0.41 0.11 -0.12 10.2 4.33
## rządowe wydatki konsumpcyjne- 27 0.487 0.79 0.83 0.39 120.5 3.20
## Eksport żywności- 27 0.442 0.40 0.13 0.14 128.7 9.11
## Inflacja 27 0.178 0.61 0.53 0.15 3.3 0.97
Wyniki na wykresach
# zamiast obrazka możemy tę informację mieć przedstawioną jako ramkę danych
summary(wynik)
## Importance of components:
## PC1 PC2 PC3 PC4 PC5
## Standard deviation 1.3070 1.1382 0.9782 0.8201 0.60572
## Proportion of Variance 0.3416 0.2591 0.1914 0.1345 0.07338
## Cumulative Proportion 0.3416 0.6007 0.7921 0.9266 1.00000
# ten wykres przedstawia ile wariancji jest wyjaśnione przez kolejne zmienne
plot(wynik, main= "Ile wariancji jest wyjaśnione przez kolejne zmienne")
# narysujmy biplot dla tych wyników
par(mfcol=c(1,2))
biplot(wynik)
#narysuj wykres osypiska
screeplot(wynik, type="line", main="Wykres osypiska")
macierz przekształcenia
## PC1 PC2 PC3 PC4
## import towarów i usług (% PKB) 0.3848558 -0.48539544 0.3503057 0.6846328
## Bezrobocie gółem 0.2085484 0.54486177 -0.5877898 0.5582011
## rządowe wydatki konsumpcyjne -0.6712979 0.08437348 0.1097780 0.2213133
## Eksport żywności -0.1054178 0.59456125 0.6842880 0.2112431
## Inflacja 0.5887612 0.32694795 0.2269092 -0.3550862
# ładunki, przystępny sposób wyboru ładunków do tworzenia głównej składowej
loadings(wynik)
## NULL