From c737796abcad1bacd649edb46cc07311c3c4f162 Mon Sep 17 00:00:00 2001
From: Arthur DANJOU <arthurdanjou@outlook.fr>
Date: Fri, 13 Feb 2026 16:28:12 +0100
Subject: [PATCH] Add tseries, separate trend/seasonal/final differencing plots
 and run KPSS/ADF stationarity tests

---
 M2/Time Series/TD2.Rmd | 22 ++++++++++++++++++++--
 1 file changed, 20 insertions(+), 2 deletions(-)

diff --git a/M2/Time Series/TD2.Rmd b/M2/Time Series/TD2.Rmd
index 0bc5e3a..8ffad5d 100644
--- a/M2/Time Series/TD2.Rmd	
+++ b/M2/Time Series/TD2.Rmd	
@@ -3,6 +3,9 @@
 ## 1. Import et visualisation
 
 ```{r}
+library(tseries)
+par(mfrow = c(1, 1))
+
 data(AirPassengers)
 plot(AirPassengers)
 ```
@@ -83,9 +86,24 @@ Différences clés :
 ## 6. Différenciation
 
 ```{r}
-plot(diff(diff(log(AirPassengers), lag = 12)))
+trend_diff <- diff(log(AirPassengers))
+saison_diff <- diff(log(AirPassengers), lag = 12)
+final_diff <- diff(trend_diff, lag = 12)
+
+plot(trend_diff, main = "Différence première (Trend)", ylab = "Différence")
+plot(saison_diff, main = "Différence saisonnière (Seasonal)", ylab = "Différence saisonnière")
+plot(final_diff, main = "Différence finale (Trend + Seasonal)", ylab = "Différence finale")
 ```
 
 Une différence première ($I−B$) pour retirer la tendance.
 
-Une différence saisonnière ($I−B^{12}$) pour retirer la saisonnalité.
\ No newline at end of file
+Une différence saisonnière ($I−B^{12}$) pour retirer la saisonnalité.
+
+
+```{r}
+kpss_test_result <- kpss.test(final_diff)
+print(kpss_test_result)
+    
+adf_test_result <- adf.test(final_diff)
+print(adf_test_result)
+```
\ No newline at end of file