Rの科学表記の代わりに4桁のみを印刷するp値を行うにはどうすればよいですか？

Question

Rの科学表記の代わりに4桁のみを印刷するp値を行うにはどうすればよいですか？あなたがここに欲しいものをここでは一例です...

>options(digits=3, scipen=12) 
>Oi <- c(A=321, B=712, C=44) 
>Ei <- c(A=203, B=28, C=6) 
>chisq.test(Oi, p=Ei,rescale.p=T)" 

Answer 1

かなりわからない ... 私はオプション（数字= 3、scipen = 12）を使用しようとしたが、それはうまくいきませんでした。再現性の例えば ありがとう：出力は、オブジェクトの構造を検査

cc <- chisq.test(Oi,p=Ei,rescale.p=TRUE) 

print(cc) 
    Chi-squared test for given probabilities 

data: Oi 
X-squared = 3090, df = 2, p-value < 0.00000000000000022 

が、この場合のp値が正確にゼロにアンダーフローしていることを明らかにした

List of 9 
$ statistic: Named num 3090 
    ..- attr(*, "names")= chr "X-squared" 
$ parameter: Named num 2 
    ..- attr(*, "names")= chr "df" 
$ p.value : num 0 
$ method : chr "Chi-squared test for given probabilities" 
$ data.name: chr "Oi" 
$ observed : Named num [1:3] 321 712 44 
    ..- attr(*, "names")= chr [1:3] "A" "B" "C" 
$ expected : Named num [1:3] 922.5 127.2 27.3 
    ..- attr(*, "names")= chr [1:3] "A" "B" "C" 
$ residuals: Named num [1:3] -19.8 51.8 3.2 
    ..- attr(*, "names")= chr [1:3] "A" "B" "C" 
$ stdres : Named num [1:3] -52.29 55.2 3.25 
    ..- attr(*, "names")= chr [1:3] "A" "B" "C" 
- attr(*, "class")= chr "htest" 

したい場合は、私が思いますあなたはあなたの方法のうちのビットを行かなければならない、このテストからのp値の正確な：

(pval <- pchisq(3090,2,lower.tail=FALSE,log.p=TRUE)) 
[1] -1545 

だから、これは約10^pval/log(10) = 10 ^である（ - 671）[R .Machine$double.xminを参照してください

Answer 2

あなたが望むものはわかりませんが、< 0.0001のようなものをSASの出力と同じように探していると思います。私はそれにはformat.pval関数を使うでしょう。あるいは、あなたが持っているテストの数によってはprintCoefmatかもしれません。 epsは許容誤差です。その下の値は< [eps]として出力されます。

Oi <- c(A=321, B=712, C=44) 
Ei <- c(A=203, B=28, C=6) 
tt <- chisq.test(Oi, p=Ei,rescale.p=T) 

format.pval(tt$p.value, eps=0.0001) 
# [1] "<0.0001" 

ttp <- data.frame(Chisq=tt$statistic, p.value=tt$p.value) 
rownames(ttp) <- "Oi vs Ei" 
printCoefmat(ttp, has.Pvalue=TRUE, eps.Pvalue=0.0001)  
#   Chisq p.value  
# Oi vs Ei 3090 <0.0001 *** 
# --- 
# Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1