similar to: How to plot circular data in the directions of 0, 0.5pi, pi and 1.5pi

This is slightly off-topic but I would be interested in whether anyone has succeeded in creating a filter expression for Google's gmail system that will select messages sent through the R-help and R-devel lists. It seems as if it should be easy to select on '[R]' or '[Rd]' in the subject line but I haven't been able to work out the exact syntax that would do this and not

Sorry, I did not make my question clear. Since I have a point theme with many points, some of them may clump together. the problems here are: 1. how to find clumps in a point theme? 2. the convex-hull extension I found only deal with all the points in a theme at each time? how to make each convex hull around each point clump automatically? Thanks. Xiaohua On 5/12/06, Bob Booth <bbooth

Hi Gurus! I have a job that is to get randomly samples from point-based GIS data (sp called shape GIS files) under the total sum resctricted. For example, I would like to take random smaples under the 1000 persons in each city. The randomly sampled persons should not be over 1000 any case. Thank you in advance, -- Kum-Hoe Hwang, Ph.D. Phone : 82-31-250-3283 Email : phdhwang@gmail.com

Hi R-profs, Maybe my question is a little off topic. Could any one tell me how to calculate a linear-circular correlation coefficient and its p-values? I had a quick look at circular and CircStats packages and did not find the related function. Thanks for any kindly help. Xiaohua -- Xiaohua Dai, Dr. Centre for Systems Research, Durban Institute of Technology P.O.Box 953, Durban 4000, South

Hi R-users: Yesterday I ran a R code for 9 hours and it did not show any sign to stop. Then I interrupted it and found it had completed 82.5%. This morning I decided to wait for another 11 hours to see what will happen. Wait a minute, I heard that transforming data.frame to matrix will make R code faster. Then I made the modification in my R code. Oooh, the new code finished within 30 minutes!!

Hi R users, Maybe the question is too simple. In a IF ... ELSE ... statement "if(cond) cons.expr else alt.expr", IF and ELSE should be at the same line? For example, if (x1==12) { y1 <- 5 }else { y1 <- 3 } is right, while if (x1==12) { y1 <- 5 } else # Error: syntax error { y1 <- 3 } is wrong? Thanks -- Xiaohua Dai

Dear R helpers, I have a largish matrix (1300 x 1300) and I wish to find the row and column numbers that identify particular elements whose values I know in advance (for example, the row and column numbers for the maximum value of the matrix). I have looked in the help manual and found the functions 'row' and 'col' for finding the indices of a matrix. But, I can't figure out

Hello, is there a package in R that would allow map drawing: coastlines, country/state boundaries, maybe topography, rivers etc? Thanks for any guidance, Mark

I would like to do the Durban-Watson test on a time series of log returns. 2 questions: 1) If I am just trying to find out if there is serial correlation, what do I do for the residuals? there is no model, so do I just use the log returns (time series) itself? 2) what is the code in R to accomplish this? Regards [[alternative HTML version deleted]]

hi, i am trying to teach myself about ARIMA models. i have followed examples from a number of sources and have more or less got the hang of how it works. i would like to compare the output from the fitted model to the original data. is this possible? or even a meaningful thing to do? to be clear, for example, having generated a fit to some data using > fit <- arima(LakeHuron, order = c(1,

Dear R users, I'm trying to specify a generalized linear mixed model in R, basically a Poisson model to describe monthly series of counts in different regions. My aim is to fit subject-specific curves, modelling a non-linear trend for each region through random effects for linear splines components (see Durban et al, Stat Med 2005, or " Semiparametric regression" by Ruppert et al,

Dear all, Does anyone know where there is R or S code for the CUSUM SQUARED structural breaks approach? (Brown, Durban and Evans, 1975 - used in Pesaran and Timmerman, 2002) The problem is that the breaks package only appears to offer the standard 'unsquared' CUSUM, even though it appears most think it is inferior to the squared version. It might appear to be a relatively simple

hi, i am sure that this is a trivial question but i have not been able to find an answer by searching the mailing lists. i want to plot points on a graph, joined by lines. the command that i am using is points(x, y, type = "b", pch = 21) this plots nice open circles at the data points and draws lines between them. however, the lines do not come all the way up to the edge of the

hello, i have been using ccf() to look at the correlation between lightning and electrogamnetic data. for the most part it has worked exactly as expected. however, i have come across something that puzzles me a bit: > x <- c(1, 0, 1, 0, 1, 0) > y <- c(0, 0, 0, 0, 0, 0) > ccf(x, x, plot = FALSE) Autocorrelations of series 'X', by lag -4 -3 -2 -1 0

Is there any routine to generate for one vector matrix. If I have X I want to generate start from zero to maximum value each vector. For example, I have a vector x = (4,2,3,1,4) I want to generate n=6 times, for 4, start 0 to 4, then 2 start 0 to 2, ect. The result something like this: generate(x,n=6) 1,1,2,1,4 1,2,3,0,3 4,0,1,1,1 3,1,0,1,4 0,0,3,0,0 4,1,3,0,4 Could anyone help me. Thanks.

Dear Dr. Bivand and R-users, I have a 5 by 5 grid, say, location[1:5,1:5], and I want to know the indices of 8 neighbours of each cell. For example, for location[2,2], its neighbour coordinates are [1,1:3], [2,1], [2,3] and [3,1:3]. Sometimes I also need to remove edge effects (torus = TRUE). I have tried "cell2nb" function in your spdep package. Here's my example: > neigh

hi, i am trying to figure out how to put a caption on the colour scale of a levelplot. there does not seem to be an option for this in levelplot(). i tried using mtext() but as soon as you put the text far out enough on the right of the plot, it goes beyond the plot boundary. so i tried to extend the margin on the right of the plot using par(mar) but this did not have any effect on the plot area.

Zulu time will (necessarily) be reckoned relative to Durban, South Africa (capital of KwaZulu-Natal province, South Africa) where the Zulu people live. Brian Brunner brian.t.brunner at gai-tronics.com (610)796-5838 >>> sam at wa4phy.net 11/07/05 07:37PM >>> Peter Farrow wrote: > For those of us who live the promised land where GMT is absolute time > all the this

Dear all I was looking at the aov documentation page and came across the following which seems like a contradiction to me: " This provides a wrapper to |lm| for fitting linear models to balanced or unbalanced experimental designs." (I presume 'This' refers to aov) and "|aov| is designed for balanced designs, and the results can be hard to interpret without

Howdy Is there any reason why I cannot mount to a Dfs share on a Windows 2000 box from Linux? I have tried: # mount -t smbfs //server/share$ /mnt/point -o username=user Which works... # mount -t smbfs //server/share$/home/pupils/user /mnt/point -o username=user Trying to mount any share within the Dfs root retruns: 23479: tree connect failed: ERRDOS - ERRnosuchshare (You specified an