Plot usando as estatísticas de um arquivo com o mesmo valor no intervalo x

1

Eu gostaria de traçar médias e intervalo de confiança com o Gnuplot. Meus dados tem várias amostras com o mesmo valor x

21  27  4   12  20  100 50000   false   0.1 "DFSA"  2   205068  31.986
24  27  4   12  20  100 50000   false   0.1 "DFSA"  2   207798  32.49448
22  27  4   12  20  100 50000   false   0.1 "DFSA"  2   207499  32.60746
23  27  4   12  20  100 50000   false   0.1 "DFSA"  2   214065  33.17844
25  27  4   12  20  100 50000   false   0.1 "DFSA"  2   198386  31.4488
26  27  4   12  20  100 50000   false   0.1 "DFSA"  2   208307  32.77412
28  27  4   12  20  100 50000   false   0.1 "DFSA"  2   208669  32.77474
27  27  4   12  20  100 50000   false   0.1 "DFSA"  2   214324  33.16506
30  27  4   12  20  100 50000   false   0.1 "DFSA"  2   209924  32.77186
29  27  4   12  20  100 50000   false   0.1 "DFSA"  2   212959  33.3218
32  27  4   12  20  100 50000   false   0.2 "DFSA"  2   486241  64.3479
31  27  4   12  20  100 50000   false   0.2 "DFSA"  2   487487  64.65076
33  27  4   12  20  100 50000   false   0.2 "DFSA"  2   494703  65.06718
34  27  4   12  20  100 50000   false   0.2 "DFSA"  2   488164  64.77828
36  27  4   12  20  100 50000   false   0.2 "DFSA"  2   476513  63.3158
35  27  4   12  20  100 50000   false   0.2 "DFSA"  2   491005  65.02426
38  27  4   12  20  100 50000   false   0.2 "DFSA"  2   487454  64.44962
37  27  4   12  20  100 50000   false   0.2 "DFSA"  2   490494  65.06572
39  27  4   12  20  100 50000   false   0.2 "DFSA"  2   472081  63.31234
40  27  4   12  20  100 50000   false   0.2 "DFSA"  2   498294  65.02114

Tem dez amostras com coluna 9 = 0.1 e o mesmo com valor 0.2. Eu quero plotar a média da coluna 13 (última coluna) com o intervalo de confiança e usar o Gnuplot yerrorlines , que precisa de dados em quatro colunas x mean min max . A fórmula dos valores min e max usa os STATS_mean, STATS_ssd e STATS_records da coluna 13 com valor igual na coluna 9.

Eu tentei gravar os valores x mean min max em um arquivo com este

system("touch lbt.dat")
stats "data.txt" u ($9==0.1?$13:1/0) name "lbt01"
system("(0.1 lbt01_mean (lbt01_mean - 2.262*lbt01_ssd/sqrt(lbt_records)) (lbt01_mean + 2.262*lbt01_ssd/sqrt(lbt_records))) >> lbt.dat")

mas recebo erro no último comando system para salvar os valores no arquivo lbt.dat. Eu também tentei

system("printf '%f\t\%f\t\%f\t\%f' 0.1 lbt01_mean (lbt01_mean - 2.262*lbt01_ssd/sqrt(lbt_records)) (lbt01_mean + 2.262*lbt01_ssd/sqrt(lbt_records))) >> lbt.dat ")

(2,262 é o valor do t-student pdf para intervalo de confiança de 95% e 9 graus de liberdade).

Parece que os resultados de stasts , por ex. lbt01_mean não pode ser gravado em um arquivo usando system .

Além disso, como o arquivo de dados tem amostras de 0.1 a 1.7 em passos de 0.1 (coluna 9), gostaria de perguntar se existe uma maneira de resumir o script do Gnuplot.

Atenciosamente

    
por Pablo 05.02.2018 / 16:44

1 resposta

3

Uma maneira fácil de imprimir em um arquivo do gnuplot é com o comando print , depois de definir o arquivo com set print . Experimente

set print "lbt.dat"
print 0.1, lbt01_mean, \
  (lbt01_mean - 2.262*lbt01_ssd/sqrt(lbt01_records)),\
  (lbt01_mean + 2.262*lbt01_ssd/sqrt(lbt01_records))

Para percorrer todos os valores de 0,1 a 1,7, você pode usar um loop for de 1 a 17 e dividir por 10.

set print "lbt.dat"
do for [i=1:17] {
  v = i/10.
  stats "data.txt" u ($9==v?$13:1/0) name "A" nooutput
  print v, A_mean,   (A_mean - 2.262*A_ssd/sqrt(A_records)),\
    (A_mean + 2.262*A_ssd/sqrt(A_records))
}
    
por 05.02.2018 / 18:20

Tags