作者w2a3n4g5 (MORRIS.W)
看板R_Language
標題[問題] 如何設定pchisq輸出的最小值
時間Sun Oct 11 16:41:00 2015
[問題敘述]:
如題!
小弟我今天在寫適合度檢定,
但檢定統計量的值非常大,我將此值
寫成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
→ celestialgod: R最多就是15~16位,要更多要用別的軟體,像是MAPLE 10/11 17:11
推 ching0629: pchisq(800,df=3,lower.tail=FALSE) 10/11 22:12
→ ching0629: 另外,如果真的很需要詳細數字,匯出log後的結果再計算 10/11 22:13
→ ching0629: pchisq(800,df=3,lower.tail=FALSE,log.p=TRUE) 10/11 22:13
→ ching0629: 至於為什麼你使用的方法無法顯示呢? 這是因為R會先匯 10/11 22:14
→ ching0629: 出9.999999*10^-1,而前面的部份位數有限 10/11 22:15