quiero hacer un loop en R que como resultado guarde un rango cuyo tamaño minimo sea 1000 e incremente su tamaño 10 puntos cada vez que la condicion no se cumple. Les mando el script que tengo pero no estoy teniendo el resultado que quiero.
set.seed(123)
newbins = seq(0,5000,1000)
binLength = 1000
stepSize = 10
blocks=NULL
end = 0
for (b in newbins[1:length(newbins)-1]){
tempBP = 0
tempBP_new = 0
start = end +1
end = binLength
while (sum(tempBP_new-tempBP) < 0.40){
tempBP_new = tempBP
tempBP_new = runif(1)*0.5
binLength = binLength + stepSize
blocks = rbind(blocks,c(start, end))
}
}
blocks
El output del script es:
[,1] [,2]
[1,] 1 1000
[2,] 1 1000
[3,] 1 1000
[4,] 1 1000
[5,] 1001 1040
[6,] 1041 1050
[7,] 1041 1050
[8,] 1041 1050
[9,] 1051 1080
[10,] 1051 1080
[11,] 1051 1080
[12,] 1081 1110
[13,] 1081 1110
[14,] 1081 1110
[15,] 1081 1110
[16,] 1081 1110
Pero yo quiero algo como por ejemplo
1 1000
1001 2040
2041 3050
3050 5000
O sea cada rango debe ser de al menos 1000 y el loop tiene que terminar cuando se llega al 5000. Tampoco entiendo por qué tengo el mismo resultado varias veces en el loop. Tienen algunas ideas?? Muchas gracias!