我有一个变量,它看起来像这样:
Var
[1] 3, 4, 5 2, 4, 5 2, 4 1, 4, 5
我需要把它拆分成一个数据帧,看起来像这样:
V1 V2 V3 V4 V5
NA NA 3 4 5
NA 2 NA 4 5
NA 2 NA 4 NA
1 NA NA 4 5
不幸的是,我找不到一个能解决我问题的职位。有人知道我是怎么做到的吗?提前非常感谢!
编辑:我根据你的回答找到了一个解决方案,并将其发布在下面。
第二版:我用Ananda的解决方案提高了代码的效率。
使用矩阵索引:
Var <- list(c(3,4,5),c(2,4,5),c(2,4),c(1,4,5))
unVar <- unlist(Var)
out <- matrix(NA, nrow=length(Var), ncol=max(unVar))
out[cbind(rep(seq_along(Var),sapply(Var,length)),unVar)] <- unVar
# and if you're using the new version of R, you can simplify a little:
out[cbind(rep(seq_along(Var),lengths(Var)),unVar)] <- unVar
# [,1] [,2] [,3] [,4] [,5]
#[1,] NA NA 3 4 5
#[2,] NA 2 NA 4 5
#[3,] NA 2 NA 4 NA
#[4,] 1 NA NA 4 5
根据OP的答案判断,"var"是一个字符串,如:
var <- c("3, 4, 5", "2, 4, 5", "2, 4", "1, 4, 5")
如果是这种情况,您可以考虑我的"splitstackshape"包中的cSplit_e
:
library(splitstackshape)
cSplit_e(data.frame(var), "var", ",", mode = "value", drop = TRUE)
# var_1 var_2 var_3 var_4 var_5
# 1 NA NA 3 4 5
# 2 NA 2 NA 4 5
# 3 NA 2 NA 4 NA
# 4 1 NA NA 4 5
如果它是list
,正如其他答案所假设的那样,您可以在"splitstackshape"中使用(未导出的(numMat
函数来为cSplit_e
供电。
var <- list(c(3,4,5), c(2,4,5), c(2,4), c(1,4,5))
splitstackshape:::numMat(var, mode = "value")
# 1 2 3 4 5
# [1,] NA NA 3 4 5
# [2,] NA 2 NA 4 5
# [3,] NA 2 NA 4 NA
# [4,] 1 NA NA 4 5
在幕后,numMat
与@thelatemail的回答中使用的方法非常相似。
如果你有-99代表NA
,并且你想排除它们,你可以尝试:
var <- c("3, 4, 5", "2, -99, 4, 5", "2, 4", "1, 4, 5, -99")
splitstackshape:::numMat(
lapply(strsplit(var, ","), function(x) as.numeric(x)[as.numeric(x) > 0]),
mode = "value")
# 1 2 3 4 5
# [1,] NA NA 3 4 5
# [2,] NA 2 NA 4 5
# [3,] NA 2 NA 4 NA
# [4,] 1 NA NA 4 5
Var <- list(c(3, 4, 5), c(2, 4, 5), c(2, 4), c(1, 4, 5))
M <- matrix(NA, nrow=length(Var), ncol=max(sapply(Var,max)))
for( L in seq(Var) ) { M [ cbind( rep( L, length(Var[[L]])), Var[[L]]) ] <- Var[[L]]}
M
[,1] [,2] [,3] [,4] [,5]
[1,] NA NA 3 4 5
[2,] NA 2 NA 4 5
[3,] NA 2 NA 4 NA
[4,] 1 NA NA 4 5
就我个人而言,我的投票建议是泰梅尔的版本,它基本上与此同构。
如果Var只是一个向量,那么我会执行以下操作:
Var = c(3,4,5,2,4,5,2,4,1,4,5)
RowIdx = c(rep(1,3),rep(2,3),rep(3,2),rep(4,3))
DF = matrix(NA,nrow=4,ncol=5)
for (idx in 1:length(Var)){
DF[RowIdx[idx],Var[idx]] = Var[idx]
}
当然,如果您有更多的数据,您可能希望找到一种以更自动化的方式生成行索引的方法
我设法根据您的回复找到了解决方案!我的最终解决方案如下:
# I had the additional problem that my variable was a factor, therefore I had to transform it first.
df <- data.frame(Var)
Var <- lapply(strsplit(as.character(df$Var), ", "), "[")
for(i in 1:length(Var)){
Var[[i]] <- as.numeric(Var[[i]])
}
# Then I created a matrix based on thelatemails and BondedDusts approach.
M <- matrix(NA, nrow=length(Var), ncol=max(sapply(Var,max)))
# Additionally, I had the problem that there were some lines with a single -99, which indicates a missing value for the complete line. I had some problems with this negative value. For this reason, I assigned NA's first.
for(i in 1:length(Var)){
Var[[i]][Var[[i]] == -99] <- NA
}
# Final assignment like suggested by BonedDust.
for( L in seq(Var) ) { M [ cbind( rep( L, length(Var[[L]])), Var[[L]]) ] <- Var[[L]]}
M
我不确定这是否是最快的解决方案,但现在一切都很好!非常感谢您快速而广泛的回答!