Here's a simple example, it might be tweakable to get want you want:

def stdplot(somedata):
    means = [N(mean(row)) for row in somedata]
    stds = [N(std(row)) for row in somedata]
    outplot = Graphics()
    for i in range(len(means)):
        outplot += point([i,means[i]])
        outplot += line([[i,means[i]-stds[i]],[i,means[i]+stds[i]]])
        outplot += line([[i-.4,means[i]-stds[i]],[i+.4,means[i]-
stds[i]]])
        outplot += line([[i-.4,means[i]+stds[i]],[i+.4,means[i]
+stds[i]]])
    return outplot


rdata = [[j + j*random()/2 for i in range(10)] for j in range(20)]
stdplot(rdata)

-Marshall Hampton

On Jun 10, 4:40 am, Marco Boretto <marco.b...@gmail.com> wrote:
>  Hello to everybody, i'm try to use sage in my physics elaboration,
> but i could not make plot point with an error bar for each point..
> there is the possibility to develop this function?
> Marco

-- 
To post to this group, send an email to sage-devel@googlegroups.com
To unsubscribe from this group, send an email to 
sage-devel+unsubscr...@googlegroups.com
For more options, visit this group at http://groups.google.com/group/sage-devel
URL: http://www.sagemath.org

Reply via email to