Mam następujące dane:

a <- c(1,1,1,1,2,2,2,2)
b <- c(2,4,6,8,2,3,4,1)
c <- factor(c("A","B","A","B","A","B","A","B"))
df <- data.frame(
    sp=a,
    length=b,
    method=c)

Mogę użyć poniższej metody, aby uzyskać liczbę próbek każdego gatunku według metody:

n <- with(df,tapply(sp,method,function(x) count(x)))

Jak uzyskać średnią długość według metody dla każdego gatunku?