Error in if (k != round(k)) { : missing value where TRUE/FALSE needed
My guess is that it has something to do with the fact that you are trying to fit 5 parameters to essentially 77 data points (with a 20 hour warmup). I think sriv might struggle with that.
I would have to dig in the code on github, where exactly the error message comes from, but maybe Joseph G knows this of the top of his head.
I think it is sriv, because running without sriv (not using line 85 in the script) and then running fitByOptim only gives a warning:
Warning messages: 1: In nseStat(aggrQ, doaggr(X), ref = aggrRef, ...) : merged time series have no data; incompatible times?
This might have to do with the fact that your last day in the time series is not a complete day? Again not sure, as I normally don't work with hourly data.
I would suggest you try to things:
1. run with a longer time series and first run without sriv
2. run with a daily time series and see if that works
Good luck
Willem
Dear Professor Willem,
I have received your reply to my hydromad hourly data challenges and I greatly appreciate your help. You suggested I try two things namely:
1. run with a longer time series and first run without sriv
2. run with a daily time series and see if that works
However, the time series data I have is flood events data and is no longer than the end of the flood of that particular year. The next available data is the flood event data of the subsequent year and thus the data is not continuous. Also, I am not in possession of reliable continuous daily time series for my area of study.
I therefore humbly would like to ask if you have any other suggestions to help with the problems I am facing with the hourly data to help me continue with my research.
I hope to hear from you soon.
Thank you.
--
Best regards,
Thelma Baddoo
--
You received this message because you are subscribed to the Google Groups "hydromad users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to hydromad-users+unsubscribe@googlegroups.com.
To post to this group, send email to hydromad-users@googlegroups.com.
Visit this group at https://groups.google.com/group/hydromad-users.
For more options, visit https://groups.google.com/d/optout.
Hi Thelma,I plan to use hydromad cwi model with hourly data and see your question. I am not sure whether this model can be used with hourly data. I am thinking one problem. Is there any differences between using hourly data and daily data. If not, may be we can change the time of hourly data to daily dates and run the model. Do you think it is feasible ?Best regards,Binru
2018-04-16 9:44 GMT+01:00 Thelma Baddoo <mzded...@gmail.com>:
Dear Binru,I have still not been able to resolve my problem with using the hydromad cwi with hourly data. I would be glad to have a sample of your R script if you have been successful with it.Thank you.
On Monday, 16 April 2018 16:25:32 UTC+8, zbr...@gmail.com wrote:Hi Thelma,I am using the hydromad cwi with hourly data. I am wondering whether you have resolved your problem and how. Thanks very much.Best regards,Binru
在 2018年1月29日星期一 UTC上午10:52:26,Thelma Baddoo写道:Hi group,I am a new user of hydromad and I am trying to use hydromad cwi on hourly data. I have followed all the instructions in the hydromad tutorial and some suggestions from this hydromad group. However, when I run the fitbyOptim function, I get this error:Error in lm.fit(x = z2[, -1, drop = FALSE], y = z2[, 1]) :0 (non-NA) casesI don't know what I have done or am doing wrong. Any help will be greatly appreciated.I have attached my raw data (2000.csv) and my hydromad R script (hydromad trial script.R) for further details.Thank you
--
You received this message because you are subscribed to the Google Groups "hydromad users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to hydromad-user...@googlegroups.com.
To post to this group, send email to hydroma...@googlegroups.com.
To unsubscribe from this group and stop receiving emails from it, send an email to hydromad-users+unsubscribe@googlegroups.com.
To post to this group, send email to hydromad-users@googlegroups.com.
--
You received this message because you are subscribed to the Google Groups "hydromad users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to hydromad-users+unsubscribe@googlegroups.com.
To post to this group, send email to hydromad-users@googlegroups.com.
To unsubscribe from this group and stop receiving emails from it, send an email to hydromad-user...@googlegroups.com.
To post to this group, send email to hydroma...@googlegroups.com.
To unsubscribe from this group and stop receiving emails from it, send an email to hydromad-users+unsubscribe@googlegroups.com.
To post to this group, send email to hydromad-users@googlegroups.com.
--
You received this message because you are subscribed to a topic in the Google Groups "hydromad users" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/hydromad-users/cnEpNwmH4Dw/unsubscribe.
To unsubscribe from this group and all its topics, send an email to hydromad-user...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/hydromad-users/90b34eb9-cf6f-4e7c-9c34-9682b167554fn%40googlegroups.com.
Hello Thelma
How should I go about calibrating the model on several events at the same time? (i.e. having a single value for each parameter of all the period of my data)
Best regards