1 Εισαγωγή δεδομένων

Στο A Comparison of Study Strategies for Passages: Rereading, Answering Questions, and Generating Questions εξετάζεται η αποτελεσματικότητα διαφόρων στρατηγικών μάθησης. Από το Statistics for the Behavioral Sciences αντλούμε τον πίνακα με τα αποτελέσματα κάποιων τεστ αξιολόγησης που εξήχθησαν.

ΔΙΑΒΑΣΜΑ ΔΙΑΒΑΣΜΑ ΚΑΙ ΑΠΑΝΤΗΣΗ ΕΡΩΤΗΣΕΩΝ ΔΙΑΒΑΣΜΑ, ΑΠΑΝΤΗΣΗ ΚΑΙ ΔΗΜΙΟΥΡΓΙΑ ΕΡΩΤΗΣΕΩΝ
2 5 8
3 9 6
8 10 12
6 13 11
5 8 11
6 9 12

Αυτόν τον εισάγουμε στην γράφοντας τα ακόλουθα:

diabasma <- c(2, 3, 8, 6, 5, 6)
diab_apantEr <- c(5, 9, 10, 13, 8, 9)
diab_apantEr_dimiourgEr <- c(8, 6, 12, 11, 11, 12)
dedomena <- data.frame(DIABASMA=diabasma,
                        DIABASMA_APANTISI_EROTISEON=diab_apantEr,
                        DIABASMA_APANTISI_EROTISEON_DIMIOURGIA_ER=diab_apantEr_dimiourgEr)

Δεχόμαστε ως μηδενική υπόθεση (\(\boldsymbol{\mathcal{H}_0}\)) την παραδοχή ότι δεν έχουν καμιά διαφορά οι τρεις αυτές στρατηγικές μάθησης. Η εναλλακτική υπόθεση, αυτό δηλαδή που περιμένουμε να συμβαίνει αν η \(\boldsymbol{\mathcal{H}_0}\) είναι εσφαλμένη, ισχυρίζεται πως, ναι, έχουν διαφορά αυτές οι τρεις στρατηγικές, παίζει, δηλαδή, ρόλο στα μαθησιακά μας αποτελέσματα το ποια υιοθετούμε.

Ακολούθως θα χρησιμοποιήσουμε τη συνάρτηση stack() που γνωρίσαμε στην ενότητα Αναδιαμόρφωση πίνακα, ώστε να μετατρέψουμε τον πίνακα σε μακράς μορφής. Ο λόγος είναι ότι θέλουμε να φτιάξουμε δύο στήλες:

  • μια στήλη με ενδείξεις τις τιμές του παράγοντα «στρατηγική μελέτης», δηλαδή τα «διάβασμα», «διάβασμα και απάντηση ερωτήσεων» και «διάβασμα, απάντηση και δημιουργία ερωτήσεων» (αυτή ονοματίστηκε αυτόματα ind) και

  • μία στήλη με τις επιδώσεις των υποκειμένων στην εκάστοτε στρατηγική μάθησης (αυτή ονοματίστηκε αυτόματα ως values).

dedomenaLong <- stack(dedomena)
values ind
2 DIABASMA
3 DIABASMA
8 DIABASMA
6 DIABASMA
5 DIABASMA
6 DIABASMA
5 DIABASMA_APANTISI_EROTISEON
9 DIABASMA_APANTISI_EROTISEON
10 DIABASMA_APANTISI_EROTISEON
13 DIABASMA_APANTISI_EROTISEON
8 DIABASMA_APANTISI_EROTISEON
9 DIABASMA_APANTISI_EROTISEON
8 DIABASMA_APANTISI_EROTISEON_DIMIOURGIA_ER
6 DIABASMA_APANTISI_EROTISEON_DIMIOURGIA_ER
12 DIABASMA_APANTISI_EROTISEON_DIMIOURGIA_ER
11 DIABASMA_APANTISI_EROTISEON_DIMIOURGIA_ER
11 DIABASMA_APANTISI_EROTISEON_DIMIOURGIA_ER
12 DIABASMA_APANTISI_EROTISEON_DIMIOURGIA_ER

2 Ανάλυση διασποράς – ANOVA

Μία πρώτη προσέγγιση στο ζήτημα θα επιτευχθεί μέσω της ANOVA. Αυτή προϋποθέτει την κανονικότητα του πληθυσμού, από τον οποίο αντλήθηκαν τα δείγματα. Οπότε είναι σημαντικό να έχει γίνει ο αντίστοιχος έλεγχος.

Θα χρησιμοποιήσουμε, λοιπόν, τη συνάρτηση aov() και αμέσως μετά με την summary() θα δούμε μια περίληψη των αποτελεσμάτων. Η συνάρτηση aov(🎻🎻🎻 ~ as.factor(🧑🧑🧑), data = 💻💻💻) εξετάζει:

  • μέσα από τον πίνακα δεδομένων 💻💻💻 (δηλαδή τον dedomenaLong)

  • την επίδραση του παράγοντα 🧑🧑🧑 (εδώ αυτός είναι οι τρεις στρατηγικές μελέτης που καταγράφονται στη στήλη ind)

  • στις τιμές 🎻🎻🎻 (εν προκειμένω, στην επίδωση των υποκειμένων όπως αυτή αναγράφεται στη στήλη values).

