R语言两正态总体方差相等检验拒绝域代码

在R语言中，对于两个正态总体方差相等的假设检验，可以使用F检验。F检验的拒绝域可以使用qf函数计算得到，具体代码如下：

n <- 100
m <- 100
alpha <- 0.05
df1 <- n - 1
df2 <- m - 1
q1 <- qf(alpha/2, df1, df2, lower.tail = FALSE)
q2 <- qf(1-alpha/2, df1, df2, lower.tail = FALSE)
rejection.region <- c(0, q1, q2, Inf)

其中，n和m分别为两个正态分布的样本大小，alpha为显著性水平，df1和df2分别为自由度，q1和q2为临界值，rejection.region为拒绝域。

如果你想进行双侧检验，那么拒绝域为$F < \frac{1}{F_{1-\alpha/2}(n-1,m-1)}$和$F > F_{1-\alpha/2}(n-1,m-1)$，其中$F_{1-\alpha/2}(n-1,m-1)$是自由度为$(n-1,m-1)$，显著性水平为$1-\alpha/2$的F分布的上分位数。如果你想进行单侧检验，那么只需要使用$1-\alpha$的显著性水平计算F分布的上分位数即可。