Dear Jon
I have to thank you so much for solving this issue. I did't notice that
there was a projection and grid problem.
I hope I would be able to move forward with my analysis with R and perhaps
still need your kind and much appreciated help.
Cheers
Anand
On Wed, Aug 5, 2015 at 12:08 PM,
Hi
Thank you Roman for your kind response.
The data is as attached which is part of the whole dataset.
My other problem is that im not able to create an appropriate prediction
grid with some data included.
The codes which failed to show colour ramps are
map - krige(logcarbseq ~ 1, locations =
DearJon
Thnaks. I will try that soon.
Cheers
Anand
On Tue, Aug 4, 2015 at 4:55 PM, Jon Skoien jon.sko...@jrc.ec.europa.eu
wrote:
Anand,
I could not open your .grd file in R, I think the .gri file is missing.
However, looking at it in a text editor, and comparing with your R-script,
it
Hi
I have kriged my data but the map shows uniform colours - no gradual
contour colour plot on map
code: print(spplot(k40, var1.pred, asp=1,
col.regions=bpy.colors(64),main=OK prediction, log-ppm Zn))
thnaks for any help
Anand Sookun
--
---
Dear,
Does anyone know how to convert multiband (5 bands) 16 bit unsigned integer
(tif file) to 8 band tif with R? Help will be appreciated.
Thanks
Zia
From: r-sig-geo-boun...@r-project.org r-sig-geo-boun...@r-project.org on
behalf of Tomislav Hengl
On 27 Mar 2015 23:50, Robert J. Hijmans r.hijm...@gmail.com wrote:
Are these really reasonable reasons? I do not think so, given that the
question had nothing to do with map navigation and the person asking
appears to live in the UK.
??? How did you figure that out? He gave some phone numbers
If you have lat-long data that crosses two UTM zones then its
generally okay to just pick *one* and transform all the points to
that. Use the one that has the most points in. Basically use the UTM
zones as guidelines to pick one UTM zone coordinate system. Unless
your data spans several zones and
Hello
I have animal movement data that I have converted from Lat/Long to UTM,
unfortunately the data falls in two UTM zones (34S and 35S). For some reason R
cannot display both of them in the same window (the 35S data is way off the
expected location).
The question is how do I convert the data
A number of field folks prefer UTM because
-it matches legacy paper USGS quad map series traditionally used for field
navigation
-units are in meters and can be used to gauge field distances from a coordinate
readout
On Mar 27, 2015, at 10:48 AM, Michael Sumner mdsum...@gmail.com wrote:
They are reasonable reasons, but traversing zones is a pain, you should see
if using one or the other is sufficient. I would check carefully the
distances you get against ellipsoidal calculations.
Cheers, Mike
On Sat, 28 Mar 2015 07:30 Andrew Duff andrewad...@gmail.com wrote:
A number of field
Are these really reasonable reasons? I do not think so, given that the
question had nothing to do with map navigation and the person asking
appears to live in the UK. Moreover, other projections have, or can
have, their units in meters as well (or feet or miles or whatever you
might fancy). UTM
I would like to know whether R gstat package (or any other) can be set up for
the following scenario.
1)Kriging/Simulation of Gaussian random fields over 128x128 grids
with values known at about 100 locations to be preserved - i.e. conditional
simulation.
2)The
I have the following problem:
I am using the reprojectHDF() function (http://r-gis.net/?q=ModisDownload
I change the input syntax:
input-list.files(path=.,.hdf,all.files=T,recursive=T)
and now my error is:
reprojectHDF(hdfName=input,filename=outname,MRTpath=~/MRT/bin,+
resample_type=NEAREST_NEIGHBOR,+
[mailto:r-sig-geo-bounces@r-
project.org] On Behalf Of A. Robles R.
Sent: 29. juli 2014 08:47
To: r-sig-geo@r-project.org
Subject: Re: [R-sig-Geo] Help me please
I change the input syntax:
input-list.files(path=.,.hdf,all.files=T,recursive=T)
and now my error is:
reprojectHDF(hdfName
The attached file contains collected soil salinity samples and extracted
pixels values from Ikonos image.
I am trying to interpolate soil salinity data and use regression models
with the axillary data
I created a grid as follows:
mydata - read.table(C://mydata.txt, header = TRUE)
-project.org r-sig-geo@r-project.org
Subject: [R-sig-Geo] Help mosiacing rasters OR replacing NA values in
   one   raster based on cell coordinates
Message-ID:
   1404304631.91596.yahoomail...@web142504.mail.bf1.yahoo.com
Content-Type: text/plain
Hi
Sorry for the complicated example, but I'm
Hi
Sorry for the complicated example, but I'm really stuck:
I am trying to merge two rasters/datasets: one for the continental USA and one
for Mexico. The dataset for Mexico was initially derived from a larger raster
of North America that was at the same extent, resolution etc. as the US
Mario,
The issue with the current version of intamap is that the anisotropy
detection method is too slow and too memory consuming for larger data
sets. The old version used a method based on the akima package, but this
was discontinued because of incompatible licenses (akima is strictly
Dear fellow user of R
I need a Little help
I am trying to estimate Anisotropy parameters, from a large data set (13
thousand data) that have Anistropy. All of these for correct the
coordinates, allowing use of a directional Variogram for estimate the
parameters that will be used for do a
I need help to connect R with Grass. In RStudio when I send the line:
loc - initGRASS(C:/Program Files/R/R-3.1.0/bin/i386, home=tempdir(),
override=TRUE)
and it says: Error in system(paste(g.dirseps.exe -g,
shQuote(Sys.getenv(WINGISRC))), :
'g.dirseps.exe' not found
Guys, can you help me?
, 10 Jan 2014 11:04:26 -0500
From: Jill Deines jillian.dei...@gmail.com
To: r-sig-geo@r-project.org
Subject: Re: [R-sig-Geo] Help with getting and processing MODIS data
Message-ID:
CAKCTE4zfgoS-1h0UkuHZO=hphoi6gjr5rqd7jknqbyciu4j...@mail.gmail.com
Content-Type: text/plain
Hi Julie,
I
Hi,
I have a distribution of artifact within an archaeological surface
(dataset: ID, DIMENSION, X, Y), and I need to verify which is the degree of
aggregation of these artifacts at different scales. I applied the L
Function (Lest in spatstat), and plotted the resulting observed values of
L(r)
Yes, it is possible and very easy. How do you extract your p-value
depends on wether you are making a pointwise or a global test. See the
help page of envelope(). You can also try a maximum absolute deviation
test with dclf.test().
Cheers,
Marcelino
El 13/01/2014 16:52, Francesco Carrer
I would like to add:
* If you perform a pointwise test, the value of r at which you conduct
the test much be chosen a priori --- before collecting, or at least
before looking at the data. Otherwise the associated p-value is
meaningless. It is *very* unlikely that you had an a priori value
,utm_zone=30,datum=WGS84,pixel_size=1000)
Good luck,
Babak
Date: Fri, 10 Jan 2014 11:04:26 -0500
From: Jill Deines jillian.dei...@gmail.com
To: r-sig-geo@r-project.org
Subject: Re: [R-sig-Geo] Help with getting and processing MODIS data
Message-ID:
CAKCTE4zfgoS-1h0UkuHZO
From: Rolf Turner [r.tur...@auckland.ac.nz]
Sent: Tuesday, 14 January 2014 3:40 AM
To: Marcelino de la Cruz
Cc: r-sig-geo@r-project.org; Adrian Baddeley
Subject: Re: [R-sig-Geo] Help L Function
I would like to add:
* If you perform a pointwise test, the value of r
.
The package has pretty good help files.
Jill
Message: 8
Date: Thu, 9 Jan 2014 11:59:14 -0800 (PST)
From: Julie Lee-Yaw jullee...@yahoo.ca
To: r-sig-geo@r-project.org r-sig-geo@r-project.org
Subject: [R-sig-Geo] Help with getting and processing MODIS data
Message-ID
'runGdal' or 'runMrt' to process with either tool.
The package has pretty good help files.
Jill
Message: 8
Date: Thu, 9 Jan 2014 11:59:14 -0800 (PST)
From: Julie Lee-Yaw jullee...@yahoo.ca
To: r-sig-geo@r-project.org r-sig-geo@r-project.org
Subject: [R-sig-Geo] Help with getting
Hi allÂ
I am fairly new to GIS and R and am trying to download, reproject, mosaic and
save MODIS NDVI and EVI data as a geotiff. I need tiles from across North
America, so a lot of data and would prefer to find an automatic solution. I've
come across modisDownload but can't seem to get it to
Dear James,
I have been able to run your code. See attached files. Not sure why
OpenBUGS crashed but I found a few strange symbols (vertical bars) when
I copied and pasted your model from your e-mail. OpenBUGS complained and
it worked fine after I had these symbols removed.
I have run it on an
!
James
From: Virgilio Gómez-Rubio [virgilio.go...@uclm.es]
Sent: 27 November 2013 20:29
To: James Rooney
Cc: r-sig-geo@r-project.org
Subject: Re: [R-sig-Geo] Help - calling OpenBugs through R2OpenBugs problem
Dear James,
I'll be happy to look
...@gmail.com]
Sent: 28 November 2013 09:11
To: James Rooney
Cc: r-sig-geo@r-project.org
Subject: Re: [R-sig-Geo] Help - calling OpenBugs through R2OpenBugs problem
Hi All
I try to follow the above example but i got stuck in reading Spatial Data Frame
any Idea or Help?
Here are my R codes with error. I am
...@gmail.com]
Sent: 28 November 2013 09:11
To: James Rooney
Cc: r-sig-geo@r-project.org
Subject: Re: [R-sig-Geo] Help - calling OpenBugs through R2OpenBugs problem
Hi All
I try to follow the above example but i got stuck in reading Spatial Data
Frame any Idea or Help?
Here are my R codes
From: Yuster Ronoh [yust...@gmail.com]
Sent: 28 November 2013 10:21
To: James Rooney
Cc: r-sig-geo@r-project.org
Subject: Re: [R-sig-Geo] Help - calling OpenBugs through R2OpenBugs problem
Hi James
I Downloaded the correct file from
http://www.cso.ie/en/census/census2011boundaryfiles
I don't know then Yuster. Is your rgdal package up to date ?
From: Yuster Ronoh [yust...@gmail.com]
Sent: 28 November 2013 11:08
To: James Rooney
Cc: r-sig-geo@r-project.org
Subject: Re: [R-sig-Geo] Help - calling OpenBugs through R2OpenBugs problem
Yes
Rooney [roone...@tcd.ie]
Sent: 23 November 2013 16:30
To: r-sig-geo@r-project.org
Subject: [R-sig-Geo] Help - caling OpenBugs thourgh R2OpenBugs problem
Hi all,
I am working on a Bayesian spatial smoothing project - and have implemented the
BYM model, similar to the description in Chapter 11
: 23 November 2013 16:30
To: r-sig-geo@r-project.org
Subject: [R-sig-Geo] Help - caling OpenBugs thourgh R2OpenBugs problem
Hi all,
I am working on a Bayesian spatial smoothing project - and have implemented
the BYM model, similar to the description in Chapter 11 of the ASDAR book
: [R-sig-Geo] Help - calling OpenBugs through R2OpenBugs problem
Dear James,
I'll be happy to look into this but I will need more information about
your system. Also, could you run the code using 'debug=TRUE' as seen in
the error message?
Best wishes,
Virgilio
On mié, 2013-11-27 at 20:13 +
Hi all,
I am working on a Bayesian spatial smoothing project - and have implemented the
BYM model, similar to the description in Chapter 11 of the ASDAR book. However,
I have been having a huge problem calling OpenBugs from within R and I cannot
figure out why. I would veyr much liek to do
Subj: Help for optim
Dear all,
I want to optimize a function using optim. I
have two problems:
1-how can i minimize function which fn value
be zero not negative. Because i want optimize parameters which lead to zero fn.
2- how can i change it to maximization?
Best,
-boun...@r-project.org [r-sig-geo-boun...@r-project.org] p#229;
vegne av Saman Monfared [samanmonfar...@gmail.com]
Sendt: 30. juli 2013 4:11
Til: r-sig-geo@r-project.org
Emne: [R-sig-Geo] help for optim
Subj: Help for optim
Dear all,
I want to optimize a function using optim. I
have two problems:
1
from the section Details of
?optim
By default ‘optim’ performs minimization, but it will maximize if
‘control$fnscale’ is negative.
Alternatively your chance the sigh of what is returned by your function
Paulo Justiniano Ribeiro Jr
LEG (Laboratorio de Estatistica e Geoinformacao)
Dear all,
I want to optimize a function using optim. I have two problems:
1-how can i minimize function which fn value be positive because it can be
negative?
2- how can i change it to maximization?
Best,
[[alternative HTML version deleted]]
Dear all,
I want to optimize a function using optim. I have two problems:
1-how can i minimize function which fn value be positive because it can be
negative?
2- how can i change it to maximization?
Best,
[[alternative HTML version deleted]]
library(gstat)
my.show.vgms = show.vgms
my.show.vgms = edit(my.show.vgms)
my.show.vgms()
what exactly you have to change during the edit phase is beyond the
context of this list, it is basic R programming.
On 07/10/2013 07:51 PM, Saman Monfared wrote:
Dear all,
How can I add text and change
Dear All,
How can I plot some different plots in one page in lattice
like par(mfrow=c(2,2)) in graphics package?
p1-bwplot(m~PRED,data=PRED)
p2-bwplot(R~Residual,data=PRED)
xyplot(obs~PRED.v|method,data=PRED.v)
histogram(~PRED.v|method,data=PRED.v)
Thanks,
Saman.
--
Saman Monfared
Msc,
Have you tried `grid.arrange` from `gridExtra` package?
http://stackoverflow.com/questions/6681145/how-can-i-arrange-an-arbitrary-number-of-ggplots-using-grid-arrange
Cheers,
Roman
On Wed, Jul 10, 2013 at 2:14 PM, Saman Monfared samanmonfar...@gmail.comwrote:
Dear All,
How can I plot some
?print.trellis
Next time use the R-help mailing list
(https://stat.ethz.ch/mailman/listinfo/r-help) for these questions.
Best,
Oscar.
Saman Monfared samanmonfar...@gmail.com writes:
Dear All,
How can I plot some different plots in one page in lattice
like par(mfrow=c(2,2)) in graphics
Dear all,
How can I add text and change color of variograms in show.vgms()?
Best,
Saman.
--
Saman Monfared
Msc, Department of Statistics, Shiraz University,
Shiraz 71454, Iran
Email: samanmonfar...@gmail.com
Tel: +98 917 5305167
___
R-sig-Geo
Thank you,
Why sigma = 10?, what function can I use for choose the best sigma?. And I
looked the info about density.ppp and it says it use the Gaussian Kernel.
But, if we look on wikipedia:
http://en.wikipedia.org/wiki/Kernel_(statistics) , it says, there are more
kernel function than the Gaussian
In general it is acknowledged that it is more important to select a
bandwith (i.e., a sigma for the Gaussian kernel) than the type of kernel
(see for example the 1st ed. of the ASDAR book, p. 167).
You can select a sensible bandwith with the functions bw.diggle in
spatstat or mse2d in splancs.
I have calculated time series NDWI (float TIFFF raster ) values from Landsat
images. In few pixels of all raster (in a raster stack) objects have values
-1.000 and 1. . I like to convert values that are -1 to -1.000 and
values that are 1 to +1.000 of all raster and save as getif (float)
Zia,
I like to use the plot function from the raster package.
?raster:::plot
To continue your example, I'm able to plot the individual layers of a stack
using the lapply() function. I'm not sure if there is a simpler way to do
it, but this way works for me:
library(raster)
library(sp)
Dear All,
I try to plot some variogram models in a panel by show.vgms in gstat.
my program is:
show.vgms(min =0, max =20, n = 50, sill =c(v.ols$psill[2],v.wls$psill[2],
v.gls$psill[2],v.rml$sigmasq,v.reml$sigmasq,v.opvme$psill[2],v.opv$psill[2],
v.opvab$psill[2],v.opvlin$psill[2]/3),range
Hi,
I need to create a weight matrix in R for point dataset of 78,000
observations from a shapefile( from ArcGIS) which has both sets of
longitudes and latitudes and X-Y coordinates (based on state plane feet).
So far I projected the dataset in ESRI coordinates using rgdal maptools.
I
Jonathan,
I think you are confusing many aspects here...
The ENFA works on two objects:
* a PCA, which provides the environmental data;
* a vector, which provides the utilization distribution.
1) The PCA is build on the data slot of a SpatialPixelsDataFrame, which
stores all environmental
Hi there
I have managed to preapre my data finally to do the enfa. However, when I
run the test I get an error and warning message shown below:
ENFA - enfa(pc, pr, scannf = FALSE)
Error in crossprod(Ze, DpZ) : non-conformable arguments
In addition: Warning messages:
1: In x * pr :
longer object
Hi there
I have been working on this over the weekend. I have managed to make some
progress.
Now I am trying to run the enfa using this code: enfa(pc, pr, scan=FALSE).
I can define pc using my text file of the locality data for my species using
this code: pc - dudi.pca(locs, scannf=FALSE).
Hi,
Typing ?enfa (after loading adehabitatHS) and searching through the
examples that are provided, you should find these lines
## We prepare the data for the ENFA
tab - slot(map, data)
pr - slot(count.points(locs, map), data)[,1]
Maybe you could start with this, i.e. understanding how the
Hi,
Typing ?enfa (after loading adehabitatHS) and searching through the
examples that are provided, you should find these lines
## We prepare the data for the ENFA
tab - slot(map, data)
pr - slot(count.points(locs, map), data)[,1]
Maybe you could start with this, i.e. understanding how the
Hi Alexandre
I have looked at the examples already and those in SP package too before I
posted my question.
When I run this line of code from the example: pr - slot(count.points(locs,
join), data)[,1], I get the following error:
Error in count.points(locs, join) :
w should inherit the class
Hi Robert and Geno,
thank you for your help i tried both solutions but none of them worked.
when i transformed the grid raster to tif in the second method, instead of
having NAs i had 128 as value in all the attribute table
On Tue, Apr 9, 2013 at 9:45 AM, Robert J. Hijmans
Hi Mathieu
Thank you greatly, I have been using your advice and it seems to be working.
I used readAsciiGrid to import my 15 variables and then cbind to create the
new data frame. The result is below:
join - cbind(x@data, y@data, a@data, b@data, c@data, d@data, e@data,
f@data, g@data, h@data,
Hi Mathieu
Thank you again. You are really helping me. Please forgive some of my
questions because I am very knew to R and I have to teach this to myself.
AdehabitatHS is using a different method to read the raster files. So when I
input this code: var1- import.asc(Data/bio_01.asc) I get the
Dear Jonathan,
To import the maps, have a look at the function readAsciiGrid from the
package maptools.
To import the points, read.csv or read.table (use your prefered way to
import data in R).
Con
HTH,
Clément Calenge.
On 04/10/2013 09:59 AM, JonathanAronson wrote:
Hi Mathieu
Thank you
Hi Clement
Thank you very much for your input. I have read the maptools information and
seen the examples for readAsciiGrid. I would like to create a
SpatailPixelDataFrame which contains the raster maps of 15 environmental
variables. Therefore, similar to the lynxjura$map
perhaps I can clarify my question a lit bit further:
With the adehabitat package I used the function as.kasc to create a single
data frame with my 15 variables as Ive shown below:
map- as.kasc(list(bio_01 = var1, bio_12 = var2, Aridity = var3, altitude =
var4, unpalatable = var5, thicket = var6,
Jonathan,
Again, this is no more an adehabitat question: the whole idea of the new
series of packages (including adehabitatHS) is to rely on spatial R classes
that are available with sp. All adehabitat functions to handle maps (such
as as.kasc) have thus been removed.
This said, this should
Eliane,
I think that is a bug (even though it should not affect the results as I
take it this does not stop the function). I think it will be gone withe
next version of raster. You can probably also avoid it by doing
a-C:/QA_soil/inttotalrain
aa-raster(a, RAT=FALSE)
On Mon, Apr 8, 2013 at 4:44
Hello there
I am really hoping someone out there can help me. I have read the various
forums but havent found information that can help me.
I am running an ENFA in R for 12 species using 15 environmental variables.
Below is my code and some questions:
library(ade4)
library(adehabitat)
test-
Dear Jonathan,
Here are a few comments:
1) I strongly suggest that you update to 'adehabitatHS', which has
superseded the good-old 'adehabitat'. See [1]. It will solve your error
with 'niche.test'.
2) For the tolerance, 'niche.test' will give you both the marginality and
the tolerance.
Hi Mathieu
Thank you so much for your quick response. I really appreciate the effort
and thank you for your answers.
I am now using the AdehabitatHS package but still have some questions.
1) I am not sure what (renfa$expvar), (renfa$pvalue) and (renfa$obs) are
actually telling me in the output.
Le 04/09/2013 03:51 PM, JonathanAronson a écrit :
Hi Mathieu
Thank you so much for your quick response. I really appreciate the effort
and thank you for your answers.
I am now using the AdehabitatHS package but still have some questions.
1) I am not sure what (renfa$expvar), (renfa$pvalue)
Hi Eliane,
Not sure that it will work, but I would recommend converting your raster layers
to .tiff files prior to importing them into R. Then try
aa-raster('C:/QA_soil/inttotalrain.tif',datatype=INT1U)
(note the file extension, and data type argument).
Good luck!
Geno
hi list,
i am using R
hi list,
i am using R 2.15.3 (64 bit)
i want to stack several raster objects to use for prediction purposes. All
my raster files are integer and has attribute tables and none has NA
values. ( i checked it in arcgis).however when i try to import the raster
to R and check the values in each
Hi Mike,
Needed to flip the y axis to get the orientation right. Levelplot() show
images that look like something I expected to see, and in the right part of
the world - all good. I will need to check the x,y resolution of the
original matlab file, but it is in the right ballpark.
Thanks for
hi list,
i am using R 2.15.3 (64 bit)
i want to stack several raster objects to use for prediction purposes. All
my raster files are integer and has attribute tables and none has NA
values. ( i checked it in arcgis).however when i try to import the raster
to R and check the values in each
Hi All,
I have been provided with a .mat file containing a time series of Sea
Surface Temperature data (50 x 42 cells with 92 time layers). It was a
Struct object in Matlab. I can happily import the file in to R using
R.matlab, creating a list as follows. My question is how then to convert
this
Hi Craig, see the array method for ?brick:
brick(x, xmn=0, xmx=1, ymn=0, ymx=1, crs=NA
So, something like
x - brick(test$data[[2], xmn = min(test$data[[5]]), ymx =
max(test$data[[5]]), ymn = min(test$data[[4]]), ymx = max(test$data[[4]]))
You might need to mess around to get the x/y ranges
Dear list,
I have 2 raster that represent exactly the same area. the 2 raster have the
same coordinate system and the same resolution but 2 different origins and
extents. They are both float. My goal is to perfectly overlap the 2
rasters, so i can use it for analysis.
I have tried to use resample
Dear Roger Bivand and list members,
I am looking for advice on how to interpret a correlogram of species
richness data collected from photographic transects of the seabed. I am
an amateur in R, and am completely new to spatial analysis, so please
forgive the naivety of my issue.
I have four
Hi
I'm work whit ggmap and I'd like to know how I can plot plot the rate of
one disease using the cols[findInterval...]
How I do that?
My script is:
stores - data.frame(name=c(Jaboatão dos Guararapes, Olinda,Paulista,
Recife, Maceió),
longitude=c(-35.00329,-34.86707,-34.88437,-34.93294,
Dear all.
I try to fit conditional auto regressive models (CAR and SAR) in package spdep.
Also, I have fitted some other models like GLM, Empirical Bayes and ...
My program to CAR and SAR is:
esar1f1 - spautolm(IMR.m~ 0+PCFP+Sup+B.F,data =data1,
listw=nb2listw(nb6, style=W), family=SAR,
Dear all.
I try to fit conditional auto regressive models (CAR and SAR) in package spdep.
Also, I have fitted some other models like GLM, Empirical Bayes and ...
My program to CAR and SAR is:
esar1f1 - spautolm(IMR.m~ 0+PCFP+Sup+B.F,data =data1,
listw=nb2listw(nb6, style=W), family=SAR,
Dear all.
I try to fit CAR and SAR model.
I have fitted these model by spatulm but I don't understand the mean
of lambda???
I couldn't find good renascences to distinct the difference between SAR and CAR.
What is the difference between BYM and CAR model??
Hw can I predict results of spatulm for
On Tue, 18 Dec 2012, Saman Monfared wrote:
Dear all.
I try to fit conditional auto regressive models (CAR and SAR) in package spdep.
Also, I have fitted some other models like GLM, Empirical Bayes and ...
My program to CAR and SAR is:
esar1f1 - spautolm(IMR.m~ 0+PCFP+Sup+B.F,data =data1,
Dear list,
I have an ESRI shape file with county names and boundaries for Florida, and an
ESRI shape file with sea grass cover for Florida. I would like to know the
total sea grass cover per county. Hence I would like to identify the
intersections of the county polygons with the sea grass
Please, does anyone know what I need to do - my call lm(regression model)
is not giving the full output result. Thanks in advance
Sam
On 9/20/12 2:13 AM, Ross Dwyer ross.dw...@uq.edu.au wrote:
Hi Tony,
I just stumbled on this link in the archive which will convert
SpatialPolygons to a
What do you mean by full output result? The ANOVA table? If so, try
summary(lm(regression model)).
Cheers,
Roman
On Thu, Sep 20, 2012 at 11:06 AM, Idowu, Oluseun S. (UMKC-Student)
oluseun.id...@mail.umkc.edu wrote:
Please, does anyone know what I need to do - my call lm(regression model)
Hello All:
I am working on parcel-level housing dataset to estimate the impact of
various variables on home sale prices (N is approx. 90,000).
I created the spatial weight metrics using sale year of four nearest houses
to assign weights. Next, I ran LM tests and the spatial lag and error
Epidemiology
Imperial College London
St. Mary's Hospital
London
W2 1PG
From: Manuel Spínola [mspinol...@gmail.com]
Sent: 19 July 2012 03:26
To: O'Hanlon, Simon J
Cc: r-sig-geo@r-project.org
Subject: Re: [R-sig-Geo] Help with gpclib
Thank you very much Simon,
I get
London
St. Mary's Hospital
London
W2 1PG
From: Manuel Sp?nola [mspinol...@gmail.com]
Sent: 19 July 2012 03:26
To: O'Hanlon, Simon J
Cc: r-sig-geo@r-project.org
Subject: Re: [R-sig-Geo] Help with gpclib
Thank you very much Simon,
I get (I don't know why):
library
Hospital
London
W2 1PG
__**__
From: Manuel Sp?nola [mspinol...@gmail.com]
Sent: 19 July 2012 03:26
To: O'Hanlon, Simon J
Cc: r-sig-geo@r-project.org
Subject: Re: [R-sig-Geo] Help with gpclib
Thank you very much Simon,
I get (I don't know why):
library(maptools
On Thu, 19 Jul 2012, hadley wickham wrote:
The chief problem is the totally unnecessary use of gpclib in fortify() in
ggplot2. gpclib should never, ever be used, as its licence only permits
hobby use to learn about polygon clipping - imposed against the author's
will by his former university.
That code is out of date. The current code is
https://github.com/hadley/ggplot2/blob/master/R/fortify-spatial.r
The released code is as quoted, which is what users see. I think github does
not provide automatic check, build, and binary build of development versions
of packages (like
-geo-boun...@r-project.org] on
behalf of Manuel Spínola [mspinol...@gmail.com]
Sent: 15 July 2012 16:02
To: r-sig-geo@r-project.org
Subject: [R-sig-Geo] Help with gpclib
Dear list members,
I am running R 2.15.1 on a MAC OSX 10.7.
I receive an error message regarding gpclib (see below):
library
From: r-sig-geo-boun...@r-project.org [r-sig-geo-boun...@r-project.org]
on behalf of Manuel Spínola [mspinol...@gmail.com]
Sent: 15 July 2012 16:02
To: r-sig-geo@r-project.org
Subject: [R-sig-Geo] Help with gpclib
Dear list members,
I am running R 2.15.1 on a MAC OSX 10.7.
I receive an error
Dear list members,
I am running R 2.15.1 on a MAC OSX 10.7.
I receive an error message regarding gpclib (see below):
library(rgdal)
library(maptools)
library(rgeos)
Loading required package: stringr
rgeos: (SVN revision 330)
GEOS runtime version: 3.3.2-CAPI-1.7.2
Polygon checking: TRUE
Dear R users,
I am trying to write a code that could efficiently determine which voxel a
point in space lies inside. For example, I have a matix with x,y,z values:
library(rgl)
x = runif(500, 0, 10)
y = runif(500, 0, 10)
z = runif(500, 0, 10)
myPoints - matrix(c(x,y,z), ncol=3)
plot3d(myPoints)
101 - 200 of 286 matches
Mail list logo