R 使用在行中存储不同p值的新列更新现有矩阵
我写了以下代码:R 使用在行中存储不同p值的新列更新现有矩阵,r,R,我写了以下代码: my.data = matrix(rnorm(20000), ncol = 20) View(my.data) empty5 = character() empty10 = character() for (variable in 1:nrow(my.data)) { nrFirst5 = my.data[variable,c(1,2,3,4,5)] nrSecond5 = my.data[variable,c(6,7,8,9,10)] nrFirst10
my.data = matrix(rnorm(20000), ncol = 20)
View(my.data)
empty5 = character()
empty10 = character()
for (variable in 1:nrow(my.data)) {
nrFirst5 = my.data[variable,c(1,2,3,4,5)]
nrSecond5 = my.data[variable,c(6,7,8,9,10)]
nrFirst10 = my.data[variable,c(1,2,3,4,5,6,7,8,9,10)]
nrSecond10 = my.data[variable,c(11,12,13,14,15,16,17,18,19,20)]
#print(nrSecord)
ttColumn5 = t.test(nrFirst5, nrSecond5, data = my.data, paired = TRUE)
pp5 = ttColumn5$p.value
#print(pp5)
ttColumn10 = t.test(nrFirst10, nrSecond10, data = my.data, paired = TRUE)
pp10 = ttColumn10$p.value
#print(tt)
#my.data$p5 = my.data[pColumn5]
#my.data$p10 = my.data[pColumn10]
pplace5 = append(empty5,pp5)
pplace10 = append(empty10,pp10)
}
my.data <- cbind(my.data, pplace5)
View(my.data)
但这给了我一个错误。像这样:
$运算符对于原子向量无效
选项1
my.data = matrix(rnorm(20000), ncol = 20)
pplace5 = character()
pplace10 = character()
for (variable in 1:nrow(my.data)) {
pp5 = t.test(my.data[variable,c(1:5)], my.data[variable,c(6:10)], paired = TRUE)$p.value
pp10 = t.test(my.data[variable,c(1:10)], my.data[variable,c(11:20)], paired = TRUE)$p.value
pplace5 = append(pplace5, pp5)
pplace10 = append(pplace10, pp10)
}
my.data <- cbind(my.data, pplace5)
my.data <- cbind(my.data, pplace10)
my.data=matrix(rnorm(20000),ncol=20)
pplace5=字符()
pplace10=字符()
for(1:nrow中的变量(my.data)){
pp5=t.test(我的数据[变量,c(1:5)],我的数据[变量,c(6:10)],配对=真)$p.value
pp10=t.test(我的数据[变量,c(1:10)],我的数据[变量,c(11:20)],配对=真)$p.value
pplace5=附加(pplace5,pp5)
pplace10=附加(pplace10,pp10)
}
我的数据
my.data = matrix(rnorm(20000), ncol = 20)
pplace5 = character()
pplace10 = character()
for (variable in 1:nrow(my.data)) {
pp5 = t.test(my.data[variable,c(1:5)], my.data[variable,c(6:10)], paired = TRUE)$p.value
pp10 = t.test(my.data[variable,c(1:10)], my.data[variable,c(11:20)], paired = TRUE)$p.value
pplace5 = append(pplace5, pp5)
pplace10 = append(pplace10, pp10)
}
my.data <- cbind(my.data, pplace5)
my.data <- cbind(my.data, pplace10)
my.data = matrix(rnorm(20000), ncol = 20)
aa <- apply(my.data, 1, function(x) {
pp5 = t.test(x[1:5], x[6:10], paired = TRUE)$p.value
pp10 = t.test(x[1:10], x[11:20], paired = TRUE)$p.value
c(pp5, pp10)
})
my.data <- cbind(my.data, t(aa))