In [1]:

    

using DataArrays, DataFrames, Gadfly, Distributions, Grid


    

In [13]:

    

df = DataFrame()
df[:EspesorPlomo] = [0,7,21,26.3,31.6,36.9,42.2,47.5]
df[:Medicion1] = [145,35,8,5,6,3,2,3]
df[:Medicion2] = [145,28,9,5,5,4,2,3]
df[:Medicion3] = [139,29,8,7,4,2,4,2]
df[:Medicion4] = [149,33,9,7,4,3,2,1]
df[:Medicion5] = [140,35,9,7,5,3,1,0]
df[:Promedio] = [(df[:Medicion1][1]+df[:Medicion2][1]+df[:Medicion3][1]+df[:Medicion4][1]+df[:Medicion5][1])/5,
    (df[:Medicion1][2]+df[:Medicion2][2]+df[:Medicion3][2]+df[:Medicion4][2]+df[:Medicion5][2])/5,
    (df[:Medicion1][3]+df[:Medicion2][3]+df[:Medicion3][3]+df[:Medicion4][3]+df[:Medicion5][3])/5,
    (df[:Medicion1][4]+df[:Medicion2][4]+df[:Medicion3][4]+df[:Medicion4][4]+df[:Medicion5][4])/5,
    (df[:Medicion1][5]+df[:Medicion2][5]+df[:Medicion3][5]+df[:Medicion4][5]+df[:Medicion5][5])/5,
    (df[:Medicion1][6]+df[:Medicion2][6]+df[:Medicion3][6]+df[:Medicion4][6]+df[:Medicion5][6])/5,
    (df[:Medicion1][7]+df[:Medicion2][7]+df[:Medicion3][7]+df[:Medicion4][7]+df[:Medicion5][7])/5,
    (df[:Medicion1][8]+df[:Medicion2][8]+df[:Medicion3][8]+df[:Medicion4][8]+df[:Medicion5][8])/5]
df[:Error] = [trunc(std([df[:Medicion1][1],df[:Medicion2][1],df[:Medicion3][1],df[:Medicion4][1],df[:Medicion5][1]]),1),
    trunc(std([df[:Medicion1][2],df[:Medicion2][2],df[:Medicion3][2],df[:Medicion4][2],df[:Medicion5][2]]),1),
    trunc(std([df[:Medicion1][3],df[:Medicion2][3],df[:Medicion3][3],df[:Medicion4][3],df[:Medicion5][3]]),1),
    trunc(std([df[:Medicion1][4],df[:Medicion2][4],df[:Medicion3][4],df[:Medicion4][4],df[:Medicion5][4]]),1),
    trunc(std([df[:Medicion1][5],df[:Medicion2][5],df[:Medicion3][5],df[:Medicion4][5],df[:Medicion5][5]]),1),
    trunc(std([df[:Medicion1][6],df[:Medicion2][6],df[:Medicion3][6],df[:Medicion4][6],df[:Medicion5][6]]),1),
    trunc(std([df[:Medicion1][7],df[:Medicion2][7],df[:Medicion3][7],df[:Medicion4][7],df[:Medicion5][7]]),1),
    trunc(std([df[:Medicion1][8],df[:Medicion2][8],df[:Medicion3][8],df[:Medicion4][8],df[:Medicion5][8]]),1)]
df


    

    
Out[13]:




EspesorPlomo
Medicion1
Medicion2
Medicion3
Medicion4
Medicion5
Promedio
Error
1
0.0
145
145
139
149
140
143.6
4.0
2
7.0
35
28
29
33
35
32.0
3.3
3
21.0
8
9
8
9
9
8.6
0.5
4
26.3
5
5
7
7
7
6.2
1.0
5
31.6
6
5
4
4
5
4.8
0.8
6
36.9
3
4
2
3
3
3.0
0.7
7
42.2
2
2
4
2
1
2.2
1.0
8
47.5
3
3
2
1
0
1.8
1.3

In [14]:

    

ymins = df[:Promedio] .- df[:Error]
ymaxs = df[:Promedio] .+ df[:Error]
plot(layer(df,x="EspesorPlomo",y="Promedio", ymin=ymins, ymax=ymaxs, Geom.point, Geom.errorbar))


    

    
Out[14]:

In [21]:

    

std([3,4,4,2,5])


    

    
Out[21]:





1.140175425099138