Este metodo es útil también para definir estadìos de desarrollo, basados por ejemplo en medidas de distancia.
Obviamente esta función también ha sido escrita para R y es realmente muy sencilla, así que manos a la obra:
>estadios<-c(18.30,18.70,19.
#Escribimos un vector con la lista de medidas que tenemos, en este caso se utilizan las medidas de longitudes estándar de individuos de una especie de Pez de la familia Loricariidae, colectados en un Rio.
> nclass.Sturges(estadios)
#Obtenemos el Nùmero de clases en los que se debe dividir la muestra, para que sea consistente, en este caso serà 9.
[1] 9
> summary (estadios)
#Obtenemos los estadísticos descriptivos básicos que serán utilizados para hallar la amplitud de rango de cada estadìo o clase.
Min. 1st Qu. Median Mean 3rd Qu. Max.
14.23 26.72 45.36 49.56 71.23 110.10
> RangoTotal<-(110.10-14.23)
#Obtenemos el rango total que tiene muestra muestra.
> RangoTotal
[1] 95.87
> amplitudtotal<-(RangoTotal+1)
#Sumamos 1 al rango total, para obtener la amplitud de rango total, necesaria para calcular la amplitud de rango para nuestras clases.
> amplitudtotal
[1] 96.87
> amplitudderango<-(
#Hallamos la amplitud de rango de nuestras clases o estadìos.
> amplitudderango
[1] 10.76333
Y listo!!!! de este modo, tendremos en este ejemplo 9 clases o estadios de amplitud de rango de 10.76333
así:
Estadio 1: 14.23 a 24.99
Estadio 2: 25.00 a 35.75
Estadio 3: 35.76 a 46.51
Estadio 4: 46.52 a 57.27
Estadio 5: 57.28 a 68.03
Estadio 6: 68.04 a 78.79
Estadio 7: 78.80 a 89.55
Estadio 8: 89.56 a 100.31
Estadio 9: 100.31 a 111.07
Que les sirva y que les rinda!!!!
Algo más automatico
ResponderEliminarnclass_e <- nclass.Sturges(estadios)
range_e <- range(estadios)
rango <- diff(range_e)
amplitud_total <- rango + 1
amplitud_rango <- amplitud_total/nclass_e
cortes <- cumsum(c(min(estadios), rep(amplitud_rango, 9)))
Este comentario ha sido eliminado por el autor.
ResponderEliminarHola Durlak o Luciano..
ResponderEliminarMuchas gracias por comentar.
Definitivamente si es mucho mas automático, (aunque permanece con el mismo numero de comandos, estos son mas sintetizados).
Casualmente, Daniel Miranda también estaba preocupado por el calculo de la amplitud de rango y configuro una pequeñisima función en Rstudio (ya que esta pegadisimo y engomado con Rstudio) que decia lo siguiente:
>amplitud_rango <- function (estadios) {amprango<-((diff(range(estadios))+1)/(nclass_e))}
He decidio unir los apoortes de los dos y formular la siguiente linea de comandos, en orden de realizar mucho mas rápido los cálculos:
> range(estadios)
> nclass_e <- nclass.Sturges(estadios)
> amplitud_rango <- function (estadios) {amprango<-((diff(range(estadios))+1)/(nclass_e))}
>valor<-amplitud_rango(estadios)
>cortes <- cumsum(c(min(estadios), rep(valor, nclass_e)))
Obviamente la función podría ser universal (y cambiar los rótulos de "estadios" por nombres fijos (como "X") del vector que entremos, de este modo aplicando la linea de comandos nos serviría para cualquier set de datos que incluyamos)
Muchas gracias a los dos por comentar!!
Hasta la Vista....
Iron Clipper
ResponderEliminarIron Clipper. ford ecosport titanium Iron Clipper is an original glass-cast stainless steel and aluminum frame does titanium have nickel in it with two sided blades titanium tube which will provide a titanium trimmer as seen on tv convenient joico titanium blade.