# (c) Karol Bonenberg - All rights reserved

cat(paste(lubienie[2,][2-3],sep=","), file='dup', append=T, fill=T)



cat(paste(lubienie[1,][3],lubienie[1,][2],sep=","), file='dup', append=T, fill=T)

vec <- numeric(0)

for (i in 1:5) { x[length(x)+1]=i }
for (i in 1:5) { x<-c(x,i) }

for(i in 1:nrow(lubienie)) {
   if 
   row <- paste(lubienie[i,][3])
   

}

test=rbind(x)


sprawdza czy w test jest linia t
nrow(merge(test,t,1,1))>0


rm(tmp)
tmp = 0

*** skopiowanie UNIKALNYCH modularity class do nowej kolumny,
for(i in 1:nrow(lubienie)) { if(nrow(merge(tmp,lubienie[i,][3],1,1))==0) { tmp = rbind (tmp,lubienie[i,][3]) } }
*** sortowanie 


GOTOWY SKRYPT
-----------------------------
tmp = 0
for(i in 1:nrow(lubienie)) { if(nrow(merge(tmp,lubienie[i,][3],1,1))==0) { tmp = rbind (tmp,lubienie[i,][3]) } }
tmp$Modularity.Class = sort(tmp$Modularity.Class)

for(i in 1:30) {
   tmp = cbind(tmp,-1)
}

for(j in 1:nrow(tmp)) {
   for(i in 1:nrow(lubienie)) {
      if(tmp[j,][,1]==lubienie[i,][,3]) {
         for(k in 1:length(tmp[j,])) {
            if(tmp[j,][,k]==-1) { 
               tmp[j,][,k]=lubienie[i,][,2] 
               break
            }
         }
      }
   }
}

total = 0
for(k in 1:nrow(tmp)) {
   for(l in 1:length(tmp[k,])) {
      if(tmp[k,][,l]>=0) {
         total = total + 1
      }
   }
   total = total - 1
} 


population=0
estimation=350000
for(k in 1:nrow(tmp)) {
   percent = 0
   for(l in 1:length(tmp[k,])) {
      if(tmp[k,][,l]>=0) {
         percent = percent + 1
      }
   }
   percent = percent - 1
   population = rbind(population,(percent/total)*estimation)
} 
population <- population[c(-1),]
cat(paste(population), file='population.txt', append=F, fill=T)




str.mu.est <- function(N.vec,y,details="no") {
   # N.vec is a vector of the stratum sizes
   # y is a list object with each component being
   # a stratum sample
   N.ttl <- sum(N.vec)
   n.vec <- unlist(lapply(y,length))
   fpc <- (N.vec-n.vec)/N.vec
   ybar <- unlist(lapply(y,mean))
   yvar <- unlist(lapply(y,var))
   mu.hat <- sum(N.vec*ybar)/N.ttl
   var.mu <- sum(N.vecˆ2*fpc*yvar/n.vec)/(N.ttlˆ2)
   Conf.Int <- 2*sqrt(var.mu)
   if(details=="no") {
      cbind(mu.hat,var.mu,Conf.Int)
   }
   else{
      cbind(mu.hat,var.mu,Conf.Int,n.vec,ybar,yvar)
   }
}



line=0
tmp$Modularity.Class <- NULL
lista <-list()
for(m in 1:nrow(tmp)) {
   line=0
   for(n in 1:length(tmp[m,])) {
      if(tmp[m,][,n]>=0) {
         line=append(line,tmp[m,][,n])
      }
   }   
   lista = append(lista,list(line))
}



wynik = str.mu.est(N.vec=population,y=lista, details="no")
summary(wynik)



for(k in 1:nrow(tmp)) {
   for(l in 1:length(tmp[k,])) {
      if(tmp[k,][,l]==-1) {
         tmp[k,][,l]=numeric(0)
      }
   }
}
cat(paste(tmp), file='tmp.txt', append=F, fill=T)