SPI指数中的Shapiro-Wilk检验计算

Question

我正在尝试通过 Shapiro Wilk 测试验证 SPI 指数是否服从正态分布，但我不知道我是否需要使用 $fitted 和来自 Shapiro 测试的 $p.value 或 $w 或者可能需要转换时间数据框的系列格式，如下所示：

install.packages("SPEI")  
install.packages("readr")
library(SPEI)
library(readr)

options(scipen=999)

df <- read_delim("https://raw.githubusercontent.com/victoresende19/testing/main/PRECPT.csv", 
                           ";", escape_double = FALSE, locale = locale(decimal_mark = ","), 
                           trim_ws = TRUE)


##### SPI #####
SPI12 <- spi(ts(df$`PRECIPITACAO`, freq = 12, start=c(1992, 1)), scale = 12)
plot(SPI12)

#Teste de normalidade
SPI_shapiro <- shapiro.test(SPI12$fitted)$p.value

SPI 指数必须返回服从正态分布的指数。但我认为我做错了什么，因为 shapiro 测试返回 p 值 0.0019，拒绝 H0。