search for: circstat

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 negative mean direction like -0.8309982. I think that this might be wrong. If it is correct, what...

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: http://www.nabble.com...

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?

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

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 points. I'm sorry if this is trivial to most of you, but I really appreciate any help. I...

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

...ust can't figure out the correct place to change the definitions. Along these same lines, what does 'strip' mean when referring to strip text? Q2. How can I move the legend defining bar height into the plot, so that it overlays the lines they refer to? Consider the same figure using Circstats: library(CircStats) dat.rad <- (dat*((2*pi)/8)) -(2*pi)/16 rose.diag(dat.rad, bins = 8) #note the origin is to the right rather than on top Q3. The key difference is that CircStats uses an area-based calculation for the size of each slice, which makes for a different presentation than ggplot...

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 R1.5.1 (under Windows 2000). Furthermore, the results were generally incorect (please see the attached r...

Hi, If I have two sets of directional data (in radians) and want to compare them with a multifactorial anova. Is it even legitimate to compare circular data with an anova? The books I've picked up from the library don't really say, but it looks unlikely. If it is allowable, is my having stored the data as circular (X = as.circular(A)) something the aov() function will take into account,

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

[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

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 Africa Tel: +27-31-2042737(O) Fax: +27-31-2042736(O) Mobile: +27-723682954 Publications: http://w...

...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 circular package work for me? > R.version.string [1] "R version 2.6.1 (2007-11-26)" > library(circular) > y <- circular(c(45, 135), units = "degrees", template = "geographics") > range.circular(y) Circular Data: Type = ang...

...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 problems. To confirm, 'library(lme4)' produced the following error: Error: package 'Matrix' required by 'lme4' could not be found What do you suggest? Spencer Graves p.s. I get the same result using several different (US)...

Hi, I develop a little pice of code to perform a watson test. I try to compare the result obtain with my code and those given by the R library CircStats and it's watson(x,alpha,dist) provided function. My trouble is that I've got a set of test data (with it's watson test result) and the result given by R's watson function is not the same. Is there anybody who knows witch function is used in the R's watson function ? Thanks --...

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: +44-(0)1334-462637 Fax: +44-(0)1334-462632 Web: http://smub.st.and.ac.uk

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

Dear All, Ulric Lund and me are developing a new package called 'circular' that as soon as possible will substitute 'CircStats' package. A pre-released is available at: http://www.dst.unive.it/~claudio/R/circular_0.1.tar.gz I kindly ask everybody is using CircStats to check this new package and send me comments, suggestions and bugs to claudio at unive.it and not to the list. As soon as the package will work fine we w...

Dear all, I would like to fit a von mises mixture to my data using the EM Algorithm. Unfornately, existing packages (including circstats) doesn't provide such method. So, does anybody have written such code ? Big thanks, Christophe, (R newbie)

Dear All, Ulric Lund and me are developing a new package called 'circular' that as soon as possible will substitute 'CircStats' package. A pre-released is available at: http://www.dst.unive.it/~claudio/R/circular_0.1.tar.gz I kindly ask everybody is using CircStats to check this new package and send me comments, suggestions and bugs to claudio at unive.it and not to the list. As soon as the package will work fine we w...