stratigikiANOVA <- aov(values ~ as.factor(ind), data = dedomenaLong)
summary(stratigikiANOVA)
##                Df Sum Sq Mean Sq F value  Pr(>F)   
## as.factor(ind)  2     84   42.00   7.159 0.00656 **
## Residuals      15     88    5.87                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Βλέπουμε λοιπόν ότι η p-τιμή είναι 0.00656 (=0.656%), άρα πολύ μικρότερη της 5% που έχουμε συνήθως ως επίπεδο σημαντικότητας. Επομένως ο παράγοντας «στρατηγική μάθησης» με τις στρατηγικές που αναφέραμε παίζει ρόλο στο τελικό μαθησιακό αποτέλεσμα.

3 Έλεγχος Kruskal-Wallis

Όπως είπαμε, η ANOVA έχει κάποιες απαιτήσεις από το δείγμα μας, που μπορεί να την καταστήσουν αναξιόπιστη. Αν δεν θέλουμε να μπλέξουμε με έλεγχο κανονικότητας, μπορούμε εναλλακτικά αυτής να εκτελέσουμε έναν άλλον έλεγχο, τον Kruskal-Wallis.

Εν προκειμένω θα χρησιμοποιήσουμε τη συνάρτηση kruskal.test(). Η συνάρτηση kruskal.test(🎻🎻🎻 ~ as.factor(🧑🧑🧑), data = 💻💻💻) εξετάζει:

  • μέσα από τον πίνακα δεδομένων 💻💻💻 (δηλαδή τον dedomenaLong)

  • την επίδραση του παράγοντα 🧑🧑🧑 (εδώ αυτός είναι οι τρεις στρατηγικές μελέτης που καταγράφονται στη στήλη ind)

  • στις τιμές 🎻🎻🎻 (εν προκειμένω, στην επίδωση των υποκειμένων όπως αυτή αναγράφεται στη στήλη values).

Γράφουμε:

kruskal.test(values ~ as.factor(ind), data = dedomenaLong)
## 
##  Kruskal-Wallis rank sum test
## 
## data:  values by as.factor(ind)
## Kruskal-Wallis chi-squared = 8.0559, df = 2, p-value = 0.01781

Και πάλι η μηδενική υπόθεση απορρίπτεται αφού η p-τιμή προέκυψε μικρότερη του 0.05. Έχουμε έναν λόγο παραπάνω, λοιπόν, να θεωρούμε πως ο παράγοντας «τρόπος διαβάσματος» παίζει ρόλο αν θα είναι το απλό διάβασμα ή το απλό διάβασμα σε συνδυασμό με απαντήσεις ερωτήσεων ή και σε συνδυασμό με την κατασκευή ερωτήσεων.

Συνολικά ο κώδικας που γράψαμε είναι ο παρακάτω:

diabasma <- c(2, 3, 8, 6, 5, 6)
diab_apantEr <- c(5, 9, 10, 13, 8, 9)
diab_apantEr_dimiourgEr <- c(8, 6, 12, 11, 11, 12)
dedomena <- data.frame(DIABASMA=diabasma,
                        DIABASMA_APANTISI_EROTISEON=diab_apantEr,
                        DIABASMA_APANTISI_EROTISEON_DIMIOURGIA_ER=diab_apantEr_dimiourgEr)
dedomenaLong <- stack(dedomena)
stratigikiANOVA <- aov(values ~ as.factor(ind), data = dedomenaLong)
summary(stratigikiANOVA)
kruskal.test(values ~ as.factor(ind), data = dedomenaLong) 

4 Έλεγχος Wilcoxon

4.1 Για αναξέρτητα δείγματα

Υπάρχουν περιπτώσεις που ο παράγοντας, του οποίου την επίδραση εξετάζουμε, παίρνει δύο τιμές. Π.χ. ο παράγοντας «φύλο» παίρνει τις τιμές «άρρεν» και «θήλυ», ο παράγοντας «είναι νόμιμη η κάνναβη» παίρνει τις τιμές «ναι» και «όχι» κ.τ.λ. Σ’ αυτές τις περιπτώσεις, αν και η ANOVA και ο έλεγχος Kruskal-Wallis αποδίδουν, θα απευθυνθούμε σ’ έναν άλλο που είναι εξειδικευμένος γι’ αυτές τις περιπτώσεις. Αυτός είναι ο έλεγχος Wilcoxon.

Εδώ αντλήσαμε δεδομένα από τον OECD σχετικά με το μορφωτικό επίπεδο των ενηλίκων και από την The World Bank – Gender Data Portal για το αν στην εκάστοτε χώρα οι γυναίκες αμείβονται εξίσου με τους άντρες για την ίδια δουλειά. Αυτά αποθηκεύτηκαν στα αρχεία epipedoEkp.csv και IsiAntimErg.csv, οι 100 πρώτες γραμμές των οποίων παρατίθενται παρακάτω.

  • epipedoEkp.csv
