SID_tcdf.R

library(devtools)

load_all('C:\\Users\\flori\\Downloads\\processx')
load_all('C:\\Users\\flori\\Downloads\\pcalg')

load_all('C:\\Users\\flori\\Downloads\\SID')


architecturelist <- c("fork", "Mediator", "Vstructure", "Diamond", "7TS", "7ts2h")

nbvarlist <- c(3,3,3,4,7,7)

sidlist <- c()

folder <- "tunningts=1000"



for(k in 1:6){
	architecture <- architecturelist[k]
	nbvar <- nbvarlist[k]
	print(architecture)


	dir <- paste('C:\\Users\\flori\\Downloads\\causal_discovery_for_time_series-master\\causal_discovery_for_time_series-master\\baselines\\scripts_python\\python_packages\\TCDF-master\\TCDF-Master\\', folder, sep="")
	listdir <- list.files(dir)


	gtfile <- paste('C:\\Users\\flori\\Desktop\\Stage_2A\\Code\\Results\\groundtruth\\', architecture, '_groundtruth.csv', sep="")
	gtgraph <- read.csv(gtfile, header=FALSE)
	gtmatrix <- matrix(0,nbvar,nbvar)

	for(i in 1:nrow(gtgraph)){
		gtmatrix[gtgraph[i,1]+1,gtgraph[i,2]+1] <- 1			
	}
	print(gtmatrix)

	for(file in listdir){
		if(endsWith(file, '.csv') & file!='Evaluation.csv' & substring(file,1,nchar(file)-5)==architecture){
			graph <- read.csv(paste(dir,'\\',file, sep=""),header= FALSE)
			adjmatrix <- matrix(0,nbvar,nbvar)
			for(i in 1:nrow(graph)){
				adjmatrix[graph[i,1]+1,graph[i,2]+1] <- 1			
			}

			for(i in 1:ncol(adjmatrix)){
				for(j in 1:nrow(adjmatrix)){
					if(adjmatrix[i,j]==1 & adjmatrix[j,i]==1 & i!=j){
						if(adjmatrix[i,j]==gtmatrix[i,j] & adjmatrix[j,i]!=gtmatrix[j,i]){
							adjmatrix[i,j] <- 0
						}
						else if(adjmatrix[j,i]==gtmatrix[j,i] & adjmatrix[i,j]!=gtmatrix[i,j]){
							adjmatrix[j,i] <- 0
						}
					}
				}
			}
			#print(adjmatrix)

			sid = tryCatch({
				structIntervDist(gtmatrix,adjmatrix)
				}, error=function(e){
				print("error in SID")
				print(file)
				}
				)
			sid <- sid[[1]]
			#print(sid)
			#sidlist <- c(sidlist, c(file,sid[[1]]))
			#print(typeof(sid))
			if(typeof(sid[1])=="double"){
				sidlist <- c(sidlist, sid[1])
			}
		}
	}
	

}

print("list of the SID")
print(sidlist)
print(length(sidlist))
print(mean(sidlist))
print(sd(sidlist))