r - R中单个统计量的Wald检验

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.

标签： rstatistics

Answer 1

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