sciruby 0.1.3 → 0.2.1

data/Rakefile

@@ -1,178 +0,0 @@
-#!/usr/bin/ruby
-# -*- ruby -*-
-# -*- coding: utf-8 -*-
-$:.unshift(File.dirname(__FILE__)+'/lib/')
-
-require 'rubygems'
-require 'rdoc/task'
-require 'sciruby'
-require 'hoe'
-
-# Hoe.plugin :compiler
-# Hoe.plugin :cucumberfeatures
-# Hoe.plugin :gem_prelude_sucks
-Hoe.plugin :git
-Hoe.plugin :gemspec
-Hoe.plugin :bundler
-# Hoe.plugin :inline
-# Hoe.plugin :manifest
-# Hoe.plugin :newgem
-# Hoe.plugin :racc
-# Hoe.plugin :rubyforge
-# Hoe.plugin :website
-
-desc "Ruby Lint"
-task :lint do
-  executable=Config::CONFIG['RUBY_INSTALL_NAME']
-  Dir.glob("lib/**/*.rb") {|f|
-    if !system %{#{executable} -w -c "#{f}"}
-        puts "Error on: #{f}"
-    end
-  }
-end
-
-desc "Open an irb session preloaded with sciruby"
-task :console do
-  sh "irb -rubygems -I lib -r sciruby.rb"
-end
-
-desc "Start the plotter without the console"
-task :plotter, [:script] => [] do |t,args|
-  if args.script.empty?
-    raise ArgumentError, "Need a script, e.g.: rake plotter[script.rb]"
-  else
-    sh "ruby -rubygems -I lib -r sciruby.rb -e 'SciRuby::Plotter.new(\"#{args.script}\")'"
-  end
-end
-
-desc "Start the plotter without the console"
-task :editor do
-  sh "ruby -rubygems -I lib -r sciruby.rb -e 'SciRuby::Editor.new'"
-end
-
-task :release do
-  system %{git push origin master}
-end
-
-h = Hoe.spec 'sciruby' do
-  self.version = SciRuby::VERSION
-  self.require_ruby_version ">=1.9"
-  self.developer('SciRuby Development Team', 'sciruby-dev@googlegroups.com')
-  self.extra_deps = {'distribution' => ">=0.4.0",
-                     'green_shoes' => ">=1.0.282",
-                     'statsample' => ">=1.1.0",
-                     'integration' => ">= 0",
-                     'minimization' => ">= 0",
-                     'gsl' => "~> 1.14.5",
-                     'rsvg2' => '~> 1.0.0',
-                     'simpler' => '>=0.1.0',
-                     'rubyvis' => '>=0.4.0'        }.to_a
-
-
-  self.extra_dev_deps = {'hoe' => "~> 2.12",
-                         'rdoc' => ">=0",
-                         'rspec' => ">=2.0",
-                         'haml' => ">=0", # for Rubyvis
-                         'coderay' => ">=0", # for Rubyvis
-                         'nokogiri' => ">=0", # for Rubyvis
-                         'RedCloth' => ">=0", # for Rubyvis
-                         'gtksourceview2' => ">=0", # for editor
-                         'shoulda' => "~> 2.11",
-                         'hoe-gemspec' => "~> 1.0",
-                         'hoe-bundler' => "~> 1.1",
-                         'minitest' => "~> 2.0" }.to_a
-
-  #self.executables = `git ls-files -- bin/*`.split("\n").map { |f| File.basename(f) }
-
-  self.post_install_message = <<-EOF
-***********************************************************
-Welcome to SciRuby: Tools for Scientific Computing in Ruby!
-
-                     *** WARNING ***
-Please be aware that SciRuby is in ALPHA status. If you're
-thinking of using SciRuby to write mission critical code,
-such as for driving a car or flying a space shuttle, you
-may wish to choose other software (for now).
-
-In order to leverage the GUI features, you will need to
-install gtk2 and optionally gtksourceview2:
-
-  $ gem install gtk2 gtksourceview2
-
-You will probably first need to install the headers for
-a number of required packages. In Ubuntu, use:
-
-  $ sudo apt-get install libgtk2.0-dev libgtksourceview2-dev \\
-      librsvg2-dev libcairo2-dev
-
-If you have trouble with Green Shoes, you should look at
-the Green Shoes wiki:
-
-* http://github.com/ashbb/green_shoes/wiki
-
-For Mac OSX Green Shoes:
-
-* https://github.com/ashbb/green_shoes/wiki/Building-Green-Shoes-on-OSX
-
-More explicit instructions for SciRuby should be available
-at our website, sciruby.com, or through our mailing list
-(which can also be found on our website).
-
-Thanks for installing SciRuby! Happy hypothesis testing!
-
-***********************************************************
-  EOF
-
-  self.need_rdoc = false
-
-end
-
-RDoc::Task.new(:docs) do |rd|
-  rd.main = h.readme_file
-  rd.options << '-d' if (`which dot` =~ /\/dot/) unless
-    ENV['NODOT'] || Hoe::WINDOZE
-  rd.rdoc_dir = 'doc'
-
-  rd.rdoc_files.include("lib/**/*.rb")
-  rd.rdoc_files += h.spec.extra_rdoc_files
-  rd.rdoc_files.reject! {|f| f=="Manifest.txt"}
-  title = h.spec.rdoc_options.grep(/^(-t|--title)=?$/).first
-  if title then
-    rd.options << title
-
-    unless title =~ /\=/ then # for ['-t', 'title here']
-    title_index = spec.rdoc_options.index(title)
-    rd.options << spec.rdoc_options[title_index + 1]
-    end
-  else
-    title = "#{h.name}-#{h.version} Documentation"
-    title = "#{h.rubyforge_name}'s " + title if h.rubyforge_name != h.name
-    rd.options << '--title' << title
-  end
-end
-
-desc 'Publish rdocs with analytics support'
-task :publish_docs => [:clean, :docs] do
-  ruby %{aggregate_adsense_to_doc.rb}
-  path = File.expand_path("~/.rubyforge/user-config.yml")
-  config = YAML.load(File.read(path))
-  host = "#{config["username"]}@rubyforge.org"
-
-  remote_dir = "/var/www/gforge-projects/#{h.rubyforge_name}/#{h.remote_rdoc_dir
-  }"
-  local_dir = h.local_rdoc_dir
-  Dir.glob(local_dir+"/**/*") {|file|
-    sh %{chmod 755 #{file}}
-  }
-  sh %{rsync #{h.rsync_args} #{local_dir}/ #{host}:#{remote_dir}}
-end
-
-require 'rspec/core/rake_task'
-namespace :spec do
-  desc "Run all specs"
-  RSpec::Core::RakeTask.new
-  # options in .rspec in package root
-end
-
-
-# vim: syntax=ruby

data/bin/sciruby-plotter

@@ -1,12 +0,0 @@
-#!/usr/bin/env ruby
-
-require 'rubygems'
-require 'sciruby'
-
-if ARGV.size == 0
-  STDERR.puts "sciruby-plotter: missing plot script"
-  STDERR.puts "Usage: sciruby-plotter myscript.rb"
-  exit false
-end
-
-SciRuby::Plotter.new ARGV.first

data/data/r/man/AirPassengers.Rd

@@ -1,51 +0,0 @@
-% File src/library/datasets/man/AirPassengers.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2007 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{AirPassengers}
-\docType{data}
-\alias{AirPassengers}
-\title{Monthly Airline Passenger Numbers 1949-1960}
-\description{
- The classic Box & Jenkins airline data.  Monthly totals of
- international airline passengers, 1949 to 1960.
-}
-\usage{AirPassengers}
-\format{
-  A monthly time series, in thousands.
-}
-\source{
-  Box, G. E. P., Jenkins, G. M. and Reinsel, G. C. (1976)
-  \emph{Time Series Analysis, Forecasting and Control.}
-  Third Edition. Holden-Day. Series G.
-}
- \examples{
-\dontrun{
-## These are quite slow and so not run by example(AirPassengers)
-
-## The classic 'airline model', by full ML
-(fit <- arima(log10(AirPassengers), c(0, 1, 1),
-              seasonal = list(order=c(0, 1 ,1), period=12)))
-update(fit, method = "CSS")
-update(fit, x=window(log10(AirPassengers), start = 1954))
-pred <- predict(fit, n.ahead = 24)
-tl <- pred$pred - 1.96 * pred$se
-tu <- pred$pred + 1.96 * pred$se
-ts.plot(AirPassengers, 10^tl, 10^tu, log = "y", lty = c(1,2,2))
-
-## full ML fit is the same if the series is reversed, CSS fit is not
-ap0 <- rev(log10(AirPassengers))
-attributes(ap0) <- attributes(AirPassengers)
-arima(ap0, c(0, 1, 1), seasonal = list(order=c(0, 1 ,1), period=12))
-arima(ap0, c(0, 1, 1), seasonal = list(order=c(0, 1 ,1), period=12),
-      method = "CSS")
-
-## Structural Time Series
-ap <- log10(AirPassengers) - 2
-(fit <- StructTS(ap, type= "BSM"))
-par(mfrow=c(1,2))
-plot(cbind(ap, fitted(fit)), plot.type = "single")
-plot(cbind(ap, tsSmooth(fit)), plot.type = "single")
-}}
-\keyword{datasets}

data/data/r/man/BJsales.Rd

@@ -1,34 +0,0 @@
-% File src/library/datasets/man/BJsales.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2007 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{BJsales}
-\docType{data}
-\alias{BJsales}
-\alias{BJsales.lead}
-\title{Sales Data with Leading Indicator}
-\description{
-  The sales time series \code{BJsales} and leading indicator
-  \code{BJsales.lead} each contain 150 observations.
-  The objects are of class \code{"ts"}.
-}
-\usage{
-BJsales
-BJsales.lead
-}
-\source{
-  The data are given in Box & Jenkins (1976).
-  Obtained from the Time Series Data Library at
-  \url{http://www-personal.buseco.monash.edu.au/~hyndman/TSDL/}
-}
-\references{
-  G. E. P. Box and G. M. Jenkins (1976):
-  \emph{Time Series Analysis, Forecasting and Control},
-  Holden-Day, San Francisco, p. 537.
-
-  P. J. Brockwell and R. A. Davis (1991):
-  \emph{Time Series: Theory and Methods},
-  Second edition, Springer Verlag, NY, pp. 414.
-}
-\keyword{datasets}

data/data/r/man/BOD.Rd

@@ -1,53 +0,0 @@
-% File src/library/datasets/man/BOD.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2010 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{BOD}
-\docType{data}
-\alias{BOD}
-\title{ Biochemical Oxygen Demand }
-\description{
-  The \code{BOD} data frame has 6 rows and 2 columns giving the
-  biochemical oxygen demand versus time in an evaluation of water
-  quality.
-}
-\usage{BOD}
-\format{
-  This data frame contains the following columns:
-  \describe{
-    \item{Time}{
-      A numeric vector giving the time of the measurement (days).
-    }
-    \item{demand}{
-      A numeric vector giving the biochemical oxygen demand (mg/l).
-    }
-  }
-}
-\source{
-  Bates, D.M. and Watts, D.G. (1988),
-  \emph{Nonlinear Regression Analysis and Its Applications},
-  Wiley, Appendix A1.4.
-
-  Originally from Marske (1967), \emph{Biochemical
-    Oxygen Demand Data Interpretation Using Sum of Squares Surface}
-  M.Sc. Thesis, University of Wisconsin -- Madison.
-}
-\examples{
-\testonly{options(show.nls.convergence=FALSE)
-old <- options(digits = 5)}
-require(stats)
-# simplest form of fitting a first-order model to these data
-fm1 <- nls(demand ~ A*(1-exp(-exp(lrc)*Time)), data = BOD,
-   start = c(A = 20, lrc = log(.35)))
-coef(fm1)
-fm1
-# using the plinear algorithm
-fm2 <- nls(demand ~ (1-exp(-exp(lrc)*Time)), data = BOD,
-   start = c(lrc = log(.35)), algorithm = "plinear", trace = TRUE)
-# using a self-starting model
-fm3 <- nls(demand ~ SSasympOrig(Time, A, lrc), data = BOD)
-summary(fm3)
-\testonly{options(old)}
-}
-\keyword{datasets}

data/data/r/man/ChickWeight.Rd

@@ -1,68 +0,0 @@
-% File src/library/datasets/man/ChickWeight.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2011 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{ChickWeight}
-\docType{data}
-\alias{ChickWeight}
-\title{Weight versus age of chicks on different diets}
-\description{
-  The \code{ChickWeight} data frame has 578 rows and 4 columns from an
-  experiment on the effect of diet on early growth of chicks.
-}
-\usage{ChickWeight}
-\format{
-  This object of class \code{c("nfnGroupedData", "nfGroupedData",
-    "groupedData", "data.frame")} containing the following columns:
-  \describe{
-    \item{weight}{
-      a numeric vector giving the body weight of the chick (gm).
-    }
-    \item{Time}{
-      a numeric vector giving the number of days since birth when
-      the measurement was made.
-    }
-    \item{Chick}{
-      an ordered factor with levels
-      \code{18} < \dots < \code{48}
-      giving a unique identifier for the chick.  The ordering of
-      the levels groups chicks on the same diet together and
-      orders them according to their final weight (lightest to
-      heaviest) within diet.
-    }
-    \item{Diet}{
-      a factor with levels 1,\dots,4 indicating which
-      experimental diet the chick received.
-    }
-  }
-}
-\details{
-  The body weights of the chicks were measured at birth and every
-  second day thereafter until day 20.  They were also measured on day
-  21.  There were four groups on chicks on different protein diets.
-  
-  This dataset was originally part of package \code{nlme}, and that has
-  methods (including for \code{[}, \code{as.data.frame}, \code{plot} and
-  \code{print}) for its grouped-data classes. 
-}
-\source{
-  Crowder, M. and Hand, D. (1990), \emph{Analysis of Repeated Measures},
-  Chapman and Hall (example 5.3)
-
-  Hand, D. and Crowder, M. (1996), \emph{Practical Longitudinal Data
-    Analysis}, Chapman and Hall (table A.2)
-
-  Pinheiro, J. C. and Bates, D. M. (2000) \emph{Mixed-effects Models in
-    S and S-PLUS}, Springer.
-}
-\seealso{
-  \code{\link{SSlogis}} for models fitted to this dataset.
-}
-
-\examples{\donttest{
-require(graphics)
-coplot(weight ~ Time | Chick, data = ChickWeight,
-       type = "b", show.given = FALSE)
-}}
-\keyword{datasets}

data/data/r/man/DNase.Rd

@@ -1,63 +0,0 @@
-% File src/library/datasets/man/DNase.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2011 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{DNase}
-\docType{data}
-\alias{DNase}
-\title{Elisa assay of DNase}
-\description{
-  The \code{DNase} data frame has 176 rows and 3 columns of data
-  obtained during development of an ELISA assay for the recombinant
-  protein DNase in rat serum.
-}
-\usage{DNase}
-\format{
-  This object of class \code{c("nfnGroupedData", "nfGroupedData",
-    "groupedData", "data.frame")} containing the following columns:
-  \describe{
-    \item{Run}{
-      an ordered factor with levels \code{10} < \dots < \code{3}
-      indicating the assay run.
-    }
-    \item{conc}{
-      a numeric vector giving the known concentration of the
-      protein. 
-    }
-    \item{density}{
-      a numeric vector giving the measured optical density
-      (dimensionless) in the assay.  Duplicate optical density
-      measurements were obtained. 
-    }
-  }
-}
-\details{
-  This dataset was originally part of package \code{nlme}, and that has
-  methods (including for \code{[}, \code{as.data.frame}, \code{plot} and
-  \code{print}) for its grouped-data classes. 
-}
-\source{
-  Davidian, M. and Giltinan, D. M. (1995) \emph{Nonlinear Models for
-    Repeated Measurement Data}, Chapman & Hall (section 5.2.4, p. 134)
-
-  Pinheiro, J. C. and Bates, D. M. (2000) \emph{Mixed-effects Models in
-    S and S-PLUS}, Springer.
-}
-\examples{
-require(stats); require(graphics)
-\testonly{options(show.nls.convergence=FALSE)}
-coplot(density ~ conc | Run, data = DNase,
-       show.given = FALSE, type = "b")
-coplot(density ~ log(conc) | Run, data = DNase,
-       show.given = FALSE, type = "b")
-## fit a representative run
-fm1 <- nls(density ~ SSlogis( log(conc), Asym, xmid, scal ),
-    data = DNase, subset = Run == 1)
-## compare with a four-parameter logistic
-fm2 <- nls(density ~ SSfpl( log(conc), A, B, xmid, scal ),
-    data = DNase, subset = Run == 1)
-summary(fm2)
-anova(fm1, fm2)
-}
-\keyword{datasets}

data/data/r/man/EuStockMarkets.Rd

@@ -1,28 +0,0 @@
-% File src/library/datasets/man/EuStockMarkets.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2007 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{EuStockMarkets}
-\docType{data}
-\alias{EuStockMarkets}
-\title{
-  Daily Closing Prices of Major European Stock Indices, 1991--1998
-}
-\description{
-  Contains the daily closing prices of major European stock indices:
-  Germany DAX (Ibis), Switzerland SMI, France CAC, and UK FTSE.  The
-  data are sampled in business time, i.e., weekends and holidays are
-  omitted.
-}
-\usage{
-EuStockMarkets
-}
-\format{
-  A multivariate time series with 1860 observations on 4 variables.
-  The object is of class \code{"mts"}.
-}
-\source{
-  The data were kindly provided by Erste Bank AG, Vienna, Austria.
-}
-\keyword{datasets}

data/data/r/man/Formaldehyde.Rd

@@ -1,44 +0,0 @@
-% File src/library/datasets/man/Formaldehyde.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2009 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{Formaldehyde}
-\docType{data}
-\alias{Formaldehyde}
-\title{Determination of Formaldehyde}
-\description{
-  These data are from a chemical experiment to prepare a standard curve
-  for the determination of formaldehyde by the addition of chromatropic
-  acid and concentrated sulphuric acid and the reading of the resulting
-  purple color on a spectrophotometer.
-}
-\usage{Formaldehyde}
-\format{
-  A data frame with 6 observations on 2 variables.
-  \tabular{rlll}{
-    [,1] \tab carb\tab numeric \tab Carbohydrate (ml) \cr
-    [,2] \tab optden \tab numeric \tab Optical Density
-  }
-}
-\source{
-  Bennett, N. A. and N. L. Franklin (1954)
-  \emph{Statistical Analysis in Chemistry and the Chemical Industry}.
-  New York: Wiley.
-}
-\references{
-  McNeil, D. R. (1977) \emph{Interactive Data Analysis.}
-  New York: Wiley.
-}
-\examples{
-require(stats); require(graphics)
-plot(optden ~ carb, data = Formaldehyde,
-     xlab = "Carbohydrate (ml)", ylab = "Optical Density",
-     main = "Formaldehyde data", col = 4, las = 1)
-abline(fm1 <- lm(optden ~ carb, data = Formaldehyde))
-summary(fm1)
-opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
-plot(fm1)
-par(opar)
-}
-\keyword{datasets}

data/data/r/man/HairEyeColor.Rd

@@ -1,77 +0,0 @@
-% File src/library/datasets/man/HairEyeColor.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2007 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{HairEyeColor}
-\docType{data}
-\alias{HairEyeColor}
-\title{Hair and Eye Color of Statistics Students}
-\description{
-  Distribution of hair and eye color and sex in 592 statistics students.
-}
-\usage{HairEyeColor}
-\format{
-  A 3-dimensional array resulting from cross-tabulating 592 observations
-  on 3 variables.  The variables and their levels are as follows:
-
-  \tabular{rll}{
-    No \tab Name \tab Levels \cr
-    1 \tab Hair \tab Black, Brown, Red, Blond \cr
-    2 \tab Eye \tab Brown, Blue, Hazel, Green \cr
-    3 \tab Sex \tab Male, Female
-  }
-}
-\details{
-  The Hair \eqn{\times}{x} Eye table comes rom a survey of students at
-  the University of Delaware reported by Snee (1974).  The split by
-  \code{Sex} was added by Friendly (1992a) for didactic purposes.
-
-  This data set is useful for illustrating various techniques for the
-  analysis of contingency tables, such as the standard chi-squared test
-  or, more generally, log-linear modelling, and graphical methods such
-  as mosaic plots, sieve diagrams or association plots.
-}
-\source{
-  \url{http://euclid.psych.yorku.ca/ftp/sas/vcd/catdata/haireye.sas}
-
-  Snee (1974) gives the two-way table aggregated over \code{Sex}.  The
-  Sex split of the \sQuote{Brown hair, Brown eye} cell was changed in
-  \R 2.6.0 to agree with that used by Friendly (2000).
-}
-\references{
-  Snee, R. D. (1974)
-  Graphical display of two-way contingency tables.
-  \emph{The American Statistician}, \bold{28}, 9--12.
-
-  Friendly, M. (1992a)
-  Graphical methods for categorical data.
-  \emph{SAS User Group International Conference Proceedings}, \bold{17},
-  190--200.
-  \url{http://www.math.yorku.ca/SCS/sugi/sugi17-paper.html}
-
-  Friendly, M. (1992b)
-  Mosaic displays for loglinear models.
-  \emph{Proceedings of the Statistical Graphics Section},
-  American Statistical Association, pp. 61--68.
-  \url{http://www.math.yorku.ca/SCS/Papers/asa92.html}
-
-  Friendly, M. (2000)
-  \emph{Visualizing Categorical Data.}
-  SAS Institute, ISBN 1-58025-660-0.
-}
-\seealso{
-  \code{\link{chisq.test}},
-  \code{\link{loglin}},
-  \code{\link{mosaicplot}}
-}
-\examples{
-require(graphics)
-## Full mosaic
-mosaicplot(HairEyeColor)
-## Aggregate over sex (as in Snee's original data)
-x <- apply(HairEyeColor, c(1, 2), sum)
-x
-mosaicplot(x, main = "Relation between hair and eye color")
-}
-\keyword{datasets}

data/data/r/man/Harman23.cor.Rd

@@ -1,25 +0,0 @@
-% File src/library/datasets/man/Harman23.cor.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2007 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{Harman23.cor}
-\docType{data}
-\alias{Harman23.cor}
-\title{Harman Example 2.3}
-\description{
-  A correlation matrix of eight physical measurements on 305 girls between
-  ages seven and seventeen.
-}
-\usage{Harman23.cor}
-\source{
-  Harman, H. H. (1976)
-  \emph{Modern Factor Analysis}, Third Edition Revised,
-  University of Chicago Press, Table 2.3.
-}
-\examples{
-require(stats)
-(Harman23.FA <- factanal(factors = 1, covmat = Harman23.cor))
-for(factors in 2:4) print(update(Harman23.FA, factors = factors))
-}
-\keyword{datasets}

data/data/r/man/Harman74.cor.Rd

@@ -1,28 +0,0 @@
-% File src/library/datasets/man/Harman74.cor.Rd
-% Part of the R package, http://www.R-project.org
-% Copyright 1995-2007 R Core Development Team
-% Distributed under GPL 2 or later
-
-\name{Harman74.cor}
-\docType{data}
-\alias{Harman74.cor}
-\title{Harman Example 7.4}
-\description{
-  A correlation matrix of 24 psychological tests given to 145 seventh and
-  eight-grade children in a Chicago suburb by Holzinger and Swineford.
-}
-\usage{Harman74.cor}
-\source{
-  Harman, H. H. (1976)
-  \emph{Modern Factor Analysis}, Third Edition Revised,
-  University of Chicago Press, Table 7.4.
-}
-\examples{
-require(stats)
-(Harman74.FA <- factanal(factors = 1, covmat = Harman74.cor))
-for(factors in 2:5) print(update(Harman74.FA, factors = factors))
-Harman74.FA <- factanal(factors = 5, covmat = Harman74.cor,
-                        rotation="promax")
-print(Harman74.FA$loadings, sort = TRUE)
-}
-\keyword{datasets}