Skip to content
Prev 219300 / 398500 Next

percent by subset

how about:

d=data.frame(ht=rnorm(20,60,20),sex=sample(0:1,20,rep=T)); d
with(d,by(ht,sex,mean))
with(d,by(ht,sex,function(x)mean(x>=60)))

hth, david freedman