Εμφάνιση αναρτήσεων με ετικέτα Παλινδρόμηση. Εμφάνιση όλων των αναρτήσεων
Εμφάνιση αναρτήσεων με ετικέτα Παλινδρόμηση. Εμφάνιση όλων των αναρτήσεων

Παρασκευή 27 Νοεμβρίου 2015

F-έλεγχος

Από Quantitative


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

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

Πλήρες Μοντέλο Y = a + b1x1 + ... + bmxm + bm + 1 + xm + 1 + ... + bkxk + v

Αναγωγικό Μοντέλο Υ = a + b1x1 + ... + bmxm + v: < math > bm + 1 = ... = bk = 0 έναντι της εναλλακτικής Α: τουλάχιστον μία από τις παραμέτρους b είναι διαφορετική του 0; Χρησιμοποιούμε την εξής στατιστική ελέγχου:

F = \frac{(SSE_A-SSE_P)/(k-m)}{SSE_P/[n-(k+1)]} = \frac{(SSE_A-SSE_P)/ w}{MSE_P}
όπου, SSEA = Άθροισμα τετραγώνων σφαλμάτων για το αναγωγικό μοντέλο
MSEP = Μέσο τετραγωνικό σφάλμα για το πλήρες μοντέλο
k-m = Πλήθος των προς έλεγχο b παραμέτρων
k+1 = Σύνολο παραμέτρων του μοντέλου μας
w = πλήθος των προς έλεγχο παραμέτρων

Πέμπτη 26 Νοεμβρίου 2015

Πρόβλεψη στο πολυδιάστατο μοντέλο παλινδρόμησης

Από Quantitative


Ας θεωρήσουμε ότι έχουμε δύο ανεξάρτητες μεταβλητές x1,x2 και θέλουμε να προβλέψουμε την τιμή της Y, με \hat{Y}, όταν x1 = x10 και x2 = x10.

Τότε μπορεί να δειχθεί οτι:

var(Y-\hat{Y})=s2\left(1+\frac{1}{n}\right)+(x_{10}-\bar{x_1})^2var(\hat{b_1})+2(x_{10}-\bar{x_1})(x_{20}-\bar{x_2})Cov(\hat{b_1},\hat{b_2})+(x_{20}-\bar{x_2})^2var\hat{b_2}

Αν s2 αντικατασταθεί με τον εκτιμητή του S2=\frac{SSE}{n-3}, λαμβάνουμε τον εκτιμητή S_{Y-\bar{Y}} για την var_{Y-\bar{Y}} και το διάστημα εμπιστοσύνης για την πρόβλεψη του y είναι το
[\bar{y}-t_{n-3;a/2} S_{Y-\bar{Y}},\bar{y}+t_{n-3;a/2} S_{Y-\bar{Y}}]

Αν είχαμε k ανεξάρτητες μεταβλητές τότε κάτω από κατάλληλες αλλαγές, όπως αντικαθιστώντας το tn − 3;a / 2 με tn − (k + 1);a / 2, έχουμε το αντίστοιχο διάστημα εμπιστοσύνης για την πρόβλεψη του yi.

Ανάλυση διακύμανσης και έλεγχοι υποθέσεων

Τετάρτη 25 Νοεμβρίου 2015

Στατιστικά συμπεράσματα στο πολυδιάστατο γραμμικό μοντέλο

Από Quantitative


Θεωρούμε οτι τα τυχαία σφάλματα vi του πολυδιάστατου γραμμικού μοντέλου Y = a + b1x1 + b2x2 + ... + bkxk + v ακολουθούν την Κανονική κατανομή Ν(0,1) και χωρίς βλάβη της γενικότητας να περιέχει μόνο τις δύο ανεξάρτητες τ.μ. X1,X2. Τότε μπορούμε εύκολα να συνοψίσουμε και να αποδείξουμε τα ακόλουθα:

1. Οι κατανομές των εκτιμητών των παραμέτρων

