[R] error message in gev

2013-07-17 Thread Roslina Zakaria
 
Hi r-users,
 
I would like to use gev and my data (annual rainfall ) is as follows:
 
 head(dat,20)  ABCDEFGHI J
1  45.1 41.5 58.5 50.1 46.0 49.1 37.7 49.1 59.8  54.0
2  50.3 39.8 49.4 56.4 49.4 48.8 42.1 49.8 49.4  58.3
3  41.7 39.3 44.6 39.1 35.7 41.5 40.8 40.8 38.5  45.6
4  50.7 33.9 48.4 28.2 35.5 39.1 61.4 17.0 30.7  38.3
5  39.3 30.6 46.9 23.8 25.8 17.1 37.9 40.1 35.4  30.1
6  50.2 43.6 55.5 26.7 36.6 20.6 35.2 46.0 52.2  30.5
7  46.0 37.0 55.3 36.3 51.9 50.6 43.2 44.9 38.2  45.4
8  46.0 15.4 35.7 27.3 40.2 21.2 41.1 65.7 35.7  28.7
9  52.0 36.5 49.6 55.0 38.3 32.0 43.8 43.8 50.3  36.5
10 57.5 51.0 47.0 63.0 43.0 46.5 41.5 82.8 50.5  38.5
11 71.0 48.0 59.5 48.5 47.0 57.5 55.8 64.5 60.0  43.4
12 57.5 52.5 58.5 42.0 62.0 44.0 47.5 59.5 45.0  51.5
13 47.0 38.0 49.9 50.5 55.0 57.5 51.5 54.5 45.0  49.7
14 66.5 44.4 60.0 99.5 54.5 41.3 57.0 69.5 33.0  58.4
15 61.5 43.6 68.5 50.7 48.0 50.5 50.5 58.8 53.2  47.0
16 55.5 58.5 55.5 45.8 48.0 57.0 52.0 48.5 58.5 437.5
17 96.7 96.7 96.7 45.5 65.1 46.5 51.0 96.7 61.3 310.5
18 71.5 55.5 63.0 52.5 59.0 55.5 55.5 51.1 47.0  49.5
19 71.5 42.5 64.0 56.5 61.0 54.0 46.5 53.7 43.0  44.0
20 57.5 51.5 51.3 34.0 55.0 61.0 55.5 51.3 51.3  44.5
 
When I apply the gev function, I have the following error message:
 
out - gev(dat, 20) Error in gev(dat, 20) : (list) object cannot be coerced to 
type 'double'
 
Can anybody explain to me what does it means.  I tried to read the package 
manual but could not understand about the block.
 
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.


Re: [R] error message in gev

2013-07-17 Thread Rui Barradas

Hello,

You should say which package is gev() coming from. I believe it's from 
evir, and what follows assumes that.
The problem is that 'dat' is a data.frame, a special type of list. It is 
a list of vectors and gev is expecting just one vector. You could try 
instead, to fit generalized extreme value distributions to each of the 
vectors in the data.frame,


out.list - lapply(dat, gev, block = 20)


to get a list of 10 objects of class gev, like the help page says.

Hope this helps,

Rui Barradas

Em 17-07-2013 10:33, Roslina Zakaria escreveu:


Hi r-users,

I would like to use gev and my data (annual rainfall ) is as follows:


head(dat,20)  ABCDEFGHI J

