1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
local({
## Compute
vars <- rk.list (test50x, test10y)
results <- data.frame ("Variable Name"=I(names (vars)), check.names=FALSE)
for (i in 1:length (vars)) {
var <- vars[[i]]
results[i, "Length"] <- length (var)
results[i, 'NAs'] <- sum (is.na(var))
try ({
test <- Box.test (var, lag = 1, type = "Box-Pierce")
results[i, "X-squared"] <- test$statistic
results[i, "degrees of freedom"] <- test$parameter
results[i, "p-value"] <- test$p.value
})
}
## Print result
rk.header ("Box-Pierce and Ljung-Box Tests", parameters=list("lag"="1",
"test to be performed"="Box-Pierce"))
rk.results (results)
})
|