LOCATION INDICATOR SUBJECT MEASURE FREQUENCY TIME Value Flag.Codes
AUS EDUADULT BUPPSRY PC_25_64 A 1989 44.65064 NA
AUS EDUADULT BUPPSRY PC_25_64 A 1991 44.12706 NA
AUS EDUADULT BUPPSRY PC_25_64 A 1993 47.15905 NA
AUS EDUADULT BUPPSRY PC_25_64 A 1994 49.80203 NA
AUS EDUADULT BUPPSRY PC_25_64 A 1995 44.93585 NA
AUS EDUADULT BUPPSRY PC_25_64 A 1997 46.69527 NA
AUS EDUADULT BUPPSRY PC_25_64 A 1998 43.95596 NA
AUS EDUADULT BUPPSRY PC_25_64 A 1999 42.61008 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2000 41.21231 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2001 41.07614 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2002 39.05032 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2003 37.51982 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2004 35.87904 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2005 34.96797 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2006 33.28984 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2007 31.81673 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2008 30.06402 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2009 28.97252 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2010 26.79274 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2011 25.92254 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2012 23.56339 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2013 24.28443 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2014 22.89584 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2015 20.98090 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2016 20.06735 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2017 19.01266 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2018 18.10864 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2019 17.13266 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2020 16.24498 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2021 15.46888 NA
AUS EDUADULT BUPPSRY PC_25_64 A 2022 14.76060 NA
AUS EDUADULT TRY PC_25_64 A 1989 30.58148 NA
AUS EDUADULT TRY PC_25_64 A 1991 31.15881 NA
AUS EDUADULT TRY PC_25_64 A 1993 22.46919 NA
AUS EDUADULT TRY PC_25_64 A 1994 23.07237 NA
AUS EDUADULT TRY PC_25_64 A 1995 24.31545 NA
AUS EDUADULT TRY PC_25_64 A 1997 24.30127 NA
AUS EDUADULT TRY PC_25_64 A 1998 25.42071 NA
AUS EDUADULT TRY PC_25_64 A 1999 26.69340 NA
AUS EDUADULT TRY PC_25_64 A 2000 27.47575 NA
AUS EDUADULT TRY PC_25_64 A 2001 28.97732 NA
AUS EDUADULT TRY PC_25_64 A 2002 30.75979 NA
AUS EDUADULT TRY PC_25_64 A 2003 31.31415 NA
AUS EDUADULT TRY PC_25_64 A 2004 30.82308 NA
AUS EDUADULT TRY PC_25_64 A 2005 31.69115 NA
AUS EDUADULT TRY PC_25_64 A 2006 33.02221 NA
AUS EDUADULT TRY PC_25_64 A 2007 33.73671 NA
AUS EDUADULT TRY PC_25_64 A 2008 36.14999 NA
AUS EDUADULT TRY PC_25_64 A 2009 36.88135 NA
AUS EDUADULT TRY PC_25_64 A 2010 37.60431 NA
AUS EDUADULT TRY PC_25_64 A 2011 38.34207 NA
AUS EDUADULT TRY PC_25_64 A 2012 41.28236 NA
AUS EDUADULT TRY PC_25_64 A 2013 39.53993 NA
AUS EDUADULT TRY PC_25_64 A 2014 41.90185 NA
AUS EDUADULT TRY PC_25_64 A 2015 42.88876 NA
AUS EDUADULT TRY PC_25_64 A 2016 43.74390 NA
AUS EDUADULT TRY PC_25_64 A 2017 45.35567 NA
AUS EDUADULT TRY PC_25_64 A 2018 45.72748 NA
AUS EDUADULT TRY PC_25_64 A 2019 47.12998 NA
AUS EDUADULT TRY PC_25_64 A 2020 49.33745 NA
AUS EDUADULT TRY PC_25_64 A 2021 49.76787 NA
AUS EDUADULT TRY PC_25_64 A 2022 51.49590 NA
AUS EDUADULT UPPSRY PC_25_64 A 1989 24.76788 NA
AUS EDUADULT UPPSRY PC_25_64 A 1991 24.71413 NA
AUS EDUADULT UPPSRY PC_25_64 A 1993 30.37176 NA
AUS EDUADULT UPPSRY PC_25_64 A 1994 27.12561 NA
AUS EDUADULT UPPSRY PC_25_64 A 1995 30.74870 NA
AUS EDUADULT UPPSRY PC_25_64 A 1997 29.00345 NA
AUS EDUADULT UPPSRY PC_25_64 A 1998 30.62333 NA
AUS EDUADULT UPPSRY PC_25_64 A 1999 30.69652 NA
AUS EDUADULT UPPSRY PC_25_64 A 2000 31.31195 NA
AUS EDUADULT UPPSRY PC_25_64 A 2001 29.94654 NA
AUS EDUADULT UPPSRY PC_25_64 A 2002 30.18990 NA
AUS EDUADULT UPPSRY PC_25_64 A 2003 31.16603 NA
AUS EDUADULT UPPSRY PC_25_64 A 2004 33.29789 NA
AUS EDUADULT UPPSRY PC_25_64 A 2005 33.34088 NA
AUS EDUADULT UPPSRY PC_25_64 A 2006 33.68794 NA
AUS EDUADULT UPPSRY PC_25_64 A 2007 34.44656 NA
AUS EDUADULT UPPSRY PC_25_64 A 2008 33.78599 NA
AUS EDUADULT UPPSRY PC_25_64 A 2009 34.14612 NA
AUS EDUADULT UPPSRY PC_25_64 A 2010 35.60295 NA
AUS EDUADULT UPPSRY PC_25_64 A 2011 35.73539 NA
AUS EDUADULT UPPSRY PC_25_64 A 2012 35.15425 NA
AUS EDUADULT UPPSRY PC_25_64 A 2013 36.17564 NA
AUS EDUADULT UPPSRY PC_25_64 A 2014 35.20231 NA
AUS EDUADULT UPPSRY PC_25_64 A 2015 36.13035 NA
AUS EDUADULT UPPSRY PC_25_64 A 2016 36.18875 NA
AUS EDUADULT UPPSRY PC_25_64 A 2017 35.63167 NA
AUS EDUADULT UPPSRY PC_25_64 A 2018 36.16388 NA
AUS EDUADULT UPPSRY PC_25_64 A 2019 35.73737 NA
AUS EDUADULT UPPSRY PC_25_64 A 2020 34.41757 NA
AUS EDUADULT UPPSRY PC_25_64 A 2021 34.76325 NA
AUS EDUADULT UPPSRY PC_25_64 A 2022 33.74350 NA
AUT EDUADULT BUPPSRY PC_25_64 A 2004 23.66626 NA
AUT EDUADULT BUPPSRY PC_25_64 A 2005 23.05507 NA
AUT EDUADULT BUPPSRY PC_25_64 A 2006 19.87161 NA
AUT EDUADULT BUPPSRY PC_25_64 A 2007 20.06289 NA
AUT EDUADULT BUPPSRY PC_25_64 A 2008 19.12009 NA
AUT EDUADULT BUPPSRY PC_25_64 A 2009 18.35498 NA
AUT EDUADULT BUPPSRY PC_25_64 A 2010 17.64749 NA
  • IsiAntimErg.csv
