Question

我在曲线上的两个点（低点和高点）具有一系列危险率，并带有相应的标准误差。我通过将高点危险率除以低点危险率来计算危险比。这是小时数列。现在，在下一列中，我想显示通过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.

Answer 1

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

R中的单一统计量的Wald检验

1 个答案: