Frank Johannes wrote:
> HI,
> Suppose I have the following data structure.
> LRT tp
> 1 1.50654010 522
> 2 0.51793929 522
> 3 0.90340299 522
> 4 1.20293325 522
> 5 1.05578774 523
> 6 0.01617942 523
> 7 0.68183543 523
> 8 0.43820244 523
> 9 1.14123995 524
> 10 0.05809550 524
> 11 0.93061597 524
> 12 1.39739700 524
> 13 1.05220953 525
> 14 0.03471461 525
> 15 0.63168798 525
> 16 1.40592603 525
> 17 1.41884492 526
> 18 0.23388479 526
> 19 0.21881064 526
> 20 0.99710830 526
> 21 2.02054187 527
> 22 1.99872887 527
> 23 1.04187450 527
> 24 1.31556807 527
> 25 2.77775190 528
> 26 2.94778561 528
> 27 1.88800177 528
> 28 2.08249941 528
>
>
> I have succesfully used a command line such as the one below to get
> maxima for each "tp-category'
>
> data.out<-data[tapply(LRT,tp, function(x) which(LRT==max(x))),]
>
> However, when I try it on the above data, it gives me the following
> error message:
>
>>Error in "[.data.frame"(data, tapply(LRT, tp, function(x)
which(LRT == :
>
> invalid subscript type
Works for me. Look at your data structures and check whether your data
frame is OK.
Or much better easier:
tapply(LRT, tp, max)
Uwe Ligges
> I don't know what to do.
> Thanks for your help
>
> --
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html