Capture the results of the apply command into an object and then work with
that. Here is one way to do it:
> res <- apply(C, 2, fitdist, "gamma")
> out <- c( res$A$estimate["shape"], res$B$estimate["shape"],
res$A$estimate["rate"], res$B$estimate["rate"])
> names(out) <- c("A shape","B shape","A rate","B Rate")
> print(out)
# A shape B shape A rate B Rate
# 3.702253 31.300800 1.234126 3.912649
HTH,
Eric
On Mon, Jan 29, 2018 at 10:25 AM, smart hendsome via R-help <
r-help@r-project.org> wrote:
> Hi,
> Let say I have data by two columns A and B, and I have fit each column
> using the gamma distribution by 'fitdist' . I just want the result show
> only the shape and rate only.
>
> Eg:
> library(fitdistrplus)
>
> A <-c(1,2,3,4,5)
>
> B<-c(6,7,8,9,10)
>
> C <-cbind(A,B)
> apply(C, 2, fitdist, "gamma")
> Output show like this:
> $A
> Fitting of the distribution ' gamma ' by maximum likelihood
> Parameters:
> estimate Std. Error
> shape 3.702253 2.2440052
> rate 1.234126 0.8011369
>
> $B
> Fitting of the distribution ' gamma ' by maximum likelihood
> Parameters:
> estimate Std. Error
> shape 31.300800 19.69176
> rate 3.912649 2.48129
>
> I want the output to be like this:
> A B
> shape 3.702253 31.300800rate 1.234126 3.912649
> Can anyone solve my problem? Many thanks.
>
> Regards,
> Zuhri
>
>
>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/
> posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
