On 17/02/15 12:59, smart hendsome wrote:
I'm very new to r-programming. I have rainfall data. I have tried to fit gamma
into my data but there is error. Anyone can help me please.
My rainfall data as I uploaded. When I try run the coding:
library(MASS)
KLT1<-read.csv('C:/Users/User/Dropbox/PhD
Materials/Coding_PhD_Thesis/Kelantan_Average/K1.csv')
KLT<-KLT1$Amount
fd_g <- fitdistr(KLT, "gamma")
Error in stats::optim(x = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, :
initial value in 'vmmin' is not finite
Anyone can help me? Thanks.
No. No one can help. Not unless you provided a self-contained
reproducible example as the posting guide requests.
cheers,
Rolf Turner
P. S. I am not sure, but the error may be indicating that your data
begin with a long string of zeros. (In which case fitting a gamma
distribution is simply idiotic.) Rainfall data often contain many
zeroes. (There are, contrary to popular belief, many dry days.)
What you *may* want to fit is a mixture of a gamma distribution and a
point-mass ("atom") at 0.
R. T.
--
Rolf Turner
Technical Editor ANZJS
Department of Statistics
University of Auckland
Phone: +64-9-373-7599 ext. 88276
Home phone: +64-9-480-4619
______________________________________________
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.