similar to: Contour Plot on a non Rectangular Grid

Hi, I'm currently hunting a bug in the akima library, especially in the code behind the interpp.old function (bi-variate linear interpolation). It is based on a triangulation algorithm, interpolation at a given point needs to know the triangle which contains this point, then the interpolation is a straightforward calculation based on the three vertexes. The problem is: Sometimes the triangle

Bare test code: My simple Java test class source and R test code follow: public class RJavTest { public static void main(String[]args) { RJavTest rJavTest=new RJavTest(); } public final static String conStg="testString"; public final static double con0dbl=10000001; public final static double[]con1Arr=new double[] { 10001,10002,10003,10004,10005,10006 }; public final static

Dear R-Listers, I have spatial data on an equidistant rectangular grid, similar to topographic data. I know that there are quite a few R-packages or base functions that provide nice iso-contours plot, but I don't want a plot, just the smoothed isocontour line of ONE level (e.g. 10 mm). Data sets are large, so it would be preferable if the availability of regular grid data could be exploited,

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 I have a rectangular matrix of size 920 by 85. I'd like to coerce it into a square matrix such that all row/col names are present in the new matrix and the additional values are zero. As an example: A B C D A 1 2 3 4 E 5 6 7 8 F 9 10 11 12 Would be coerced to: A B C D E F A 1 2 3 4 5 9 B 2 0 0 0 6 10 C 3 0 0

Is there any function to show the points like this example of SPSS? http://biostatistic.de/temp/reg.jpg The point size should represent the number of data at this point. with regards Knut Krueger

I have a data set that I display using barplot. I don't know what you call it, but when I look at it, it looks like this: > lsu (0,0.1] (0.1,0.2] (0.2,0.3] (0.3,0.4] (0.4,0.5] (0.5,0.6] A 0.052631579 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 B 0.000000000 0.000000000 0.001007049 0.003021148 0.000000000 0.000000000 E 0.200000000 0.000000000

Hi, I'm trying to customize a window with 2 graphs. I'm able to do the first one with something like this general example par(mfrow=c(1,2),cex.axis=0.85,cex.lab=0.80,mai=c(1.3,1,0.5,0),las=3) bplot<-barplot(bar.values,names.arg=cf.names,width=0.5,ylab="% Area held") abline(h=0.3,lty=3,col="red") abline(h=0.1,lty=3,col="blue")

Hello! I have created a bplot-figure using this code: *file <- "2dcali_red.ttt" ux<-as.matrix(read.table(file, dec = ",")) mode(ux)<-'numeric' vel<-ux[,1] ang<-ux[,2] x<-ux[,3] y<-ux[,4] dat<- data.frame(ang=ang, x=x,y=y) require(rms) ddist2 <- datadist(dat) options(datadist="ddist2") fitn <- lrm(ang ~ rcs(x,4) +

I have a plot produced by function bplot (package = rms) that is really a lattice plot (class="trellis"). It is similar to this plot produced by a very minor modification of the first example on the bplot help page: requiere(rms) n <- 1000 # define sample size set.seed(17) # so can reproduce the results age <- rnorm(n, 50, 10) blood.pressure <- rnorm(n, 120,

Hello everybody, I'm a new user of R and I'm working right now with the ACE function from the acepack library. I Have a question: Is there a way to predict new responses using ACE? What I mean is doing something similar to the following code that uses PPR (Projection Pursuit Regression): library(MASS) x <- runif(20, 0, 1) xnew <- runif(2000, 0, 1) y <- sin(x) a <- ppr(x, y,

Hello, Is there any other bivariate pointwise interpolation command besides akima's interpp? I tried to search through the J. Baron's page without luck. The problem is that I have got regularly spaced data (in x and y) what is not acceptable for interpp. I am not very much interested in the method of interpolation as the data are dense and the error would not be to high. Thanks in

... finding myself doing this by hand one too many times ... maybe someone else has found a nice general solution ... I have a data frame containing n rows on p+q variables which I will call p.1, p2, ..., pp and q1, q2, ... qq. For simplicity all of these variables are continuous numeric. I want to create the following arrangement of plots: A scatterplot matrix which contains p times q panels

Hi. I'm trying to take a data set with two independent and one dependent variable and enter a x,y value to predict the dependent with a nonparametric technique. I've been using interpp in the akima package, (windows xp, R 2.4.1), but get values that are orders of magnitude off when the predictors are slightly out of the range of the data set. Can you recommend a function for me?

Hi Duncan, Thanks for your reply. Yes, it does seem to be specific to the CTYPE setting to Chinese on Windows. If I set it to English using Sys.setlocale() there is no problem, then back to Chinese and the authors disappear: Sys.setlocale("LC_ALL","English") citation("readr") #' To cite package ?readr? in publications use: #' #' Hadley Wickham, Jim

Dear all, I am running a code using bootstraps for estimating standard errors but the mac crashes. When I use small number of bootstraps (100) it works fine but if I increase that number it crashes. Thanks in advance dimitris The code, the error and my mac characteristics are the following: ##############code##################### qr.1<- rq(y~factor(year)+factor(state)+x1+I(x^2)+I(x^3),

Hello, One of the less friendly lists I am on is the Postfix list, but one of the features that Postfix has, and the use of which is required when asking for help on their list, is a command called 'postconf', which simply dumps the current configuration parameters into a neat, orderly output that you can then cut-n-paste into an email request for help. The postfix guys won't

Hi, I''m currently trying to set up a solution involving multiple servers using the same index over nfs. The problem is that from what I have seen, ferret doesn''t support multiple processes writing to the same index. Using a DRb service is not an option since this would create a single point of failure. I tried using Ferret::Store::FSDirectory to create a write lock on the the

Dear All, The following code save my graphs as pdf: pdf("j:/mix.pdf", width = 18, height = 16) grid.newpage() pushViewport(viewport(layout = grid.layout(3,1))) vplayout <- function(x, y) viewport(layout.pos.row = x, layout.pos.col = y) print(Aplot, vp = vplayout(1, 1)) print(Bplot, vp = vplayout(2, 1)) print(Cplot, vp = vplayout(3, 1)) dev.off() How can I save it in

Recently I was trying to cite a package where the authors have ? and ? in their names. I found that on Windows the citation() function did not return the authors' names at all, but on Linux there was no problem (sessionInfos at the bottom): On Windows, no author names are returned: #--------------- > citation("readr") To cite package ?readr? in publications use: (2017).

Hi Everyone, I''m using Ferret 0.10.11 with acts_as_ferret from SVN (same results with 0.10.10 and 0.10.9 though). I''m running into an odd problem where the scores of my top-ranking search results are ridiculously small - even when the query is one that should match at least one document with a decent score. To give an example, I have just the names of 5 businesses indexed