[英]R: one missing value ruins my Jarque-Bera test
此代碼不省略 N/As 並正確執行 Jarque-Bera 測試。 什么代碼呢? 謝謝你的建議。
library(moments)
jarque.test(mydata$item1, na.rm = TRUE)
我的資料
structure(list(item1 = c(6, 7, NA, 7, 6, 2, 7, 3, 4, 7, 4, 5,
4, 7, 7, 6, 5, 7, 6, 7, 4, 7, 5, 6, 5, 4, 7, 5, 4, 6, 7, 5, 5,
7, 7, 5, 7, 7, 7, 4, 5, 7, 7, 7, 5, 7, 6, 7, 7, 5), item2 = c(6,
7, 6, 6, 6, 3, 7, 3, 3, 4, 5, 6, 4, 7, 6, 6, 4, 6, 6, 7, 6, 3,
5, 5, 3, 2, 7, 5, 6, 6, 7, 3, 5, 7, 6, 5, 6, 5, 6, 4, 6, 7, 7,
7, 7, 7, 6, 7, 4, 7)), row.names = c(NA, -50L), class = c("tbl_df",
"tbl", "data.frame"))
如評論 state 所示,由於moments::jarque.test
中沒有na.rm
參數,並且它只接受一個參數,因此解決方法是僅對數據進行子集化,不包括NA
s。
moments::jarque.test(mydata$item1[!is.na(mydata$item1)])
# Jarque-Bera Normality Test
#
# data: mydata$item1[!is.na(mydata$item1)]
# JB = 4.9442, p-value = 0.08441
# alternative hypothesis: greater
我們可以使用JarqueBeraTest
形式的DescTools
,與其他的相比,它有一個na.rm
參數:
library(DescTools)
JarqueBeraTest(df$item1, robust = TRUE, na.rm = TRUE)
data: structure(c(6, 7, 7, 6, 2, 7, 3, 4, 7, 4, 5, 4, 7, 7, 6, 5, 7, 6, 7, 4, 7, 5, 6, 5, 4, 7, 5, 4, 6, 7, 5, 5, 7, 7, 5, 7, 7, 7, 4, 5, 7, 7, 7, 5, 7, 6, 7, 7, 5), na.action = structure(3L, class = "omit"))
X-squared = 3.9623, df = 2, p-value = 0.1379
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.