MSE-R/import.R at master · tatlchri/MSE-R · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
import<-function(filename){

	#Importing data
	data <- read.delim(filename, header = TRUE, sep="", as.is=TRUE)

	#Sort is m,u,d
	data <- data[order(data[1],data[2],data[3]),]

	#Header names
	header<-colnames(data)

	#Calculating number of attributes
	noAttr<-length(data)-3-1

	#Calculating number of markets
	noM=nrow(unique(data[1]))

	#Creating a 3 dimensional list of variable length to access attributes by {market,upstream,downstream}
 	data<-setNames(data, rep(" ", length(data)))
	data<-split(data, data[[1]])
	data<-lapply(data, function(x) split(x, x[[2]]))
	#data<-lapply(data, function(x) lapply(x, function(y) lapply(y, function(z) list(z))))

	#Calculating number of upstreams and down streams in each market
	x<-(lapply(data,function(x) lapply(x,'[',2)))
	y<-lapply(x, function(y) do.call(rbind, y))
	noU<-lapply(y, function(z) lapply(z,function(k) length(unique(k))))

	x<-(lapply(data,function(x) lapply(x,'[',3)))
	y<-lapply(x, function(y) do.call(rbind, y))
	noD<-lapply(y, function(z) lapply(z,function(k) length(unique(k))))

	#Calculating distance matrices
	distanceMatrices<-lapply(data,function(x) lapply(x,'[',(3+1):(3+noAttr)))

	#Calculating matchMatrix
	matchMatrix<-lapply(data,function(x) as.vector(lapply(x,'[',3+1+noAttr)))
	matchMatrix<-lapply(seq_along(matchMatrix), function(i) lapply(seq_along(matchMatrix[[i]]), function(j) as.vector(as.matrix(matchMatrix[[i]][[j]]))))

	#Calculating mate
	mate<-lapply(matchMatrix, function(x) lapply(x,function(y) which(y==1)))
	index<-lapply(matchMatrix, function(x) 1:length(x))
	mate <- lapply(seq_along(mate), function(i){
		mate$line <-  list(index[i],mate[[i]])})

	return(list("header"=header,"noM"=noM,"noU"=noU,"noD"=noD,"noAttr"=noAttr,"distanceMatrices"=distanceMatrices,"matchMatrix"=matchMatrix,"mate"=mate))
}