【发布时间】:2020-12-02 21:19:28
【问题描述】:
我正在使用 bootstrap 创建一个新样本,以生成一个控制年龄和性别的对照组。每个变量在 2 到 8 个值之间的数字长度不同
control_40_50_HRV_female2 <- abs(parametric_bootstrap_boot2(control_40_50_HRV_female_SDNN))
control_50_60_HRV_male2 <- abs(parametric_bootstrap_boot2(control_50_60_HRV_male_SDNN))
control_50_60_HRV_female2 <-abs(parametric_bootstrap_boot2(control_50_60_HRV_female_SDNN))
control_60_70_HRV_male2 <- abs(parametric_bootstrap_boot2(control_60_70_HRV_male_SDNN))
control_60_70_HRV_female2 <-abs(parametric_bootstrap_boot2(control_60_70_HRV_female_SDNN))
control_70_80_HRV_male2 <-abs(parametric_bootstrap_boot2(control_70_80_HRV_male_SDNN))
如何将它们全部放在一组(一个变量)?所以我可以开始做 t 检验?我希望这很清楚
这是变量的示例输出
> control_40_50_HRV_female2
[1] 29.08388 102.49869
> control_50_60_HRV_male2
[1] 36.81686 127.47986 13.40681
> control_50_60_HRV_female2
[1] 25.50313
> control_60_70_HRV_male2
[1] 39.93050 140.75967 13.89545 316.45988 158.91477
> control_60_70_HRV_female2
[1] 26.27908 106.40483
当我运行这个命令时
out <- stack(mget(ls(pattern = '^control_\\d{2}_\\d{2}_\\w+_')))[2:1]
dim(out)
我明白了..这是我从脚本开头创建的所有变量的列表..
> dim(out)
[1] 96683 2
> head(out)
ind values
1 control_40_50_HRV_female_BPM 63.48
2 control_40_50_HRV_female_BPM 52.67
3 control_40_50_HRV_female_BPM 88.92
4 control_40_50_HRV_female_BPM 69.04
5 control_40_50_HRV_female_BPM 53.46
6 control_40_50_HRV_female_BPM 63.64
虽然我只需要这些变量的列表 control_40_50_HRV_female2, control_50_60_HRV_male2 .... control_70_80_HRV_male2
【问题讨论】:
-
您的评论
while I only need a list of these variable control_40_50_HRV_female2, control_50_60_HRV_male2是否需要数据子集,即subset(out, ind %in% c("control_40_50_HRV_female2", "control_50_60_HRV_male2", "control_70_80_HRV_male2"))
标签: r statistics statistics-bootstrap