Rose Diagrams for Geology

Hello everyone,

In geology we often do rose diagrams showing the number of features along a
certain compass direction within a given range (bin) of angle (0-180
degrees).  I was wondering if anybody has had experience with this in R and
if they could recommend a package.

I looked at the circular package but it seems to deal only in radian and we
normally use degrees.

I've also looked a little at openair being rose diagrams are often used for
wind directions.

Any suggestions / guidance would be greatly appreciated.

findFn("rose diagram")

unique( findFn("rose diagram")$Package )

On Nov 18, 2014, at 1:42 PM, David Doyle wrote:

Hello everyone,

In geology we often do rose diagrams showing the number of features

along a

certain compass direction within a given range (bin) of angle (0-180
degrees).  I was wondering if anybody has had experience with this in R

and

if they could recommend a package.

I looked at the circular package but it seems to deal only in radian and

we

normally use degrees.

I've also looked a little at openair being rose diagrams are often used

for

wind directions.

Any suggestions / guidance would be greatly appreciated.

Learn to search, Grasshopper.

Choose one of:
-----------------
install.packages("sos")
library(sos)

findFn("rose diagram")

found 19 matches;  retrieving 1 page

Downloaded 14 links in 9 packages.

unique( findFn("rose diagram")$Package )

found 19 matches;  retrieving 1 page

Downloaded 14 links in 9 packages.
[1] "circular"  "GEOmap"    "oce"       "CircStats" "plotrix"   "HistData"
[7] "climatol"  "openair"   "xergm"
--------------
OR use Rseek.org
------------
OR  http://stackoverflow.com/search?q=[r]+rose+diagram

--

Thank you for your time.
David Doyle

      [[alternative HTML version deleted]]

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

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

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

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