similar to: A bug in circ.plot of CircStats package

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?

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

Hi, I try to compute circular means for a matrix with NAs, but the function "circ.mean" return only means for lines with complete values and do not accept "na.omit=T" or "na.rm=T", or "na.action=na.omit", or "na.fail=T". Also, I try to use "polar.plot" of the package poltrix with the same matrix, but angles are not display because of

Hello, I am quite new to R and although I have read some material (apparently not enough) I'm getting very frustrated. I'm trying to do some simple circular statistics in R (CircStats, primarily circ.mean) and I'm confused about the results. What is a negative radian or degree mean in R? I thought it would be a coterminal angle, but that doesn't make sense with my data

Hello, Once again, I posted a message without a subject line. Sorry.... here is the question again. Is there a simple way to modify the circ.mean function in the CircStats package to include a vector of weights to obtain a weighted average angle? Thanks! Martin -- Martin Biuw Sea Mammal Research Unit Gatty Marine Laboratory, University of St Andrews St Andrews, Fife KY16 8PA Scotland Ph:

I have a vector of dimension 67, but CircStat circ.summary and circ.disp seems to conclude that the dimension is 1. I am new to R so maybe I'm just missing something. > dim(data) [1] 67 1 > circ.mean(data) [1] 0.2200071 > circ.summary(data) n mean.dir rho 1 1 0.2200071 43.12534 > circ.disp(data) n r rbar var 1 1 43.12534 43.12534 -42.12534 >

I have a problem that deals with correlating wind velocity to seed collection data. The problem lies in that I have a wind data set that is contains 2000+ data points and weed collection data on the order of a couple hundred. Both data sets were collected for the same time period, but there is not a one-to-one wind velocity->seed location match. My understanding of correlation is that you

Hola a todos, es la primera vez que escribo a la lista, espero explicarme correctamente. Trabajo en cronobiología de mamíferos y quiero comenzar a analizar mis datos con R (haciendo análisis de varianza ya llevo unos años en el entorno). Me interesa saber si cada una de mis variables podría ajustarse a una función sinusoidal; tengo medidas cada dos horas para cada una. Querría extraer el mesor y

I want to get the minimum arc (in degrees) needed to include a set of compass directions. I would like to use the range.circular() function of the package circular, because that package understands a compass-type of angle convention, but it gives results I don't understand. Howver, I can get the correct answer in the example below, 90 degrees, using the CircStats package. How can I make the

Dear R community, I tried to post a query about the use of the rose.diag function from the CircStats package to the r-help forum a few days ago. It gave all appearances of having been successful. It would have been my first post; I am correctly recieving all the R-help emails, but  but the message has not appeared anywhere on the forum. Am I held up by the moderator somewhere? any help would be

Hello I will make an interpolation of data which represents azimuth direction ( angle from north in clockwise direction) values. But there is a problem. Say, for instance, while 1 and 359 indicate somewhat same direction, interpolation puts values in the range from 1 to 359. What can I do to solve the problem ? Anything you offer ? thanks in advance Ahmet Temiz General Directory of

Hello all, I am trying to use xyplot to give a six panel plot and to highlight only points (in any panel) that meet a certain criterion. With the plot command I would do something like: plot.default(filein$Site,filein$circ.conc) points(filein$Site,filein$circ.conc,type="p", pch=ifelse(filein$p.value<5e-02,19,21)) I had thought I could just stick in the pch line from above into

[message needed manual approvement: it had a single line "help" -- MM ] I have a table containing two columns and some 600 rows. The first column is a time series in the sense contains the amount of time the model that i am studying stays in a particular state(ie sojourn time). I extract the contents of the table using read.table. and extract the first column from the table

Hello everyone, I'm stumped. I'd like to create a scatterplot matrix with circular reference lines. Here is an example in 2d: library(ellipse) set.seed(1) dat <- matrix(rnorm(300), ncol = 3) colnames(dat) <- c("X1", "X2", "X3") dat <- as.data.frame(dat) grps <- factor(rep(letters[1:4], 25)) panel.circ <- function(x, y, ...) { circ1

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

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

Dear Sir/Madam, I am running a R program for my SNP data. There are some errors when I run glm model in Hapassoc software, sometimes it is over the memory and sometimes the matrix is singular. I want to ingore these errors and excute the next statement. Otherwise, I have a big big trouble. Do you have some idea about this problem of ingore errors. Wish to get your help assp. thanks. -- Wei

Hi guys Here is my problem... Using Ajax to update a div with link_to_remote. Problem is I have a layout defined in my class, so, the div gets the layout when I call the method, how di I tell a method to ingore the layout? Cheers. Pieter. -------------- next part -------------- An HTML attachment was scrubbed... URL:

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