@@ -382,11 +382,9 @@ runSilhouetteIMG <- function(data, k) {

       sil.c$clus.avg.silwidths=summary(sil.w)$clus.avg.widths

       sil.c$avg.silwidths=summary(sil.w)$avg.width

-

       estable[[i.datos]]$kmk.dynamic.bs.or.numeric=part.onto

       estable[[i.datos]]$sil.width=sil.w

-      #DESDE AQUI new

       g.main=paste(metric.name,sep="")

       plot(sil.w, col=colores, main=g.main, border=NULL,

@@ -400,12 +398,14 @@ runSilhouetteIMG <- function(data, k) {

       legend(x=x.leyenda+0.1,y=sil.c$n+1, legend=expression(' s'['i']), col="black",

              xjust=0, yjust=0, bty="n", xpd=TRUE, inset=c(-0.1,0), cex=escalax)

       xleyenda=rep(x.leyenda,k.cl)

+      #yleyenda=(sil.c$cluster.size==1)*0.6*(sil.c$n-cumsum(sil.c$cluster.size))+

+      #  (sil.c$n-cumsum(sil.c$cluster.size))+sil.c$cluster.size*3/k.cl+2

       yleyenda=(sil.c$cluster.size==1)*0.6*(sil.c$n-cumsum(sil.c$cluster.size))+

-        (sil.c$n-cumsum(sil.c$cluster.size))+sil.c$cluster.size*3/k.cl+2

+        (sil.c$n-cumsum(sil.c$cluster.size))+sil.c$cluster.size/k.cl+2

       leyenda=paste(names(sil.c$clus.avg.silwidths),rep(": ",k.cl),

                     sil.c$cluster.size,"|",round(sil.c$clus.avg.silwidths,digits=2),sec="")

       for (i.leyenda in 1:k.cl){

-        legend(list(x=xleyenda[i.leyenda],y=yleyenda[i.leyenda]-0.7), legend=leyenda[i.leyenda], col="black",

+        legend(list(x=xleyenda[i.leyenda],y=yleyenda[i.leyenda]), legend=leyenda[i.leyenda], col="black",

                xjust=0, yjust=1, bty="n", xpd=TRUE, inset=c(-0.1,0), cex=escalal)

         }

     }

@@ -276,11 +276,12 @@ runStabilityIndexK_IMG <- function(bs, k.min, k.max) {

   for (j.k in i.min:i.max) {

     cur.data = e.stab.global[[j.k]]

     cur.data = cur.data[!is.na(cur.data)]

+    ymin = min(cur.data)

     xnames=as.character(names.metr)

     ynames="Global Stability Indices"

     g.main=paste(" St. Indices of the metrics for k=", j.k,sep="")

     plot(cur.data, main=g.main, axes=TRUE, col.axis="white",

-         xlim=c(0.75,length(xnames)+0.25), xlab="", ylim=c(0.3,1),

+         xlim=c(0.75,length(xnames)+0.25), xlab="", ylim=c(ymin,1),

          ylab=ynames, col=colores[1],type="o", lwd=1, lty=ltype[1])

     axis(1,at=1:length(xnames),labels=xnames,las=2,cex.axis=0.75)

     axis(2,las=3,cex.axis=0.85)

@@ -316,13 +317,13 @@ runStabilityIndexMetric_IMG <- function(bs, k.min, k.max) {

   for (i.metr in 1:length(names.metr)) {

     cur.data = m.stab.global[[i.metr]]

     cur.data = cur.data[!is.na(cur.data)]

-

+    ymin = min(cur.data)

     xnames=c(k.min:k.max)

     ynames="Global Stability Indices"

     g.main=paste(" St. Indices of '", names.metr[i.metr], "' for k in [",

                  k.min, ",", k.max,"]",sep="")

     plot(cur.data, main=g.main, axes=TRUE, col.axis="white",

-         xlim=c(0.75,length(k.min:k.max)+0.25), xlab="", ylim=c(0.3,1),

+         xlim=c(0.75,length(k.min:k.max)+0.25), xlab="", ylim=c(ymin,1),

          ylab=ynames, col=colores[1],type="o", lwd=1, lty=ltype[1])

     axis(1,at=1:length(k.min:k.max),labels=xnames,las=1,cex.axis=escalax)

     axis(2,las=3,cex.axis=0.85)

@@ -2,7 +2,6 @@ library(evaluomeR)

 data("rnaMetrics")

-

 dataFrame <- stability(data=rnaMetrics, k=2, bs=20, getImages = FALSE)

 dataFrame <- stabilityRange(data=rnaMetrics, k.range=c(2,5), bs=20, getImages = FALSE)