setwd("C:/sepp/work/fabia/experiments")


library(fabia)
library(biclust)
library(truecluster)
library(BicARE)

m <- 100


indrorig_s<-list()
indrfabia_s<-list()
indrfabias_s<-list()
indrfabiap_s<-list()
indrmfsc_s<-list()
indrplaid_ss_s<-list()
indrplaid_ms_s<-list()
indrplaid_ms_5_s<-list()
indrisa_s<-list()
indrisa_1_s<-list()
indropsm_s<-list()
indrsamba_s<-list()
indrxmotif_s<-list()
indrbimax_s<-list()
indrcc_s<-list()
indrplaid_t_ab_s<-list()
indrplaid_t_a_s<-list()
indrspec_s<-list()
indrfloc_s<-list()

sindrorig_s<- as.vector(rep(0,12))
sindrfabia_s<- as.vector(rep(0,12))
sindrfabias_s<- as.vector(rep(0,12))
sindrfabiap_s<- as.vector(rep(0,12))
sindrmfsc_s<-  as.vector(rep(0,12))
sindrplaid_ss_s<- as.vector(rep(0,12))
sindrplaid_ms_s<-  as.vector(rep(0,12))
sindrplaid_ms_5_s<- as.vector(rep(0,12))
sindrisa_s<-  as.vector(rep(0,12))
sindrisa_1_s<- as.vector(rep(0,12))
sindropsm_s<-  as.vector(rep(0,12))
sindrsamba_s<- as.vector(rep(0,12))
sindrxmotif_s<- as.vector(rep(0,12))
sindrbimax_s<-  as.vector(rep(0,12))
sindrcc_s<-  as.vector(rep(0,12))
sindrplaid_t_ab_s<- as.vector(rep(0,12))
sindrplaid_t_a_s<-  as.vector(rep(0,12))
sindrspec_s<-  as.vector(rep(0,12))
sindrfloc_s<-  as.vector(rep(0,12))



