table <-read.table(file="/u/home/praktikum/RNA_SandraGerstl/Tabelle_Daten", sep="\t", header=T)

#Histogramm erstellen, in dem die Laengen der Sequenzen aufgetragen sind, abspeichern der Daten

#pdf("Histogramm_Laengenverteilung.pdf")
#hist(table$Sequenzlaenge, main="Histogramm ueber die Sequenzlaenge", xlab="Sequenzlaenge")
#dev.off()

#Daten_Hist <-hist(table$Sequenzlaenge)
#Daten_Tab <- data.frame(X=Daten_Hist$mids, Y=Daten_Hist$counts) #mids = Mitte der Kategorie
Steigung_list<-vector()

   for (i in 1:1000)

  {
    #Erste Kategorie (0 - einschliessl. 500)
    A<-table[table[2]>0 & table[2] <=500,]
    B<-A[sample(1:length(A[,1]), 10, replace=FALSE),]
  

    #Zweite Kategorie (501 - einschliessl. 1000)
    C<-table[table[2]>500 & table[2] <=1000,]
    D<-C[sample(1:length(C[,1]), 10, replace=FALSE),]
   
    #Dritte Kategorie (1001 - einschliessl. 1500)
    E<-table[table[2]>1000 & table[2] <=1500,]
    F<-E[sample(1:length(E[,1]), 10, replace=FALSE),]
  

    #Vierte Kategorie (1501 - einschliessl. 2000)
    G<-table[table[2]>1500 & table[2] <=2000,]
    H<-G[sample(1:length(G[,1]), 10, replace=FALSE),]
   

    #Fuenfte Kategorie (2001 - einschliessl. 2500)
    I<-table[table[2]>2000 & table[2] <=2500,]
    J<-I[sample(1:length(I[,1]), 10, replace=FALSE),]
    
    #Sechste Kategorie (2501 - einschliessl. 3000)
    K<-table[table[2]>2500 & table[2] <=3000,]
    L<-K[sample(1:length(K[,1]), 10, replace=FALSE),]
  

    #Siebente Kategorie (3001 - einschliessl. 3500)
    M<-table[table[2]>3000 & table[2] <=3500,]
    N<-M[sample(1:length(M[,1]), 10, replace=FALSE),]
  

    Zufall_gesamt<-rbind(B,D,F,H,J,L,N)
    #plot(B$Sequenzlaenge, B$MCC)
    Steigung_ges<-lm(Zufall_gesamt$MCC~Zufall_gesamt$Sequenzlaenge)$coefficients[2]

    Steigung_list<-c(Steigung_list, as.numeric(Steigung_ges))

}

pdf("MCC_Steigung.pdf")
plot(Steigung_list,ylim=c(-.001,.001), pch='.', xlab="Distanz", ylab="Steigung der Regressionsgeraden")
dev.off()

savePlot(filename="Bucket.png", type="png")