Sha256: 781d558d94f2903902a09a4214b0473b8eec87e381842e75cb9cfd81c6040fc5
Contents?: true
Size: 1009 Bytes
Versions: 289
Compression:
Stored size: 1009 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));rug(jitter(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));rug(jitter(values))')
EOF
end
`op '#{plot || tmp}'`
end
Version data entries
289 entries across 289 versions & 1 rubygems