积分。

Question

我正在使用生存软件包运行survdiff，p值为0.02。我希望看到它具有更高的精度（即0.02xxxx）。我可以传递一个参数来指定pvalue的长度吗？我阅读了生存程序包的文档，但没有找到有关如何指定生存程序的任何信息。

survdiff(surv_object~access_sam2$Area_mTLSHL)

Answer 1

积分。

const a;类对象的p值的计算并不完全明显。我必须查看该类对象的"survdiff"方法中发生了什么，才能了解计算自由度的方式。

下面的代码是print的简化代码，因此功劳归功于

print.survdiff

代码本身可以在源代码中查看或通过运行

citation("survival")
#
#Therneau T (2015). _A Package for Survival Analysis
#in S_. version 2.38, <URL:
#https://CRAN.R-project.org/package=survival>.
#
#Terry M. Therneau, Patricia M. Grambsch (2000).
#_Modeling Survival Data: Extending the Cox Model_.
#Springer, New York. ISBN 0-387-98784-3.
#
#To see these entries in BibTeX format, use
#'print(<citation>, bibtex=TRUE)', 'toBibtex(.)', or
#set 'options(citation.bibtex.max=999)'.

现在是问题的问题。

我已经编写了一个通用的getAnywhere("print.survdiff")函数，以便更轻松地为函数pvalue返回的类的对象调用方法。该示例摘自该功能的帮助页面。
返回值是一个包含3个成员的命名列表，这些名称不言自明。其中之一，survdiff是chisq返回的值的重复。为了完整起见，我将其包括在内。

survdiff

Answer 2

我不确定生存包，您没有提供可复制的代码（请下次提供）。但总的来说，如果您想查看更多数字，您需要做的是

print(value, digits= n)
# n is the number of digits you want to see

您的情况是

 print(survdiff(surv_object~access_sam2$Area_mTLSHL), 6)

在Survdiff中获得更高的pvalue精度

2 个答案:

积分。

现在是问题的问题。