Performs Q-mode and R-mode factor analysis simultaneously on spatial and non-spatial data. A single function, qrfactor(), carries out principal component analysis, R-mode factor analysis, Q-mode factor analysis, simultaneous R- and Q-mode factor analysis, principal coordinate analysis and multidimensional scaling. 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.