- LB
- Posts : 22
Join date : 2010-02-11
Cook's distance for linear mixed effects models
Fri 29 Jun 2012 - 0:43
Καλησπέρα σε όλους και όλες.
Αυτή τη περίοδο ασχολούμαι με linear mixed effects models πάνω σε ιατρικά δεδομένα. Χρησιμοποιώ την library nlme για την εφαρμογη του μοντέλου lme(.) Σύμφωνα με τα diagnostic plots, υπάρχουν κάποιες παρατηρήσεις που φαίνεται να είναι outliers.
Έχω την πεποίθηση ότι αφορούν συγκεκριμένο άτομο, του οποίου τη συμπεριφορά την έχω δει μέσα από profile plots. Ωστόσο, ψάχνω να αποδείξω ότι το συγκεκριμένο άτομο είναι outlier μέσω του Cook's distance.
Η απορία μου είναι λοιπόν η εξής: Η αντίστοιχη εντολή που χρησιμοποιούμε στα απλά lm(.) δεν δουλεύει σε αυτή την περίπτωση. Όσο έψαχνα βρήκα ότι μπορώ να υπολογίσω το Cook's distance αλλά για ΜΟΝΟ μοντέλα με τη χρήση του πακέτου lme4. Υπάρχει τρόπος να βρώ το cook's d για το πακέτο nlme? Ή απλά επικαλούμαι τα profile plots από την exploratory analysis που έκανα αρχικά για να αιτιολογήσω την απόφασή μου να βγάλω τις παρατηρήσεις του ατόμου αυτού από το dataset?
Ευχαριστώ.
Αυτή τη περίοδο ασχολούμαι με linear mixed effects models πάνω σε ιατρικά δεδομένα. Χρησιμοποιώ την library nlme για την εφαρμογη του μοντέλου lme(.) Σύμφωνα με τα diagnostic plots, υπάρχουν κάποιες παρατηρήσεις που φαίνεται να είναι outliers.
Έχω την πεποίθηση ότι αφορούν συγκεκριμένο άτομο, του οποίου τη συμπεριφορά την έχω δει μέσα από profile plots. Ωστόσο, ψάχνω να αποδείξω ότι το συγκεκριμένο άτομο είναι outlier μέσω του Cook's distance.
Η απορία μου είναι λοιπόν η εξής: Η αντίστοιχη εντολή που χρησιμοποιούμε στα απλά lm(.) δεν δουλεύει σε αυτή την περίπτωση. Όσο έψαχνα βρήκα ότι μπορώ να υπολογίσω το Cook's distance αλλά για ΜΟΝΟ μοντέλα με τη χρήση του πακέτου lme4. Υπάρχει τρόπος να βρώ το cook's d για το πακέτο nlme? Ή απλά επικαλούμαι τα profile plots από την exploratory analysis που έκανα αρχικά για να αιτιολογήσω την απόφασή μου να βγάλω τις παρατηρήσεις του ατόμου αυτού από το dataset?
Ευχαριστώ.
- PhD opportunity at UEA: NETWORK META-ANALYSIS OF DIAGNOSTIC ACCURACY STUDIES VIA FACTOR COPULA MIXED MODELS
- AUEB STATS SEMINARS 26/1/2017: Analyzing supersaturated designs for discrete responses via generalized linear models
- AUEB SEMINARS MAY-JUNE 2014 (ΕΝΗΜΕΡΩΜΕΝΟ 21/5/2014)
- NTUA Seminar Talk: CATHERINE HUBER-CAROL - Comparing prediction ability of models and sub-models
- AUEB Stats Seminars 2/11/2022: Non-linear Network Autoregression by K. Fokianos (University of Cyprus)
Permissions in this forum:
You cannot reply to topics in this forum