similar to: polynomial fitting

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]

Hello R users, Can someone point me to a solution to the problem below? The surf.ls() call works fine in other contexts, e.g. the example in help(surf.ls). I'm using R 1.3.1, and get the same error regardless of the platform (debian Linux, MS Windows 98). Many thanks, Stewart > library(spatial) > topo.kr <- surf.ls(2,midxy) Error in surf.ls(2, midxy) : NA/NaN/Inf in foreign

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

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,

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

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

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,

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

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

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

Hi guys, After all the discussion yesterday about persp and color, I decided to have a more closer look at demo(persp), and decided to write a function to generate 'topo-like' colours to plot perspectives (Thanks a lot to Uwe Ligges for his enlightning comments regarding the code in the demo). Here it goes, I believe that this function will be pretty useful to a lot of people:

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

Am Dienstag, 25. April 2017, 11:21:31 schrieb Dirk Eddelbuettel: > On 25 April 2017 at 16:11, Johannes Ranke wrote: > | This looks similar to what I got this morning when I tested my > | (unreleased) > | backport of R 3.4.0 to Debian jessie. My test was > | > | library(MASS) > | example(rlm) > | > | and there was an object that was not found. I am on a train on the way

Hi R-List, I have a question regarding R-language formats, I think. I am producing a series of graphs (using plot, barplot, barchart, and bwplot, using either text or mtext to place values on the graphs) and tables for a Francophone country. In fact, I have already done so. However, while they are pleased with the results they've requested I convert all of my decimal points into the French

Hello, (maybe a quite basic statistical question, but I'm just struggling with it) I'm doing an anova: > Res1<-aov(H2O~location+topo+horizon+pF+Error(location:topo:horizon)) (water retention of soils) which gives a significant difference at factor "location". Which function should I use to now which locations (there are 3: A, B and C) differ significantly and which do

With snarls and curses and gnashing of teeth, I still keep one laptop and one dual-bootable PC hard drive with XPProSP2 installed. (Everything else is Fedora 10.) I do so purely for one class of apps, proprietary map software. Wine can now install and run at least one, and the most important, of the four suites of such apps that I have. (That one is Garmin MapSource, including Topo US

Dear R people, This is probably a FAQ, and a silly question besides, but... a lot of R packages have datasets contained within them. cf. spatial has "topo". Once I have loded the appropriate library, how do I access the datasets contained within? I know this must be very simple, but I have spent some (fruitless, frustrated time) trying to do this, so am asking for help.

Am Dienstag, 25. April 2017, 08:50:34 schrieb Dirk Eddelbuettel: > On 25 April 2017 at 14:58, G?ran Brostr?m wrote: > | hello, > | > | I just installed R-3.4.0 from scratch: > | > | $ sudo apt install r-base > | > | but when I try > | > | > library(survival, lib.loc = "/usr/lib/R/library") > | > fit <- coxph(Surv(exit, event) ~ x, data =