3 
df1=data.frame(c(2,1,2),c(1,2,3,4,5,6),seq(141,170)) #create data.frame
names(df1) = c("gender","age","height") #column names
df1$gender <- factor(df1$gender,
levels=c(1,2),
labels=c("female","male")) #gives levels and labels to gender
df1$age <- factor(df1$age,
levels=c(1,2,3,4,5,6),
labels=c("16-24","25-34","35-44","45-54","55-64","65+")) # gives levels and labels to age groups

性別と年齢でサブセット化された高さの値の要約を作成しようとしています。

subsetと関数を as として使用するとby、必要な出力が得られます。

females<-subset(df1,df1$gender==1) #subsetting by gender
males<-subset(df1,df1$gender==2)

foutput=by(females$height,females$age,summary) #producing summary subsetted by age
moutput=by(males$height,males$age,summary)

ただし、XLconnect を使用してこれらの結果を度数分布表と一緒にエクスポートできるように、data.frame に入れる必要があります。

おそらくplyrを使用して、出力をdata.frameまたはエレガントな代替手段に変換する方法はありますか?

4

2 に答える 2

4

plyrを使用する1つのアプローチは次のとおりです。

> ddply(df1, c("gender", "age"), function(x) summary(x$height))
  gender   age Min. 1st Qu. Median Mean 3rd Qu. Max.
1 female 25-34  142     148    154  154     160  166
2 female 55-64  145     151    157  157     163  169
3   male 16-24  141     147    153  153     159  165
4   male 35-44  143     149    155  155     161  167
5   male 45-54  144     150    156  156     162  168
6   male   65+  146     152    158  158     164  170
于 2012-02-23T12:31:56.883 に答える
2

からの出力byは実際には ですが、メソッドlistのために異なって見えます。print.by

したがって、要素を a に使用do.callして、それを呼び出すことができます。rbindmatrixdata.frame

data.frame(do.call(rbind,by(mtcars$hp,mtcars$cyl,summary)),check.names=FALSE)
  Min. 1st Qu. Median   Mean 3rd Qu. Max.
4   52    65.5   91.0  82.64    96.0  113
6  105   110.0  110.0 122.30   123.0  175
8  150   176.2  192.5 209.20   241.2  335

check.names列名のサニタイズを避けるために引数を使用することに注意してください。

于 2012-02-23T12:59:25.450 に答える