search for: topo

...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, span = 0.25, normalize=F) topo.mar <- list(x = seq(0, 6.5, 0.1), y=seq(0, 6.5, 0.1)) topo.lo <- predict(topo.loess, expand.grid(topo.mar), se=T) lo.1 <- as.data.frame(topo.lo) lo.2 <- as.matrix...

...#39;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) | topo.mar<-list(x=seq(0,6.5,0.2),y=seq(0,6.5,0.2)) | topo.lo=predict(topo.loess,expand.grid(topo.mar)) | par(pty="s") | contour(topo.mar$x,topo.mar$y,topo.lo,xlab="",ylab=""...

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, topo.mar$y, topo.lo$fit, levels = seq(700, : no proper `z' matrix specified I've checked the sizes...

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: usr/src/cmd/fm/topo/files/sparc/SUNW,Sun-Fi...

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] 6.0000 357.5059 # and if the AIC values differ, shouldn't they do so by an additive constant? > (extractAIC(surf.ls(2, topo))-extractAIC(lm(z ~...

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 mus...

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 instead of topo maps, because I can't upload the good ones from linux (FC1). I'm getting alm...

...FF", "cyan", "#7FFF7F", "yellow", "#FF7F00", "red", "#7F0000")) Other predefined functions for creating and managing color palettes that I know of, are: * R provides functions for creating palettes of continuous colors: rainbow, topo.colors, heat.colors, terrain.colors.colors, gray * tim.colors in package fields contains palette similar to Matlab's jet palette (see examples for simpler implementation) * rich.colors in package gplots contains two palettes of continuous colors. * Functions brewer.pal from RColorBr...

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 set up an intermediate goal of run...

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]]

...(0, 10, length=50) > y <- x > f <- function (x, y) {x^2 + y} > library(spatial) > test <- data.frame(x=x, y=y, z=f(x, y)) > test.kr <- surf.ls(2, test) > test.kr$beta [1] 0 0 0 0 0 0 When I try the example from the help I get: > library(MASS) > data(topo, package="MASS") > topo.kr <- surf.ls(2, topo) > topo.kr$beta [1] 801.217617 -11.018887 68.229148 -73.992968 3.343573 8.342717 Why is my test data causing problems? Also it seems that the beta attribute from the object returned the surf.ls correspond with the terms of...

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 this message in context: http://r.789695.n4.nabble.com/Colors-vector-based-on-group-tp4361425p4361425.html Sent from th...

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 im...

...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, density.inf="none",trace="none", scale="none") ------- The outputs are attached here. Could anyone help me figure out why? Thanks a lot:)

...et 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 platform: R version 4.2.2 Patched (2022-11-10 r83330) -- "Innocent and Trusting" Copyright (C) 2022 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) Bug or feature? Sigbert -- https://hu.berl...

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', &...

...R 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 biolord bioinf 6 Feb 5 2003 PHYLIP -> EMBOSS/ ... lrwxrwxrwx 1 biolord bioinf 6 Feb 5 2003 TOPO -> EMBOSS/ ... software/EMBOSS: total 87942 -rw-r--r-- 1 biolord bioinf 551090 Jul 15 18:54 DOMAINATRIX-1.0.0.tar.gz -rw-r--r-- 1 biolord bioinf 12100223 Jan 21 2003 EMBOSS-2.6.0.tar.gz -rw-r--r-- 1 biolord bioinf 14024961 Nov 30 2003 EMBOSS-2.8.0.tar.gz -rw-r--r-- 1 biolor...