Hi R-users,
I'm not sure what's wrong and how do I check the error?
I have 3 sets of data
> aug_dt1
[1] 22.8 49.6 77.9 67.6 91.2 48.9 124.6 54.6 85.8 62.8 63.3 174.7
32.3 22.7 99.9 123.2 4.8 149.0 66.0
[20] 73.9 37.4 22.4 69.5 87.3 66.9 87.2 81.5 109.1 47.4 22.8 92.4
60.9 77.5 91.3 91.9 78.2 39.4 107.5
[39] 57.9 65.5 76.0 47.2 145.5 95.9 106.0 120.4 108.0 96.8 30.2 16.6
54.0 91.0 63.0 120.6 10.0 144.2 126.4
[58] 160.6
> dec_dt1
[1] 3.3 44.2 157.5 1.0 22.1 179.7 42.0 37.1 184.2 31.5 1.5 3.8
54.7 19.3 19.9 0.5 45.5 15.8 42.7
[20] 104.7 52.6 15.4 67.4 21.7 66.5 11.8 124.4 39.3 17.3 77.5 11.4
18.0 68.2 77.0 41.3 16.3 61.6 24.6
[39] 197.3 4.6 62.1 24.4 37.9 97.7 1.3 69.8 13.9 98.4 30.6 8.6
43.4 1.0 102.0 35.0 21.8 55.4 18.2
[58] 67.6
> jan_dt1
[1] 49.7 21.8 5.6 24.2 0.8 32.0 88.9 61.6 36.8 74.0 16.3 2.6
7.6 149.7 13.0 65.2 3.8 72.1 49.0
[20] 3.8 24.1 4.3 19.6 32.3 5.6 54.6 94.2 40.7 101.3 0.5 11.7
9.9 86.4 24.3 91.3 66.5 13.2 0.8
[39] 16.0 24.7 18.1 13.5 126.5 18.1 62.5 98.9 212.6 19.2 42.4 41.6
73.2 8.6 15.0 43.2 43.0 1.4 167.4
[58] 11.4
This is my code to estimate alpha and beta:
library(MASS)
alpha.beta <- function(dt)
{ tol <- 1E-6
fit <- list()
x <- dt
a <- x [x > tol]
fit <- fitdistr(a,dgamma, list(shape = 1, rate = 0.01), lower = 0.01)
fit
}
When I run I got this message:
> alpha.beta(dec_dt1)
Error in fitdistr(a, dgamma, list(shape = 1, rate = 0.01), lower = 0.01) :
optimization failed
> alpha.beta(aug_dt1)
Error in fitdistr(a, dgamma, list(shape = 1, rate = 0.01), lower = 0.01) :
optimization failed
but it works for january data
> alpha.beta(jan_dt1)
shape rate
0.89517560 0.02063325
(0.14460597) (0.00437990)
Thank you so much for any help given.
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
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.
