similar to: function arrows.circular not working

I'm wondering if I am dealing with a limitation in lattice. It's probably due to my own limitations though. I'm working with a lattice dotplot. The x-axis is set at "free". In one panel, there are only two data points that have equal x-axis values. Since no range exists, lattice isn't plotting those points. It reports: "Error Packet using 6 'x' and

I'm trying to use fine axis controls in lattice for each panel. Specifically, within each panel, I want to set the limits for x and y equal to each other since it is paired data (using the max value of the two). Of course, I have no problems setting the limits for the entire plot but I am having trouble setting them for each specific panel. Could someone please provide me some guidance?

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,

Hi all, I'm curious if anyone has coded an Expectation-Maximization algorithm that could help me model some circular data I have. I'd like to model it as a mixture of uniform and Von Mises centered on 0, so the only free parameters is the mixing proportion and the kappa of the Von Mises. I couldn't find anything in the contributed packages that seemed to suit this purpose. Any

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

Hi, I am trying to figure out the observed acceptance rate and M, using generalised rejection sampling to generate a sample from the posterior distribution for p. I have been told my code doesn't work because I need to "take the log of the expression for M, evaluate it and then exponentiate the result." This is because R is unable to calculate high powers such as 545.501. As

Dear all, I am using plot.circular(x, stack=TRUE) to plot a "histogram" from a list of angle. I would also like to draw a line from the origin at the angle of the mean (mean.circular), preferably with the resultant's length (rho.circular) as length. How do I achieve this on the circular plot, please? Thanks a lot, Tim -- -- Tim Gruene Institut fuer anorganische Chemie

Hello, Consider MCMC sampling with metropolis / metropolis hastings proposals and a density function with a given valid parameter space. How are MCMC proposals performed if the parameter could be located at the very extreme of the parameter space, or even 'beyond that' ? Example to express it and my very nontechnical 'beyond that': The von Mises distribution is a circular

Hi, I've noticed I get different results fitting a function to some data on my laptop to when I do it on my computer at work. Here's a code snippet of what I do: ##------------------------------------------------------------------ require(circular) ## for Bessel function I.0 ## Data: dd <- c(0.9975948929787, 0.9093316197395, 0.7838819026947, 0.9096108675003, 0.8901804089546,

Hi, I have a density that I need to get MLEs from, which includes definite integrals both in the denominator and in the numerator of the density function. It looks like the outcome depends on the initial values given. My program is shown below: library(circular) ######################################## 4 parameters ######################################## z<-rvonmises(100,0,1)

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

Hello Friends!!! I would want draw a circular histogram, and I would like draw a polygon with a circular side. This is easy if I use the functions polygon and arc, but I want that the polygon with a circular side have background colour. The polygon created with function polygon can have background colour, but the surface created with function arc can?t have background colour. How I could create a

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

Hi, I installed "circular" package and I wanted to load it automatically when R starts up, so I added the following lines in file Rprofile.site, .First <- function() library(circular) When R starts up, it gives me the following error: ---------- Loading required package: boot Loading required package: graphics Loading required package: stats Error : .onAttach failed in

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

require(circular) c <- circular(rep(0, 20), zero = pi/2, rotation = 'clock') plot(c, stack = TRUE, shrink = 1.5) Can anyone tell me why the stack is offset from 0? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick

Hi, I'm new to R. I'm trying to plot my data into a circle. my data sort of looks like 12,12,4,5,6,5,11,10,3,9,9,9,12,12,2 total of 15 numbers. I'm trying to add all the same numbers, such that, there are 4 of 12s,1 of 11, 1 of 10, 3 of 9s, and such.... so the circle plot would have 4 parts of 12, 1 part of 11, 1 part of 10, 3 part of 9, and such... I tried

In conducting studies of animal orientation and displacement, I need to produce circular histograms of angles (bearings in radians 0-2pi) where the centre of the circle indicates very few observations for a given bin of angles and outwardly concentric circles indicate greater frequencies of observations for a given bin of angles. I'd like not to have to write the function myself but I

I have a little image gallery application in rails enough for small galleries (20< photos). But as images increase and the user needs to scroll down the page it starts getting anoying to see all pictures... I thought pagination but then I saw this AJAX goodness http://www.apple.com/getamac/ads/ As you can see the videos scroll right/left when the arrows are pressed in a circular manner

Dear r-helpers, I want to show that time is flowing CCW in the following: require(circular) len <- 8 labl <- as.character(c(0, 1, 1, 1, 0, 0, 1, 0)) r <- circular(2*pi* (rep(c(1, 3, 6), each = 200)/len + rnorm(600, 0, 0.025))) r.dens <- density(r, bw = 25, adjust = 4, kernel = 'vonmises') plot(r, shrink = 2.5, axes = FALSE, ticks = FALSE, pch = 1, col =