r - 2 つの要因による変数の平均

Question

次のデータがあります。

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)

以下を使用して、メソッドごとに各種のサンプル数を取得できます。

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

種ごとにメソッドごとに平均長を取得するにはどうすればよいですか?