1  45.1 41.5 58.5 50.1 46.0 49.1 37.7 49.1 59.8  54.0
2  50.3 39.8 49.4 56.4 49.4 48.8 42.1 49.8 49.4  58.3
3  41.7 39.3 44.6 39.1 35.7 41.5 40.8 40.8 38.5  45.6
4  50.7 33.9 48.4 28.2 35.5 39.1 61.4 17.0 30.7  38.3
5  39.3 30.6 46.9 23.8 25.8 17.1 37.9 40.1 35.4  30.1
6  50.2 43.6 55.5 26.7 36.6 20.6 35.2 46.0 52.2  30.5
7  46.0 37.0 55.3 36.3 51.9 50.6 43.2 44.9 38.2  45.4
8  46.0 15.4 35.7 27.3 40.2 21.2 41.1 65.7 35.7  28.7
9  52.0 36.5 49.6 55.0 38.3 32.0 43.8 43.8 50.3  36.5
10 57.5 51.0 47.0 63.0 43.0 46.5 41.5 82.8 50.5  38.5
11 71.0 48.0 59.5 48.5 47.0 57.5 55.8 64.5 60.0  43.4
12 57.5 52.5 58.5 42.0 62.0 44.0 47.5 59.5 45.0  51.5
13 47.0 38.0 49.9 50.5 55.0 57.5 51.5 54.5 45.0  49.7
14 66.5 44.4 60.0 99.5 54.5 41.3 57.0 69.5 33.0  58.4
15 61.5 43.6 68.5 50.7 48.0 50.5 50.5 58.8 53.2  47.0
16 55.5 58.5 55.5 45.8 48.0 57.0 52.0 48.5 58.5 437.5
17 96.7 96.7 96.7 45.5 65.1 46.5 51.0 96.7 61.3 310.5
18 71.5 55.5 63.0 52.5 59.0 55.5 55.5 51.1 47.0  49.5
19 71.5 42.5 64.0 56.5 61.0 54.0 46.5 53.7 43.0  44.0
20 57.5 51.5 51.3 34.0 55.0 61.0 55.5 51.3 51.3  44.5

When I apply the gev function, I have the following error message:

out - gev(dat, 20) Error in gev(dat, 20) : (list) object cannot be coerced to 
type 'double'

Can anybody explain to me what does it means.  I tried to read the package 
manual but could not understand about the block.

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.



__
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.


Re: [R] error message in gev

2013-07-17 Thread Rui Barradas

Hello,

It's better if you Cc the list, the odds of getting more and better 
answers are greater.


Have you tried with different block values?

Rui Barradas

Em 17-07-2013 13:33, Roslina Zakaria escreveu:

Hi Rui,
Yes, I am using gev from 'evir'
install.packages(evir)
library(evir)
I have tried your suggestion and follow are what I got.


out.list - lapply(dat, gev, block = 20)

Error in optim(theta, negloglik, hessian = TRUE, ..., tmp = data) :
   non-finite finite-difference value [1]
In addition:Warning messages:
1: In FUN(X[[7L]], ...) : optimization may not have succeeded
2: In FUN(X[[7L]], ...) : optimization may not have succeeded
3: In sqrt(diag(varcov)) : NaNs produced
4: In FUN(X[[7L]], ...) : optimization may not have succeeded
5: In FUN(X[[7L]], ...) : optimization may not have succeeded
6: In FUN(X[[7L]], ...) : optimization may not have succeeded
7: In FUN(X[[7L]], ...) : optimization may not have succeeded



Thank you.


*From:* Rui Barradas ruipbarra...@sapo.pt
*To:* Roslina Zakaria zrosl...@yahoo.com
*Cc:* r-help@r-project.org r-help@r-project.org
*Sent:* Wednesday, July 17, 2013 5:51 PM
*Subject:* Re: [R] error message in gev

Hello,

You should say which package is gev() coming from. I believe it's from
evir, and what follows assumes that.
The problem is that 'dat' is a data.frame, a special type of list. It is
a list of vectors and gev is expecting just one vector. You could try
instead, to fit generalized extreme value distributions to each of the
vectors in the data.frame,

out.list - lapply(dat, gev, block = 20)


to get a list of 10 objects of class gev, like the help page says.

Hope this helps,

Rui Barradas

