Hotelling 多元控制图

Question

在我的工作中，我们做了很多 control charts for n variables for monitoring multivariate measurement processes. I am trying to implement Hotelling multivariate control charts so we can consider the correlation between variables and assess when a sample has gone out-of-control. The package MSQC 具有函数 mult.chart 来轻松实现这一点。函数 return 提供有关数据的一些信息，控制图的绘图包括限制，此外，当样本失控时，它会对其进行分解，以便确定哪个变量（或多个变量）是轮班的负责人，信息以列表的形式给出。我无法确定如何提取该信息，因为 return 的结构不遵循 name_of_data_frame$name_of_variable_of_interest 的传统形式，至少对于分解矩阵而言是这样。

library(MSQC)

data("carbon1")
Xmv <- mult.chart(carbon1, type = "t2") $Xmv
S <- mult.chart(carbon1, type = "t2") $covariance
colm<-nrow(carbon1)
#Phase II
data("carbon2")
Hot<-mult.chart(carbon2, type = "t2", Xmv = Xmv, S = S, colm = colm)

The following(s) point(s) fall outside of the control limits[1] 4

$`Decomposition of`
[1] 4

我试过str(Hot)但是$`Decomposition of`的部分没有出现。我怎样才能得到这类信息？

Answer 1

mult.chart 的输出中不包含 ATM 测试统计量 (T2) 的分解。我会在下一次更新中包含它。 解决它的一种方法是使用：

library(MSQC)
data("carbon1") # dataset used in Phase I
Xmv <- mult.chart(carbon1, type = "t2") $Xmv
S <- mult.chart(carbon1, type = "t2") $covariance
colm<-nrow(carbon1)

#Phase II
data("carbon2")

co <- capture.output(Hot<-mult.chart(carbon2, type = "t2", Xmv = Xmv, S = S, colm = colm))
write(co[5:length(co)], "C:\1\decomp.txt") # saving it 
df <- read.table("C:\1\decomp.txt", sep = "", header = F)

希望对您有所帮助。