Aquí les dejo el programita en C y el ejecutable del juego del caos.
Como siempre, bajan todos los archivos y los ponen en un mismo directorio. El archivo b.seq es la secuencia del cromosoma 30 del caballo domestico, el archivo jdc2.c es el código fuente.
Teniendo un compilador de C instalado (GCC para Unix, o DJGPP para Windows), compilan con:
gcc -lm jdc2.c
Tras correr el ejecutable anteriormente generado, se genera un archivo puntos.dat que contiene las coordenadas de los puntos. Pueden gráficarlos con gnuplot usando el siguiente comando (dentro de gnuplot): plot "puntos.dat" using 1:2 with dots 1 14
La gráfica la pueden hacer con cualquier otra cosa, incluso en excel.
Porfavor experimenten con otras secuencias, pueden bajar de esta página, o este archivo.
Archivos:
Codigo Fuente.
Secuencia caballo.
Ejecutable para Windows 32 bits y para linux 2.x 32 bits.
Un artículo relacionado.
Adicionalmente les aviso que el Proximo VIERNES la clase será en el observatorio virtual Ixtli en DGSCA, donde Omar Suarez nos mostrará gráficas generadas a partir de secuencias de DNA.
Saludos.
2 comentarios:
Saludos y antes que nada, una disculpa para todos aquellos que se quedaron con dudas con respecto a la pequeña exposición que les hice acerca de redundancia cíclica.
Lamentablemente despues ví que tuve un par de errores en cuanto lo expuesto, pero en fin ,les dejo un par de páginas donde pueden ver mejor este tema, si es que se quedaron con dudas todavía.
http://www.emezeta.com/articulos/
codigos-de-redundancia-ciclica
http://www.argo.es/~jcea/artic/
ecc-crc.htm
No los copien y peguen, sino que escribanlos directo, para que les funcione. Tuve que usar un salto de linea para que no se saliese del espacio.
De nuevo, mil disculpas si quedó muy chafa la exposición, pero con esto creo qu ya no quedarán mas dudas.
Saludos de nuevo y nos vemos.
saludos
Baje el programita, y compilo bien pero a la hora de correr el paso de gnuplot obtuve esto:
cheko@isis:~/tests/teoin> plot "puntos.dat" using 1:2 with dots 1 14
#PLOT 1
o000 0$780000
0,780000
0.465000 0.890000
392�|[1��564@0$(D�4�116<�4k<01(D�4�1*6�4�1E063 0$�VL6R(1,�VL6R(15�1��g�R2�e�5�X26E5U{�5�266
0.3 0.038203
0.3768$02
0.1881��/��702
plot: unrecognized command `0xa' encountered in input
plot: could not parse input file `puntos.dat'
0.188408 0.00E$408 \.59)4204(/�0Excheko@isis:~/tests/teoin>
aunque ahora me dispongo a usarlo con las demas secuencias.
Publicar un comentario