2

我已经在 R 上尝试了一段时间，但我似乎无法提取 p 值（“Pr(>f) 的值）”用于 levene 测试。对 R 进行统计测试的常用方法是在测试命令末尾加上 $p.value。然而，这似乎不适用于 Levene 测试，如下所示：

> leveneTest(all.vec,factors)
Levene's Test for Homogeneity of Variance (center = median)
        Df F value    Pr(>F)    
group    3  8.9261 6.982e-06 ***
      2607                      
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> leveneTest(all.vec,factors)$p.value
NULL

此外，对于 jarque bera 测试，我无法以同样的方式提取 p.value

> jarque.bera.test(lg.day.ret.vec)

    Jarque Bera Test

data:  lg.day.ret.vec
X-squared = 63087.83, df = 2, p-value < 2.2e-16

> jarque.bera.test(lg.day.ret.vec)$p.value
X-squared 
        0

谢谢您的帮助

4

3 回答 3

3

有几种方法。这是一个可重现的示例：

library("car")
test <- with(Moore, leveneTest(conformity, fcategory))

首先，查看返回对象的结构，因为它通常会告诉您正在玩什么：

str(test)

这给出了：

> str(test)
Classes ‘anova’ and 'data.frame':   2 obs. of  3 variables:
 $ Df     : int  2 42
 $ F value: num  0.046 NA
 $ Pr(>F) : num  0.955 NA
 - attr(*, "heading")= chr "Levene's Test for Homogeneity of Variance (center = median)"

我们看到对象是一个数据框，p值在第三列。因此，以下任何一项都将提取数据

test[,3]          # pull out the entire 3rd column
test[1,3]         # pull out only the none NA p-value
test$`Pr(>F)`     # pull out the P-value column by name
test$`Pr(>F)`[1]  # as above, but then take only the 1st element

对于上面的示例，这些给出：

> test[,3]
[1] 0.9550975        NA
> test[1,3]
[1] 0.9550975
> test$`Pr(>F)`
[1] 0.9550975        NA
> test$`Pr(>F)`[1]
[1] 0.9550975
于 2015-02-19T02:12:56.560 回答
0

您还可以使用包levene.test中的功能lawstat：

> df<-data.frame(group=c(rep(LETTERS[1],20),rep(LETTERS[2],20)),
               value=c(sample(1:10,size=20,replace=T),sample(30:40,size=20,replace=T)))

> df
   group value
1      A     8
2      A     4
3      A     8
4      A     3
5      A     2
6      A     3
7      A     4
8      A     5
9      A     8
10     A     6
11     A     1
12     A     4
13     A     9
14     A     7
15     A     8
16     A     4
17     A     2
18     A    10
19     A     8
20     A     7
21     B    40
22     B    37
23     B    39
24     B    30
25     B    30
26     B    30
27     B    34
28     B    39
29     B    40
30     B    34
31     B    30
32     B    33
33     B    32
34     B    39
35     B    36
36     B    37
37     B    35
38     B    39
39     B    34
40     B    34

> library(lawstat)
> levene.test(df$value,df$group,location='median')

modified robust Brown-Forsythe Levene-type test based on the absolute deviations from the median

data:  df$value
Test Statistic = 1.7537, p-value = 0.1933

> levene.test(df$value,df$group,location='median')$p.value
[1] 0.1933243
于 2015-02-19T03:24:29.753 回答
0

你也可以写 leeneTest(all.vec,factors)$F value

于 2020-12-18T00:35:17.323 回答