similar to: Multiple comparison for circular data

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-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, 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

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

Hello, I have found a problem with using the glm function with a gamma family. I have a vector of data, assumed to be generated by a gamma distribution. The parameters of this gamma distribution are estimated in two ways (i) using the glm() function, (ii) "by hand", using the optim() function. I find that the -2*likelihood at the maximum found by (i) is substantially larger than that

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, 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

Hi R users, Sorry for a simple question: I found different people use different i in set.seed(i), are there any rules to choose an i or one can choose as he likes? Thanks Xiaohua

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

Hi all, >From late April, everytime when I update my packages, the package "gap" is always on the list of "packages to be updated". I try to download and install the gap zip file from CRAN. However, the same problem still remains. No such problem for other packages. > Sys.info() sysname release

Assume a C program invoked by .Call, that returns a list. Near the top of the program we allocate space for all the list elements. (It is my habit to use "xyz2" for the name of the R object and "xyz" for the pointer to its contents.) PROTECT(means2 = allocVector(REALSXP, nvar)); means = REAL(means2); PROTECT(u2 = allocVector(REALSXP, nvar)); u =

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

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

Dear R users, I have been trying to solve for mle's of a function that involves an integral and I keep getting an error. I created an example to work on first and even the simple example doesn't give me the mle's. I am getting the error "Error in integrate(integrand, 0, Inf) : non-finite function value". I divided my likelihood function into two parts, one part involves

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

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.

The Partek package (www.partek.com) allows only two selections for Multiple Test Correction: Bonferroni and Dunn-Sidak. Can anyone suggest why Partek implemented Dunn-Sidak and not the other methods that R has? Is there any particular advantage to the Dunn-Sidak method? R knows about these methods (in R 2.1.1): > p.adjust.methods [1] "holm" "hochberg" "hommel"

Matt, Here's a Bonferroni-corrected multiple one-sample t-test that I wrote some years ago. It took a while to get it into R, as na.omit doesn't seem to handle vectors and I had to write a quick kludge (na.remove). Another more general point was that I discovered that the help page for t.test gives the name "parameters" for the degrees of freedom, as in S. However, the name

Hi! I need to run a wilcoxon (Mann-whitly, in fact) test with bonferroni correction, as I am running 10 consecutive wilcoxon test not independent, and I know that bonferroni will partially correct for this problem, but I have no idea how to do it with R, I have been looking in the archive but couldn't understand how to do it. The format I am using at the moment is r4_o <- [1]

Hi- Is anyone working on multiple comparisons of means or post-hoc tests (ie: Tukey, Bonferroni) for R? I saw in the winter 98' archives of the R mailing lists that these tests had not been implemented yet, I was just wondering if I could look foward to having them. ;-) I also looked through the contributed packages and didn't see anything that offered such tests. I guess I could check