similar to: add non-linear line

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

I am trying to fit a curve to a cumulative mortality curve (logistic) where y is the cumulative proportion of mortalities, and t is the time in hours (see below). Asym. at 0 and 1 > y [1] 0.00000000 0.04853859 0.08303777 0.15201970 0.40995074 0.46444992 0.62862069 0.95885057 1.00000000 [10] 1.00000000 1.00000000 > t [1] 0 13 20 24 37 42 48 61 72 86 90 I tried to find starting values for

Red Hat Enterprise Linux Server release 6.5 (Santiago) selinux is disabled. The following commands were all run on the RedHat Server on which I am running samba. *The following ports are open* 5 ACCEPT tcp -- 0.0.0.0/0 0.0.0.0/0 state NEW tcp dpt:137 6 ACCEPT tcp -- 0.0.0.0/0 0.0.0.0/0 state NEW tcp dpt:138 7 ACCEPT tcp --

Hello all   I have a geology map that has three level, bellow   <-geology lithology    landscape   landform     landform level is used as covariate (with codes=1,2,3,4,5) for training of neural network, but this level has missing data as NA. I want to replace the missing data of landform level with 0 (zero). Finally, landform will have codes  0,1,2,3,4,5.   please help me   Thanks alot.

Hi, I am interested in ggplot2 and I found this lattice code very interesting (http://addictedtor.free.fr/graphiques/graphcode.php?graph=48). Code: library(lattice) lattice.options(default.theme = canonical.theme(color = FALSE)) tmp <- expand.grid(geology = c("Sand","Clay","Silt","Rock"), species = c("ArisDiff",

Dear R users, If I manage to sort out this graph, it is certainly a candidate for the new R graph gallery (http://addictedtor.free.fr/graphiques/displayGallery.php)! I created the following lattice graph: library(lattice) tmp <- expand.grid(geology = c("Sand","Clay","Silt","Rock"), species =

Dear All   I want to import the vector file (   .shp) to R. I could import the file by rgdal package before, by following:   geology<-readOGR('C:/geology//saga/geo.geom','finalgeology')   but now there is an error:   Error in ogrInfo(dsn = dsn, layer = layer, input_field_name_encoding = input_field_name_encoding) :           GDAL Error 4: .shx file is unreadable, or corrupt.

Can someone point me to an implementation of the "runs up and runs down" test, or does such beast exist in R? From web searches the "runs up and runs down" test is commonly used for testing pseudo-random number generators and in simulations. John C. Davis describes its use in geology in his book on geological statistics (2002). Having searched the full site, I can locate only

Hi, While doing the ADF test in R using the following command I am getting the error and the result.."> x.ct=ur.df(rev$REVENUE,start=1,end=length(rev$REVENUE),frequency=1) Error in ur.df(rev$REVENUE, start = 1, end = length(rev$REVENUE), frequency = 1) : unused argument(s) (start = 1, end = 4, frequency = 1) >

Hi, I have a 3-d geology problem. I have an xyz fault data set. I am able to view the data set in R using plot3-d but cannot create a raster from it. I need to create a raster and export so that it can be read in arcscene. Can anybody help? Thanks -- Shane [[alternative HTML version deleted]]

Is there any way to maintain spaces, slashes, and parentheses in variable names when reading these into R? Of course, read.table converts these to periods. However, I know that it's not strictly illegal to have these characters in variable names as I am able to add them using the "variable editor" portion of the "data editor." I need to batch produce dozens of histograms

Dear R graphics gurus, Another question about lattice graphics. This time I would like to plot means and confidence intervals by group factor in a lattice graph. I can not find any working lattice examples. Maybe a custom panel function is the answer, but that is a bit beyond me for now. The individual plots within the lattice graph could look like this: # Example with confidence intervals

Hi, I would like to put the sample number beside each lable in a boxplot. How do I do this? Essentially, I need to count the sample size for each factor, see below: Thanks boxplot(DATA$K_Merge~factor(DATA$UnitName_1),axes=FALSE,col=colours) title(main=list("Tukey Boxplot by Geology:\n K(%)",cex=cexlb)) axis(1, 1:21, labels=FALSE, las=2) text(seq(1, 21, by=1), par("usr")[3],

Dear R users, I face to a nested pattern and despite the numerous examples in the help I am still confused. I sampled bugs in different habitats within sites which were within rivers themselves within different regions. The habitat correspond to different substrata (not systematically present in all sites). For rivers and sites, I have environemental variables (e.g. altitude and slope of

Dear list users How is it possible to visualise both a linear trend line and a quadratic trend line on a plot of two variables? Here my almost working exsample. data(Duncan) attach(Duncan) plot(prestige ~ income) abline(lm(prestige ~ income), col=2, lwd=2) Now I would like to add yet another trend line, but this time a quadratic one. So I have two trend lines. One linear trend line

Dear R-project users I would like to calculate a linear trend versus time taking into account a first order autoregressive process of a single time series (e.g. data$S80 in the following example) using th gls function. gls(S80 ~ tt,data=data,corAR1(value, form, fixed)) My question is what number to set in the position of value within corAR1? Should it be the acf at lag 1? I look forward for

Hi I am struggling to fit a non-linear trend using the likfit function in geoR. Specifically I want a sigmoidal function, something like SSfpl in the nls package to fit the trend. But it seems trend.spatial in geoR only works with lm or glm type models. Any ideas how I can specify the model to calculate the kriging parameters using REML, including the parameters of a sigmoidal trend function

I am trying to use the levelplot function from a command file. Here is the code: library(sp) library(gstat) library(lattice) gatherData <- read.table("~/gather.txt", header = TRUE) grd = makegrid(gatherData$x, gatherData$y, cell.size = 5) k <- krige(z~x+y, ~x+y, data = gatherData, newdata = grd, nmax = 5) levelplot(var1.pred~x+y, k, aspect = mapasp(k), main = "Predicted

I have come across three problems in the past few days, in spell of heavy R ( version R0.50-a4/Sun Solaris2.5.1 ) use. 1. I was using lwd=2 to get thicker lines on plots for printing, but although the 'lwd' parameter works fine with x11(), the thickened lines do not print with print.plot, or by using postscript() directly. 2. Try the following, plot(1:10, -(1:10))

Setup: Sparcstation 4 with Solaris 2.5.1 and Sun compilers After downloading R-0.16.1 and runnning configure, make fails as follows cd main; make cc -o R.binary arithmetic.o array.o attrib.o bind.o builtin.o character.o coerce.o colors.o context.o cov.o cum.o dataframe.o debug.o deparse.o dotcode.o dstruct.o duplicate.o edit.o envir.o errors.o eval.o format.o gram.o