R Delta方法与聚集标准误差

我有一个问题，当我有聚集的标准错误时，如何应用delta方法。考虑下面的数据集和（简单）回归（（请注意，这个问题不一定是围绕着“我们”的聚类或这个回归的正确性/有用性）是否有意义。

#使用软件包
图书馆（多路VCOV）
图书馆（三明治）
图书馆（lmtest）
图书馆（msm）
#加载数据
数据（mtcars）
#运行回归
model1只是编辑
deltamethod
调用以输出一个答案-我不知道这个答案是否对您想要做的事情有意义

deltamethod(
  g = formula('~x2/x1'), 
  mean = model1$coefficients, 
  cov = vcov, 
  ses = TRUE)

你是如何应用这个方法的？你得到了什么错误？如果没有一个可重复的例子，很难看出你做错了什么。你好，奥利弗，很抱歉我的问题不完整。希望这有帮助。感谢你研究它。看看帮助，
？deltamethod
。第一个参数应该是一个公式，你提供了两个比率数字。你好，伪自旋，你是对的。我在我的R中正确地运行了它，但在我的问题中键入了错误。当我插入“g”而不是分数时，错误保持不变。所以我仍然不确定我应该在这里更改什么。你知道吗？
deltamethod(
  g = formula('~x2/x1'), 
  mean = model1$coefficients, 
  cov = vcov, 
  ses = TRUE)