在成对的调查年之后过滤并返回 Fischer-方差比等式的函数

如何解决在成对的调查年之后过滤并返回 Fischer-方差比等式的函数

我有一个数据集，其中包含从 1980 年到 2020 年这几年收集的信息。该数据集是横截面的、具有国家代表性的 (n > 1,500)。

总的来说，有 150 个感兴趣的变量（例如，生活满意度、态度等）。我想知道相邻调查年份（例如 1980-1982、1982-1984、1984-1986 等）之间 150 个感兴趣变量中的每一个的方差是否相等。

最终，我希望能够绘制方差估计和置信区间的相等性。

year  var1  var2  var3  var4 ...

1980   1     2      5     3
1980   2     3      2     3
1980   1     Na     2     3
1982   3     1      2     5
1982   4     2      2     4
1982   5     3      Na    2
1984   2     1      3     2
1984   1     4      5     1
...

这适用于单个 var ~ year 组合。但总体而言，有 >15 年对 x 150 个变量，因此对于每个组合几乎不可能做到。

df %>%
filter(year==1980 | year==1982) %>%
var.test(ep01 ~ year,.) -> ftst # calculates the Fischer F-test for equality of variances

下面从上面的 ftst 中选择 F-test 估计。这就是我想要为每个 >15 年对 x 150 个变量保存到一个单独的数据框中（用于以后绘图）的内容。

estimate <- { 
    f.estm <- (ftst$estimate)
   f.estm[[1]] 
  }
ratio <- print(paste(estimate))