man/listinfo/r-sig-geo
Assuming the result of 'fit.variogram' is in the variable 'spam', this
will get you the SSErr:
attr(spam, "SSErr")
cheers,
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaan 1
fitting.
>
> http://en.wikipedia.org/wiki/Universal_Transverse_Mercator_coordinate_system
>
> might point you to the appropriate zone.
...in addition, for reprojection of data you can use the spTransform
function from the rgdal package.
regards,
Paul
> With best regards,
>>
>&g
__
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
You can use rgdal with the hdf5 driver installed.
good luck,
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI
oord_equal() + scale_x_continuous("", breaks = NULL) +
scale_y_continuous("", breaks = NULL))
# ...and couple the color to the ffreq:
gg = ggplot(aes(x = x, y = y, size = om, color = ffreq), data = meuse) +
geom_point()
print(gg + coord_equal() + scale_x_continuous(
ailing list
>>> R-sig-Geo@r-project.org
>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>>
>>>
>>
>> --
>> In God we trust, all others bring data.
>>
>
>
>
> ___
> R-sig-Geo mailing
to your .bashrc and restart a terminal. Now it should be possible to
install rgdal without manually specifying the location of the
executables and libraries, because they are in the search path.
cheers,
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute
On 02/27/2012 10:22 AM, Paul Hiemstra wrote:
> On 02/26/2012 08:24 PM, Hodgess, Erin wrote:
>> Hi R Sig Geo People:
>>
>> May I ask an ArcGIS question here, or should I try to find an ArcGIS group,
>> please?
>>
>>
>>
>> Erin M. Hodgess, PhD
&
r, I have a hard time helping you.
Some notes:
- Try fortify(NUTS1, region = "id") on the unaltered NUTS1 object, i.e.
just after calling readShapeSpatial. If this works, the problem is not
in fortify, but in your code.
- There is no need to call fortify.SpatialPolygonsDataFrame e
nd there are a lot of ArcGIS questions there (and also
a lot of know how I assume).
good luck,
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaan 10 | 3732 GK | De Bilt | Kamer B 3.39
P.O. Box 201 | 3730 AE | De Bilt
tel: +31 30 2
the cov. matrix of a 2 GB data file with
almost no memory usage.
cheers,
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaan 10 | 3732 GK | De Bilt | Kamer B 3.39
P.O. Box 201 | 3730 AE | De Bilt
tel: +31 30 2206 494
h
todetected)
>
> I suspect something is going on with the environment variables between the
> Stat-ET launching of R, and the command-line launching. A Sys.getenv()
> reveals that the environment variables are different, but I'm unclear which
> variable(s?) I should be lo
L version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminala
nInfo() command.
Give us more to work with, or it will be very hard to help you.
regards,
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaan 10 | 3732 GK | De Bilt | Kamer B 3.39
P.O. Box 201 | 3730 AE | De Bilt
tel: +31 30 2206
ituto di Fisiologia Clinica
> Pisa Italia
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorolo
units (m/s)^2. For the x-axis they are most probably the distance
unit in which the coordinate system of your data set is, e.g. kilometer.
hope this helps,
Paul
On 11/16/2011 03:08 PM, carolina monmany wrote:
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute
Hi,
Please provide a reprodicble example which shows this problem. This
makes it much easier for us to understand what your problem exactly is
and help you solve it.
good luck,
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaan 10 |
y success.
>
> Thanks
> Javier
> ...
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological In
gt;
> Please advice
>
>
> I want the datframe in a csv file where i can edit this
>
>
> ,also i have to batch process this so cant do many changes
>
>
> manually.
>
>
> Thanks in advance
>
>
> Sandeep
>
> [[alternative HTML version deleted]]
>
please give me more instructions on how to fix this?
> Thanks a lot.
>
> On Mon, Oct 17, 2011 at 3:42 PM, Paul Hiemstra <mailto:paul.hiems...@knmi.nl>> wrote:
>
> On 10/16/2011 02:10 PM, zhijie zhang wrote:
> > cygwin warning:
> > MS-DOS style path detect
gt;
> Sandeep
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
?spDists, ?zerodist
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
nt variable option "nodosfilewarning" turns off this
> warning.
Have you tried this suggestion?
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaan 10 | 3732 GK | De Bilt | Kamer B 3.39
P.O. Box 201 | 3730 AE | De
berg.html
>>
>> ___
>> R-sig-Geo mailing list
>> R-sig-Geo@r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
> [[alternative HTML version deleted]]
>
> ___
> R-
nell. I think the book is ok, especially when teaching to under
graduates. It covers a wide range of topics and does not explicitely use
one type of software (e.g. ArcGIS).
regards,
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaan
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institut
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaan 10 | 3732 GK | De Bilt | Kamer B 3.39
P.
e welcome
>>
>> Thanks
>>
>> Sandeep
>>
>> [[alternative HTML version deleted]]
>>
>> _______
>> R-sig-Geo mailing list
>> R-sig-Geo@r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
P
th programming, but determined to figure this out.
> What exactly do you mean by rebuild? I use the OSGeo4W package to
> update gdal every few months, so I think I have the most recent
> version.
>
> Thanks,
> Ezra
>
>
>
> On Mon, Oct 3, 2011 at 3:21 AM, Paul Hiemstra wrote:
&g
ce will they be willing to exchange a few emails?
>
>
> Thanks
>
> --
> Sandeep R Patil
> Graduate Student
> M.S (Hydrological Sciences)
>
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaan 10 | 3732 GK | De Bilt | K
icular, I need to get forecast data
>> into a format where I can overlay a shapefile, but the wqgrib2 program
>> won't work because hybrid coordinate system.
>>
>> Thanks,
>> Ezra
>
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands
esirable scenario, insert the
> needed "#include " in the given *.c file yourself and
> compile.
> *
> Can anyone make out anything from this , i want to tweak the configure
> command but do not know how to proceed..
>
> Sandeep
>
> [[alternative HTML versi
ntDataFrame coming from grass point vector data,
> who have an irregular boundary and 5 km spaced.
>
> Any hits are very appreced,
> bast regards
>
> Eugenio
>
> [[alternative HTML version deleted]]
>
> _
archive at Nabble.com.
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaa
ight
>>
>> Gianni
>>
>>[[alternative HTML version deleted]]
>>
>> ___
>> R-sig-Geo mailing list
>> R-sig-Geo@r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
> [[alter
provding us with a small piece of code that reproduces your problem
would be of great help here.
regards,
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaan 10 | 3732 GK | De Bilt | Kamer B 3.39
P.O. Box 201 | 3730 AE | De Bilt
il/r-sig-geo/2010-December/010363.html
I sent you the scripts off-list, which includes example code.
good luck,
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaan 10 | 3732 GK | De Bilt | Kamer B 3.39
P.O. Box 201 | 3730 AE | De B
courses on spatial statistics.
They are, however, often in Europe. But as you are willing to travel all
across the North-Americas, a trip to Europe might also be worthwhile...
good luck,
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelmi
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Hi,
Saving the file as a comma separated file (csv) should be fine for
reading in Excel. See the write.csv function.
good luck,
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaan 10
Evolution
>> Stony Brook University
>>
>>[[alternative HTML version deleted]]
>>
>> ___
>> R-sig-Geo mailing list
>> R-sig-Geo@r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
&g
1] ""
>
> $major
> [1] "2"
>
> $minor
> [1] "13.1"
>
> $year
> [1] "2011"
>
> $month
> [1] "07"
>
> $day
> [1] "08"
>
> $`svn rev`
> [1] "56322"
>
> $language
> [1] "R"
upport. Alternatively, find a
(command)-line tool which can translate HDF5 to geotiff or netCDF (e.g.
gdal_translate or ) and then use readGDAL or the ncdf package.
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaan 10 | 3732
orks perfectlyAny ideas?
>
>
>
> regards, m
>
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Paul Hiemstra, Ph.D.
Glo
ite? I can guess the
> query looks like: SELECT * FROM table1 AS a, the_geom FROM table2 AS b
> WHERE a.id = b.id
>
> Dan
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Dear Sweta J.P.,
The function autofitvariogram does not support fitting anisotropic
variograms. The function does do something however. The example below
illustrates what happens (under the hood automap uses fit.variogram to
fit the variogram):
library(automap)
loadMeuse()
# No anisotropy
vgm1
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Hi,
As I understand your (brief) question, take a look at the residuals
generic function.
lm.cars = lm(speed~dist, cars)
residuals(lm.cars)
cheers,
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorological Institute (KN
ongitude-latitude.html>?
>
>
> Thanks in advance, any help would be appreciated
>
> Mathieu
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz
om the R-sig-geo mailing list archive at Nabble.com.
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Meteorol
alues)?
> The color table is the same, but the scaling is independent for each
> raster object and cannot be compared.
>
> Thanks
>
> Agus
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailm
Take a look at any of the plotting packages in R: the standard plotting,
lattice or ggplot (I prefer the latter). Using the standard plotting
facilities:
library(gstat)
loadMeuse()
plot(meuse.grid$dist, meuse.grid$soil)
cheers,
Paul
--
Paul Hiemstra, Ph.D.
Global Climate Division
Royal Netherlands Me
try:
spplot(meuse.grid, sp.layout=list("sp.polygons", meuse.riv,
first = FALSE))
I added this now to the spplot documentation.
Thanks!
Tom
--
Paul Hiemstra, Ph.D.
Department of Physical Geography
Faculty of Geosciences
University of Utrecht
Heidelberglaan 2
P.O. Box 8
[[alternative HTML version deleted]]
>
>
>
> _______
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Paul Hiemstra, MSc
Global Climate Division
Royal Netherlands Meteorologica
On 04/08/2011 08:34 AM, Paul Hiemstra wrote:
> On 04/07/2011 11:26 PM, Narayani Barve wrote:
>> Hi All,
>>
>> I am using raster package for various operations. It is a useful package,
>> but manual / pdf contains lines which are broken and not complete. This
>> ha
ou can read the documentation online using
R itself, which does show all the text.
cheers,
Paul
--
Paul Hiemstra, MSc
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaan 10 | 3732 GK | De Bilt | Kamer B 3.39
P.O. Box 201 | 3730 AE | De Bilt
tel: +31 30 2206
l
>
> I see it's possible to calculate the residuals with
>
> residuals(lm(zinc~sqrt(dist), meuse))
>
> and pass them directly. But I wonder whether a more convenient direct
> passing of the trend coefficients is possible.
>
> Thanks!
> Tom
>
--
Paul Hiemstra, MS
On 03/28/2011 10:47 AM, Mauricio Zambrano wrote:
2011/3/26 Paul Hiemstra:
Hi,
The miscFitOptions parameter does not include the 'variogram' paramters. I
will look into adding them, together with the option to set the boundaries
manually. However, autofitVariogram at this stage
27;cutoff'
passed to the 'autofitVariogram' function.
In any case, I think that even after solving the issue related to
'cutoff' and 'boundaries', it will still be impossible to pass the
'cutoff' argument to the 'autoKrige' function, because the
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Paul Hiemstra, MSc
Global Climate Division
Royal Netherlands Meteorological Institute (KNMI)
Wilhelminalaan 10 | 3732 GK | De Bilt | Kamer B 3.39
P.O. Box
I think I know
>
> If I don't know I know
> I think I don't know
>
> -- Knots
> R.D. Laing
>
> [[alternative HTML version deleted]]
>
> __
or even faster.
cheers,
Paul
--
Paul Hiemstra, MSc
Department of Physical Geography
Faculty of Geosciences
University of Utrecht
Heidelberglaan 2
P.O. Box 80.115
3508 TC Utrecht
Phone: +3130 253 5773
http://intamap.geo.uu.nl/~paul
http://nl.linkedin.com/pub/paul-hiemstra/20/30b/770
c
Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Hi TM,
See the geoR package.
Paul
--
Paul Hiemstra, MSc
Department of Physical Geography
Faculty of Geosciences
University of Utrecht
Heidelberglaan 2
P.O. Box 80.115
3508 TC Utrecht
Phone: +3130 253 5773
http://intamap.geo.uu.nl/
x27;x'
and 'y' axis. Since the projection I'm trying to use dictates the ratio
between the x and y axis, I don't understand how this can be.
Any suggestion on how to proceed would be greatly appreciated..
Thanks,
Elena
--
Paul Hiemstra, MSc
Department of Physical Geogra
ger Bivand:
On Wed, 15 Dec 2010, Paul Hiemstra wrote:
Hi people,
I posted a similar question to the ggplot2 mailing list and with their
help and a lot of tinkering I got a well working function to add a scalebar
to a ggplot plot. I could add the function to automap, but is there anoth
On 12/15/2010 09:53 AM, Roger Bivand wrote:
On Wed, 15 Dec 2010, Paul Hiemstra wrote:
Hi people,
I posted a similar question to the ggplot2 mailing list and with
their help and a lot of tinkering I got a well working function to
add a scalebar to a ggplot plot. I could add the function to
, size = 6, data =
textScaleBar, hjust = 0.5, vjust = 1.2, legend = FALSE))
}
library(ggplot2)
library(sp)
data(meuse)
data(meuse.grid)
ggobj = ggplot(aes(x = x, y = y, color = zinc), data = meuse) + geom_point()
# Make sure to increase the graphic device a bit
addScaleBar(ggobj, meuse, &quo
tcome of the krige function.
#' @author Paul Hiemstra, \email{p.h.hiems...@gmail.com}
#' @export
#' @examples
#' data(meuse)
#' coordinates(meuse) = ~x+y
#' data(meuse.grid)
#' gridded(meuse.grid) = ~x+y
#'
#' meuse_tps = doTps(zinc~dist, meuse, meuse.gr
listinfo/r-sig-geo
--
Paul Hiemstra, MSc
Department of Physical Geography
Faculty of Geosciences
University of Utrecht
Heidelberglaan 2
P.O. Box 80.115
3508 TC Utrecht
Phone: +3130 253 5773
http://intamap.geo.uu.nl/~paul
http://nl.linkedin.com/pub/paul-hiemstra/20/30b/770
65 matches
Mail list logo