Οι εκτιμητές \hat{a},\hat{b_1},\hat{b_2} ακολουθούν Κανονικές κατανομές με μέσους a,b1 και b2.
Αν με r12 συμβολίσουμε τον συντελεστή συσχέτισης μεταξύ των ανεξάρτητων τυχαίων μεταβλητών X1 και X2, με

S11 = \sum_{i=1}^{n}x^2_{1i}-n\bar{x_1}^2 και S22 = \sum_{i=1}^{n}x^2_{2i}-n\bar{x_2}^2,

τότε οι διακυμάνσεις δίνονται από τις εξής σχέσεις

var(\hat{b_1}) = \frac{s^2}{S_{11}(1-r^2_{12})},

var(\hat{b_2}) = \frac{s^2}{S_{22}(1-r^2_{12})},

var(\hat{a}) = \frac{s^2}{n} + \bar {x_1}^2var(\hat{b_1} )+ 2\bar{x_1}\bar{x_2}Cov(\hat{b_1},\hat{b_2}) + \bar{x_2}^2var(\hat{b_2}),

Cov (\hat{b_1},\hat{b_2}) = \frac{-s^2r_{12}^2}{S_{12}(1-r_{12}^2)},
Cov (\hat{a},\hat{b_1}) = -[\bar{x_1}var(\hat{b_1}) + \bar{x_2}Cov(\hat{b_1},\hat{b_2})],

Cov (\hat{a},\hat{b_2}) = -[\bar{x_1}Cov(\hat{b_1},\hat{b_2}) + \bar{x_2}var(\hat{b_2})].

Τρίτη 24 Νοεμβρίου 2015

Εκτίμηση διακύμανσης της πολλαπλής παλινδρόμησης και πολυδιάστατος συντελεστής συσχέτισης

Από Quantitative


Γενικά

Έχουμε συναντήσει οτι με τη βοήθεια της Μεθόδου των Ελαχίστων Τετραγώνων το πολυδιάστατο γραμμικό μοντέλο εκτιμάται από την:
\hat{Y} = \hat{a} + \hat{b_1}x_1+...+\hat{b_k}x_k,
όπου τα \hat{a},\hat{b_1},...,\hat{b_k} επιλέγονται έτσι ώστε να ελαχιστοποιούν την:
SSE = \sum_{i=1}^{n}(Y_i-\hat{Y_i})^2
Η τελευταία ποσότητα, όπως και στην απλή παλινδρόμηση, χρησιμοποιείται για να εκτιμήσουμε τη διακύμανση s2. Έτσι ένας εκτιμητής για το s2 της Πολυδιάστατης Γραμμικής Παλινδρόμησης με k ανεξάρτητες μεταβλητές ορίζεται από τη σχέση:
S2 = (S^2_{Y(1...k)}) = \frac{SSE}{n- (plithos~parametrwn)} = \frac{SSE}{n-(k+1)}

Γράψαμε το συμβολισμό μέσα στην παρένθεση για να υπενθυμίσουμε στον αναγνώστη οτι έχουμε τη διακύμανση των τιμών της Y σε σχέση με τις k ανεξάρτητες μεταβλητές.
Η μέση μεταβλητότητα είναι αμερόληπτος εκτιμητής της διακύμανσης,S2, των σφαλμάτων του πληθυσμού , e , και συμβολίζεται : MSE=\frac {SSE} {n-(k-1)}=\frac {\sum_(Y-\hat{Y})^2}{n-(k-1)}

Προσεγγιστικά : S=\sqrt MSE<math>    Ο δειγματοληπτικός συντελεστής συσχέτισης <math>r^2, αποτελεί ένα μέτρο κατά πόσο τα δεδομένα μας είναι προσαρμοσμένα γύρω από την ευθεία παλινδρόμησης. Για να μετρήσουμε πόσο καλά προσαρμόζονται γύρω από το υπερεπίπεδο παλινδρόμησης χρησιμοποιούμε τον πολυδιάστατο συντελεστή παλινδρόμησης που το συμβολίζουμε με R2.