for (i in 1:m) {

    load(paste("exp_",as.character(i),"_ind_res.RData",sep = ""))


    indrfabia_s[[i]]<-indrfabia
    indrfabias_s[[i]]<-indrfabias
    indrfabiap_s[[i]]<-indrfabiap
    indrmfsc_s[[i]]<-indrmfsc
    indrplaid_ss_s[[i]]<-indrplaid_ss
    indrplaid_ms_s[[i]]<-indrplaid_ms
    indrplaid_ms_5_s[[i]]<-indrplaid_ms_5
    indrisa_s[[i]]<-indrisa
    indrisa_1_s[[i]]<-indrisa_1
    indropsm_s[[i]]<-indropsm
    indrsamba_s[[i]]<-indrsamba
    indrxmotif_s[[i]]<-indrxmotif
    indrbimax_s[[i]]<-indrbimax
    indrcc_s[[i]]<-indrcc
    indrplaid_t_ab_s[[i]]<-indrplaid_t_ab
    indrplaid_t_a_s[[i]]<-indrplaid_t_a
    indrspec_s[[i]]<-indrspec
    indrfloc_s[[i]]<-indrfloc

    sindrfabia_s<-sindrfabia_s  + indrfabia
    sindrfabias_s<- sindrfabias_s + indrfabias
    sindrfabiap_s<-sindrfabiap_s  + indrfabiap
    sindrmfsc_s<- sindrmfsc_s + indrmfsc
    sindrplaid_ss_s<- sindrplaid_ss_s + indrplaid_ss
    sindrplaid_ms_s<- sindrplaid_ms_s + indrplaid_ms
    sindrplaid_ms_5_s<-sindrplaid_ms_5_s  + indrplaid_ms_5
    sindrisa_s<- sindrisa_s + indrisa
    sindrisa_1_s<-sindrisa_1_s  + indrisa_1
    sindropsm_s<- sindropsm_s + indropsm
    sindrsamba_s<-sindrsamba_s  + indrsamba
    sindrxmotif_s<-sindrxmotif_s  + indrxmotif
    sindrbimax_s<- sindrbimax_s + indrbimax
    sindrcc_s<- sindrcc_s + indrcc
    sindrplaid_t_ab_s<-sindrplaid_t_ab_s  + indrplaid_t_ab
    sindrplaid_t_a_s<- sindrplaid_t_a_s + indrplaid_t_a
    sindrspec_s<- sindrspec_s + indrspec
    sindrfloc_s<- sindrfloc_s + indrfloc
}

    sindrfabia_s<-sindrfabia_s/m
    sindrfabias_s<- sindrfabias_s/m
    sindrfabiap_s<-sindrfabiap_s/m
    sindrmfsc_s<- sindrmfsc_s/m
    sindrplaid_ss_s<- sindrplaid_ss_s/m
    sindrplaid_ms_s<- sindrplaid_ms_s/m
    sindrplaid_ms_5_s<-sindrplaid_ms_5_s/m
    sindrisa_s<- sindrisa_s/m
    sindrisa_1_s<-sindrisa_1_s/m
    sindropsm_s<- sindropsm_s/m
    sindrsamba_s<-sindrsamba_s/m
    sindrxmotif_s<-sindrxmotif_s/m
    sindrbimax_s<- sindrbimax_s/m
    sindrcc_s<- sindrcc_s/m
    sindrplaid_t_ab_s<-sindrplaid_t_ab_s/m
    sindrplaid_t_a_s<- sindrplaid_t_a_s/m
    sindrspec_s<- sindrspec_s/m
    sindrfloc_s<- sindrfloc_s/m

    sindrfabia_v<- as.vector(rep(0,12))
    sindrfabias_v<- as.vector(rep(0,12))
    sindrfabiap_v<- as.vector(rep(0,12))
    sindrmfsc_v<-  as.vector(rep(0,12))
    sindrplaid_ss_v<- as.vector(rep(0,12))
    sindrplaid_ms_v<-  as.vector(rep(0,12))
    sindrplaid_ms_5_v<- as.vector(rep(0,12))
    sindrisa_v<-  as.vector(rep(0,12))
    sindrisa_1_v<- as.vector(rep(0,12))
    sindropsm_v<-  as.vector(rep(0,12))
    sindrsamba_v<- as.vector(rep(0,12))
    sindrxmotif_v<- as.vector(rep(0,12))
    sindrbimax_v<-  as.vector(rep(0,12))
    sindrcc_v<-  as.vector(rep(0,12))
    sindrplaid_t_ab_v<- as.vector(rep(0,12))
    sindrplaid_t_a_v<-  as.vector(rep(0,12))
    sindrspec_v<-  as.vector(rep(0,12))
    sindrfloc_v<-  as.vector(rep(0,12))

    for (i in 1:m) {


    sindrfabia_v<-(indrfabia_s[[i]]-sindrfabia_s)^2
    sindrfabias_v<-(indrfabias_s[[i]]-sindrfabias_s)^2
    sindrfabiap_v<-(indrfabiap_s[[i]]-sindrfabiap_s)^2
    sindrmfsc_v<-(indrmfsc_s[[i]]-sindrmfsc_s)^2
    sindrplaid_ss_v<-(indrplaid_ss_s[[i]]-sindrplaid_ss_s)^2
    sindrplaid_ms_v<-(indrplaid_ms_s[[i]]-sindrplaid_ms_s)^2
    sindrplaid_ms_5_v<-(indrplaid_ms_5_s[[i]]-sindrplaid_ms_5_s)^2
    sindrisa_v<-(indrisa_s[[i]]-sindrisa_s)^2
    sindrisa_1_v<-(indrisa_1_s[[i]]-sindrisa_1_s)^2
    sindropsm_v<-(indropsm_s[[i]]-sindropsm_s)^2
    sindrsamba_v<-(indrsamba_s[[i]]-sindrsamba_s)^2
    sindrxmotif_v<-(indrxmotif_s[[i]]-sindrxmotif_s)^2
    sindrbimax_v<-(indrbimax_s[[i]]-sindrbimax_s)^2
    sindrcc_v<-(indrcc_s[[i]]-sindrcc_s)^2
    sindrplaid_t_ab_v<-(indrplaid_t_ab_s[[i]]-sindrplaid_t_ab_s)^2
    sindrplaid_t_a_v<-(indrplaid_t_a_s[[i]]-sindrplaid_t_a_s)^2
    sindrspec_v<-(indrspec_s[[i]]-sindrspec_s)^2
    sindrfloc_v<-(indrfloc_s[[i]]-sindrfloc_s)^2


}

    sindrfabia_v<-sindrfabia_v/(m-1)
    sindrfabias_v<- sindrfabias_v/(m-1)
    sindrfabiap_v<-sindrfabiap_v/(m-1)
    sindrmfsc_v<- sindrmfsc_v/(m-1)
    sindrplaid_ss_v<- sindrplaid_ss_v/(m-1)
    sindrplaid_ms_v<- sindrplaid_ms_v/(m-1)
    sindrplaid_ms_5_v<-sindrplaid_ms_5_v/(m-1)
    sindrisa_v<- sindrisa_v/(m-1)
    sindrisa_1_v<-sindrisa_1_v/(m-1)
    sindropsm_v<- sindropsm_v/(m-1)
    sindrsamba_v<-sindrsamba_v/(m-1)
    sindrxmotif_v<-sindrxmotif_v/(m-1)
    sindrbimax_v<- sindrbimax_v/(m-1)
    sindrcc_v<- sindrcc_v/(m-1)
    sindrplaid_t_ab_v<-sindrplaid_t_ab_v/(m-1)
    sindrplaid_t_a_v<- sindrplaid_t_a_v/(m-1)
    sindrspec_v<- sindrspec_v/(m-1)
    sindrfloc_v<- sindrfloc_v/(m-1)

     write.table(t(as.vector(sindrfabia_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =FALSE)
     write.table(t(as.vector(sindrfabias_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrfabiap_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrmfsc_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrplaid_ss_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrplaid_ms_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrplaid_ms_5_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrisa_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrisa_1_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindropsm_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrsamba_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrxmotif_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrbimax_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrcc_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrplaid_t_ab_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrplaid_t_a_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrspec_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrfloc_s)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)



     write.table(t(as.vector(sindrfabia_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrfabias_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrfabiap_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrmfsc_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrplaid_ss_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrplaid_ms_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrplaid_ms_5_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrisa_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrisa_1_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindropsm_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrsamba_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrxmotif_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrbimax_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrcc_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrplaid_t_ab_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrplaid_t_a_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrspec_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)
     write.table(t(as.vector(sindrfloc_v)), file ="final_result.txt", quote = FALSE, sep = "\t",row.names = FALSE,col.names = FALSE,append =TRUE)