Re: [R-sig-Geo] Alternate statistical test to linear regression?

2019-10-24 Thread Greg Snow
> mildly skewed in distribution? I have large sample sizes for all of my > datasets, and the variables are continuous. > > That would pretty much cover all of my questions concerning this! > > Thank you, once again, for your time! > > -Original Message- > Fro

Re: [R-sig-Geo] Alternate statistical test to linear regression?

2019-10-23 Thread Greg Snow
gt; Pearson correlation coefficient, since both x and y variables are > quantitative. Does that make sense? > > Thanks again, > > > -Original Message- > From: Greg Snow <538...@gmail.com> > To: rain1290 > Cc: r-sig-geo > Sent: Wed, Oct 23, 2019 1

Re: [R-sig-Geo] Alternate statistical test to linear regression?

2019-10-23 Thread Greg Snow
Note that the normality assumptions are about the residuals (or about y conditional on x), not on the x variable(s) or all of y (non-conditional). If x is highly skewed and the residuals are normal then diagnostics just on y will also show skewness (if there is a relationship between x and y).

Re: [R-sig-Geo] Question on histograms from Raster and RasterVIS packages

2014-11-18 Thread Greg Snow
Nuno, Does rasterVIS use base or grid graphics? The place I usually see the error that plot.new has not been called yet is when trying to use a base graphics function to augment a plot created using grid graphics. I suspect that the plot created used grid graphics and that is why abline is not

Re: [R-sig-Geo] maps showing post-Soviet states?

