Wald测试两个分数微分参数是否相等

Question

我有两个ARFIMA估计，如下所示：

short_rate = ts(yields$`1M`)
y_s = short_rate - mean(short_rate)
short_rate_FD = arfima(y_s)
summary(short_rate_FD)
d_r = summary(short_rate_FD)$coef[[1]][1]

intermediate_rate = ts(yields$`2Y`)
y_m = intermediate_rate - mean(intermediate_rate)
summary(intermediate_rate_FD)
intermediate_rate_FD = arfima(y_m)
d_m = summary(intermediate_rate_FD)$coef[[1]][1]

我要测试d_r = d_m的相等性，如下所示：

library(systemfit)
m <- matrix(c(d_r,d_l), nrow = 1)
R <- matrix(c(1,-1), nrow = 1)
linearHypothesis(m, R, test = "Chisq")

但是由于出现错误错误，代码无法正常工作：“ $运算符对原子向量无效”。任何人都可以帮助我解决问题或向我展示如何对两个分数微分参数的相等性进行Wald检验。谢谢！