Tengo los siguientes datos:

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)

Puedo usar lo siguiente para obtener un recuento del número de muestras de cada especie por método:

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

¿Cómo obtengo también la longitud media por método para cada especie?