search for: topos

R- Users: Can someone indicate what I am during wrong? This is a script essentially from Venerable & Ripley's text on interpolating a surface with loess function but I can not get it to run. Thanks. John Lewis Professor McGill University Montreal library(MASS) library(modreg) data(topo) par(mfcol=c(2,2), pty="s") topo.loess <- loess(z ~ x * y, topo, degree=2,

I'm reading Ripley-Venables "Modern Applied Statistics with S - Fourth edition" , at the same time trying the examples proposed in the book using R 1.8.1 under linux. Now I'm trying the following code from the book (example code of spatial statistics at page 76) with R : | data(topo) library("spatial") topo.loess<-loess(z ~ x * y, topo, | degree= 2,span=0.25) |

Hi Folks, I'm probably overlooking some small point, but can't see it. Trying the "Local trend surfaces" example from p.437 of MASS (3rd edn) by V&R, all goes well until the line contour(topo.mar$x, topo.mar$y, topo.lo$fit, levels=seq(700,1000,25), xlab="fi t", ylab="") which produces the response: Error in contour.default(topo.mar$x,

Hi, I was wondering if anyone knows of a way by which one can estimate values from a contour plot created by using the loess function? I am hoping to use the loess contour plot as a means of interpolation to identify the loess created values at points at pre-defined (x,y) locations. Could anyone point me in the right direction please? Thanks. Laura Quinn Institute of Atmospheric Science School

Author: mikechr Repository: /hg/zfs-crypto/gate Revision: 3a473014cdbb72a0ecad997583eaabfe2d32283c Log message: 6345621 FMA topo support for Boston/Seattle Files: create: usr/src/cmd/fm/topo/files/sparc/SUNW,Sun-Fire-V215/Makefile create: usr/src/cmd/fm/topo/files/sparc/SUNW,Sun-Fire-V245/Makefile create: usr/src/cmd/fm/topo/files/sparc/SUNW,Sun-Fire-V245/platform.topo create:

Although the 'spatial' documentation doesn't mention that extractAIC works, it does seem to give an output. I may have misunderstood, but shouldn't the following give at least the same d.f.? > library(spatial) > data(topo, package="MASS") > extractAIC(surf.ls(2, topo)) [1] 46.0000 437.5059 > extractAIC(lm(z ~ x+I(x^2)+y+I(y^2)+x:y, topo)) [1]

Dear useRs, I got stuck trying to generate a palette of topographic colors that would satisfy these two requirements: - the pallete must be 'anchored' at 0 (just like on a map), with light blue/lawn green corresponding to data values close to 0 (dark blue to light blue for negative values, green-yellow-brown for positive values) - the brown must get darker for higher positive

Hello all, I have two questions about anova (one is probably VERY basic...) 1 - when one asks for a summary of a trend surface created with surf.ls, he/she gets: > summary(g3r) Analysis of Variance Table Model: surf.ls(np = 3, x = gradiente$east, y = gradiente$north, z = gradiente$num1) Sum Sq Df Mean Sq F value Pr(>F) Regression 215.7182 9 23.968693976

Can anyone please hazard a guess at an estimate for a poor old squirrel hunter with a Very Very Dumb Question? Is wine likely to be able to install and run any of the proprietary GPS-interfacing topo map programs (DeLorme, Garmin, Mapquest, et al.) at any roughly predictable future date? The season opens in three weeks, my GPS (Garmin's etrex vista) is full -- and it still contains road maps

I will second Roger's suggestion, colorRampPalette is a great function for creating your own palettes. For example, Matlab's jet palette (also available in fields package under peculiar name 'tim.colors') can be defined by: jet.colors = colorRampPalette(c("#00007F", "blue", "#007FFF", "cyan", "#7FFF7F", "yellow",

Hello all, The subject line is a little ambitious, although that is my ultimate goal, to be able to run the National Geographic TOPO! program to look at maps of Arizona and California on my laptop using only Linux. I don't even own a copy of Windows, nor do I desire to, so hopefully Wine will be my friend. So far it isn't though. I haven't even got as far as trying to run TOPO! I

If a run a LOESS model and then produce a smoothed surface: Is there any way to determine the coordinates of the local maxima on the surface? [[alternative HTML version deleted]]

I'm trying to find a way to fit a polynomial of degree n in x and y to a set of x, y, and z data that I have and obtain the coefficients for the terms of the fitted polynomial. However, when I try to use the surf.ls function I'm getting odd results. > x <- seq(0, 10, length=50) > y <- x > f <- function (x, y) {x^2 + y} > library(spatial) > test <-

Hi everyone. I have a vector like this: x = c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2) I would like to associate a colour to each of these entry. For example, I could use /topo.colors(3)/ (since I have 3 groups). I know it is easy to do this with /if/else/ statements, but I would like to find a proper way to do this. Thanks in advance, Phil -- View

Dear R-helpers! What it's means "Rank failure in Choleski decomposition" by using of function surf.gls {spatial}? Sincerely Yours Branimir K. Hackenberger

Can any of it talk to a GPS yet?? I have suites from Garmin, Maptech, Topo.com, and Topo (as sold through National Geographic -- the least satisfactory of all). Last year at about this time, it still took massive effort, time, and help just to get any of them to launch. Now a couple at least launch, and even run, quite nicely, just in plain wine -- except in the most important respect of all.

I have a huge dataset(15k X 18) and tried to use the heatmap in R to examine the patterns. However, I found that heatmap and heatmap.2 gave me completely different outputs. Here are the codes: ------------ > dim(as.matrix(data.dcpm)) [1] 15462??? 18 > > heatmap(as.matrix(data.dcpm), col=topo.colors(100)) > heatmap.2(as.matrix(data.dcpm), col=topo.colors(100), key=TRUE,

Hi, I would have expected that I get always 3 colors as result which is not true: hcl.colors(3, alpha=c(0, 0.5, 1)) # 3 colors rainbow(3, alpha=c(0, 0.5, 1)) # 3 colors heat.colors(3, alpha=c(0, 0.5, 1)) # 3 colors terrain.colors(3, alpha=c(0, 0.5, 1)) # 6 colors cm.colors(3, alpha=c(0, 0.5, 1)) # 6 colors topo.colors(3, alpha=c(0, 0.5, 1)) # 9 colors R-Version and

Dear list members I'm currently working on some topographic (elevation) data, and was somewhat surprised that the 'topo.colors' and 'terrain.colors' are of little to no use here. The problem is that these functions only return a palette of colours; they don't map depth values to colours. So if I plot (using 'image', 'persp' or similar functions) and

Hi, this is possibly a bug report (I'm not sure if this is a feature). It's related to the --keep-dirlinks option, when combined with --delete . I have the following directory structure on server A: ls -lR software software: total 238 drwxr-xr-x 2 biolord bioinf 1024 Sep 20 10:49 EMBOSS/ lrwxrwxrwx 1 biolord bioinf 6 Feb 5 2003 MSE -> EMBOSS/ lrwxrwxrwx 1