Displaying 20 results from an estimated 500 matches similar to: "plot two variograms on a same graph"
2006 Nov 24
2
plot() and Jpeg() increase font size and resolution
Dear list,
I am having troubles increasing the fontize when plotting a
variogram{gstat} and its model (vgm) with plot and using jpeg(). Also
the resolution in the jpeg call does not work. I am using R2.4 on
Mandriva 10.2 linux.
I can change fontsize with cex.axis in a normal plot, so I presume it
has to do with plotting the variogram model. Any help on how to increase
the font size and
2006 Aug 20
3
plot problem
Hello.
I'm pretty much new to R and I'm trying to produce some figures.
It seems to me, that R has some asynchronous way of plotting figures.
When I run this code:
#constructs the semivariogram of SC1929
vgm1 <- variogram(SC1929~1,~U+V,puerto.map$att.data)
# trying to make new plot
dev.set(which=dev.next())
plot(vgm1)
title(main="Semivariogram",font.main=4)
2009 Mar 24
1
Variogram with Gstat
Dear all can you help me? i have this problem: i have a dataset in a text file
in a matrix of 3 columns: x, y, z where x and y are the coordinates and Z are
the mesurements. How can i do a variogram with R?
2008 Apr 01
1
spatial cross-correlation
Hi;
I cannot find in the R html documentation a way to evaluate
cross-correlation in 2D data sets.
I would like to evaluate cross-correlation in a series of moving windows
between two maps.
i,e, specify several windows inside the complete 2D spatial matrixes and
for each one ofthese windows evaluate the 2D cross-correlation (commonly
conducted in the spectral domain).
Thanks in advance and best
2009 Sep 08
1
gstat---2 basic plot questions
Hi all--
I'm new to R, statistics and programming, so sorry if this is a really basic
question!
I have plotted a directional variogram, and I want to
a. overlay the omni-directional line over each directional panel
b. display the directional variograms in a single panel with a legend that
associated each line to each degree measurement.
The line I'm using is
2009 Aug 06
2
random between two values
Hi,
I would obtain a random value between two (for example between 40.15 and 56.58 I would have only one value).
I'm looking for a package/a function which could do this.
Could anybody help me please?
Cordialement
Damien Landais
2007 Jan 05
1
gstat package. "singular" attibute
Hello,
I'm using the gstat package within R for an automated procedure that
uses ordinary kriging.
I can see that there is a logical ("singular") atrtibute of some
adjusted model semivariograms:
.- attr(*, "singular")= logi TRUE
I cannot find documentation about the exact meaning and the implications
of this attribute, and I dont know anything about the inner calculations
2009 Jun 26
3
to raise in a loop more than 1
I would raise x,y and z in a loop but I won't raise of 1. I tried this but it doesn't work
mydata=matrix(nrow=1500,ncol=3)
i=1
for(x in 0:10){
for(y in 0:20){
for(z in 0:10){
mydata[i,]=c(x,y,z)
i=i+1
z=z+2}
y=y+4}
x=x+2}
And I would have something like that
x y z
0 0 0
0 0 2
0 0 4
0 0 6
0 0 8
0 0 10
0 4 0
0 4 2
...
Could anybody help me?
I
2006 Apr 25
3
persp plot increasing 'x' and 'y' values expected
hello,
i do the following in order to get an persp-plot
x<-c(2,2,2,2,2,2,3,3,3,3)
y<-c(41,41,83,83,124,166,208,208,208,208)
z<-c(90366,90366,92240,92240,92240,96473,100995,100995,100995,100995)
x<-data$x
y<-data$y
z<-matrix(data$z,length(y),length(x))
persp(x,y,z, col="gray")
but i always get the error message increasing 'x' and 'y' values
expected,
2012 Sep 12
2
valgrind crashing
I am trying to do a classic
R -d valgrind --vanilla < mypkg-Ex.R
as described in
http://cs.swan.ac.uk/~csoliver/ok-sat-library/internet_html/doc/doc/R/2.9.1/doc/manual/R-exts.html#Using-valgrind
The problem is valgrind crashes imediately.
I am using Ubuntu 10.04 LST (actually Biolinux) on a 64bit Dell
Precision laptop with 4 core-i5 Intel processors. I have R-2-15-0
installed from source
2015 Aug 04
2
Duda interpolación (package ' gstat ')
Hola,
# Hacemos el KED. Ver función "krige()":
KED.rad <- krige(
formula=pluvPcp~layer, # covariable -> radar
locations=lluvia.rad.pluv.spdf,
newdata=radarGrid, # podría ser cualquier objeto Spatial
model=v.fit, # modelo de semivariograma.
maxdist=Inf
2017 Jul 26
5
[PATCH 0/2] daemon: Reimplement handling of lvm.conf and filters.
Simplify how we handle lvm.conf.
2012 May 02
3
Consulta gráfica
Hola,
Por favor, ¿podríais indicarme qué recursos (librerías o ideas) pueden resultar de utilidad para crear un gráfico del estilo del de la figura 3.8 del siguiente link?
http://www.tsc.uvigo.es/BIO/Bioing/ChrLDoc3.html#3.5
Actualmente estoy utilizando funciones muy básicas y la verdad es que no me encuentro muy satisfecha con el resultado.
Muchas gracias.
Eva
[[alternative HTML
2011 Mar 04
2
make an own (different) color legend with spplot()
Hi!
Is there a way to manually costumize the color legend in an spplot() -
especially where to draw ticks and labels for the ticks?
The reason I'm asking: Usually spplot() automatically divides the data
into fitting slices and makes a color legend (also automatically).
I want to assign the slices myself and have a fixed scale instead of an
automatic/dynamic scale.
I think what I want
2003 Mar 07
1
REML option in gstat
Hi, please help!!
I've been trying to fit variogram models using the REML method in the gstat
package. Every time the Windows GUI crashes. For example
library(gstat)
data(meuse)
x <- variogram(zinc ~ 1, ~x + y, meuse)
v <- vgm(140000, "Sph", 800, nug = 10000)
plot(x, model = fit.variogram(x, model = v, fit.method=5))
Other fit methods are non problematic (eg. fit.method=7
2011 Mar 05
3
Change panel background color in spplot()
Hi!
How does one change the background color of the map-panel in spplot()?
Example:
library(sp)
data(meuse.grid)
gridded(meuse.grid) = ~x+y
spplot(meuse.grid, "part.a")
How would I get another background-color for the map-panel (but not for
the whole plot) here?
Thank you!
Marcel
2014 Oct 09
1
live migration to localhost with lvm based storage
I have lvm based (thin pool) storage on local disks. I need to move
vps from one vg on one disk to another. Does it possible to migrate to
localhost with blockcopy migration to another vg?
I'm understand that i can move lv from one vg to another, but i don't need that.
--
Vasiliy Tolstov,
e-mail: v.tolstov@selfip.ru
jabber: vase@selfip.ru
2002 Apr 15
2
krige and polygon limit problem
Dear all,
I'm new on R and this mailing list. We work on spatial rainfall estimation
with R and Grass.
We have a problem with the krige function from the sgeostat package. We would
like to limit the estimated area with a polygon limit.
I use a 50 points polygon to describe my work area. The krige function work
quiet well without limit. But if I use this option I have the following error
2003 Jun 26
1
krige in gstat() package
HI,
I wonder does anyone have experience with doing sequential gaussian
simulation with krige() function in gstat?
I find it VERY slow compared to use krige() to achieve kriging function
itself.. I wonder why, is that because it has to model the variogram, and
do the kriging separately for each point to be simulated?
so it would be N times slower to achieve the simulation than the kriging
if
2002 Dec 17
2
Cross-correlograms or cross-variograms in R?
Hello group,
For my PhD I'm working on a spatial sampling grid. I do have two data sets
which I'd like to compare using cross-correlograms or cross-variograms.
Is this an option in one of the R-packages? I've been searching the R-help
archive and the available package-documentations, but I can't find how to do
this.
Thanks in advance,
Ren?.