jueves, 26 de septiembre de 2013

Hoy tratarémos de hacer plots en julia. Existen varios paquetes para hacer plots, pero esta muy temprano en su desarrollo por el momento.  Primero veremos a Winston, que implemanta plots en el estili de octave (o Matlab...).

julia> using Winston
Warning: using Winston.plot in module Main conflicts with an existing identifier.

julia> x = [-pi:pi/60:pi];

julia> length(x)
121

julia> y=sin(x);

Terminar la line de comando con ;  inhibe que se imprime todo el resultado (como en octave). Pero en julia no es necesario usar ; en las definiciones de funciones, solo en uso interactivo.

julia> p = FramedPlot()
FramedPlot(...)

julia> add(p, Curve(x,y))

julia> setattr(p, "xlabel", "x")
"x"

julia> setattr(p, "ylabel", "\\sin(x)")
"\\sin(x)"

julia> setattr(p, "title", "curva de sinus")
"curva de sinus"

julia> file(p,"plot_ejemplo_1.png")


Por el momento parece que no hay manera de construir el plot interactivamente, es necesrio escribirlo a un fichero como hemos hecho arriba.

Para uso estadístico, sin embargo, el paquete Gadfly parece más interesante, como implementa el "Grammar of Graphics". Pero por el momento solo me da mensajes de error ...



No hay comentarios:

Publicar un comentario