【问题标题】:Wald-test for single statistic in RR中单个统计量的Wald检验
【发布时间】:2019-11-15 09:39:19
【问题描述】:

我在曲线中的两个点(低点和高点)有一系列风险率,并带有相应的标准误差。我通过将高点风险率除以低点风险率来计算风险比。这是比率列。现在在下一列中,我想使用 Wald 检验显示该比率与 1 显着不同的概率(p 值)。

我尝试使用 aods3 包中的 wald.test() 执行此操作,但我不断收到错误消息。似乎代码只允许比较两个相关的回归模型。

你会怎么做呢?

> wald
     fit.low      se.low   fit.high     se.high    hratio
1 0.09387638 0.002597817 0.09530283 0.002800329 0.9850324
2 0.10941588 0.002870383 0.10831292 0.003061924 1.0101831
3 0.02549611 0.001054303 0.02857411 0.001368525 0.8922802
4 0.02818208 0.000917136 0.02871669 0.000936373 0.9813833
5 0.04857652 0.000554676 0.04897211 0.000568229 0.9919222
6 0.05121328 0.000565592 0.05142951 0.000554893 0.9957956
> library(aods3)
> wald$pv <- wald.test(b=wald$hratio)
Error in wald.test(b = wald$hratio) : 
  One of the arguments Terms or L must be used.

【问题讨论】:

    标签: r statistics


    【解决方案1】:

    定义L=NULL, Terms=NULL, Sigma = vcov(b)

    【讨论】:

    • 我收到同样的错误信息:wald$pv &lt;- wald.test(b=wald$hratio, L=NULL, Terms=NULL, Sigma = vcov(b)) Error in wald.test(b = wald$hratio, L = NULL, Terms = NULL, Sigma = vcov(b)) : One of the arguments Terms or L must be used.
    猜你喜欢
    • 1970-01-01
    • 1970-01-01
    • 2020-06-06
    • 1970-01-01
    • 2015-09-16
    • 2012-11-28
    • 1970-01-01
    • 1970-01-01
    • 2021-06-18
    相关资源
    最近更新 更多