如何通过从 Excel 导入数据来运行函数 multivari (plotrix)

Question

我对 R 比较陌生，并使用 RStudio 编写了以下脚本来创建 MultiVari 图表。当我想从 Excel 工作表导入数据作为数据集时，我收到一条错误消息（= 数据集 MVData）：

dat_split[[i]][fac1][!is.na(data[[var]]), ]

中的错误： ！无法使用

!is.na(data[[var]])

对行进行子集化。 ✖ 逻辑下标

!is.na(data[[var]])

的大小必须为 1 或 24，而不是 48。

但是，当我使用 data.frame（= 数据集 RRData）构建数据集时，脚本运行没有问题。

library(SixSigma)
library(dplyr)
library(plotrix)
library(readxl)

# Create data set
waarde = c(0.68, -2.09, 1.74, 1.84, -1.7, 0.24, -0.6, -0.23, -0.63, -0.3, 0.75, -0.34, 0.1, 0.32, 0.35, 0.11, 0.31, 0.19, 0.05, -0.01, -0.06, 0.3, 0.33, 0.16, 5.06, 4.18, 5.14, 4.34, 4.29, 4.06, 6.07, 3.7, 6.31, 5.45, 4.72, 5.26, 3.31, 3.47, 3.12, 3.34, 3.27, 2.78, 3.59, 3.57, 3.88, 3.54, 3.67, 3.82)
proces = c('Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'Blank', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA', 'PSA')
type = c('Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Vierkant', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat', 'Plaat')
operator = c('BL', 'BL', 'BL', 'BL', 'BL', 'BL', 'JM', 'JM', 'JM', 'JM', 'JM', 'JM', 'BL', 'BL', 'BL', 'BL', 'BL', 'BL', 'JM', 'JM', 'JM', 'JM', 'JM', 'JM', 'BL', 'BL', 'BL', 'BL', 'BL', 'BL', 'JM', 'JM', 'JM', 'JM', 'JM', 'JM', 'BL', 'BL', 'BL', 'BL', 'BL', 'BL', 'JM', 'JM', 'JM', 'JM', 'JM', 'JM')
run = c('R1', 'R1', 'R1', 'R2', 'R2', 'R2', 'R1', 'R1', 'R1', 'R2', 'R2', 'R2', 'R1', 'R1', 'R1', 'R2', 'R2', 'R2', 'R1', 'R1', 'R1', 'R2', 'R2', 'R2', 'R1', 'R1', 'R1', 'R2', 'R2', 'R2', 'R1', 'R1', 'R1', 'R2', 'R2', 'R2', 'R1', 'R1', 'R1', 'R2', 'R2', 'R2', 'R1', 'R1', 'R1', 'R2', 'R2', 'R2')
ID = c('1', '2', '3', '1', '2', '3', '1', '2', '3', '1', '2', '3', '1', '2', '3', '1', '2', '3', '1', '2', '3', '1', '2', '3', '1', '2', '3', '1', '2', '3', '1', '2', '3', '1', '2', '3', '1', '2', '3', '1', '2', '3', '1', '2', '3', '1', '2', '3')
meting = c('1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1', '1')
RRData = data.frame(waarde,proces, type, operator, run, ID, meting)
print(RRData)

# Import data from excel sheet
MVData <- read_excel("InputMultiVari.xlsx", 
          col_types = c("numeric", "text", "text","text", "text", "text", "text"))
print(MVData)

# Multi Vari chart
multivari("waarde", "ID", "operator", "run", "type",
          data = MVData,
          fac.cex = 1.5,
          main = "Multi-Vari Chart"
)

有人知道如何解决这个问题并使用 Excel 文件作为导入吗？

我尝试更改不同列的变量类型，但这并没有改变结果。

谢谢

结果

dput(MVData)

> dput(MVData)
structure(list(waarde = c(0.68, -2.09, 1.74, 1.84, -1.7, 0.24, 
-0.6, -0.23, -0.63, -0.3, 0.75, -0.34, 0.1, 0.32, 0.35, 0.11, 
0.31, 0.19, 0.05, -0.01, -0.06, 0.3, 0.33, 0.16, 5.06, 4.18, 
5.14, 4.34, 4.29, 4.06, 6.07, 3.7, 6.31, 5.45, 4.72, 5.26, 3.31, 
3.47, 3.12, 3.34, 3.27, 2.78, 3.59, 3.57, 3.88, 3.54, 3.67, 3.82
), proces = c("Blank", "Blank", "Blank", "Blank", "Blank", "Blank", 
"Blank", "Blank", "Blank", "Blank", "Blank", "Blank", "Blank", 
"Blank", "Blank", "Blank", "Blank", "Blank", "Blank", "Blank", 
"Blank", "Blank", "Blank", "Blank", "PSA", "PSA", "PSA", "PSA", 
"PSA", "PSA", "PSA", "PSA", "PSA", "PSA", "PSA", "PSA", "PSA", 
"PSA", "PSA", "PSA", "PSA", "PSA", "PSA", "PSA", "PSA", "PSA", 
"PSA", "PSA"), type = c("Vierkant", "Vierkant", "Vierkant", "Vierkant", 
"Vierkant", "Vierkant", "Vierkant", "Vierkant", "Vierkant", "Vierkant", 
"Vierkant", "Vierkant", "Plaat", "Plaat", "Plaat", "Plaat", "Plaat", 
"Plaat", "Plaat", "Plaat", "Plaat", "Plaat", "Plaat", "Plaat", 
"Vierkant", "Vierkant", "Vierkant", "Vierkant", "Vierkant", "Vierkant", 
"Vierkant", "Vierkant", "Vierkant", "Vierkant", "Vierkant", "Vierkant", 
"Plaat", "Plaat", "Plaat", "Plaat", "Plaat", "Plaat", "Plaat", 
"Plaat", "Plaat", "Plaat", "Plaat", "Plaat"), operator = c("BL", 
"BL", "BL", "BL", "BL", "BL", "JM", "JM", "JM", "JM", "JM", "JM", 
"BL", "BL", "BL", "BL", "BL", "BL", "JM", "JM", "JM", "JM", "JM", 
"JM", "BL", "BL", "BL", "BL", "BL", "BL", "JM", "JM", "JM", "JM", 
"JM", "JM", "BL", "BL", "BL", "BL", "BL", "BL", "JM", "JM", "JM", 
"JM", "JM", "JM"), run = c("R1", "R1", "R1", "R2", "R2", "R2", 
"R1", "R1", "R1", "R2", "R2", "R2", "R1", "R1", "R1", "R2", "R2", 
"R2", "R1", "R1", "R1", "R2", "R2", "R2", "R1", "R1", "R1", "R2", 
"R2", "R2", "R1", "R1", "R1", "R2", "R2", "R2", "R1", "R1", "R1", 
"R2", "R2", "R2", "R1", "R1", "R1", "R2", "R2", "R2"), ID = c("1", 
"2", "3", "1", "2", "3", "1", "2", "3", "1", "2", "3", "1", "2", 
"3", "1", "2", "3", "1", "2", "3", "1", "2", "3", "1", "2", "3", 
"1", "2", "3", "1", "2", "3", "1", "2", "3", "1", "2", "3", "1", 
"2", "3", "1", "2", "3", "1", "2", "3"), meting = c("1", "1", 
"1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
"1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
"1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
"1", "1", "1", "1", "1", "1", "1")), class = c("tbl_df", "tbl", 
"data.frame"), row.names = c(NA, -48L))

结果

identical(MVData, RRData)

> identical(MVData, RRData)
[1] FALSE

Answer 1

我确实相信这是因为 RR- 和 MVData 的类不一致。

[...]

class(RRData)
#> [1] "data.frame"
class(MVData)
#> [1] "tbl_df"     "tbl"        "data.frame"

MVData <- MVData |> 
  as.data.frame()

multivari("waarde", "ID", "operator", "run", "type",
          data = MVData,
          fac.cex = 1.5,
          main = "Multi-Vari Chart"
)

^{创建于 2024 年 10 月 26 日，使用 reprex v2.1.0}