2014-07-01 Thread Greg Snow
The combination of the sp and maptools packages will read in shapefiles and plot them, you just need to find a shapefile that has the countries that you are interested in (google for shape files, there are a lot of free ones available, I have not checked for the countries you specified, but I

Re: [R-sig-Geo] Creating a map with metropolitan areas

2013-12-26 Thread Greg Snow
What have you tried so far? Do you want base graphics? lattice graphics? gglot2 graphics? or any of them? I would probably start with the sp and maptools packages, but there are other tools listed in the spatial task view that may help as well. You are more likely to receive useful help if you

Re: [R-sig-Geo] random sampling along contour

2013-10-09 Thread Greg Snow
Look at the spsample function in the sp package. If you can convert the information on your contour(s) into a spatial line object, then spsample can help with choosing random points along that line/contour. That may accomplish what you want. On Wed, Oct 9, 2013 at 5:32 AM, Seth Bigelow

Re: [R-sig-Geo] Batch conver from img to rst

2013-09-16 Thread Greg Snow
Use the list.files function to get a list of your file names, save that in a variable, then use lapply or a loop to run your above code for each filename. The sub function is one option for changing the input filename to the output filename. On Mon, Sep 16, 2013 at 11:06 AM, Eddie Smith

Re: [R-sig-Geo] Batch conver from img to rst

2013-09-16 Thread Greg Snow
? On Mon, Sep 16, 2013 at 6:09 PM, Greg Snow 538...@gmail.com wrote: Use the list.files function to get a list of your file names, save that in a variable, then use lapply or a loop to run your above code for each filename. The sub function is one option for changing the input filename

Re: [R-sig-Geo] making an inset in a map figure

2013-09-13 Thread Greg Snow
Here is an example using the subplot function: library(TeachingDemos) library(maptools) data(wrld_simpl) plot(wrld_simpl,border=NA,col='blue',axes=TRUE, xlim=c(100,130), ylim=c(-40,30)) tmp - subplot( plot(wrld_simpl,border=NA,col='green',bg='white',axes=FALSE), 'bottomleft', size=c(2,1),

Re: [R-sig-Geo] R or R studio 3d manipulate

2013-07-17 Thread Greg Snow
step is the same but by visualizing a surface any suggestion? best paolo -- *Da:* Greg Snow [538...@gmail.com] *Inviato:* giovedì 11 luglio 2013 22.27 *A:* Paolo Piras *Cc:* r-sig-geo@r-project.org *Oggetto:* Re: [R-sig-Geo] R or R studio 3d manipulate

Re: [R-sig-Geo] R or R studio 3d manipulate

2013-07-11 Thread Greg Snow
? Anyway...thanks a lot!! Awaiting an answer from you I will try to play with TeachinDemos more accurately best regards Paolo -- *Da:* Greg Snow [538...@gmail.com] *Inviato:* lunedì 8 luglio 2013 17.53 *A:* Paolo Piras *Cc:* r-sig-geo@r-project.org *Oggetto

Re: [R-sig-Geo] R or R studio 3d manipulate

2013-07-08 Thread Greg Snow
If you can write a function that takes the arguments that you want to change and creates the graph you want to see, then there are several tools that can be used to create controls that will change those arguments. The `tkexamp` function in the TeachingDemos package is one such tool, but there

Re: [R-sig-Geo] Finding the Circumscribing Circle of a polygon

2013-01-26 Thread Greg Snow
The dists variable represents the distance between the points and the circle measured along the radius of the circle. Positive distances mean that the point/vertex is inside the circle and negative distances mean that the point is outside of the circle. Without some type of penalty the optimal

Re: [R-sig-Geo] background colour when add=T

2012-09-13 Thread Greg Snow
One possibility would be to define a spatial polygon that represents the entire plotting area, then use the gDifference function from the rgeos package to create a map that excludes the land areas (or whatever you don't want to cover) and add the resulting polygon(s). On Thu, Sep 13, 2012 at 8:04

Re: [R-sig-Geo] interacting with plots

2012-08-16 Thread Greg Snow
You could try something like: library(rgdal) library(sp) scot_BNG - readOGR(dsn=dsn, layer=scot_BNG) plot(scot_BNG) tmp - SpatialPoints(locator(1)) proj4string(tmp) - proj4string(scot_BNG) over(tmp, scot_BNG) # or tmp2 - over(scot_BNG, tmp) plot( scot_BNG[!is.na(tmp2),], add=TRUE, col='green')

Re: [R-sig-Geo] interacting with plots

2012-08-16 Thread Greg Snow
maptools? Thanks, Ty On Aug 16, 2012, at 6:33 PM, Greg Snow wrote: You could try something like: library(rgdal) library(sp) scot_BNG - readOGR(dsn=dsn, layer=scot_BNG) plot(scot_BNG) tmp - SpatialPoints(locator(1)) proj4string(tmp) - proj4string(scot_BNG) over(tmp, scot_BNG) # or tmp2

Re: [R-sig-Geo] interacting with plots

2012-08-16 Thread Greg Snow
, the output plots to Quartz, which doesn't seem to have a pan function and the zoom function is very very limited. Thanks! Ty On Aug 16, 2012, at 7:30 PM, Greg Snow wrote: My example used locator(1) to get one click, but if you just use locator() (without the 1) then you can click in as many

Re: [R-sig-Geo] looping over files

2012-03-03 Thread Greg Snow
A couple of options: Use the paste or sprintf functions to create your file names vector. This works if you want all combinations, or a very specific set of names. Use list.files to get the list of all the file names in the folder/directory, then use the grep function to extract just the files

Re: [R-sig-Geo] Overlaying rasters in spplot

2011-11-29 Thread Greg Snow
This is not an answer to the spplot question, but if your main goal is to plot your own data on top of google map tiles, then you may want to consider using the RgoogleMaps package which does exactly that. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare

Re: [R-sig-Geo] small data set for an example

2011-09-22 Thread Greg Snow
In the TeachingDemos package there are 3 datasets 'ccc', 'h2h', and 'towork' which are real, not overly large (but you can decide if they are small or not) which contain lat, long, and other variables. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare

Re: [R-sig-Geo] combining different plot types

2011-02-18 Thread Greg Snow
You need the gridBase package (for combining grid and base graphics). -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original Message- From: r-sig-geo-boun...@r-project.org [mailto:r-sig-geo-bounces@r- project.org]