First check the length of your variable, e.g. variable length is 12 and if you are trying to call (i+1) when i is 12 The R will give you this error as it cant find 13th id in your variable.

Hope this will be helpful

Thank you

Dear Chintan thank you very much, it is now working without any error message

I am not sure from your comment how the problem was resolved. Anyway, for future reference may be I would like to point out that RClimDex explicitly requires an extra record at the beginning if the start time matches with the start of the data set. For instance, if you are interested to find out the indices from 1979 to 2010 and your data set also ranges from 1979 to 2010 then just add a dummy record of the last day of 1978. e.g. 1978 12 31 -99.9 -99.9 -99.9.

Hope it helps in future.

Exactly Mr. Jamal, you are right that you have to add dummy data or nan at the start and end of the data to overcome this problem.

Dear Jamal Uddin Khan , so I just need to add one line of data at the start as you indicated 1978 12 31 -99.9 -99.9 -99.9. Just one line?. I thanked Chintan Jethva since the Rclimdex runned only for data quality check, but when I run the Indices calculation I am still getting errors like:

Error in file(file, "rt") : cannot open the connection In addition: Warning message: In file(file, "rt") : cannot open file 'P0_Station_List.txt': No such file or directory

Mr. Philbert, please refer to this video as it will give you some idea about the Error, also check second link for more details.

https://www.youtube.com/watch?v=xao_G3c5LBU

http://etccdi.pacificclimate.org/RClimDex/RClimDexUserManual.doc

Thank you

Dear Chintan Jethva , thank you very much.

would you help to organize my file. attached the text file, so that I can run it in RClimdex?

Dear Mr. Philbert, The necessary modification has been applied in attached file which runs on my PC without any error.

I have just added dummy values at the start and end of the file and replaced zeros('0') with -99.99 as per its rule. No other changes have been made.

Also please run Indices calculation from your original year not from the dummy value.

A screen shot of output indices is attached to verify it with your output.

Hope now this will resolve your problem.

Dear Chintan Jethva , thank you very much for your help.

Dear Mr. Philbert, I hope every thing works fine now. However if you find any difficulty please inform.

Thank you.

Dear Chintan Jethva, Thank you very much for your detailed answer. I believe we do not need to replace the zero (0) rainfall in the input dataset with -99.9. This will mark those records as "missing" and many indices will not be calculated. I believe the correct approach will be to add only one line above if you have to start from the very beginning and one line below if you have to continue till the end.

I am quoting from the manual here about the impact of missing data - "Monthly indices are calculated if no more than 3 days are missing in a month, while annual values are calculated if no more than 15 days are missing in a year. No annual value will be calculated if any one month’s data are missing. For threshold indices, a threshold is calculated if at least 70% of data are present."

Could you please make a cross-check? I believe you might see a number of gaps in rainfall indices. As you have removed all the days without rainfall, you should see no Consecutive Dry Days (CDD) in your result.

For the dataset Mr. Luhunga provided I have run rclimdex taking baseline from 1989 to 2039, and added one line before and after the datarange. Then the rclimdex is run. It worked OK. Please find attached the results.

Thank you very much.

Yes Md Jamal Uddin Khan, I tried yesterday to run using the file provided by Chintan Jethva and it was computing the quality control part. But I tried the file you provided, it compute the indices without any error. Thank you. This could also be justified scientifically missing and zero are different things.

Another question if I have 100 stations, how can I simplify the computation of indices, To compute at one by one station takes time, is there a simple way to compute indices to these stations at once or by commanding one line.?

Mr. Md Jamal, I considered zero(0) values as missing one but in regarding to your suggestion if on that particular day if there is no rainfall than you have to use zero(0) instead of missing (-99.9), because this will change the entire statistical calculations, so you are right if that zeros shows No Rainfall than you have to keep it as zero for correct statistics.

So I would suggest Mr. Philbert to consider those zeros in data file and also to figure out and confirm that weather that zeros are indicates No Data or No rainfall.

Let's say you want to compute indices from Jan 1 1981 to Dec. 31 , 2010.

As first line of record of your text file, add this:

1980 12 31 -99.9 -99.9 -99.9

As last line of record add this:

2011 01 01 -99.9 -99.9 -99.9

Of course, in the user-defined parameters window, the first and last years of your base period will be 1981 and 2010, respectively.

Thanks.

Thank you very much, but if I have 100 stations, how to compute all these stations and save time. To compute at each station one by one takes time

This discussion has helped me a lot with errors too, thank you guys.

Thanks a lot, i have got my answer or find my solution here by changing first and last line with default parameter.

Please am also trying to use Rclimdex and this is what am seeing:

Error in nast[, , , i] : subscript out of bounds

can someone suggest a possible solution?,my data is as attached below.@ Chintan Jethva

Dear Christian Anayo Njoku ,

Please find the modified copy of your file.

As mentioned in previous answer; I just added dummy entries above and below of your data set to get rid of this error.

Write 1971 in the first and last year for indices calculation.

Hope this will help you.

The threshold values depends to locations and therefore study need to be carried out to establish threshold for specific location. However, you can use the default thresholds as they do not affect the computation of many indices such as (TXx, TNn etc).

I agree with Philbert here. In practice, user defined thresholds allows to look into the statistics that are not readily available. For instance, in our recently published paper (Article Observed trends in climate extremes over Bangladesh from 1981 to 2010

However, the objective of the predefined indices was to derive the same statistics over the whole globe, and many useful climate statistics is readily available from them. It provides a means to distinguish one region in the context of global climate.

Hi everyone, I am trying to use ClimDex Ver 1.3 (excel) to compute the indices with using different set of example datasets. End up, error pop up. Pls someone can guide me how to solve this problem

Have you checked the length of the data?

The length/dim might not match

**sometimes that we're using multiple stations/grids/GCMs/etc, the length of different inputs might differ and it might be one of the reasons.

Can you provide more detail about your code/data?

I had this kind of problem, please check how people helped me.

Thank you guys for this interactive and helpful scientific chat.. problem solved... Please can anyone suggest me a good article/material to read on the use, application and interpretation of this Rclimdex data

Hello everyone, I'm using ClimPACT2 for calculating climate indices.

R calculated temperature indices but did not calculate precipitation indices. Are there anyone who had the same problem? Or any idea how can I get the precipitation indices?

I am attaching my data, if anyone can check it and help me, I would be grateful.

Thank you.

Hi İzem Çetinkaya

Do you mind letting me know what the 4th, 5th and 6th columns are?

If you have tasmin, tasmax and precipitation data, you can directly use "climdex.pcic" package

https://cran.r-project.org/web/packages/climdex.pcic/climdex.pcic.pdf

Hi Jason Mahmoudi

Those columns are precipitation, tmax and tmin, respectively.

Thank you, I'll check the link.

Do you have any idea about the error which I attached the screenshot of it?

There are many missing data for precipitation , could it be the reason?

climdexdata

İzem Çetinkaya TBH I haven't used ClimPACT2 before.

I usually code in R

Thank you so much. Jason Mahmoudi

hi, is there any standard to select threshold value for precipitation in RClimdex while calculating indices?

Thanks a lot guys for this discussion, helped me to solve my problem.

Hammad Salahuddin

I am unable to get plots of indices. the GUI is only creating indices and log folder but no plot folder when using the RclimDex package (1.9-3). can anyone help.