check_model_stanarm: Diagnostiche

In masspastore/ADati: Analisi dei dati in Psicologia

Diagnostiche

Description

Produce i grafici per le diagnostiche di un modello ottenuto con il pacchetto rstanarm.

Usage

check_model_stanarm( fit, all = FALSE, 
      cook_levels = c(.5,1), cex = 3, 
      whichplot = 1:6  ) 

Arguments

fit	oggetto ottenuto con la funzione stan_glm() o stan_glmer().
all	logico, se posto a TRUE (il default) produce tutti i grafici insieme.
cook_levels	vettore numerico, per default c(.5,1), indica le soglie per rappresentare la distanza di Cook. Deve avere al massimo due valori.
cex	valore numerico che indica la grandezza dei punti.
whichplot	vettore numerico tra 1 e 6, che indica quali grafici produrre, per default tutti.

Value

Restituisce sei grafici: 1) residui del modello rispetto ai valori attesi, 2) Quantile-Quantile rispetto alla normale, 3) Scale-Location: radice quadrata dei residui standardizzati in valore assoluto rispetto ai valori attesi, 4) residui del modello rispetto ai leverages con le curve relative alla distanza di Cook, 5) k di Pareto per le osservazioni, 6) Posterior predictive check.

I primi quattro grafici sono gli stessi che si ottengono con il comando plot() a partire da un oggetto ottenuto con la funzione lm(), il quinto grafico, è quello che si ottiene con un oggetto ricavato dalla funzione loo() ed il sesto si ottiene con la funzione pp_check().

Note

La funzione richiede il pacchetto grafico ggplot2.

Author(s)

Massimiliano Pastore

Examples

require(rstanarm)
data(sherifdat)
fit <- stan_glmer(y~time+condition+(1|group),data=sherifdat)
check_model_stanarm(fit)

## per ottenere un grafico unico
require(cowplot)
PLOTS <- check_model_stanarm(fit, TRUE)
plot_grid( plotlist = PLOTS )

masspastore/ADati documentation built on June 12, 2025, 7:28 p.m.