Мне нужно повторить этот код несколько раз (44 раза) и получить 44 разные матрицы mat. Данные загружаю отсюда https://intl-atlas-downloads.s3.amazonaws.com/index.html догадываюсь, что это можно сделать при помощи "for-loop", но пока до конца не понимаю, как именно. Буду признателен за подсказки.
library(readstata13)
library(diverse)
library(plyr)
setwd('C:/Users/Admin/Desktop/work/world HH')
filenames <- paste0('S2_final_',1962:2016,'.dta')
filenames <- c(filenames,paste0('S2_final_',1962:2016,'.dta'))
matr=NULL
g=23
for (j in 1:g) {
dat <- read.dta13(filenames[j]) #pickup
data = aggregate(dat$export_value, by = list(dat$exporter,dat$commoditycode), FUN = sum)
colnames(data) = c("land","product","value")
dt = split(data, f = data$product)
land = as.data.frame(sort(unique(data[, 1])))
nds = seq(1, nrow(land), by = 1)
texmat = cbind(nds, land)
colnames(texmat) = c("num", "land")
for (i in 1:length(unique(data[, 2]))) {
(join(texmat, dt[[i]], by = "land", type = "left")$value)
}
mt = sapply(1:length(unique(data[, 2])), function(i) join(texmat, dt[[i]], by = "land", type = "left")$value)
colnames(mt) = unique(data[, 2])
rownames(mt) = sort(unique(data[, 1]))
mt[is.na(mt)] = 0
rcamat=values(mt, category_row = FALSE, norm = "rca",filter = 1, binary = TRUE)
rcamat[is.na(rcamat)] = 0
tmat = rcamat[rowSums(rcamat) != 0, , drop = TRUE]
mat = t(tmat)
matr[[j]] <- mat
}
state=NULL
for (n in 1:g) {
state[[n]]=matr[[n]][,"FIN"]
}
stcom=sapply(state, '[', seq(max(sapply(state,length))))
stc <- stcom[complete.cases(stcom), ]
colnames(stc)=paste("FIN", 1:g, sep = "")