search for: circ

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

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

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, ...) {...

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 either xyplot or xYplot with groups set to Site like: xyplot(circ.conc~Day|Site,data=allsites.byday,groups=Site, pch=ifelse...

Dear all, After going through the Lattice doc and R-help list and google, I got the feeling that there is no function in lattice or other package to compute a pie chart object of class "trellis". Although pie charts are obviously not considered optimal even in the pie() doc ;-) , pie chart trellis objects would be easy positioned e.g. over a map drawn with the grids package. Can

OK, it seems to be an authentication issue. At least w.r.t. the member server. When I add the user name, it does ask me for the password and succeeds. I don't see any difference when using IP numbers instead of host names. root at circe:~# samba-tool time ulysses -U peter Password for [HOME\peter]: Fri Feb 27 14:31:46 2015 CET

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

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

Hai, ? Can somebody try to run : ? samba-tool time and samba-tool time servername.fqdn ? when i run it, im getting : ? ERROR(runtime): uncaught exception - Connection to SRVSVC pipe of server 'rotterdam.bazuin.nl' failed: NT_STATUS_NONE_MAPPED ? File "/usr/lib/python2.7/dist-packages/samba/netcmd/__init__.py", line 175, in _run ??? return self.run(*args, **kwargs) ? File

I've just gotten through translating my notes for the course in basic statistics for health researchers. I've placed a copy on ftp://blueberry.kubism.ku.dk/priv/R-intro.tgz This is essentially a translation of the Danish notes that I wrote during the Autumn of 1997 for use with Rsept, i.e. roughly R v.50. They lean heavily on the textbook, Altman: Practical Statistics for Medical

I've got a netbook, circe 2009. When I got it, not that I was wild about ubuntu, but there was specifically an ubuntu netbook remix. Well, it's a few years later - has anyone put CentOS (6, preferably) on a netbook, and were there any problems? mark

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

...chrieb Peter Serbe: >> OK, it seems to be an authentication issue. At least w.r.t. the member server. >> When I add the user name, it does ask me for the password and succeeds. >> I don't see any difference when using IP numbers instead of host names. >> >> root at circe:~# samba-tool time ulysses -U peter >> Password for [HOME\peter]: >> Fri Feb 27 14:31:46 2015 CET >> > -- > To unsubscribe from this list go to the following URL and read the > instructions: https://lists.samba.org/mailman/options/samba >

Hello everybody. I am trying to cluster circular data (data points which are angles), thus i can not use the "dist" function in "mclust" to generate my distance matrix, I am using the function " Dij = 0.5*( 1 - cos(theta_i - theta_j)). The thing is "hclust" will not accept this distance matrix, i tried to pu...

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

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

...ta 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 need equal data sets. Is there a way to compare means and standard deviations (or kappa in circular stats) that doesn't require equal data sets? Thank you in advance. Joe Example (simplified): library(CircStats) wind<-c(2.1, 2.3, 2.6, 2.1, 2.6, 2.4, 1.2, 1.5) # in radians seed<-c(1.7, 1.3, 1.5, 1.9, 1.8, 2.3) circ.cor(wind, seed, test=T) Joseph Dauer Master's...

...to create a heatmap using 2 matrices I have: z and v. Both matrices represent different correlations for the same independent variables. The problem I have is that I wish to have the values from matrix z to be represented by color intensity while having the values from matrix v to be represented by circle size. I currently have the following in front of me and an unsure of what to add or change in order to achieve that goal. panel.corrgram.2 = function(x, y, z, subscripts, at = pretty(z), scale = 0.8, ...) { require("grid", quietly = TRUE) x <- as.numeric(x)[subscripts] y <-...

...is identify the 3 break/inflection points in the population trend. That way, I can make an argument that the break points corresponded to shifts in government policy with respect to land use management. I've been looking at the segmented package, and initially I looked at change.pt test in the circ.stats package (which is inappropriate b/c my data is not amenable to circular statistical analysis). Any ideas on what I could do would be appreciated! Thank you! -Charlotte