I would like to use ETAS package for R (Prof. Abdollah Jalilian), but I find an error when I try to use my catalog.
I have a catalog with the following header: Datenum, Lon, Lat, Mag and Depth.
I use the following comand in R to prepare ppx object and run the package:
> data Newdata timedomain = c(floor(min(Newdata$time)),ceiling(max(Newdata$time)))
> longdomain = c(floor(min(Newdata$long)),ceiling(max(Newdata$long)))
> latdomain = c(floor(min(Newdata$lat)),ceiling(max(Newdata$lat)))
> ppxdomain = boxx(t=timedomain,lon=longdomain,lat=latdomain)
> iceland.quakes = ppx(data=Newdata, domain=ppxdomain, coord.type =c("t", "s", "s", "m", "m"))
> check.temporal(iceland.quakes)
> plot(iceland.quakes)
> win plot(iceland.quakes$data[, 2:3])
> plot(win, add=TRUE)
> proj tperiod param01 res Error in if (tstart2 < tstart || tstart2 > tend || tlength + tstart > : missing value where TRUE/FALSE needed
But as you saw at the end I recieve an Error message.
Is there somebody that can help me to understand where is the error in the code?
Thank you very much
Francesco
Does your data have missing time values? If so, try creating a copy of the data frame with those rows removed and see if etas runs.
Thank you very much. Is there a specific comand in R to verify if some missing data are in the time column?
easiest thing to do is summary (Newdata) and see which variables report any NA values.
Your error message states:
Error in if (tstart2 < tstart || tstart2 > tend || tlength + tstart > : missing value where TRUE/FALSE needed
That means the code within the if statement doesn't give the expected T/F output. Debug by printing out each of the variables within the if statement or split the code to parts to find out if the output of each part gives a T/F like so:
tstart2 < tstart
tstart2 > tend
- Badges
- Science topic
- Similar topics
- Mathematics
- Statistics
- Statistical Software
- R Statistical Package