hat_values_stanarm: Hat values

In masspastore/ADati: Analisi dei dati in Psicologia

Hat values

Description

Calcola gli hat values di un modello ottenuto con il pacchetto rstanarm.

Usage

hat_values_stanarm( fit ) 

Arguments

fit	oggetto ottenuto con la funzione stan_glm() o stan_glmer().

Details

Per i modelli lineari semplici i valori sono gli stessi che si ottengono con la funzione hatvalues su un modello di classe lm; nel caso di modelli con effetti random, al posto degli hat values sono calcolati i k di Pareto ricavati da un oggetto ottenuto con la funzione loo().

Value

Restituisce un vettore con gli hat values (per i modelli lineari semplici) o i k di Pareto (per i modelli multilivello) per ciascuna osservazione.

Author(s)

Massimiliano Pastore

Examples

require(rstanarm)
data(sherifdat)

# modello lineare semplice
fit_lm <- lm(y~time+condition,data=sherifdat)
hat_lm <- hatvalues( fit_lm )

fit_stan <- stan_glm(y~time+condition,data=sherifdat)
hat_stan <- hat_values_stanarm( fit_stan )

all.equal(hat_lm,hat_stan) # verifica

# ---------------------------
# modello multilivello 
fit <- stan_glmer(y~time+condition+(1|group),data=sherifdat)
HV <- hat_values_stanarm(fit)
plot( HV, ylab = "hat values", xlab = "data" )

# rappresentazione alternativa
PLOTS <- check_model_stanarm(fit,all=TRUE)
PLOTS[[5]]

masspastore/ADati documentation built on Jan. 31, 2025, 4:49 p.m.