qrfactor

Simultaneous Q-mode and R-mode factor analysis for spatial data.

qrfactor performs Q-mode and R-mode factor analysis simultaneously on spatial and non-spatial data. A single function, qrfactor(), carries out

• principal component analysis (PCA),
• R-mode factor analysis,
• Q-mode factor analysis,
• simultaneous R- and Q-mode factor analysis,
• principal coordinate analysis (PCoA), and
• multidimensional scaling (MDS).

Loadings and scores are returned from the fitted object, and the plot() method provides annotated biplots for combinations of eigenvectors, loadings and scores. Input may be supplied as an ESRI shapefile, a delimited text file or a data frame.

Installation

Install the development version from GitHub:

# install.packages("remotes")
# remotes::install_github("gowusu/qrfactor")
install.packages("qrfactor", repos = c('https://gowusu.r-universe.dev', 'https://cloud.r-project.org'))

For the spatial (shapefile) workflow you will also need the sf and sp packages:

install.packages(c("sf", "sp"))

Documentation website

Full documentation — the complete manual plus a function reference — is online at https://gowusu.github.io/qrfactor/.

Manual

A full, worked manual — every input mode, all six analyses, the plot method, the spatial maps and diagnostics, with every number and figure produced by running the package on the bundled African freshwater data:

• Read it online — on the documentation website
• Source on GitHub (Markdown) — renders in the repo
• manual/qrfactor_manual.docx — Word version (Save As → PDF for a printable copy)
• manual/reproduce_manual.R — regenerates every table and figure

Quick start

Non-spatial (data frame) example:

library(qrfactor)
data(UScereal, package = "MASS")
v <- c("calories", "protein", "sodium", "carbo", "sugars", "potassium")
m <- qrfactor(UScereal[v], scale = "n")
plot(m, plot = "r", type = "loading")

Spatial (shapefile) example, using the bundled African freshwater data:

p <- system.file("external", package = "qrfactor")
m <- qrfactor(p, layer = "Africanfreshwater",
              var = c("Domestic", "Industry", "Agricultur"))
plot(m, plot = "map")

Status

This is the 1.5 restoration of the package (originally archived from CRAN in 2018). The archival cause — an unregistered S3 method — has been fixed, and the retired rgdal / maptools / mgraph dependencies have been replaced with an sf/sp compatibility layer.

Author

George Owusu (gowusu@gmail.com)

License

GPL-2.