Majid Javanmard
2016-Feb-20 10:33 UTC
[R] Is there any simple guide for spatialVx R package !?
I was searching a windows-based software to execute CRA(contiguous rain area) method, I failed, but I found a R package and its example to execute this method but I encountered some ambiguities !! here is the example code : data(pert000) data(pert004) data(ICPg240Locs) hold <- make.SpatialVx(pert000, pert004, loc=ICPg240Locs, projection=TRUE, map=TRUE, loc.byrow = TRUE, field.type="Precipitation", units="mm/h", data.name=c("ICP Perturbed Cases", "pert000", "pert004")) look <- FeatureFinder(hold, smoothpar=10.5, thresh = 5) plot(look) look2 <- minboundmatch(look, verbose = TRUE) plot(look2) craer( look2 ) 1) What kind of dataset should I import here !!? I have an observation and forecast files contains Lat/Lon and Value columns(observation 9 rows, forecast 16 rows). 2) Should I interpolate these files to the same bound by Arcmap ?! I really cant understand the tutorial of SpatialVx package, I appreciate if someone could response on this issue !! Many Thanks Majid Javanmard [[alternative HTML version deleted]]
Sarah Goslee
2016-Feb-20 13:12 UTC
[R] Is there any simple guide for spatialVx R package !?
The data commands in the example load sample data. You can run those lines in R and look at what the sample data are set up like, if the explanations aren't clear to you. There may also be help files specifically for the sample data. There's some information at https://www.ral.ucar.edu/projects/icp/SpatialVx/ but sadly no tutorial. Still, those folks are probably better to contact for package-specific questions. Sarah On Saturday, February 20, 2016, Majid Javanmard <javanmard.majid at gmail.com> wrote:> I was searching a windows-based software to execute CRA(contiguous rain > area) method, I failed, but I found a R package and its example to execute > this method but I encountered some ambiguities !! > > here is the example code : > > data(pert000) > data(pert004) > data(ICPg240Locs) > > hold <- make.SpatialVx(pert000, pert004, > loc=ICPg240Locs, projection=TRUE, map=TRUE, loc.byrow = TRUE, > field.type="Precipitation", units="mm/h", > data.name=c("ICP Perturbed Cases", "pert000", "pert004")) > > look <- FeatureFinder(hold, smoothpar=10.5, thresh = 5) > plot(look) > > look2 <- minboundmatch(look, verbose = TRUE) > plot(look2) > > craer( look2 ) > > 1) > > What kind of dataset should I import here !!? I have an observation and > forecast files contains Lat/Lon and Value columns(observation 9 rows, > forecast 16 rows). > > 2) > > Should I interpolate these files to the same bound by Arcmap ?! > > > I really cant understand the tutorial of SpatialVx package, I appreciate > if someone could response on this issue !! > > Many Thanks > Majid Javanmard > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at r-project.org <javascript:;> mailing list -- To UNSUBSCRIBE and > more, see > 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. >-- Sarah Goslee http://www.stringpage.com http://www.sarahgoslee.com http://www.functionaldiversity.org [[alternative HTML version deleted]]
Majid Javanmard
2016-Feb-20 17:11 UTC
[R] Is there any simple guide for spatialVx R package !?
I checked datasets,however, it is my ambiguities, I mean why there is no Lat/Lon in pert000 and pert004 datasets, and how conncet data to each other !!? another question : does this package match the two entities which RMSE is the lowest or correlation coefficient is the highest one ?! On Sat, Feb 20, 2016 at 4:42 PM, Sarah Goslee <sarah.goslee at gmail.com> wrote:> The data commands in the example load sample data. You can run those lines > in R and look at what the sample data are set up like, if the explanations > aren't clear to you. There may also be help files specifically for the > sample data. > > There's some information at > https://www.ral.ucar.edu/projects/icp/SpatialVx/ > but sadly no tutorial. Still, those folks are probably better to contact > for package-specific questions. > > Sarah > > > On Saturday, February 20, 2016, Majid Javanmard <javanmard.majid at gmail.com> > wrote: > >> I was searching a windows-based software to execute CRA(contiguous rain >> area) method, I failed, but I found a R package and its example to execute >> this method but I encountered some ambiguities !! >> >> here is the example code : >> >> data(pert000) >> data(pert004) >> data(ICPg240Locs) >> >> hold <- make.SpatialVx(pert000, pert004, >> loc=ICPg240Locs, projection=TRUE, map=TRUE, loc.byrow = TRUE, >> field.type="Precipitation", units="mm/h", >> data.name=c("ICP Perturbed Cases", "pert000", "pert004")) >> >> look <- FeatureFinder(hold, smoothpar=10.5, thresh = 5) >> plot(look) >> >> look2 <- minboundmatch(look, verbose = TRUE) >> plot(look2) >> >> craer( look2 ) >> >> 1) >> >> What kind of dataset should I import here !!? I have an observation and >> forecast files contains Lat/Lon and Value columns(observation 9 rows, >> forecast 16 rows). >> >> 2) >> >> Should I interpolate these files to the same bound by Arcmap ?! >> >> >> I really cant understand the tutorial of SpatialVx package, I appreciate >> if someone could response on this issue !! >> >> Many Thanks >> Majid Javanmard >> >> [[alternative HTML version deleted]] >> >> ______________________________________________ >> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see >> 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. >> > > > -- > Sarah Goslee > http://www.stringpage.com > http://www.sarahgoslee.com > http://www.functionaldiversity.org >[[alternative HTML version deleted]]