Em 17-07-2013 10:33, Roslina Zakaria escreveu:
 
  Hi r-users,
 
  I would like to use gev and my data (annual rainfall ) is as follows:
 
  head(dat,20)  ABCDEFGHIJ
  1  45.1 41.5 58.5 50.1 46.0 49.1 37.7 49.1 59.8 54.0
  2  50.3 39.8 49.4 56.4 49.4 48.8 42.1 49.8 49.4  58.3
  3  41.7 39.3 44.6 39.1 35.7 41.5 40.8 40.8 38.5  45.6
  4  50.7 33.9 48.4 28.2 35.5 39.1 61.4 17.0 30.7  38.3
  5  39.3 30.6 46.9 23.8 25.8 17.1 37.9 40.1 35.4  30.1
  6  50.2 43.6 55.5 26.7 36.6 20.6 35.2 46.0 52.2  30.5
  7  46.0 37.0 55.3 36.3 51.9 50.6 43.2 44.9 38.2  45.4
  8  46.0 15.4 35.7 27.3 40.2 21.2 41.1 65.7 35.7  28.7
  9  52.0 36.5 49.6 55.0 38.3 32.0 43.8 43.8 50.3  36.5
  10 57.5 51.0 47.0 63.0 43.0 46.5 41.5 82.8 50.5  38.5
  11 71.0 48.0 59.5 48.5 47.0 57.5 55.8 64.5 60.0  43.4
  12 57.5 52.5 58.5 42.0 62.0 44.0 47.5 59.5 45.0  51.5
  13 47.0 38.0 49.9 50.5 55.0 57.5 51.5 54.5 45.0  49.7
  14 66.5 44.4 60.0 99.5 54.5 41.3 57.0 69.5 33.0  58.4
  15 61.5 43.6 68.5 50.7 48.0 50.5 50.5 58.8 53.2  47.0
  16 55.5 58.5 55.5 45.8 48.0 57.0 52.0 48.5 58.5 437.5
  17 96.7 96.7 96.7 45.5 65.1 46.5 51.0 96.7 61.3 310.5
  18 71.5 55.5 63.0 52.5 59.0 55.5 55.5 51.1 47.0  49.5
  19 71.5 42.5 64.0 56.5 61.0 54.0 46.5 53.7 43.0  44.0
  20 57.5 51.5 51.3 34.0 55.0 61.0 55.5 51.3 51.3  44.5
 
  When I apply the gev function, I have the following error message:
 
  out - gev(dat, 20) Error in gev(dat, 20) : (list) object cannot be
coerced to type 'double'
 
  Can anybody explain to me what does it means.  I tried to read the
package manual but could not understand about the block.
 
  Thank you so much for any help given.
  [[alternative HTML version deleted]]
 
 
 
  __
  R-help@r-project.org mailto: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.
 




__
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.


Re: [R] error message in gev

2013-07-17 Thread Roslina Zakaria
Hi Rui,
 
I have tried different block values and below are the error message I got:
  out - gev(dat, A) Error in as.POSIXlt.default(attributes(data)$times) : 
do not know how to convert 'attributes(data)$times' to class “POSIXlt”
 
out.list - lapply(dat, gev, block = 10) Error in optim(theta, negloglik, 
hessian = TRUE, ..., tmp = data) :  non-finite finite-difference value [2]
 


 From: Rui Barradas ruipbarra...@sapo.pt

Cc: 'r-help' r-help@r-project.org 
Sent: Wednesday, July 17, 2013 8:49 PM
Subject: Re: [R] error message in gev
  

Hello,

It's better if you Cc the list, the odds of getting more and better 
answers are greater.

Have you tried with different block values?

Rui Barradas

Em 17-07-2013 13:33, Roslina Zakaria escreveu:
 Hi Rui,
 Yes, I am using gev from 'evir'
 install.packages(evir)
 library(evir)
 I have tried your suggestion and follow are what I got.

