将数据从向量添加到矩阵

Question

我有一个国家向量 p 和一个国家矩阵 Z。 p中的一些国家不包括在Z中。 我想要一个论点，如果它们尚未包含在 Z 中，则将国家从 p 添加到 Z，但前提是它们的值 > 0。

以下是数据示例。在预期输出中，应将澳大利亚添加到 Z，所有矩阵列中都为 0，因为它没有要报告的任何矩阵数据。矩阵尺寸应为 5x5。

p <- c(Afghanistan = 2769, Albania = 93893300, Algeria = 0, Argentina=4550430, Australia=50, Ban = 2)
p <- stack(p)
names(p) <- c("Value", "Area")
#vector

Z <- matrix(c(0,138201.333333333,0,0,0,0,1162.33333333333,0,38.3333333333333,0,0,0,0,0,0,0,300238,0,9675,0), nrow = 4, ncol = 4, byrow=TRUE)
dimnames(Z) = list(c("Afghanistan","Albania","Algeria", "Argentina"),c("Afghanistan","Albania","Algeria", "Argentina")) #matrix
#matrix

Answer 1

cbind 和 rbind 的组合可以解决问题。我在 p 中添加了一个附加元素，只是为了表明它甚至可以使用更多名称。

# your data EDITED BASED ON THE DPUT you posted!!
p <- c(Afghanistan = 2769, Albania = 93893300, Algeria = 0, Argentina=4550430, Australia=50, Ban = 2)
p <- stack(p)
names(p) <- c("Value", "Area")

Z <- matrix(c(0,138201.333333333,0,0,0,0,1162.33333333333,0,38.3333333333333,0,0,0,0,0,0,0,300238,0,9675,0), nrow = 4, ncol = 4, byrow=TRUE)
dimnames(Z) = list(c("Afghanistan","Albania","Algeria", "Argentina"),c("Afghanistan","Albania","Algeria", "Argentina")) #matrix

# solution
to_add <- setdiff(p[p$Value>0, "Area"], unlist(dimnames(Z)))

Z <- cbind(Z, matrix(0, nrow(Z), length(to_add), dimnames = list(NULL, to_add)))
Z <- rbind(Z, matrix(0, length(to_add), ncol(Z), dimnames = list(to_add, NULL)))

Z
#>             Afghanistan  Albania  Algeria Argentina Australia Ban
#> Afghanistan     0.00000 138201.3    0.000         0         0   0
#> Albania         0.00000      0.0 1162.333         0         0   0
#> Algeria        38.33333      0.0    0.000         0         0   0
#> Argentina       0.00000      0.0    0.000         0         0   0
#> Australia       0.00000      0.0    0.000         0         0   0
#> Ban             0.00000      0.0    0.000         0         0   0

Answer 2

在您的评论中有新信息后的编辑版本：检查：

# your given matrix
p <- as.matrix(c(Afghanistan = 2769, 
                 Albania = 93893300, 
                 Algeria = 0, 
                 Argentina=4550430, 
                 Australia=50)) # a matrix of 5 values

p_transpose <- t(p)   # transpose matrix p

for(j in seq(ncol(p)))
  if(! all(p[, j] == p_transpose[j, ])) stop("wrong transpose")  # https://stat.ethz.ch/R-manual/R-devel/library/base/html/t.html


p_transpose_dataframe <- as.data.frame(p_transpose)  # p_transpose to dataframe


# your given main matrix
Z <- matrix(c(0,138201.333333333,0,0,
                  0,0,1162.33333333333,0,
                  38.3333333333333,0,0,0,
                  0,0,0,0,
                  300238,0,9675,0), nrow = 5, ncol = 4, byrow=TRUE) # matrix Z
colnames(Z) <- c("Afghanistan","Albania","Algeria", "Argentina")    # name columnes

Z <- as.data.frame(Z)   # make a dataframe Z

# helper function: code from here: # reference: https://stackoverflow.com/questions/45857787/adding-column-if-it-does-not-exist
fncols <- function(data, cname) {
  add <-cname[!cname%in%names(data)]
  
  if(length(add)!=0) data[add] <- 0
  data
} 

newZ <- fncols(Z, colnames(p_transpose_dataframe))   # helper functions to add not existing columns to Z and assign to new dataframe


result_Z <- rbind(newZ, p_transpose_dataframe)     # result dataframe
View(result_Z)                                     # View result dataframe

请提供反馈。