作者nearfly (细枝小叶)
看板NCU97Stat
标题[转录][统算] chisq.test
时间Sat Jun 6 16:24:45 2009
※ [本文转录自 NCU96Stat 看板]
作者: Chris7462 (～烤焦面包～) 看板: NCU96Stat
标题: Re: [统算] chisq.gof()
时间: Wed Jun 4 00:48:56 2008
# chisq.test.r
- rescale.p = FALSE, simulate.p.value = FALSE, B = 2000)
+ rescale.p = FALSE, simulate.p.value = FALSE, B = 2000, unspec.m = 0)
else {
- PARAMETER <- length(x) - 1
+ if(unspec.m){
+ PARAMETER <- length(x) - unspec.m - 1
+ }
+ else {
+ PARAMETER <- length(x) - 1
+ }
PVAL <- pchisq(STATISTIC, PARAMETER, lower.tail = FALSE)
}
## Example 9c ##
> N <- c(6, 2, 1, 9, 7, 5)
> prob <- c(0.0550, 0.1596, 0.2312, 0.2237, 0.1622, 0.1682)
> chisq.test(N, p=prob, unspec.m=1)
Chi-squared test for given probabilities
data: N
X-squared = 19.887, df = 4, p-value = 0.0005257
Warning message:
In chisq.test(N, p = p, unspec.m = 1) :
Chi-squared approximation may be incorrect
--
改一下 chisq.test 的函式，加入一个 unspec.m（unspecified parameter 的个数）
这样跑出来的结果就跟 Example 9c 一样了，也就是我们要的 goodness of fit :p
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.115.221.212
※ 编辑: Chris7462 来自: 140.115.221.212 (06/04 01:25)
1F：推 anniecs:好厉害 06/04 12:47
2F：推 yayacy:真的超强的!!!感恩感恩 06/05 01:25
3F：推 luckyb:感激Q__Q 06/05 14:58
--
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 118.161.197.102
4F：→ nearfly:学长3q~ 06/06 16:34
5F：推 Chris7462:不客气。突然看到自己 po 的文还满奇怪的 ＝ ＝ 06/06 21:05
6F：推 ricindy:太感谢了~!! 06/07 23:17