out.list - lapply(dat, gev, block = 20)
 Error in optim(theta, negloglik, hessian = TRUE, ..., tmp = data) :
    non-finite finite-difference value [1]
 In addition:Warning messages:
 1: In FUN(X[[7L]], ...) : optimization may not have succeeded
 2: In FUN(X[[7L]], ...) : optimization may not have succeeded
 3: In sqrt(diag(varcov)) : NaNs produced
 4: In FUN(X[[7L]], ...) : optimization may not have succeeded
 5: In FUN(X[[7L]], ...) : optimization may not have succeeded
 6: In FUN(X[[7L]], ...) : optimization may not have succeeded
 7: In FUN(X[[7L]], ...) : optimization may not have succeeded



 Thank you.


 *From:* Rui Barradas ruipbarra...@sapo.pt

 *Cc:* r-help@r-project.org r-help@r-project.org
 *Sent:* Wednesday, July 17, 2013 5:51 PM
 *Subject:* Re: [R] error message in gev

 Hello,

 You should say which package is gev() coming from. I believe it's from
 evir, and what follows assumes that.
 The problem is that 'dat' is a data.frame, a special type of list. It is
 a list of vectors and gev is expecting just one vector. You could try
 instead, to fit generalized extreme value distributions to each of the
 vectors in the data.frame,

 out.list - lapply(dat, gev, block = 20)


 to get a list of 10 objects of class gev, like the help page says.

 Hope this helps,

 Rui Barradas

 Em 17-07-2013 10:33, Roslina Zakaria escreveu:
  
   Hi r-users,
  
   I would like to use gev and my data (annual rainfall ) is as follows:
  
   head(dat,20)  A    B    C    D    E    F    G    H    
I    J
   1  45.1 41.5 58.5 50.1 46.0 49.1 37.7 49.1 59.8 54.0
   2  50.3 39.8 49.4 56.4 49.4 48.8 42.1 49.8 49.4  58.3
   3  41.7 39.3 44.6 39.1 35.7 41.5 40.8 40.8 38.5  45.6
   4  50.7 33.9 48.4 28.2 35.5 39.1 61.4 17.0 30.7  38.3
   5  39.3 30.6 46.9 23.8 25.8 17.1 37.9 40.1 35.4  30.1
   6  50.2 43.6 55.5 26.7 36.6 20.6 35.2 46.0 52.2  30.5
   7  46.0 37.0 55.3 36.3 51.9 50.6 43.2 44.9 38.2  45.4
   8  46.0 15.4 35.7 27.3 40.2 21.2 41.1 65.7 35.7  28.7
   9  52.0 36.5 49.6 55.0 38.3 32.0 43.8 43.8 50.3  36.5
   10 57.5 51.0 47.0 63.0 43.0 46.5 41.5 82.8 50.5  38.5
   11 71.0 48.0 59.5 48.5 47.0 57.5 55.8 64.5 60.0  43.4
   12 57.5 52.5 58.5 42.0 62.0 44.0 47.5 59.5 45.0  51.5
   13 47.0 38.0 49.9 50.5 55.0 57.5 51.5 54.5 45.0  49.7
   14 66.5 44.4 60.0 99.5 54.5 41.3 57.0 69.5 33.0  58.4
   15 61.5 43.6 68.5 50.7 48.0 50.5 50.5 58.8 53.2  47.0
   16 55.5 58.5 55.5 45.8 48.0 57.0 52.0 48.5 58.5 437.5
   17 96.7 96.7 96.7 45.5 65.1 46.5 51.0 96.7 61.3 310.5
   18 71.5 55.5 63.0 52.5 59.0 55.5 55.5 51.1 47.0  49.5
   19 71.5 42.5 64.0 56.5 61.0 54.0 46.5 53.7 43.0  44.0
   20 57.5 51.5 51.3 34.0 55.0 61.0 55.5 51.3 51.3  44.5
  
   When I apply the gev function, I have the following error message:
  
   out - gev(dat, 20) Error in gev(dat, 20) : (list) object cannot be
 coerced to type 'double'
  
   Can anybody explain to me what does it means.  I tried to read the
 package manual but could not understand about the block.
  
   Thank you so much for any help given.
       [[alternative HTML version deleted]]
  
  
  
   __
   R-help@r-project.org mailto: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.
  


