intervalo.plot<- function(li, ls, parametro,semilla) {
  y1 <- li ;y2<-ls
  n <- length(y1)
  no.condicion <-  (li >parametro) + (ls <parametro)
  suma<-sum(no.condicion)
  plot(y1, type = "n", ylim=c(min(li),max(ls)),xlab = " ",ylab = " ",
  main=paste(c("semilla=","nš intervalos="),c(semilla,suma)))
  segments((1:n)[y1<parametro & y2 >parametro], y1[y1<parametro & y2 >parametro],
  (1:n)[y1<parametro & y2 >parametro], y2[y1<parametro & y2 >parametro])
  if(sum(y1 > parametro)>0) segments((1:n)[y1>parametro], y1[y1 >parametro],(1:n)[y1>parametro],
   y2[y1>parametro],col=3,lwd=6)
  if(sum(y2 < parametro)>0) segments((1:n)[y2<parametro], y1[y2 <parametro],(1:n)[y2<parametro],
  y2[y2<parametro],col=6, lwd=6)
  abline(h=parametro)
  list("intervalos que no contienen al parametro"=suma)
  }