yyh
/
BioStaInfer
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

更新 exam02_example.Rmd

This commit is contained in:
yyh 2025-03-20 22:55:23 +08:00
parent f73e7ffb46
commit 8c42c72cd9
1 changed files with 4 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
exam02_example.Rmd
View File

@ -212,8 +212,11 @@ cat(group_levels[1], "组均值与", group_levels[2], "组均值的差异:",
# "greater" - 单侧检验，假设第一组大于第二组 # "greater" - 单侧检验，假设第一组大于第二组
# "less" - 单侧检验，假设第一组小于第二组 # "less" - 单侧检验，假设第一组小于第二组
test_alternative <- "two.sided" # 可根据研究假设修改 test_alternative <- "two.sided" # 可根据研究假设修改
test_paired <- 'F'
```
### 检查统计检验的假设 ### 检查统计检验的假设
```{r} ```{r}
if(normality_assumption_met) { if(normality_assumption_met) {
cat("数据满足正态性假设，使用t检验\n") cat("数据满足正态性假设，使用t检验\n")
@ -223,6 +226,7 @@ if(normality_assumption_met) {
cat("数据满足方差同质性假设，使用等方差t检验(var.equal = TRUE)\n") cat("数据满足方差同质性假设，使用等方差t检验(var.equal = TRUE)\n")
t_test <- t.test(value ~ group, data = data, t_test <- t.test(value ~ group, data = data,
alternative = test_alternative, alternative = test_alternative,
paired = test_paired,
var.equal = TRUE) var.equal = TRUE)
} else { } else {
cat("数据不满足方差同质性假设，使用Welch's t检验(var.equal = FALSE)\n") cat("数据不满足方差同质性假设，使用Welch's t检验(var.equal = FALSE)\n")