Indicator.Name Indicator.Code Country.Name Country.Code Year Value Disaggregation
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2022 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2021 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2020 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2019 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2018 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2017 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2016 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2015 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2014 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2013 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2012 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2011 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2010 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2009 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2008 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2007 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2006 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2005 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2004 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2003 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2002 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2001 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 2000 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1999 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1998 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1997 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1996 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1995 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1994 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1993 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1992 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1991 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1990 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1989 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1988 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1987 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1986 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1985 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1984 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1983 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1982 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1981 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1980 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1979 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1978 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1977 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1976 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1975 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1974 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1973 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1972 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1971 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Afghanistan AFG 1970 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2022 1 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2021 1 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2020 1 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2019 1 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2018 1 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2017 1 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2016 1 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2015 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2014 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2013 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2012 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2011 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2010 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2009 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2008 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2007 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2006 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2005 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2004 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2003 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2002 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2001 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 2000 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1999 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1998 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1997 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1996 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1995 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1994 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1993 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1992 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1991 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1990 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1989 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1988 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1987 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1986 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1985 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1984 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1983 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1982 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1981 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1980 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1979 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1978 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1977 0 NA
Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Albania ALB 1976 0 NA

Θα εξετάσουμε κατά πόσον ο παράγοντας «οι γυναίκες αμείβονται εξίσου με τους άντρες για την ίδια δουλειά» (με τιμές 1=«ναι» και 0=«όχι») επιδρά στο ποσοστό των γυναικών που έχουν τελειώσει την τριτοβάθμια εκπαίδευση.

Αρχικά τροποποιούμε και ενώνουμε τους δύο πίνακες σε έναν, γράφοντας:

