Hello,
I believe the following is simpler.
It changes the OP's code a bit and uses lapply, not apply.
res2 <- lapply(C, fitdist, "gamma")
do.call(rbind, lapply(res2, `[[`, "estimate"))
# shape rate
#A 3.702253 1.234126
#B 31.300800 3.912649
Hope this helps,
Rui Barradas
On 1/29/2018 10:25 AM, Eric Berger wrote:
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.
______________________________________________
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.
