**FabricPlotR** *See [change log here][1] for recent updates to scripts and packages.* This R script and supporting package allows the plotting and manipulation of c-axis orientation data imported as azimuth and inclination pairs (trend and plunge - see example datafile for format). Once loaded the data can be output as a scatter plot of *c*-axis orientations or a kernel density contour plot (filled or unfilled). The script will automatically determine the opening angle of crossed-girdle fabrics based on user-defined input. The script allows the rotation of data about any axis and, though it is focused on *c*-axis data, the script can be used with any line data. The FabricPlotR package automatically calculates opening angle temperature and fabric strength parameters (see below) and provides control on contour plotting parameters. In addition, there is a separate batch script to process entire folders of data files at once. The scripts require the following R libraries to be installed and loaded in order to run: - ggplot2 - *for plotting* - viridis - *for colour palette* - Directional - *for kernel density estimations and contours* - mapproj - *to project spherical data in 2D* - colorspace - *simple colour gradient implementation* - FabricPlotR - *custom functions called in the script* *Note that the [FabricPlotR package][2] must be downloaded from the [***Files***][3] herein and installed as a local library through your R app.* In addition if you are running R on MacOS, you will need to download and install [XQuartz][4] to enable all required functions. Plots generated are lower hemisphere approximations in an equal area (i.e. Schmidt) projection by default. The user can also choose to plot the data can be plotted in an equal angle (Wulff) projections. Density calculations use a von Mises-Fisher kernel on spherical data, as implemented by the vmf.kerncontour function in the [Directional R package][5]. Fabric opening angle determination is informed by the approach of [Hunter et al. (2018)][6] but is entirely original script. Quartz CPO opening angle temperature is calculated using the non-linear thermometer defined in Equation 3 of [Faleiros et al. (2016)][7]. Cylindricity and associated values are calculated as in [Vollmer (1990)][8]. *Intensity* is calculated as in [Lisle (1985)][9]. Uniformity Statistic is as proposed by [Mardia (1972)][10]. Density norm is calculated as the L2 norm of a von Mises-Fisher KDE for spherical data, as implemented by the vmf.kde function in the [Directional R package][11], using a user-input half-width normalized to a homogeneous distribution. Below are example output plots. Quartz *c*-axis opening angle plot (below) ![Quartz *c*-axis opening angle plot][12] Quartz c-axis scatter plot (below) !["Quartz *c*-axis scatter plot"][13] Quartz *c*-axis kernel density contour plot (note - some aspects of this plot have changed in recent versions) ![Quartz c-axis kernel density contour plot][14] [1]: https://osf.io/mjtxv/wiki/Change%20Log/ [2]: https://osf.io/g6xw4/ [3]: https://osf.io/mjtxv/files/ [4]: https://www.xquartz.org/ [5]: https://rdrr.io/cran/Directional/ [6]: http://dx.doi.org/10.1016/j.jsg.2018.04.006 [7]: http://dx.doi.org/10.1016/j.tecto.2016.01.014 [8]: https://doi.org/10.1130/0016-7606(1990)102%3C0786:AAOEMT%3E2.3.CO;2 [9]: https://doi.org/10.1016/0191-8141%2885%2990119-1 [10]: https://www.sciencedirect.com/book/9780124711501/statistics-of-directional-data [11]: https://rdrr.io/cran/Directional/ [12]: https://upload.wikimedia.org/wikipedia/commons/thumb/b/bd/Quartz_c-axis_opening_angle_plot.png/600px-Quartz_c-axis_opening_angle_plot.png [13]: https://upload.wikimedia.org/wikipedia/commons/thumb/9/93/Quartz_c-axis_scatter_plot.png/600px-Quartz_c-axis_scatter_plot.png "Quartz c-axis scatter plot" [14]: https://upload.wikimedia.org/wikipedia/commons/thumb/5/53/Quartz_c-axis_contour_plot.png/600px-Quartz_c-axis_contour_plot.png