5
を使用して新しいデータフレーム内の複数の出力を強制:は、私はこの機能を持っているddply
> λ.est <- function(x){
mle.optim <- mle2(paretoNLL,start=list(λ=-0.7),data=list(x=x),trace=TRUE)
return(summary(mle.optim)@coef[1,1:4])
}
分布をフィットし、パラメータ推定値、STDをretunsことは。エラー、z値、およびp。 私はddply機能を使用している要因pond,habitat,treatment,date
の組み合わせによって定義された私の元のデータフレームsize
の異なるサブセットにこの機能を適用する必要があり、これを行うには:
> mle.λ <- ddply(size, .(pond,habitat,treatment,date),
summarise, λ=λ.est(x=mass.wei))
問題が行うことによって、ということです
> mle.λ
pond habitat treatment date estimate std. error z value Pr(z)
- - - - - - - -
- - - - - - - -
- - - - - - - -
- - - - - - - -
- - - - - - - -
...
:これは、私が唯一の新しいデータフレーム
mle.λ
に1列の時間を追加することができ、私は次のようになり
mle.λ
4新しい列、
λ.est
基本的には何かの出力ごとに1つに追加する必要がありwereas
これまでのところ私必要な出力ごとに異なる関数を書いていますが、私はいくつかのコード経済をしたいと思います。
おかげ マッテオ
を、私も...まだに多くのことを行進しなければならないよりもはるかに簡単いつものように:何の例のデータがないので、ここではこれを行う方法を明確に示さなければならないものですRルート.. – matteo
素晴らしい!私はそれが助けてうれしいです。 –