[R-sig-eco] SPEI index calculation

[Drew Tyre]

Hi Manuel

Without a little more information about how you calculated PET it is hard to know for sure how to help. But I would take a close look at the PET object as a next step to make sure there are no NA values in there. 

Which package and version are you using for spei() function?

You can see how my collaborators and I did this at 
https://datadryad.org/stash/dataset/doi:10.5061/dryad.f2mp4


-- 
Drew Tyre
 
He/him/his
 
While I may be sending this email outside my normal office hours, 
I have no expectation to receive a reply outside of yours.
 
School of Natural Resources
University of Nebraska-Lincoln
416 Hardin Hall, East Campus
3310 Holdrege Street
Lincoln, NE 68583-0974
 
phone: +1 402 472 4054
fax: +1 402 472 2946
email: atyre2 using unl.edu
http://snr.unl.edu/tyre
http://drewtyre.rbind.io
 
 

On 1/27/20, 10:41 AM, "R-sig-ecology on behalf of Manuel Esteban Lucas Borja" <r-sig-ecology-bounces using r-project.org on behalf of ManuelEsteban.Lucas using uclm.es> wrote:

    Dear all, 
    
    I am trying to calculate SPEI index using Rstudio (please see data used below this email). First I calculated PET (latitude 40,29), and then I failed when asking for SPEI index: 
    
     spei(SPEI$PREC-PET,1). 
    
    The only result I got is NA. 
    
    Any idea about what´s wrong? 
    
    Thanks for helping and for your time, 
    
    Manuel 
    
    Year	Month 	TMED	PREC	
    2010	1	0.9	59.8	
    2010	2	1.43	66.7	
    2010	3	3.3	35	
    2010	4	7.9	23.5	
    2010	5	9.58	22.9	
    2010	6	14.3	23.1	
    2010	7	20.6	0.5	
    2010	8	18.95	7.6	
    2010	9	13.94	25.6	
    2010	10	8.29	60	
    2010	11	3.14	71.2	
    2010	12	2.05	95.8	
    2011	1	1.5	35	
    2011	2	2.8	66	
    2011	3	4.8	59.8	
    2011	4	10.1	68.6	
    2011	5	12.7	24	
    2011	6	16.7	18.2	
    2011	7	18.2	12	
    2011	8	19.9	10.6	
    2011	9	15.8	1.7	
    2011	10	10.9	8.6	
    2011	11	5.9	78.6	
    2011	12	2.7	18.5	
    2012	1	1.86	10.3	
    2012	2	-0.89	4.2	
    2012	3	4.9	4.5	
    2012	4	5.4	119.8	
    2012	5	12.8	7.7	
    2012	6	17.6	21.7	
    2012	7	19.6	1.8	
    2012	8	20.8	1.6	
    2012	9	14.7	49.5	
    2012	10	9.64	61	
    2012	11	5.55	95.5	
    2012	12	1	46	
    				
    				
    				
    				
    				
    				
    				
    				
    				
    				
    				
    _______________________________________________
    R-sig-ecology mailing list
    R-sig-ecology using r-project.org
    https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_mailman_listinfo_r-2Dsig-2Decology&d=DwIGaQ&c=Cu5g146wZdoqVuKpTNsYHeFX_rg6kWhlkLF8Eft-wwo&r=aLEaSryyUcERVqcKVZl7lQ&m=7X13dMuFRoBufpNmVGN6OG2Ex1TlLBoogThOfS0djQE&s=sl2G4mgZNLLVruSfeIjLRUROuw-rpqMxi56CNk6-prA&e=