Sha256: db2411b7805a010432d71964175eb1c47fcb7efb310c2f130061232fd2dd03fb
Contents?: true
Size: 969 Bytes
Versions: 21
Compression:
Stored size: 969 Bytes
Contents
#!/usr/bin/env ruby
require 'rbbt/util/R'
$0 = "rbbt #{$previous_commands*""} #{ File.basename(__FILE__) }" if $previous_commands
options = SOPT.setup <<EOF
Calculate density
$ rbbt stat density <file> <plot>
-h--help Help
-m--max* Upper limit
EOF
SOPT.usage if options[:help]
file = ARGV.shift
plot = ARGV.shift
file, plot = ["-", file] if plot.nil? and file.include? '.png'
if file == '-' or file.nil?
file = STDIN
else
file = Open.open file
end
values = file.read.split("\n").collect{|v| v.to_f}
max = options[:max]
res = TmpFile.with_file nil, :extension => 'png' do |tmp|
if max
R.run <<-EOF
values = #{R.ruby2R values}
d = density(values)
rbbt.png_plot(filename = '#{plot || tmp}', 'plot(d, xlim=c(0,#{max})); abline(v=mean(values))')
EOF
else
R.run <<-EOF
values = #{R.ruby2R values}
d = density(values)
rbbt.png_plot(filename = '#{plot || tmp}', 'plot(d); abline(v=mean(values))')
EOF
end
`op '#{plot || tmp}'`
end
Version data entries
21 entries across 21 versions & 1 rubygems