[[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.


Re: [R] error message in gev

2013-07-17 Thread Rui Barradas

Hello,

And what if you try it without a block value?

out.list - lapply(dat, gev)

Rui Barradas

Em 17-07-2013 14:53, Roslina Zakaria escreveu:

Hi Rui,
I have tried different block values and below are the error message I got:


out - gev(dat, A)

Error in as.POSIXlt.default(attributes(data)$times) :

   do not know how to convert 'attributes(data)$times' to class “POSIXlt”



out.list - lapply(dat, gev, block = 10)
Error in optim(theta, negloglik, hessian = TRUE, ..., tmp = data) :
   non-finite finite-difference value [2]


*From:* Rui Barradas ruipbarra...@sapo.pt
*To:* Roslina Zakaria zrosl...@yahoo.com
*Cc:* 'r-help' r-help@r-project.org
*Sent:* Wednesday, July 17, 2013 8:49 PM
*Subject:* Re: [R] error message in gev

Hello,

It's better if you Cc the list, the odds of getting more and better
answers are greater.

Have you tried with different block values?

Rui Barradas

Em 17-07-2013 13:33, Roslina Zakaria escreveu:
  Hi Rui,
  Yes, I am using gev from 'evir'
  install.packages(evir)
  library(evir)
  I have tried your suggestion and follow are what I got.
 
 out.list - lapply(dat, gev, block = 20)
  Error in optim(theta, negloglik, hessian = TRUE, ..., tmp = data) :
 non-finite finite-difference value [1]
  In addition:Warning messages:
  1: In FUN(X[[7L]], ...) : optimization may not have succeeded
  2: In FUN(X[[7L]], ...) : optimization may not have succeeded
  3: In sqrt(diag(varcov)) : NaNs produced
  4: In FUN(X[[7L]], ...) : optimization may not have succeeded
  5: In FUN(X[[7L]], ...) : optimization may not have succeeded
  6: In FUN(X[[7L]], ...) : optimization may not have succeeded
  7: In FUN(X[[7L]], ...) : optimization may not have succeeded
 
 
 
  Thank you.
 
 
  *From:* Rui Barradas ruipbarra...@sapo.pt mailto:ruipbarra...@sapo.pt
  *To:* Roslina Zakaria zrosl...@yahoo.com mailto:zrosl...@yahoo.com
  *Cc:* r-help@r-project.org mailto:r-help@r-project.org
r-help@r-project.org mailto:r-help@r-project.org
  *Sent:* Wednesday, July 17, 2013 5:51 PM
  *Subject:* Re: [R] error message in gev
 
  Hello,
 
  You should say which package is gev() coming from. I believe it's from
  evir, and what follows assumes that.
  The problem is that 'dat' is a data.frame, a special type of list. It is
  a list of vectors and gev is expecting just one vector. You could try
  instead, to fit generalized extreme value distributions to each of the
  vectors in the data.frame,
 
  out.list - lapply(dat, gev, block = 20)
 
 
  to get a list of 10 objects of class gev, like the help page says.
 
  Hope this helps,
 
  Rui Barradas
 
  Em 17-07-2013 10:33, Roslina Zakaria escreveu:
   
Hi r-users,
   
I would like to use gev and my data (annual rainfall ) is as follows:
   
head(dat,20)  ABCDEFGHI J
1  45.1 41.5 58.5 50.1 46.0 49.1 37.7 49.1 59.8 54.0
2  50.3 39.8 49.4 56.4 49.4 48.8 42.1 49.8 49.4  58.3
3  41.7 39.3 44.6 39.1 35.7 41.5 40.8 40.8 38.5  45.6
4  50.7 33.9 48.4 28.2 35.5 39.1 61.4 17.0 30.7  38.3
5  39.3 30.6 46.9 23.8 25.8 17.1 37.9 40.1 35.4  30.1
6  50.2 43.6 55.5 26.7 36.6 20.6 35.2 46.0 52.2  30.5
7  46.0 37.0 55.3 36.3 51.9 50.6 43.2 44.9 38.2  45.4
8  46.0 15.4 35.7 27.3 40.2 21.2 41.1 65.7 35.7  28.7
9  52.0 36.5 49.6 55.0 38.3 32.0 43.8 43.8 50.3  36.5
10 57.5 51.0 47.0 63.0 43.0 46.5 41.5 82.8 50.5  38.5
11 71.0 48.0 59.5 48.5 47.0 57.5 55.8 64.5 60.0  43.4
12 57.5 52.5 58.5 42.0 62.0 44.0 47.5 59.5 45.0  51.5
13 47.0 38.0 49.9 50.5 55.0 57.5 51.5 54.5 45.0  49.7
14 66.5 44.4 60.0 99.5 54.5 41.3 57.0 69.5 33.0  58.4
15 61.5 43.6 68.5 50.7 48.0 50.5 50.5 58.8 53.2  47.0
16 55.5 58.5 55.5 45.8 48.0 57.0 52.0 48.5 58.5 437.5
17 96.7 96.7 96.7 45.5 65.1 46.5 51.0 96.7 61.3 310.5
18 71.5 55.5 63.0 52.5 59.0 55.5 55.5 51.1 47.0  49.5
19 71.5 42.5 64.0 56.5 61.0 54.0 46.5 53.7 43.0  44.0
20 57.5 51.5 51.3 34.0 55.0 61.0 55.5 51.3 51.3  44.5
   
When I apply the gev function, I have the following error message:
   
out - gev(dat, 20) Error in gev(dat, 20) : (list) object cannot be
  coerced to type 'double'
   
Can anybody explain to me what does it means. I tried to read the
  package manual but could not understand about the block.
   
Thank you so much for any help given.
   [[alternative HTML version deleted]]
   
   
   
__
R-help@r-project.org mailto:R-help@r-project.org
mailto:R-help@r-project.org mailto: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
http://www.r-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
   
 
 




__
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

Re: [R] error message in gev

2013-07-17 Thread Roslina Zakaria
Hi Rui,
 
Yes, it works beatifully this time.
  
Thank you so much for your help.
 


 From: Rui Barradas ruipbarra...@sapo.pt

Cc: R help forum R-help@r-project.org 
Sent: Wednesday, July 17, 2013 11:28 PM
Subject: Re: [R] error message in gev
  

Hello,

And what if you try it without a block value?

out.list - lapply(dat, gev)

Rui Barradas

Em 17-07-2013 14:53, Roslina Zakaria escreveu:
 Hi Rui,
 I have tried different block values and below are the error message I got:

out - gev(dat, A)
 Error in as.POSIXlt.default(attributes(data)$times) :

    do not know how to convert 'attributes(data)$times' to class 
“POSIXlt”



 out.list - lapply(dat, gev, block = 10)
 Error in optim(theta, negloglik, hessian = TRUE, ..., tmp = data) :
    non-finite finite-difference value [2]


 *From:* Rui Barradas ruipbarra...@sapo.pt

 *Cc:* 'r-help' r-help@r-project.org
 *Sent:* Wednesday, July 17, 2013 8:49 PM
 *Subject:* Re: [R] error message in gev

 Hello,

 It's better if you Cc the list, the odds of getting more and better
 answers are greater.

 Have you tried with different block values?

 Rui Barradas

 Em 17-07-2013 13:33, Roslina Zakaria escreveu:
   Hi Rui,
   Yes, I am using gev from 'evir'
   install.packages(evir)
   library(evir)
   I have tried your suggestion and follow are what I got.
  
  out.list - lapply(dat, gev, block = 20)
   Error in optim(theta, negloglik, hessian = TRUE, ..., tmp = data) :
      non-finite finite-difference value [1]
   In addition:Warning messages:
   1: In FUN(X[[7L]], ...) : optimization may not have succeeded
   2: In FUN(X[[7L]], ...) : optimization may not have succeeded
   3: In sqrt(diag(varcov)) : NaNs produced
   4: In FUN(X[[7L]], ...) : optimization may not have succeeded
   5: In FUN(X[[7L]], ...) : optimization may not have succeeded
   6: In FUN(X[[7L]], ...) : optimization may not have succeeded
   7: In FUN(X[[7L]], ...) : optimization may not have succeeded
  
  
  
   Thank you.
  
  
   *From:* Rui Barradas ruipbarra...@sapo.pt mailto:ruipbarra...@sapo.pt

com
   *Cc:* r-help@r-project.org mailto:r-help@r-project.org
 r-help@r-project.org mailto:r-help@r-project.org
   *Sent:* Wednesday, July 17, 2013 5:51 PM
   *Subject:* Re: [R] error message in gev
  
   Hello,
  
   You should say which package is gev() coming from. I believe it's from
   evir, and what follows assumes that.
   The problem is that 'dat' is a data.frame, a special type of list. It is
   a list of vectors and gev is expecting just one vector. You could try
   instead, to fit generalized extreme value distributions to each of the
   vectors in the data.frame,
  
   out.list - lapply(dat, gev, block = 20)
  
  
   to get a list of 10 objects of class gev, like the help page says.
  
   Hope this helps,
  
   Rui Barradas
  
   Em 17-07-2013 10:33, Roslina Zakaria escreveu:
    
     Hi r-users,
    
     I would like to use gev and my data (annual rainfall ) is as follows:
    
     head(dat,20)  A    B    C    D    E    F    G    H  
  I J
     1  45.1 41.5 58.5 50.1 46.0 49.1 37.7 49.1 59.8 54.0
     2  50.3 39.8 49.4 56.4 49.4 48.8 42.1 49.8 49.4  58.3
     3  41.7 39.3 44.6 39.1 35.7 41.5 40.8 40.8 38.5  45.6
     4  50.7 33.9 48.4 28.2 35.5 39.1 61.4 17.0 30.7  38.3
     5  39.3 30.6 46.9 23.8 25.8 17.1 37.9 40.1 35.4  30.1
     6  50.2 43.6 55.5 26.7 36.6 20.6 35.2 46.0 52.2  30.5
     7  46.0 37.0 55.3 36.3 51.9 50.6 43.2 44.9 38.2  45.4
     8  46.0 15.4 35.7 27.3 40.2 21.2 41.1 65.7 35.7  28.7
     9  52.0 36.5 49.6 55.0 38.3 32.0 43.8 43.8 50.3  36.5
     10 57.5 51.0 47.0 63.0 43.0 46.5 41.5 82.8 50.5  38.5
     11 71.0 48.0 59.5 48.5 47.0 57.5 55.8 64.5 60.0  43.4
     12 57.5 52.5 58.5 42.0 62.0 44.0 47.5 59.5 45.0  51.5
     13 47.0 38.0 49.9 50.5 55.0 57.5 51.5 54.5 45.0  49.7
     14 66.5 44.4 60.0 99.5 54.5 41.3 57.0 69.5 33.0  58.4
     15 61.5 43.6 68.5 50.7 48.0 50.5 50.5 58.8 53.2  47.0
     16 55.5 58.5 55.5 45.8 48.0 57.0 52.0 48.5 58.5 437.5
     17 96.7 96.7 96.7 45.5 65.1 46.5 51.0 96.7 61.3 310.5
     18 71.5 55.5 63.0 52.5 59.0 55.5 55.5 51.1 47.0  49.5
     19 71.5 42.5 64.0 56.5 61.0 54.0 46.5 53.7 43.0  44.0
     20 57.5 51.5 51.3 34.0 55.0 61.0 55.5 51.3 51.3  44.5
    
     When I apply the gev function, I have the following error message:
    
     out - gev(dat, 20) Error in gev(dat, 20) : (list) object cannot be
   coerced to type 'double'
    
     Can anybody explain to me what does it means. I tried to read the
   package manual but could not understand about the block.
    
     Thank you so much for any help given.
        [[alternative HTML version deleted]]
    
    
    
     __
     R-help@r-project.org mailto:R-help@r-project.org
 mailto:R-help@r-project.org mailto:R-help@r-project.org