R软件结果P值的显示问题。

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

大虾们，我想问下额。我运行一个R程序，结果里面的P值，精度只到0.0099。小于这个值，P值都显示为<0.0099，
因为要做多重检验的修正，所以我想知道具体的P值大小，R里面有什么办法可以让这个P值不要显示小于多少多少，而是显示具体的值吗？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：r软件 R程序 软件

一般可以用"objectname$p.value"的格式提取出p-value, 这时的p-value应该是精确的。确切起见，你可以把你具体的程序贴上来看看。

hugebear 发表于 2015-4-21 21:35
一般可以用"objectname$p.value"的格式提取出p-value, 这时的p-value应该是精确的。确切起见，你可以把你具 ...

大虾。好像不行欸。。我打“print(aa$result$Perm.P)”。。。。出来的是NULL。
原句如下：

> library(mbmdr)
Loading required package: logistf
Loading required package: mice
Loading required package: Rcpp
Loading required package: lattice
mice 2.22 2014-06-10
Loading required package: mgcv
Loading required package: nlme
This is mgcv 1.8-4. For overview type 'help("mgcv-package")'.
mbmdr 2.6
>
>
> fan<-read.csv("ldlc-data.csv",header=TRUE)
> attach(fan)
> pp <- mbmdr(y=ldlcj_class414, data=fan[1:10], order=3, covar=fan[11:13], exclude=NA, risk.threshold=0.1,
+ output=NULL, adjust="covariates",family=binomial(link=logit),use.logistf=TRUE,
+ first.model=NULL, list.models=NULL,printStep1=FALSE)
>
>
> #Permutation test for all models with MIN.P <= 0.05
> order <- 3
> models <- subset(pp$result, select = 1:order)
> aa<-mbmdr.PermTest(pp,100,models)
> print(aa)