FinTS包没有了怎么做ArchTest呢？或者说怎么做LM检验呢？

复制粘贴这段代码：

1. ArchTest <- function (x, lags=12,
   
2.                 demean = FALSE)
   
3. {
   
4. # Capture name of x for documentation in the output  
   
5.   xName <- deparse(substitute(x))
   
6. #
   
7.   x <- as.vector(x)
   
8.   if(demean) x <- scale(x, center = TRUE, scale = FALSE)
   
9. #  
   
10.   lags <- lags + 1
    
11.   mat <- stats::embed(x^2, lags)
    
12.   arch.lm <- summary(stats::lm(
    
13.         mat[, 1] ~ mat[, -1]))
    
14.   STATISTIC <- arch.lm$r.squared *
    
15.     length(stats::resid(arch.lm))
    
16.   names(STATISTIC) <- "Chi-squared"
    
17.   PARAMETER <- lags - 1
    
18.   names(PARAMETER) <- "df"
    
19.   PVAL <- stats::pchisq(STATISTIC,
    
20.         df = PARAMETER, lower.tail=FALSE)
    
21.   METHOD <- paste("ARCH LM-test; ",
    
22.       "Null hypothesis:  no ARCH effects")
    
23.   result <- list(statistic = STATISTIC,
    
24.         parameter = PARAMETER,
    
25.         p.value = PVAL, method = METHOD,
    
26.         data.name = xName)
    
27.   class(result) <- "htest"
    
28.   return(result)
    
29. }

复制代码