names(IsiAntimErg)[names(IsiAntimErg) == "Country.Code"] <- "LOCATION"
names(IsiAntimErg)[names(IsiAntimErg) == "Year"] <- "TIME"
epipedoEkpALL <- epipedoEkp[epipedoEkp$SUBJECT == "TRY",]
EkpErgALL <- merge(IsiAntimErg, epipedoEkpALL, by = c("LOCATION","TIME"))
LOCATION TIME Indicator.Name Indicator.Code Country.Name Value.x Disaggregation INDICATOR SUBJECT MEASURE FREQUENCY Value.y Flag.Codes
ARG 2004 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Argentina 1 NA EDUADULT TRY PC_25_64 A 16.711861 NA
ARG 2005 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Argentina 1 NA EDUADULT TRY PC_25_64 A 17.471054 NA
ARG 2006 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Argentina 1 NA EDUADULT TRY PC_25_64 A 17.830963 NA
ARG 2008 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Argentina 1 NA EDUADULT TRY PC_25_64 A 18.973812 NA
ARG 2009 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Argentina 1 NA EDUADULT TRY PC_25_64 A 19.518599 NA
ARG 2010 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Argentina 1 NA EDUADULT TRY PC_25_64 A 20.111586 NA
ARG 2011 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Argentina 1 NA EDUADULT TRY PC_25_64 A 20.237337 NA
ARG 2012 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Argentina 1 NA EDUADULT TRY PC_25_64 A 20.556126 NA
ARG 2013 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Argentina 1 NA EDUADULT TRY PC_25_64 A 21.056189 NA
ARG 2014 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Argentina 1 NA EDUADULT TRY PC_25_64 A 20.534544 NA
ARG 2017 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Argentina 1 NA EDUADULT TRY PC_25_64 A 21.406713 NA
ARG 2018 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Argentina 1 NA EDUADULT TRY PC_25_64 A 21.924601 NA
ARG 2019 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Argentina 1 NA EDUADULT TRY PC_25_64 A 22.811752 NA
ARG 2020 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Argentina 1 NA EDUADULT TRY PC_25_64 A 22.400209 NA
ARG 2021 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Argentina 1 NA EDUADULT TRY PC_25_64 A 24.765341 NA
AUS 1989 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 0 NA EDUADULT TRY PC_25_64 A 30.581482 NA
AUS 1991 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 0 NA EDUADULT TRY PC_25_64 A 31.158810 NA
AUS 1993 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 0 NA EDUADULT TRY PC_25_64 A 22.469191 NA
AUS 1994 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 0 NA EDUADULT TRY PC_25_64 A 23.072365 NA
AUS 1995 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 0 NA EDUADULT TRY PC_25_64 A 24.315445 NA
AUS 1997 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 24.301270 NA
AUS 1998 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 25.420708 NA
AUS 1999 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 26.693399 NA
AUS 2000 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 27.475746 NA
AUS 2001 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 28.977322 NA
AUS 2002 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 30.759785 NA
AUS 2003 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 31.314152 NA
AUS 2004 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 30.823080 NA
AUS 2005 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 31.691153 NA
AUS 2006 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 33.022213 NA
AUS 2007 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 33.736710 NA
AUS 2008 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 36.149986 NA
AUS 2009 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 36.881355 NA
AUS 2010 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 37.604309 NA
AUS 2011 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 38.342072 NA
AUS 2012 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 41.282364 NA
AUS 2013 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 39.539928 NA
AUS 2014 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 41.901855 NA
AUS 2015 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 42.888756 NA
AUS 2016 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 43.743900 NA
AUS 2017 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 45.355675 NA
AUS 2018 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 45.727478 NA
AUS 2019 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 47.129978 NA
AUS 2020 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 49.337452 NA
AUS 2021 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 49.767868 NA
AUS 2022 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Australia 1 NA EDUADULT TRY PC_25_64 A 51.495899 NA
AUT 2004 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 24.783438 NA
AUT 2005 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 24.639351 NA
AUT 2006 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 25.356432 NA
AUT 2007 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 25.092672 NA
AUT 2008 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 26.191858 NA
AUT 2009 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 26.895824 NA
AUT 2010 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 27.743721 NA
AUT 2011 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 27.956152 NA
AUT 2012 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 28.739943 NA
AUT 2013 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 29.738844 NA
AUT 2014 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 29.904896 NA
AUT 2015 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 30.550726 NA
AUT 2016 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 31.383963 NA
AUT 2017 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 32.394390 NA
AUT 2018 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 32.711426 NA
AUT 2019 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 33.773777 NA
AUT 2020 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 34.205894 NA
AUT 2021 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 34.604504 NA
AUT 2022 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Austria 1 NA EDUADULT TRY PC_25_64 A 35.553291 NA
BEL 1989 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 0 NA EDUADULT TRY PC_25_64 A 16.864670 NA
BEL 1991 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 0 NA EDUADULT TRY PC_25_64 A 19.620007 NA
BEL 1992 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 0 NA EDUADULT TRY PC_25_64 A 20.226023 NA
BEL 1994 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 0 NA EDUADULT TRY PC_25_64 A 22.345255 NA
BEL 1995 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 0 NA EDUADULT TRY PC_25_64 A 24.604177 NA
BEL 2000 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 0 NA EDUADULT TRY PC_25_64 A 27.084993 NA
BEL 2001 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 0 NA EDUADULT TRY PC_25_64 A 27.626753 NA
BEL 2002 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 0 NA EDUADULT TRY PC_25_64 A 28.134819 NA
BEL 2003 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 0 NA EDUADULT TRY PC_25_64 A 29.031809 NA
BEL 2004 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 0 NA EDUADULT TRY PC_25_64 A 30.400860 NA
BEL 2005 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 0 NA EDUADULT TRY PC_25_64 A 31.045088 NA
BEL 2006 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 0 NA EDUADULT TRY PC_25_64 A 31.754845 NA
BEL 2007 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 0 NA EDUADULT TRY PC_25_64 A 32.094864 NA
BEL 2008 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 0 NA EDUADULT TRY PC_25_64 A 32.310043 NA
BEL 2009 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 1 NA EDUADULT TRY PC_25_64 A 33.397053 NA
BEL 2010 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 1 NA EDUADULT TRY PC_25_64 A 34.989502 NA
BEL 2011 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 1 NA EDUADULT TRY PC_25_64 A 34.612984 NA
BEL 2012 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 1 NA EDUADULT TRY PC_25_64 A 35.316551 NA
BEL 2013 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 1 NA EDUADULT TRY PC_25_64 A 35.535332 NA
BEL 2014 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 1 NA EDUADULT TRY PC_25_64 A 36.855022 NA
BEL 2015 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 1 NA EDUADULT TRY PC_25_64 A 36.857422 NA
BEL 2016 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 1 NA EDUADULT TRY PC_25_64 A 37.503956 NA
BEL 2017 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 1 NA EDUADULT TRY PC_25_64 A 40.257053 NA
BEL 2018 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 1 NA EDUADULT TRY PC_25_64 A 40.638546 NA
BEL 2019 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 1 NA EDUADULT TRY PC_25_64 A 40.670090 NA
BEL 2020 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 1 NA EDUADULT TRY PC_25_64 A 42.442284 NA
BEL 2021 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 1 NA EDUADULT TRY PC_25_64 A 44.873066 NA
BEL 2022 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Belgium 1 NA EDUADULT TRY PC_25_64 A 45.769352 NA
BGR 2015 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Bulgaria 1 NA EDUADULT TRY PC_25_64 A 27.498056 NA
BGR 2022 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Bulgaria 1 NA EDUADULT TRY PC_25_64 A 29.780050 NA
BRA 2007 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Brazil 0 NA EDUADULT TRY PC_25_64 A 9.623828 NA
BRA 2008 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Brazil 0 NA EDUADULT TRY PC_25_64 A 10.795494 NA
BRA 2009 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Brazil 0 NA EDUADULT TRY PC_25_64 A 10.863071 NA
BRA 2011 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Brazil 0 NA EDUADULT TRY PC_25_64 A 11.614742 NA
BRA 2012 Law mandates equal remuneration for females and males for work of equal value (1=yes; 0=no) SG.LAW.EQRM.WK Brazil 0 NA EDUADULT TRY PC_25_64 A 12.993286 NA

