similar to: How to add direction of time to plot.circular()

Thanks Dimitris!!! That's much clearer now. Still have a lot of work to do this weekend to understand every bit but your code will prove very useful. Cheers, Aziz -----Original Message----- From: Dimitrios Rizopoulos [mailto:Dimitris.Rizopoulos at med.kuleuven.be] Sent: May 12, 2006 4:35 PM To: Chaouch, Aziz Subject: RE: [R] Maximum likelihood estimate of bivariate

Dear Friends, Here is an algorithm for finding concave hulls: http://get.dsi.uminho.pt/local/ Has anyone implemented such an algorithm in R? RSiteSearch('concave hull') didn't reveal one (I think). _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology Postal Address: P.O.Box 400400, Charlottesville, VA 22904-4400 Express Parcels

Dear Friends, require(cluster) x <- rbind(cbind(rnorm(10, 0, 0.5), rnorm(10, 0, 0.5)), cbind(rnorm(15, 5, 0.5), rnorm(15, 5, 0.5))) plot(pp <- pam(x, 2), which.plots = 1) How can I extract the coordinates used in the plot? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology Postal Address: P.O.Box 400400, Charlottesville, VA 22904-4400

Kind Experts, Many thanks for your guide. I have tried to figure out something that can help me plot my own data using the examples you referred me to. I copied part of the code as: set.seed(44) N=500 events <- as.POSIXct("2011-01-01", tz="GMT") + days(floor(365*runif(N))) + hours(floor(24*rnorm(N))) + # using rnorm here

Dear R-helpers, I wonder if you can give me advice about the best way to use help(). (1) If I type ?normal because I forgot the name dnorm() I get a long list of relevant pages. Getting to right page is laborious. (2) If I remember dnorm() and want to be reminded of the call, I also get a list of pages. Advice? It seems to me that if the output of help() listed base functions first, it

Dear R-helpers, I run R version 2.8.1 (2008-12-22) on i386-apple-darwin8.11.1 I would like to call R from Java. I downloaded SJava_0.69-0.tar.gz, and then issued the command suggested on http://www.omegahat.org/ RSJava/ : % R INSTALL -c SJava_0.69-0.tar.gz ARGUMENT 'INSTALL' __ignored__ WARNING: unknown option '-c' ARGUMENT 'SJava_0.69-0.tar.gz' __ignored__ I then

Looks like you forgot to load the lubridate package library(lubridate) You are calling functions days(), hours(), minutes(), seconds(), and hour() which all come from that package. ------------------------------------- David L Carlson Department of Anthropology Texas A&M University College Station, TX 77840-4352 -----Original Message----- From: R-help [mailto:r-help-bounces at

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 dealing with wind data and I'd like to model their distribution in order to simulate data to fill-in missing values. Wind direction are typically following a vonmises distribution and wind speeds follow a weibull distribution. I'd like to build a joint distribution of directions and speeds as a VonMises-Weibull bivariate distribution. First is this a stupid question? I'm

Dear All, I am trying to generate a circular/radial plot. The script below has a result I am looking for: testlen<-rnorm(24)*2+5 testpos<-0:23+rnorm(24)/4 clock24.plot(testlen,testpos,main="Test Clock24 (lines)",show.grid=FALSE, line.col="green",lwd=3) if(dev.interactive()) par(ask=TRUE) # now do a 'daylight' plot

I use ggplot2 for all my plotting needs where you can make plots circular with the coord_polar. Maybe this will help you along: http://rstudio-pubs-static.s3.amazonaws.com/3369_998f8b2d788e4a0384ae565c4280aa47.html On Fri, 22 Apr 2016 at 08:31 Ogbos Okike <giftedlife2014 at gmail.com> wrote: > Dear All, > I am trying to generate a circular/radial plot. The script below has a >

Hey R Users I like to control the ticks and labels in a boxplot as described for a xyplot below (thread in maillinglist in may 2003). Does anybody knows how it works? Thanks in advance Thomas Thread from May 2003 (http://tolstoy.newcastle.edu.au/R/help/03a/5604.html) Hello R Users! I'm using lattice to produce some graphs with logaritmic y-scales. I use the command xyplot(hits ~

Hello R Users! I'm using lattice to produce some graphs with logaritmic y-scales. I use the command xyplot(hits ~ c(1:1024), data=eichData, type="S", scales=list(y = list(log=10))) to create the plot. This is fine, except for the automatically choosen tick marks. I'd like to have a major tick at the 10^n location and minor ticks in between which correspond with the native

Dear R-helpers, I have two data frames, op and em4: > str(op) 'data.frame': 37 obs. of 5 variables: $ m : num 0.202 0.336 0.122 0.139 0.14 ... $ lln : num 0.798 0.643 0.863 0.835 0.823 ... $ rrn : num 0.789 0.702 0.894 0.895 0.923 ... $ asym2: num 0.177 0.304 0.108 0.187 0.274 ... $ asym3: num 0.0755 0.0975 0.0818 0.0651 0.13 ... > str(rownames(op)) chr

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

Dear group, I tried to generate labels for every second tick in lattice (qqmath). Version: 0.14-16 Date: 2006/12/01 R version 2.4.1 (2006-12-18) An example: library(lattice) numy=100 y=runif(numy,min=0,max=1) sig=0.05 numsig=length(which(y<sig)) tck.no=11 # number of ticks tcks=1:tck.no labl=as.character(0.1*tcks-0.1) # label for all ticks labl[seq(2,tck.no-1,2)]="" # delete

splom() doesn't complain here, but writes no names in the diagonal boxes. What am I missing? I believe that I need to add something like grid.text(x, ...) to the diagonal panel, but I don't know how to get it cycle through the column labels. And should varname.col = 'blue', varname.cex = 1 be inside the diag.panel() function? splom(szw[, n], pscales = 0, diag.panel

Dear List, Following the link below ( http://rgm2.lab.nig.ac.jp/RGM2/func.php?rd_id=plotrix:clock24.plot) I got an interesting polar plots which displayed my data and the time of observation. Thank you very much for providing such details. However, I have two set of data which I wish to display in the same polar plot. I tried using points to add the second data but could not succeed. That is,

Dear All, This problem is over. Clock24.plot did the job. Thanks to all those who assisted me. Ogbos On Apr 22, 2016 8:34 PM, "Ogbos Okike" <giftedlife2014 at gmail.com> wrote: > Dear All, > One hand. Many thanks!! The code run as soon as I loaded lubridate. > > Please can you guide me on how to relate this code to my actual data. > My actual data is looking like:

Dear list, I am using the ?clock24.plot command in this excellent package to plot animal activity data. Does anyone know if both symbols and a line can be plotted on the same plot to show both raw data (symbols) and a line (describing a statistical model of the pattern) ? Or if more than one line etc can be plotted? Thanks Paul