[問題] 如何設定pchisq輸出的最小值

看板R_Language作者 (MORRIS.W)時間8年前 (2015/10/11 16:41), 編輯推噓1(1推 0噓 6→)
留言7則, 2人參與, 最新討論串1/1
[問題敘述]: 如題！ 小弟我今天在寫適合度檢定， 但檢定統計量的值非常大，我將此值 寫成pvalue<-1-pchisq(800,df=3) 但我得出來的直永遠都是0 我測試pchisq這個函數後發現他最多只會顯示到1E-16 再小就出不來了，用format或options都無法！ 所以想知道怎樣可以讓它顯示出來 比如說顯示pchisq(100,3) 出來就是零！可是我想知道他到底是多小就算是1E-100也好！ 麻煩各位了 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.117.247.2 ※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1444552863.A.DFE.html
→
10/11 17:11, , 1F
10/11 17:11, 1F
→
10/11 17:11, , 2F
R最多就是15~16位，要更多要用別的軟體，像是MAPLE
10/11 17:11, 2F
推
10/11 22:12, , 3F
pchisq(800,df=3,lower.tail=FALSE)
10/11 22:12, 3F
→
10/11 22:13, , 4F
另外，如果真的很需要詳細數字，匯出log後的結果再計算
10/11 22:13, 4F
→
10/11 22:13, , 5F
pchisq(800,df=3,lower.tail=FALSE,log.p=TRUE)
10/11 22:13, 5F
→
10/11 22:14, , 6F
至於為什麼你使用的方法無法顯示呢? 這是因為R會先匯
10/11 22:14, 6F
→
10/11 22:15, , 7F
出9.999999*10^-1,而前面的部份位數有限
10/11 22:15, 7F
‣ 更多 w2a3n4g5 的文章
文章代碼(AID): #1M6Y2Vt- (R_Language)