Ακολούθως γράφουμε:

wilcox.test(Value.y ~ Value.x, data = EkpErgALL)
## 
##  Wilcoxon rank sum test with continuity correction
## 
## data:  Value.y by Value.x
## W = 81579, p-value < 2.2e-16
## alternative hypothesis: true location shift is not equal to 0

Ως μηδενική υπόθεση πήραμε το ότι ο παράγοντας δεν επιδρά, δηλαδή ότι είτε οι γυναίκες αμείβονται εξίσου με τους άντρες είτε όχι, το ποσοστό γυναικών στην τριτοβάθμια δεν αλλάζει.

Δεδομένου αυτού, αυτή η μικροσκοπική p-τιμή που βρήκαμε κάνει εντελώς απίθανη αυτή την παραδοχή. Άρα όντως ο παράγοντας «οι γυναίκες αμείβονται εξίσου με τους άντρες για την ίδια δουλειά» επιδρά στο ποσοστό των γυναικών που έχουν τελειώσει την τριτοβάθμια εκπαίδευση.

Εδώ θα πρέπει να πούμε πως ο παράγοντας «οι γυναίκες αμείβονται εξίσου με τους άντρες για την ίδια δουλειά» και το ποσοστό των γυναικών που έχουν τελειώσει την τριτοβάθμια εκπαίδευση δεν έχουν σχέση αιτίου-αιτιατού. Αυτό που στην πραγματικότητα βρήκαμε είναι η αλληλεξάρτησή τους. Πολύ πιθανόν ένας άλλος παράγοντας που καθορίζει το ένα να διαμορφώνει και το άλλο.

4.2 Για ζευγαρωτά δεδομένα

Στο προηγούμενο παράδειγμα συγκρίναμε δύο δείγματα. Ένα δείγμα χωρών που οι γυναίκες αμείβονται εξίσου με τους άντρες για την ίδια δουλειά με ένα άλλο που οι γυναίκες δεν αμείβονται εξίσου με τους άντρες για την ίδια δουλειά. Προσπαθήσαμε να δούμε αν παίζει ρόλο αυτό το γεγονός στο ποσοστό των γυναικών που έχουν τελειώσει την τριτοβάθμια εκπαίδευση.

Τι γίνεται όμως αν δεν είχαμε δύο ξεχωριστά δείγματα χωρών, αλλά για την κάθε χώρα είχαμε για κάποια περίοδο ότι οι γυναίκες δεν αμείβονται εξίσου με τους άντρες για την ίδια δουλειά, ενώ για κάποια άλλη ότι αμείβονται, μαζί με τα αντίστοιχα ποσοστά των γυναικών που έχουν τελειώσει την τριτοβάθμια εκπαίδευση; Τότε δεν θα μιλάγαμε για δύο ξεχωριστά δείγματα, αλλά για ζευγάρια χωρών πριν και μετά από αυτή την αλλαγή πολιτικής προς τις γυναίκες.

