# multi03.rb.txt: ruby script (encoding: Windows-31J)
  # z_multi03.csv を出力
  # 複数回答だけに着目して、その相互の選択状況をみる
  # data01.csv: 「ID,年齢層,新聞,雑誌,テレビ,ラジオ」の98人分のデータ
  # 主な情報源(新聞,雑誌,テレビ,ラジオ)の複数回答の選択状況として、
  # 新聞と雑誌の両方を選択した人が何人いるかなどをみる
  # 各々の選択肢の組合せを一通りみる。
  # 「基本表, 該当数つきの表, 横軸パーセント, 対総数のパーセント」を出力
require "rrx"

    # Rプログラム
rpro = <<EOS
dtf <- read.csv("data01.csv", header=T, na.strings="")
dtf <- ma.refactor(dtf, 3:6, c(1,NA))
xx <- ma2(dtf, 3:6)
xx <- lround(xx, 1)
count <- ma.count(dtf, 3:6, "○")  # 複数回答の選択個数
count <- sprintf("%d個", count)
dtf <- cbind(dtf, 選択個数=count)
zz <- sa1(dtf, "選択個数")
xx$count <- join.tbl(zz$tbl2, zz$pct1, "%d(%.1f)")
robj(xx)
EOS

     # Rプログラムの実行と、その結果の出力
hs = Rrx::rexec(rpro)
xx = hs["robj1"]
csv_str = "※ 主な情報源(複数回答)に関する相互の選択状況をみる\n"
csv_str += "新聞・雑誌・テレビ・ラジオは複数選択可\n"
csv_str += "新聞と雑誌の両方を選択したのは何人かなどをみる。\n\n"
count = xx.delete("count")  # 複数回答の選択個数の状況を得る
count[0] = [nil, "選択個数", "該当者数"]
count = Rrx::table_turn(count)[1..-1]
csv_str += "** 複数回答の選択個数の分布\n"
csv_str = csv_str + Rrx::ary2str(count, ",") + "\n\n"
keys = %w(tbl1 tbl2 pct1 pct2)
titles = %w(基本表 「該当数」付きの表 横軸パーセント 対総数のパーセント)
while key = keys.shift
  title = titles.shift
  val = xx[key]
  next  unless val
  csv_str = csv_str + "** #{title}\n"
  csv_str = csv_str + Rrx::ary2str(val, ",") + "\n\n"
end
File.open("z_multi03.csv", "w") {|ff|  ff.write csv_str}

__END__
○ ma2    tjsf.rで定義されているR用の関数
・利用例：xx <- ma2(dtf, 3:6, 1, FALSE)
・機能：複数回答の選択肢相互の関係をみるための表を生成する。
・戻り値：tbl1, tbl2, pct1, pct2 の４つの表からなるリストを返す。
・引数：第１引数はデータフレーム、第２引数は複数回答の列の指定。