質問: Rを使用して生存率解析(log-rank検定)を行っています。 その際に出てくるp値の有効数字が1桁になってしまっている様で、それを3桁程度に変更したいと考えています。 書いたcodeは以下の様になっており、設定はdefaultからいじっていません。参考にしたページと同じデータを用いて実験したのですが、参考ページではp=0.033と出ていたところこちらではp=0.03と有効数字1桁に丸められてしまっています。 ``` survdiff(Surv(time,event) ~ cancer, rho=0, data=Datafile) # cancerの有無でlog-rank検定をしてます # 大元の"Datafile"にtime(観察期間),event(1:死亡、0:生存/lost)が格納されています ``` 調べたところformmat.pvalというoptionで変更できるのかと思ったのですが、こちらもlog-rankをする際に適用させることがうまくできなさそうです。。。 お分かりになる方がいらっしゃいましたらご教授のほどよろしくお願い申し上げます。
回答: result$chisqでカイ二乗値を取り出して、pchisqでp-valueに変換 ``` result <- survdiff(Surv(time,event) ~ cancer, rho=0, data=Datafile) pchisq(result$chisq, 1, lower.tail=FALSE) #自由度1 ```