Làm cách nào để thêm cột để mở rộng ma trận trong R

Question

Rất tiếc, tôi không thể nghĩ ra một tiêu đề thông tin hơn, nhưng đây là thách thức của tôi. Tôi có một ma trận và tôi cần phải thêm các cột ở những vị trí cụ thể dựa trên các tham số được mô tả bởi một vectơ. Ví dụ: nếu tôi có ma trận sau:

1, 0, 1, 2, 0 
0, 0, 1, 1, 1 
1, 1, 0, 0, 0 
2, 0, 1, 0, 2 

nhưng đối với gói R cụ thể (không được đánh dấu), tôi cần thêm cột NA vào vị trí cụ thể. Tôi có một vector liên quan đến các cột trong ma trận:

1, 1, 1, 2, 3 

Điều này cho thấy các cột 1-3 từ cùng một khoảng thời gian lấy mẫu và cột 4 và 5 là từ các giai đoạn lấy mẫu khác nhau. Tôi cần phải làm cho số cột trong ma trận bằng số lượng tối đa từ cùng một khoảng thời gian lấy mẫu lần số lần lấy mẫu. Trong trường hợp này có ba số 1 (số lượng tối đa của bất kỳ giá trị duy nhất nào trong vectơ) và tổng cộng ba khoảng thời gian lấy mẫu (số tối đa trong vectơ). Vì vậy, tôi cần một ma trận với 9 cột (3 x 3). Cụ thể, tôi cần phải thêm các cột mới của NA sau cột thứ 4 và thứ 5. Về cơ bản, tôi chỉ cần các cột của NA là phần giữ chỗ để có một ma trận trong đó số lượng quan sát (mỗi cột) là giống nhau (= 3) cho mỗi giai đoạn mẫu (được biểu thị bằng số trong vectơ). Đây là khó khăn để mô tả nhưng trong ví dụ này tưởng tượng tôi muốn kết thúc với:

1, 0, 1, 2, NA, NA, 0, NA, NA 
0, 0, 1, 1, NA, NA, 1, NA, NA 
1, 1, 0, 0, NA, NA, 0, NA, NA 
2, 0, 1, 0, NA, NA, 2, NA, NA 

này sẽ được mô tả bởi một vector mà trông giống như:

1, 1, 1, 2, 2, 2, 3, 3, 3 

mặc dù tôi không thực sự cần để tạo ra vectơ đó, chỉ là ma trận. Rõ ràng, nó rất dễ dàng để thêm những cột trong trường hợp này, nhưng đối với dữ liệu của tôi tôi có một ma trận lớn hơn nhiều mà sẽ kết thúc với ~ 200 cột. Plus tôi có thể sẽ phải làm điều này cho nhiều tập hợp dữ liệu.

Bất cứ ai có thể giúp tôi với cách để mã này trong R để tôi có thể tự động hóa quá trình mở rộng ma trận?

Cảm ơn bạn đã được tư vấn hoặc đề xuất!

---

EDIT: để làm những điều tương tự hơn một chút để dữ liệu thực tế của tôi ở đây là một ma trận tái sản xuất và vector tương tự như những người hiện tại của tôi:

m <- matrix(rpois(120*26, 1), nrow = 120, ncol = 26) 
    v <- c(1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 4, 5, 5, 6, 6, 6, 6, 7) 

Answer 1

Giả sử m là ma trận và v là vector, bạn có thể sử dụng giống như

t = table(v) 
size = dim(m)[1] * max(t) # size of each block based on the longest 
matrix(unlist(lapply(names(t), function(i) { 
       x = m[, v == i]     # get the short block 
       c(x, rep(NA, size - length(x))) # extend it to size 
     })), dim(m)[1]) 

Answer 2

Để sửa đổi ma trận cũng giống như bạn hỏi giả định ma trận là mat:

nr <- nrow(mat) 
nas <- rep(NA, nr) 
l <- lapply(4:ncol(mat), function(x) matrix(c(mat[,x],nas,nas), nrow = nr)) 
cbind(mat[,1:3], do.call(cbind,l))