Επειδή δεν μπορέσαμε να βρούμε δεδομένα σχετικά με αυτό που συζητάμε, ψάξαμε να βρούμε κάποια άλλα. Παρουσιάζουμε κάποια δεδομένα από τη σελίδα Υλικό Φυσικής-Χημείας του κυρίου Π. Μάλλιου. Ο εν λόγω καθηγητής κατέθεσε μια ερώτηση προς το ακροατήριο (https://forumning.wordpress.com/2016/10/21/πρέπει-αυτός-ο-μαθητής-να-προαχθεί-2/) «Θα έπρεπε ο μαθητής με αυτούς τους βαθμούς να προαχθεί;» παραθέτοντας τους βαθμούς του μαθητή.

Αρχικά τους περνάμε στην R, ως ακολούθως:

elegxos <- data.frame(
    mathimata = c("Θρησκευτικά", "Αρχαία Γενικής", "Νεοελληνική Γλώσσα", "Νεοελληνική Λογοτεχνία", "Άλγεβρα", "Γεωμετρία", "Φυσική", "Χημεία", "Βιολογία", "Η/Υ", "Ιστορία", "Φιλοσοφία", "Πολιτική Παιδεία", "Αρχαία Προσανατολισμού", "Κοινωνιολογία"),
    a_tetramino = c(13,12,10,11,11,11,10,11,12,13,11,10,14,13,14),
    b_tetramino = c(13,13,12,13,13,13,12,12,12,16,12,12,14,13,14),
    exetaseis = c(4, 6.6, 7, 12, 6, 3, 1, 3, 3, 5.6, 5, 6, 1, 10, 10)
  )
elegxos
##                 mathimata a_tetramino b_tetramino exetaseis
## 1             Θρησκευτικά          13          13       4.0
## 2          Αρχαία Γενικής          12          13       6.6
## 3      Νεοελληνική Γλώσσα          10          12       7.0
## 4  Νεοελληνική Λογοτεχνία          11          13      12.0
## 5                 Άλγεβρα          11          13       6.0
## 6               Γεωμετρία          11          13       3.0
## 7                  Φυσική          10          12       1.0
## 8                  Χημεία          11          12       3.0
## 9                Βιολογία          12          12       3.0
## 10                    Η/Υ          13          16       5.6
## 11                Ιστορία          11          12       5.0
## 12              Φιλοσοφία          10          12       6.0
## 13       Πολιτική Παιδεία          14          14       1.0
## 14 Αρχαία Προσανατολισμού          13          13      10.0
## 15          Κοινωνιολογία          14          14      10.0

Αυτό που θέλουμε εν προκειμένω να δούμε είναι το αν ο παράγοντας «εξετάσεις» έχει επίδραση στη βαθμολογία του μαθητή. Δηλαδή, το τάδε μάθημα στις εξετάσεις έχει διαφορετική βαθμολογία σε σύγκριση με τα τετράμηνα;

Καταρχάς ας συμπυκνώσουμε τη βαθμολογία των δύο τετραμήνων στο μέσο όρο τους.

elegxos$tetramina <- (elegxos$a_tetramino+elegxos$b_tetramino)/2
elegxos
##                 mathimata a_tetramino b_tetramino exetaseis tetramina
## 1             Θρησκευτικά          13          13       4.0      13.0
## 2          Αρχαία Γενικής          12          13       6.6      12.5
## 3      Νεοελληνική Γλώσσα          10          12       7.0      11.0
## 4  Νεοελληνική Λογοτεχνία          11          13      12.0      12.0
## 5                 Άλγεβρα          11          13       6.0      12.0
## 6               Γεωμετρία          11          13       3.0      12.0
## 7                  Φυσική          10          12       1.0      11.0
## 8                  Χημεία          11          12       3.0      11.5
## 9                Βιολογία          12          12       3.0      12.0
## 10                    Η/Υ          13          16       5.6      14.5
## 11                Ιστορία          11          12       5.0      11.5
## 12              Φιλοσοφία          10          12       6.0      11.0
## 13       Πολιτική Παιδεία          14          14       1.0      14.0
## 14 Αρχαία Προσανατολισμού          13          13      10.0      13.0
## 15          Κοινωνιολογία          14          14      10.0      14.0

Έτσι, εκτελούμε τον κάτωθι έλεγχο Wilcoxon με την προσθήκη της παραμέτρου paired = TRUE.

wilcox.test(elegxos$tetramina, elegxos$exetaseis, paired = TRUE)
## 
##  Wilcoxon signed rank test with continuity correction
## 
## data:  elegxos$tetramina and elegxos$exetaseis
## V = 105, p-value = 0.001082
## alternative hypothesis: true location shift is not equal to 0

Προκύπτει p-τιμή 0.0010816, δηλαδή αρκετά μικρότερη από το 5%, από την πιθανότητα που θα θεωρούσαμε ότι μπορεί να δικαιολογήσει την απόκλιση που βλέπουμε, αν οι εξετάσεις δεν διαφοροποιούνταν από τα τετράμηνα. Με άλλα λόγια, οι εξετάσεις παίζουν ρόλο στη βαθμολογία, διότι στην αντίθετη περίπτωση θα ήταν εντελώς απίθανο να συναντήσουμε αυτή τη διαφοροποίηση που βλέπουμε.

Συνολικά σε αυτή την ενότητα γράψαμε τον κώδικα:

elegxos <- data.frame(
    mathimata = c("Θρησκευτικά", "Αρχαία Γενικής", "Νεοελληνική Γλώσσα", "Νεοελληνική Λογοτεχνία", "Άλγεβρα", "Γεωμετρία", "Φυσική", "Χημεία", "Βιολογία", "Η/Υ", "Ιστορία", "Φιλοσοφία", "Πολιτική Παιδεία", "Αρχαία Προσανατολισμού", "Κοινωνιολογία"),
    a_tetramino = c(13,12,10,11,11,11,10,11,12,13,11,10,14,13,14),
    b_tetramino = c(13,13,12,13,13,13,12,12,12,16,12,12,14,13,14),
    exetaseis = c(4, 6.6, 7, 12, 6, 3, 1, 3, 3, 5.6, 5, 6, 1, 10, 10)
  )
elegxos
elegxos$tetramina <- (elegxos$a_tetramino+elegxos$b_tetramino)/2
elegxos
wilcox.test(elegxos$tetramina, elegxos$exetaseis, paired = TRUE)

5 Μέγεθος επίδρασης

5.1 Cohen d

Ας επιστρέψουμε τώρα στο παράδειγμα σχετικά με το μορφωτικό επίπεδο των ενηλίκων. Ένα ζήτημα που ανακύπτει κατόπιν της διαπίστωσης της επίδρασης ενός παράγοντα είναι το πόσο πολύ αυτός επιδρά. Ένας τρόπος, όπως έχουμε πει (βλ. επίσης υποενότητα Μέγεθος επίδρασης), να το ελέγξουμε αυτό είναι υπολογίζοντας το d του Cohen. Αυτό για να υπολογιστεί χρειάζεται το πακέτο effsize, οπότε γράφουμε:

if(!require(effsize)){
    install.packages("effsize")
    library(effsize)
}

Ακολούθως υπολογίζουμε το d του Cohen χρησιμοποιώντας τη συνάρτηση cohen.d(🎻🎻🎻 ~ as.factor(🧑🧑🧑), data = 💻💻💻), όπου:

  • ο πίνακας δεδομένων 💻💻💻 (δηλαδή ο EkpErgALL)

  • ο προς εξέταση παράγοντας 🧑🧑🧑 (εδώ αυτός είναι οι τιμές Value.x ως προς τα αν οι γυναίκες αμοίβονται το ίδιο καλά με τους άντρες)

  • στις τιμές 🎻🎻🎻 (εν προκειμένω, το ποσοστό των γυναικών που έχουν τελειώσει γ-βάθμια εκπαίδευση, ήτοι η στήλη Value.y). `

Γράφουμε:

cohen.d(Value.y ~ as.factor(Value.x), data = EkpErgALL)
## 
## Cohen's d
## 
## d estimate: -0.713341 (medium)
## 95 percent confidence interval:
##      lower      upper 
## -0.8399722 -0.5867098

Οπότε μας εξάγεται το αποτέλεσμα -0.713341. Αυτό υποδηλώνει μεγάλο προς μεσαίο μέγεθος επίδρασης, όπως φαίνεται και από τις κάτωθι οδηγίες.

  • Αν το d του Cohen είναι 0, ±0.1, τότε έχουμε ασήμαντο μέγεθος επίδρασης.

  • Αν το d του Cohen είναι ±0.2, ±0.3, ±0.4 τότε έχουμε μικρό μέγεθος επίδρασης.

  • Αν το d του Cohen είναι ±0.5, ±0.6, ±0.7 τότε έχουμε μεσαίο μέγεθος επίδρασης.

  • Αν το d του Cohen είναι ±0.8, ±0.9, ±1, ±1.1 κ.τ.λ. τότε έχουμε μεγάλο μέγεθος επίδρασης.

5.2 Wilcoxon r

Γενικώς το d του Cohen είναι καλό να υπολογίζεται όταν έχουμε κανονικούς πληθυσμούς. Εναλλακτικά από αυτό μπορούμε να υπολογίζουμε το r του Wilcoxon. Αυτό με τη σειρά του απαιτεί το πακέτο rstatix, οπότε γράφουμε:

if(!require(rstatix)){
    install.packages("rstatix")
    library(rstatix)
}

Τώρα θα χρησιμοποιήσουμε τη συνάρτηση wilcox_effsize(🎻🎻🎻 ~ 🧑🧑🧑, data = 💻💻💻), όπου:

  • ο πίνακας δεδομένων 💻💻💻 (δηλαδή ο EkpErgALL)

  • ο προς εξέταση παράγοντας 🧑🧑🧑 (εδώ αυτός είναι οι τιμές Value.x ως προς τα αν οι γυναίκες αμοίβονται το ίδιο καλά με τους άντρες)

  • στις τιμές 🎻🎻🎻 (εν προκειμένω, το ποσοστό των γυναικών που έχουν τελειώσει γ-βάθμια εκπαίδευση, ήτοι η στήλη Value.y). `

Γράφουμε λοιπόν:

wilcox_effsize(Value.y ~ Value.x, data = EkpErgALL)
## # A tibble: 1 × 7
##   .y.     group1 group2 effsize    n1    n2 magnitude
## * <chr>   <chr>  <chr>    <dbl> <int> <int> <ord>    
## 1 Value.y 0      1        0.343   408   676 moderate

Έτσι, προκύπτει η τιμή 0.343, η οποία υποδηλώνει μέτριο μέγεθος επίδρασης.

Γενικότερα για το r του Wilcoxon έχουμε:

  • Αν το r του Wilcoxon είναι ±0.1, ±0.2 τότε έχουμε μικρό μέγεθος επίδρασης.

  • Αν το r του Wilcoxon είναι ±0.3, ±0.4 τότε έχουμε μεσαίο μέγεθος επίδρασης.

  • Αν το r του Wilcoxon είναι ±0.5, ±0.6 κ.τ.λ. τότε έχουμε μεγάλο μέγεθος επίδρασης.

Συνολικά ο κώδικας αυτής ενότητας που αφορούσε τις αμοιβές και το εκπαιδευτικό επίπεδο των γυναικών (όχι αυτή με τη βαθμολογία του μαθητή) ήταν:

names(IsiAntimErg)[names(IsiAntimErg) == "Country Code"] <- "LOCATION"
names(IsiAntimErg)[names(IsiAntimErg) == "Year"] <- "TIME"
epipedoEkpALL <- epipedoEkp[epipedoEkp$SUBJECT == "TRY",]
EkpErgALL <- merge(IsiAntimErg, epipedoEkpALL, by = c("LOCATION","TIME"))
wilcox.test(Value.y ~ as.factor(Value.x), data = EkpErgALL)
install.packages("effsize")
library(effsize)
cohen.d(Value.y ~ as.factor(Value.x), data = EkpErgALL)
install.packages("rstatix")
library(rstatix)
wilcox_effsize(Value.y ~ Value.x, data = EkpErgALL)