[問題] 如何計算p value?

看板Statistics作者 (阿水~*)時間14年前 (2009/09/03 22:22), 編輯推噓1(1推 0噓 27→)
留言28則, 7人參與, 5年前最新討論串1/1
這個網頁 http://www.danielsoper.com/statcalc/calc11.aspx 能夠用 Chi-Square Value Degrees of Freedom 去計算p value 請問有沒有人知道計算的公式 我查了好久.... 蠻急的,感謝@@! 我接到一個案子要寫這個功能,讀檔案大量計算,所以必須知道公式並實現 不能用excel之類的軟體去算 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 219.70.226.11
→
09/04 00:40, , 1F
網路查不到？
09/04 00:40, 1F
推
09/04 01:03, , 2F
wikipedia..
09/04 01:03, 2F
→
09/04 03:03, , 3F
你的連結看起來很像色情網頁...
09/04 03:03, 3F
→
09/04 07:03, , 4F
就是卡方分配的截尾機率 對該分配做定積分阿 用數值方法
09/04 07:03, 4F
→
09/04 07:06, , 5F
是說你說的案件可能搜一搜就有code了...
09/04 07:06, 5F
→
09/04 07:10, , 7F
卡方是長上面那樣 你有自由度(k) 樣本值(x)
09/04 07:10, 7F
→
09/04 07:11, , 8F
然後算 P(χ^2(k)>x) 就是p-value了...
09/04 07:11, 8F
→
09/04 10:34, , 9F
我有查 維基 google 知識+ IEEE SCOPUS ACM
09/04 10:34, 9F
→
09/04 10:34, , 10F
雖然直覺上看起來很簡單,好像找一下就有,但實際做才知道
09/04 10:34, 10F
→
09/04 10:34, , 11F
09/04 10:34, 11F
→
09/04 10:35, , 12F
上面這網頁就有公式 另外gamma function的部份也研究了
09/04 10:35, 12F
→
09/04 10:36, , 13F
我的df一定會是偶數,所以gamma function的輸入一定整數
09/04 10:36, 13F
→
09/04 10:37, , 14F
可以直接用階乘實現 但是! 就是都照公式做了 值卻不同
09/04 10:37, 14F
→
09/04 10:37, , 15F
我不清楚問題再哪裡 所以發問
09/04 10:37, 15F
→
09/04 10:37, , 16F
我很討厭被懷疑沒做功課 >"<
09/04 10:37, 16F
→
09/04 12:37, , 17F
抱歉喔！你上面的寫法意思是查不到，可不是說有查到做
09/04 12:37, 17F
→
09/04 12:38, , 18F
不出來.....我猜測你也許是gamma function弄錯
09/04 12:38, 18F
→
09/04 12:40, , 19F
Γ(n)=(n-1)! 當n是整數時
09/04 12:40, 19F
→
09/04 15:03, , 20F
謝謝你 有板友丟水球給我 真溫暖
09/04 15:03, 20F
→
09/04 15:05, , 21F
我處理的陣列row一定是3 col=2~20
09/04 15:05, 21F
→
09/04 15:08, , 22F
所以自由度一定是偶數　因此gamma(k/2)的輸入 一定整數
09/04 15:08, 22F
→
09/04 15:11, , 23F
網友解答我的疑惑了 我的公式是pdf 我需要的是cdf
09/04 15:11, 23F
→
09/04 15:39, , 24F
pgamma() in C, 在R的nmath有原始碼.
09/04 15:39, 24F
→
09/04 15:48, , 25F
直覺上, 並不簡單.......
09/04 15:48, 25F
→
09/04 15:50, , 26F
如果量不太大,用R就行了.......
09/04 15:50, 26F
→
09/05 10:08, , 27F
問題解決了! 感謝大家~~~YA
09/05 10:08, 27F
→
01/02 14:58, 5年前 , 28F
上面這網頁就有公式 另 https://muxiv.com
01/02 14:58, 28F
‣ 更多 pizza0117 的文章
文章代碼(AID): #1Adz51XS (Statistics)