Rコード:データフレームのグループごと（セルごと）の平均を出す。

Tweet

グループごとの平均値などをまとめて算出する方法。

Rを用いて統計解析をする際に、データフレームに格納されているデータの各グループごとに平均値や合計などを計算したいことがあります。

例えば以下のデータフレームにある会社の営業部、サービス部、総務部の従業員の性別と年収が格納されています。

 
このデータフレームは以下のコードで作成することができます。

data=data.frame( 部署=c(rep("営業部",10),
                        rep("サービス部",5),
                        rep("総務部",5)),
                 性別=c(rep("男性",8),rep("女性",2),
                        rep("男性",3),rep("女性",2),
                        rep("男性",1),rep("女性",4)),
                 年収=c(550,600,600,700,800,500,600,750,
                        700,800,
                        500,400,400,
                        350,350,
                        600,
                        350,400,400,350)
                 )

このデータの部署ごと、性別ごと、部署x性別ごとの年収の平均値を出す方法を紹介します。

データフレームの各グループごとの平均値を出したいときは、aggregate()関数を使用します。

部署ごとの平均を出すには以下のコードを両します。

aggregate(年収~部署, # 部署ごとの年収
				data=data, # データフレームdataを利用
				FUN=mean) # 平均値を指定

このコードで以下のように部署ごとの平均年収が出力されます。

 
同じように性別ごとの年収の平均値は以下の通りです。

aggregate(年収~性別,data=data,FUN=mean)

 
次に部署x性別ごとの平均値の出し方です。

aggregate(年収~部署+性別, # 部署x性別ごとの年収を指定
				data=data, # データフレームdataを指定
				FUN=mean) # 平均値を算出

これにより、以下の通り部署x性別ごとの年収の平均値を出すことができます。

 
このaggregate()関数においては平均値だけでなく様々な関数の値を出すことができます。

例えば、各部署x性別ごとの年収の最大値を出したい場合にはFUN=meanの代わりにFUN=maxを指定します。

> aggregate(年収~部署+性別,
					data=data,
					FUN=max # 最大値maxを指定)

Tweet