similar to: help for linear-circular correlation

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, I use plot.circular(rad, stack=T,bins=4) and could just obtain the four stacks in the directions of 45, 135, 225 and 315 degrees. But I want them in 0, 90, 180 and 270 degrees. Is there any parameter or any other function to make it? Thanks for helping a R beginner. -- Xiaohua Dai, Dr. -------------------------------------------------------------------------------- Centre for

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 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 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 All, Does anyone know how to compare the means of several circular-data samplings? Any related websites, references and softwares? Thanks Xiaohua -- Xiaohua Dai, Dr. [[alternative HTML version deleted]]

Se ha borrado un texto insertado con un juego de caracteres sin especificar... Nombre: no disponible Url: https://stat.ethz.ch/pipermail/r-help/attachments/20060628/32226ee4/attachment.pl

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

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

Dear R-Group: Recently, I ran into a problem. I was using a function called "I.1", which evaluates the first-order modified Bessel function of the first kind, in the package "CircStats". This function is not vectorized, since it uses a couple of "if" conditions. However, when I called this function with a vector argument, I got no error/warning messages in

Hi, I've currently to work with some circular data. Unfortunately I'm not very familiar with circular statistics and would really appreciate if I could get some help concerning the CircStats package this way. My data lies in the range 0 to 2*pi, and is transformed to radians (as expected by the CircStats methods). Calculating the mean direction (circ.mean) results for some datasets in a

Dear R gurus - consider this plot: library(ggplot2) dat <- sample(1:8,100,replace=TRUE) smp <- ggplot(data.frame(dat), aes(x=factor(dat),fill=factor(dat))) + geom_bar(width=1) smp + coord_polar() Q1. How do I change the font size and weight of bar labels (1,2,3...)? I've been wallowing in the 'Themes' structure and I just can't figure out the correct place to change

Dear all I was writing some code that needed functions from packages 'MASS' and 'CircStats', and I received a warning saying that the function 'eqscplot' in one of the packages was masked by another 'eqscplot' from the other package (in fact 'eqscplot' from 'CircStats' seems a short version of 'eqscplot' from 'MASS'). This masking of

Hi, i found a bug in circ.plot of CircStats package, namely that the plot ingore the datasets with value = 0. it's simple to fix it, but i cannot seem to figure out where to submit the fix. is there a official procedure of submitting bug report/patch to the contributed packages? -- View this message in context:

circ.plot in CircStats under R 1.8.0 yields incorrect plots (easily compare hist() and rose.diag()). However, the code in circ.plot has not changed. Has anyone tracked this down yet?

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.

I have data for several rings of a left heart chamber, and which I would like to display in concentric rings, with color-encoding of the values. Each ring corresponds to one slice through the heart, and the rings correspond to positions from the base to the apex of the heart as you move from the outermost ring to the innermost one. The data have a circular pattern. These types of displays are

Yesterday, I installed R2.2.0 for Windows [Version 2.2.0 (2005-10-06 r35749)]. Unfortunately, 'install.packages("Matrix")' produced the following message: Warning in download.packages(pkgs, destdir = tmpd, available = available, : no package 'Matrix' at the repositories I installed lme4, maps, mapproj, CircStats, scatterplot3d, gregmisc, Hmisc without

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 R users, Hope the question is not too simple: How to use a for loop to generate two lists as below: testlist <- list(test1, test2, ..., test1000) stringlist <- list("test1","test2",...,"test1000") Thanks Xiaohua