similar to: H-scatter plot in geostatistics

Hi, I am trying to find m breakpoints of a linear regression model. I used the segmented package. It works fine for small number of predicators and breakpoints.(3 r.v. 3 points). However, my model has 14 variables it even would not work even for just one breakpoints!. The error message is always estimated breakpoints are out of range. Since my problem is time related problem. So I

The majority of the functionality present in the gstat stand-alone program (http://www.gstat.org/) is now available as a package/library for the S language (R, S-Plus), again called gstat. The package provides multivariable geostatistical modelling, prediction and simulation, as well as several visualisation functions. Gstat was started 10 years ago and was released under the GPL in 1996; the

The majority of the functionality present in the gstat stand-alone program (http://www.gstat.org/) is now available as a package/library for the S language (R, S-Plus), again called gstat. The package provides multivariable geostatistical modelling, prediction and simulation, as well as several visualisation functions. Gstat was started 10 years ago and was released under the GPL in 1996; the

Hello everyone! I am using GRASS with spgrass6 for my work. I will be using variograms in the process of landsat image classification. I am quite ok with GRASS but am finding R really tough. I understand that spgrass6 is a link between GRASS and R which can read and write raster/vector layers. Out of really many packages in R, for generating variograms out of landsat images which packages of R

Dear All: Pleas help me to increase the memory in R. I am trying to make euclidean distance matrix. The number of low in data is 500,000. Therefore, the dimension of euclidean distance matrix is 500,000*500,000. When I run the data in R. R could not make distance matrix because of memory allocation problem. In order increase memory, I read the FAQ and follow the instruction as below:

Hi, I read somewhere that when carrying out geostatistical analysis in R you should not use latitude and longitude...can anyone expand on this a little for me, and what would be the best coordinate system to use? I have my data in a geographic coordinate system, WGS84, decimal degrees....is this the wrong format for such analyses? I have also converted my data in the UTM projection and so have

> If you use the spatial objects provided by the > sp-package (http://cran.r-project.org/web/packages/sp/vignettes/sp.pdf) > you transform your data to other projections using the spTransform package. Thus you will need the rgdal package in complement (it actually includes spTransform). This function is extremely convenient: you can manage coordinate transformations extremely easily

I tried toy problems and there doesn't seem to be a basic problem between lm() and ts objects: X = data.frame(x=c(1,2,7,9), y=c(7,2,3,1)) lm(y ~ x, X) X <- lapply(X, function(x) ts(x, frequency=12, start=c(1994,7))) lm(y ~ x, X) and this works fine - whether you do an lm() before or after making ts objects, it's okay. But I have a situation where things aren't okay.

Dear list, In using the scatter.smooth() function (modreg) on a small data set (100 obs) the following error was produced: > scatter.smooth(Na, S) Warning message: k-d tree limited by memory. ncmax= 200 I haven't used scatter.smooth much but when I have, I haven't seen this message before. gc() returns > gc() used (Mb) gc trigger (Mb) Ncells 417693 11.2 667722

Dear List I hope you can help me: I?ve got a dataframe (df) within which I am looking for Peak Over Threshold values as well as the length of the events. An event starts when walevel equals 5.8 and it should end when walevel equals the lower threshold value (5.35). I tried ?clusters (?)? from ?evd package?, and varied r (see example) but it did not work for all events (again

On Wed, May 29, 2013 at 03:25:30PM +0200, Duncan Sands wrote: > Hi Jack, I pulled the loop vectorizer and fast math changes into the 3.3 branch, > so hopefully they will be part of 3.3 rc3 (and 3.3 final!). It would be great > if you could redo the benchmarks rc3. > Duncan, As requested, appended are the updated Polyhedron 2005 benchmark results with both RC1 and RC3 llvm 3.3

Hi,   I would like to draw a box at each corner of lower 10% and upper 10% in the scatter plot on(0,1)*(0,1)  to indicate the lower and upper tail.  I hope somebody can help me.   Here is my code: ## scatter plot plot(hume_uni[,2],beec_uni[,2], pch=19, xlab="Hume_uni", ylab="Beec_uni", col= "blue")  title(main="Scatter plot of uniformised data of Hume and

Hello, I have been using scatter.smooth and been trying to format my plots and been unable to adjust the line type or color of the smooth curve. When I specify col in scatter.smooth, I can adjust the color of the scatter plot points, but the line always stays black and solid. lty has no effect. When I specify par prior to calling scatter.smooth, col will change the color of the points,

Does anybody know how to adding colors to data points in scatter plot in R? I have the following R data frame - Voter_Name, Education_years, Income, Registered Ricky, 18, 40000, Y Lisa, 20, 34000, N Jon, 10, 35000, Y ... I am doing a scatter plot in R for those data points - Income vs. EducationYears. I would also like to using different colors to differentiate whether that data point is a

Dear List, Having had my appetite sufficiently whetted by Prof. Ripley's email about the new graphics capabilities in Unixes, I wanted to try them out. I updated to svn r44608, configured with the following options: R is now configured for x86_64-unknown-linux-gnu Source directory: .. Installation directory: /usr/local C compiler: gcc -O3 -g -std=gnu99

Hi R Users, I have a dataframe like this: id x rho A 1 0.1 B 20 0.5 C 2 0.9 ... I want to do a scatter plot of "x" versus "rho" but for each point on the scatter plot I want the corresponding entry for "id" instead of points. In STATA I can do so by twoway (scatter x rho, mlabel(id)) How can I do the same in R? I am sure there is some simple way to do

Dear All, I run the goodness of fit test using goodfit() in vcd package. The result is as follow: Goodness-of-fit test for poisson distribution X^2 df P(> X^2) Pearson 1.053348 2 0.5905661 Warning message: In summary.goodfit(gf) : Chi-squared approximation may be incorrect I want to save the the test statistics(X^2), df, and p-value. How can I save the result.

On 5/3/19 10:19 AM, Harsh Gondaliya wrote: > In my Linux VM, I am unable to change scatter-gather / tx-checksumming / > TCP-segmentation offload features by ethtool command. I am using virtio nic > with vhostuser port backend. > It shows me an error: > Cannot change tx-checksumming > Cannot change scatter-gather > Cannot change tcp-segmentation offload > Could not change

Dear useRs, I'd like to produce some scatter plots where N units on the X axis are equal to N units on the Y axis (as measured with a ruler, on screen or paper). This approach x <- sample(10:200,40) ; y <- sample(20:100,40) windows(width=max(x),height=max(y)) plot(x,y) is better than plot(x,y) but doesn't solve the problem because of the other parameters (margins etc). Is

Hi I'm trying to scatter plot 3 data sets on the one graph. I plot the first data set by calling: plot and the next two data sets by calling points points and coloring them different colors to distinguish between the data sets I've also set the xlim and ylim values appropiately The problem is points for one of the data sets will not appear when points is called and only appears when