r - 如何在循环中修复“下标越界”错误
问题描述
我正在为我的学士论文编写规范曲线。为了构建一个图表,显示哪些规范组合显示出显着的结果,哪些不是,我试图创建一个循环,导致 R 测试所有可能的组合。即使循环很短,也有是一些错误消息,不明白如何解决问题。
首先,我创建了变量来定义可能的规格:
outlier <- c("none", "40", "33,40,63", "33,40,54,63")
ineligble <- c("28", "28,38,45,48,56")
gender <- c("together", "male", "female")
regression <- c("ANOVA")
control <- c("none", "premanipulation_mean", "all hormones", "all")
data.frame“规范”定义如下:
specifications <- expand.grid(outlier = outlier, ineligble = ineligble, gender = gender, regression = regression,
control = control)
specifications <- data.frame(specifications, p_value = rep(NA, nrow(specifications)), f_value = rep(NA, nrow(specifications)),
partial_eta_square = rep(NA, nrow(specifications)), r = rep(NA, nrow(specifications)))
如您所见,我从可能的组合中构建了一个数据框,并为 p 值、f 值、部分 eta 平方和效果大小添加了列。为了用值填充额外的列,我在循环结束时添加了所需的命令。循环如下所示:
for(i in 1:nrow(specifications)){
dat <- ccy
if (specifications$outlier[i] == "none") {
dat <- dat
} else {
if (specifications$outlier[i] == "40") {
dat <- dat[-11,]
} else {
if (specifications$outlier[i] == "33,40,63") {
dat <- dat[-c(5,11,31),]
} else {
if (specifications$outlier[i] == "33,40,54,63") {
dat <- dat[-c(5,11,23,31),]
}
}
}
}
if(specifications$ineligble == "28") {
dat <- ccySC
} else {
if(specifications$ineligble == "28,38,45,48,56") {
dat <- ccy
}
}
if (specifications$gender == "together") {
dat <- dat
} else {
if(specifications$gender == "male"){
dat <- dat[which(dat$gender == "male"),]
} else {
if(specifications$gender == "female") {
dat <- dat[which(dat$gender == "female"),]
}
}
}
if (specifications$regression == "ANOVA") {
if (specifications$control == "none") {
anova <- aov(T_time2_mean ~ posecondition, data = dat)
} else {
if (specifications$control == "premanipulation_mean") {
anova <- aov(T_time2_mean ~ T_time1_mean + posecondition, data = dat)
} else {
if(specifications$control == "all hormones") {
anova <- aov(T_time2_mean ~ T_time1_mean + posecondition + C_time1_mean + C_time2_mean)
} else {
if (specifications$control == "all") {
anova <- aov(T_time2_mean ~ T_time1_mean + posecondition + C_time1_mean + C_time2_mean + sex)
}
}
}
}
specifications$p_value[i] <- drop1(anova, test = "F")$"Pr(>F)"[[3]]
specifications$f_value[i] <- drop1(anova, test = "F")$"F value"[[3]]
specifications$partial_eta_square[i] <- etaSquared(anova, type = 2, anova = F)$"eta.sq.part"[[2]]
specifications$r[i] <- sqrt(specifications$partial_eta_square[i])
specifications$k[i] <- nrow(specifications)
}
}
所以我想要它做的是用 p 值、f 值等填充四个额外的列。但是我收到以下行的错误消息“下标超出范围”:
specifications$p_value[i] <- drop1(anova, test = "F")$"Pr(>F)"[[3]]
我知道错误消息的含义,但我不知道如何解决它。在循环外的随机 anova 上测试相同的命令时,它可以工作。
使用的数据可以通过以下链接找到,称为 ccy-source-data: https ://dataverse.harvard.edu/dataset.xhtml?persistentId=doi:10.7910/DVN/FMEGS6
解决方案
[[3]]
表示您正在获取 anova 表第三行的信息。没有任何控件的模型的方差分析表只有两行。
例如,如果您对变量的结果感兴趣posecondition
,那么通过行名直接对该变量的相关信息进行子集化比假设它始终位于位置 3 更安全。
下面是一个如何为 p 值和 F 值做的例子:
results <- drop1(anova, test = "F")
specifications[i, c("p_value", "F_value")] <-
results["posecondition", c("Pr(>F)", "F value")]
推荐阅读
- jsx - 在 Oak 中使用 JSX 作为视图引擎
- php - 我没有获得所需的 URLRewrite
- python - 更新到 VS Code 1.55.0 但没有改进 python notebook
- q - Merlin - 使用 dnsmasq 将 *.localdomain 等同于 *.mydomain.com
- windows - Azure DevOps Windows 与 Ubuntu 托管的跨平台应用程序代理?
- java - 如何在项目中使用一些 libGDX 库
- python - Django:带有用于创建相关对象的字段的表单
- angular - 在 Angular 中使用守卫的更好方法
- r - 如何使用一个变量的值作为另一个变量的名称
- python-3.x - 龙卷风中的多处理