Question

我有一个90行和4列的数据帧，我想在软件R中生成矩阵（4 * 90）x4，第一列是90数据，另一个是从91到360的0。第二列是从1到0在图90中，从91到180填充数据，然后行181到360是0.依此类推。

我已经成为这样的语法

psi10<-c(c(data[,1]),rep(0,3*90))

psi20<-c(rep(0,90),c(data[,2]),rep(0,2*90))

psi30<-c(rep(0,90),rep(0,90),c(data[,3]),rep(0,90))

psi40<-c(rep(0,3*90),c(data[,4]))

Mpsi<-as.matrix(data.frame(psi10,psi20,psi30,psi40))

但是我认为它太长了，如果添加更多列，会更加困难请帮助我解决我的问题。感谢

Answer 1

你可以尝试

library(Matrix)
m1 <- as.matrix(bdiag(lapply(data, matrix, 1,90)))

或使用sapply

中的base R

n <- 4
m2 <- t(sapply(seq_len(n), function(i) {c(rep(0,90*(i-1)), data[,i],
                                              rep(0, 90*(n-i)))}))

all.equal(m1, m2, check.attributes=FALSE)
#[1] TRUE

更新

似乎OP想要“转置”，在这种情况下

m1 <- as.matrix(bdiag(data))

与

相同

 m1N <- as.matrix(t(bdiag(lapply(data, matrix, 1,90))))
 identical(m1, m1N)
 #[1] TRUE

数据

set.seed(24)
 data <- as.data.frame(matrix(sample(1:4, 90*4, replace=TRUE), ncol=4))

Answer 2

另一个想法可能是：

DF = data.frame(c1 = 1:3, c2 = 4:6, c3 = 7:9)
DF
#  c1 c2 c3
#1  1  4  7
#2  2  5  8
#3  3  6  9

mat = matrix(0, nrow = prod(dim(DF)), ncol = ncol(DF))
mat[cbind(seq_len(prod(dim(DF))), rep(seq_along(DF), each = nrow(DF)))] = unlist(DF)
mat
#      [,1] [,2] [,3]
# [1,]    1    0    0
# [2,]    2    0    0
# [3,]    3    0    0
# [4,]    0    4    0
# [5,]    0    5    0
# [6,]    0    6    0
# [7,]    0    0    7
# [8,]    0    0    8
# [9,]    0    0    9

Answer 3

像这样的工作吗？

df<-data.frame(a=100:189,b=100:189,c=100:189,d=100:189)
df<-as.matrix(df)

df<-rbind(df, matrix(0,nrow=(ncol(df)-1)*90,ncol=ncol(df)))
for (x in 1:ncol(df)){
  df[,x]<-c(rep(0,90*(x-1)),df[1:90,x],rep(0,90*(ncol(df)-x)))
}

如何用对角线制作矩阵是数据而另一个是0？

3 个答案:

更新

数据