美文网首页
【R>>scales】pvalues

【R>>scales】pvalues

作者: 高大石头 | 来源:发表于2021-04-21 07:19 被阅读0次

在实际科研中,常常会出现pvalue = 0.000000000123类似的形式,甚至出现p = 0,遇到如此尴尬的情况怎么办呢?
其实scales::pvalue()函数已经帮咱们设计好了,来看具体代码:


> lm(hwy~displ+drv, data=mpg) %>% 
+   tidy()
# A tibble: 4 x 5
  term        estimate std.error statistic  p.value
  <chr>          <dbl>     <dbl>     <dbl>    <dbl>
1 (Intercept)    30.8      0.924     33.4  4.21e-90
2 displ          -2.91     0.218    -13.4  1.73e-30
3 drvf            4.79     0.530      9.05 6.40e-17
4 drvr            5.26     0.734      7.17 1.03e-11
> lm(hwy~displ+drv, data=mpg) %>% 
+   tidy() %>% 
+   mutate(p.value=pvalue(p.value))
# A tibble: 4 x 5
  term        estimate std.error statistic p.value
  <chr>          <dbl>     <dbl>     <dbl> <chr>  
1 (Intercept)    30.8      0.924     33.4  <0.001 
2 displ          -2.91     0.218    -13.4  <0.001 
3 drvf            4.79     0.530      9.05 <0.001 
4 drvr            5.26     0.734      7.17 <0.001 

聪明的你,会发现,pvalue直接被处理成了能直接上图的格式。
参考链接:
Andrew Heiss的官方推特

相关文章

网友评论

      本文标题:【R>>scales】pvalues

      本文链接:https://www.haomeiwen.com/subject/rvholltx.html