Result show the values of fitting gamma parameter

Previous Topic Next Topic
 
classic Classic list List threaded Threaded
3 messages Options
Reply | Threaded
Open this post in threaded view
|

Result show the values of fitting gamma parameter

R help mailing list-2
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]]

______________________________________________
[hidden email] 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.
Reply | Threaded
Open this post in threaded view
|

Re: Result show the values of fitting gamma parameter

Eric Berger
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 <
[hidden email]> 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]]
>
> ______________________________________________
> [hidden email] 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]]

______________________________________________
[hidden email] 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.
Reply | Threaded
Open this post in threaded view
|

Re: Result show the values of fitting gamma parameter

Rui Barradas
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 <
> [hidden email]> 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]]
>>
>> ______________________________________________
>> [hidden email] 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]]
>
> ______________________________________________
> [hidden email] 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.
>

______________________________________________
[hidden email] 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.