Displaying 10 results from an estimated 10 matches for "grid2".
Did you mean:
grid
2002 Apr 16
0
still have problem with krige and border option
...nd) :
user>border_limite<-read.table("/home/lionel/rwork/poly_test.txt",header=FALSE)
user>polygone<-list(x=border_limite[,1],y=border_limite[,2])
So now I have a list of x,y coordinates to limit my work area.
If I call the function krige with extrap=FALSE like this
user>grid2.krige<-krige(grid2.point,logtab.point,'lz',
logtab.Exp,maxdist=NULL,extrap=FALSE)
[ where grid2.point is a point object to store prediction, logtab.point a
point object containing points data, 'lz' is the variable name in
logtab.point, logtab.Exp is a variogram object ]
all...
2009 Nov 15
1
R crashing
...URL = getURL(paste("http://maps.google.com/maps?saddr=
",x,sep=""))
y = strsplit(URL, "<div><b>")
y = strsplit(y[[1]][2], " mi</b>" )[[1]][1]
as.numeric(y)
}
}
grid2=grid[!is.na(grid)]
n = length(grid2)
for (i in 1:n) {
temp = Distances(grid2[i])
write.table(temp,"distances.csv",col.names=F,row.names=F,append=T)
}
##########################################
But R still crashes after 2 hours (all I got was around 20.000 distances).
It doesn't real...
2005 Apr 07
3
4D Plot ??
...lplot, so I duplicated all rows of the data frame and changed the
z data for the second half, in order to plot one half at a time.
#----------------------------------------------------------
# Try at a 4D contourplot:
y = x = 1:50
grid <- expand.grid(x=x, y=y)
grid$z = sqrt(x*y)
n1 = nrow(grid)
grid2 = rbind(grid,grid)
grid2$z[(n1+1):(n1*2)] = log(grid2$x[1:n1] * grid2$y[1:n1] + 10)
panel.4d <- function(x,y,z,subscripts) {
n1 = 1; n2 = length(x)/2
panel.levelplot(x[n1:n2],y[n1:n2],z[n1:n2],subscripts,region=TRUE)
n1 = n2 + 1 ; n2 = length(x)
panel.levelplot(x[n1:n2],y[n1:n2],z[n1...
2002 Apr 17
0
still have problem with krige and border option [end]
> both "_ "and "<-" are assignment operators
well, I've changed it in my script, but the result is the same.
> what do you get if you call
>
> in.polygon(grid2.point$x,grid2.point$y,polygon$x,polygon$y)
That's it !
In fact all my points in grid2.point where xy inverted. And I had no points
in polygon.
It was really a stupid mistake. Thanks for all !
Lionel
--
CNERTA-ENESAD
4 rue champs-prevois
batiment grand-champs
21000 Dijon
tel:03.80.77.28....
2007 Dec 04
1
2/3d interpolation from a regular grid to another regular grid
Hello R users,
I have numerical data sampled on two grids, each one shifted by 0.5
from the other.
For example:
grid1 = expand.grid(x=0:3, y=0.5:2.5)
grid2 = expand.grid(x=0.5:2.5, y=0:3)
gridFinal = expand.grid(x=0.5:2.5, y=0.5:2.5)
plot(gridFinal, xlim=c(0,3), ylim=c(0,3), col="black", pch=19)
points(grid1, xlim=c(0,3), ylim=c(0,3), col="red", pch=19)
points(grid2, xlim=c(0,3), ylim=c(0,3), col="blue", pch=19)
I would...
2010 Jul 30
2
Layering multiple persp plots
Hi all,
I've got two persp plots with Identical X and Y's, and I'd like to plot them
on the same graph, so that it is obvious where one plot is above the other.
I can't find any mention of this anywhere. Do I need to use wireframe?
Any help is appreciated.
Thanks,
--
Ian Bentley
M.Sc. Candidate
Queen's University
Kingston, Ontario
[[alternative HTML version deleted]]
2012 Dec 27
0
Suggestion: 'method' slot for expand.grid() (incl. diffs)
...ry
little changes (100% backward compatible) and which allows one to control this
order. Please find attached diffs against R-devel.
Cheers,
Marius
### ./src/library/base/R/expand.grid.R #########################################
--- expand.grid.R 2012-12-27 22:37:29.000000000 +0100
+++ expand.grid2.R 2012-12-27 22:41:00.331979950 +0100
@@ -16,7 +16,8 @@
# A copy of the GNU General Public License is available at
# http://www.r-project.org/Licenses/
-expand.grid <- function(..., KEEP.OUT.ATTRS = TRUE, stringsAsFactors = TRUE)
+expand.grid <- function(..., KEEP.OUT.ATTRS = TRUE, strin...
2011 Jun 22
1
lme convergence failure within a loop
...ross a spatial grid with data points collected at various
time points within each grid cell. I'm trying to use a 'for' loop to try
the model fit on each grid cell. In some cells lme does not converge,
giving me the error:
Error message: In lme.formula(logarea ~ year + summ_d, data = grid2,
random = ~year + :
nlminb problem, convergence error code = 1
message = iteration limit reached without convergence (9)
When I get the error, the program aborts, and my loop stops.
I'm not too worried about the error, as a generic mixed model structure
may not be the best fit for ev...
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
2009 Sep 23
2
[PATCH node] Introduces the virtual network administration functions.
...ecting To Physical Network"),
+ grid]
+
+ def get_summary_page(self, screen):
+ grid1 = Grid(2, 1)
+ grid1.setField(Label("Network name:"), 0, 0, anchorRight = 1)
+ grid1.setField(Label(self.__config.get_name()), 1, 0, anchorLeft = 1)
+
+ grid2 = Grid(2, 3)
+ grid2.setField(Label("Network:"), 0, 0, anchorRight = 1)
+ grid2.setField(Label(self.__config.get_ipv4_address()), 1, 0, anchorLeft = 1)
+ grid2.setField(Label("Gateway:"), 0, 1, anchorRight = 1)
+ grid2.setField(Label(self.__config.get...