[R] Rose Diagrams for Geology

[William Dunlap]

> Most field mapping data I have collected has been
> either in quadrant form (rarely) or more commonly in azimuthal form
> (0-360 degrees order clockwise from the top).

You can specify that when making the 'circular' object, by using the
zero and roation arguments.  Compare the plots made by the following:
  par(mfrow=c(1,2))
  A <- circular(rep(seq(0, 180, len=11), c(1:11)), units="degrees",
zero=pi/2, rotation="clock")
  rose.diag(A)
  B <- circular(rep(seq(0, 180, len=11), c(1:11)), units="degrees")
  rose.diag(B)



Bill Dunlap
TIBCO Software
wdunlap tibco.com

On Thu, Nov 20, 2014 at 11:57 PM, jwd <jwd at surewest.net> wrote:

> On Tue, 18 Nov 2014 22:06:03 -0600
> David Doyle <kydaviddoyle at gmail.com> wrote:
>
> > Thank you to David and David for their help.  The code below
> > generated what I needed.
> >
> >
> > library(circular)
> > mydata <- read.table("http://doylesdartden.com/R/Joints.csv",
> > header=TRUE, sep=",",)
> > x <- circular(mydata$JointsRad)
> > rose.diag(x,
> >
> >           #Set point character to use
> >           pch = 20,
> >           #sets font size
> >           cex = 1,
> >           #parameter that controls the size of the circle.
> >           #1= default <1 makes it larger > makes it smaller
> >           shrink = 1,
> >           #the color for filling the rose diagram.
> >           col=2,
> >           prop = 2,
> >           # number of bins.  36 = 10 degrees each.  18 = 20 degree
> > each bins=36,
> >           # Ticks showing bins
> >           ticks=TRUE,
> >           # Unites.
> >           units="degrees",
> >           # list main title
> >           main="Rose Diagram of XXX")
> > # for more info see
> > http://www.inside-r.org/packages/cran/circular/docs/rose.diag
> >
> I've been following this thread with some interest.  One problem that I
> might have with the code above is that as it is, the plot is labeled
> with 0-deg to the left, and numbered counter clockwise (standard
> trigonometric format). Most field mapping data I have collected has been
> either in quadrant form (rarely) or more commonly in azimuthal form
> (0-360 degrees order clockwise from the top).  Is that an issue?
>
> jwdougherty
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>

	[[alternative HTML version deleted]]