====== Analysing European data from CIA ======

[[ru:7iss:labs:cl|Clustering]]

wdir <- "C:/Users/batagelj/work/clamix/relC"

===== Analyses =====


[[notes:da:euahc|Hierarchical clustering without constraints]]

[[notes:da:eurdc|Reading constraints in R]]

[[notes:da:eudis|Clustering with relational constraint based on the class dist]]

<code>
> idx <- function(i,j) {if (i<j) return(n*(i-1) - i*(i-1)/2 + j-i) else
+   return(n*(j-1) - j*(j-1)/2 + i-j)}
>     
> d <- vector("list",length(R))
> names(d)<-names(R)
> for(i in 1:length(R)) for(j in R[[i]])
+   d[[i]] <- append(d[[i]],list(ind=j,dis=DD[idx(i,j)]))
> sd <- d
</code>


<code>
> source("C:\\Users\\batagelj\\work\\R\\RelCon\\relConH.R")
> Ri <- sRi; Ro <- sRo
> n <- nrow(D); np <- n+1; hD <- new.env()
> for(i in 1:length(R)) for(j in R[[i]])  assign(key(i,j),DD[idx(i,j)],envir=hD)
> attr(hD,"Size") <- n; attr(hD,"Labels") <- names(R)
> res <- relConH(strategy="tolerant")
Clustering with relational constraint based on a dictionary
by Vladimir Batagelj, March 2018
Method: max   Strategy: tolerant 
[1] "Started: 2018-04-06 23:24:09"
Error in if (node[q] == 0) m[k, 2] <- -q else m[k, 2] <- node[q] : 
  missing value where TRUE/FALSE needed
</code>

[[notes:da:euanhc|Hierarchical clustering without constraints]]

===== To do =====

  * make read_Pajek in [[eurdc]] more robust (initial space in each line) 

\\

====== ======
[[notes:da:ciaeu|Back to CIA Europe data]]