You want:
tapply( Outcome, predictor, mean )
Bendix Carstensen
--
Bendix Carstensen
Senior Statistician
Steno Diabetes Center
Niels Steensens Vej 2
DK-2820 Gentofte
Denmark
tel: +45 44 43 87 38
mob: +45 30 75 87 38
fax: +45 44 43 07 06
[EMAIL PROTECTED]
www.biostat.ku.dk/~bxc
You can use 'cut' to create the breaks.. Actually there are 8 in the 3-4
range:
Outcome predictor
10 1
21 2
31 2
40 3
50 3
60 2
71 3
81 4
91 4
10
I was doing something very similar.
I found it tricky to work out how to find a confidence interval for the
'percentage' of the outcome (I call it proportion).
Some of my bins had all zeroes or all ones, so I couldn't work out how to
make a variance that was sensible. Also some bins had few