У меня есть эта функция:
> λ.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. ошибка, значение z и p для моей модели.
Я должен применить эту функцию к различным подмножествам моего исходного фрейма данных , размер
, определенный комбинацией фактора водоем, среда обитания, обработка, дата
, и для этого Я использую функцию ddply:
> mle.λ <- ddply(size, .(pond,habitat,treatment,date),
summarise, λ=λ.est(x=mass.wei))
проблема в том, что, делая это, я могу добавлять только один столбец за раз в новый фрейм данных mle.λ
, хотя мне нужно добавить в mle.λ
четыре новых столбца, по одному для каждого из выходов λ.est
в основном примерно так:
> mle.λ
pond habitat treatment date estimate std. error z value Pr(z)
- - - - - - - -
- - - - - - - -
- - - - - - - -
- - - - - - - -
- - - - - - - -
...
Пока что я писал требуется отдельная функция для каждого вывода, но я хотел бы немного сэкономить на коде ... есть ли способ сделать все за один раз?
спасибо matteo