R version 2.4.1 (2006-12-18)
Copyright (C) 2006 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

  Natural language support but running in an English locale

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> columns <- c(
+   "Expressiveness",
+   "Structure",
+   "Granularity",
+   "Intended Use",
+   "Automated Reasoning",
+   "Prescriptive vs. Descriptive",
+   "Governance"
+ )
> ontologies <- c(
+   "MeSH",
+   "Gene Ontology",
+   "Basic Formal Ontology",
+   "SUMO",
+   "Engineering Math",
+   "COSMO",
+   "PSL",
+   "Environmental Data Coding Specification",
+   "Geospatial ML",
+   "International Classification of Diseases",
+   "Unified Medical Language System Metathesaurus",
+   "World Bank MetaThesaurus",
+   "Spase Ontology",
+   "FOAF",
+   "AKTiveSAOntology",
+   "CMMI Smallscale Ontology",
+   "Beer Ontology",
+   "BioMoby Services",
+   "Bio-Zen",
+   "UNSPCS",
+   "IMDB Mapping Movie Ontology",
+   "Federal Enterprise Architecture Reference Model Ontology (FEA-RMO)",
+   "Cyc",
+   "Encyclopedia of Life",
+   "Ontology for Biomedical Investigation",
+   "Functional Genomics Investigation Ontology",
+   "ProMed Mail",
+   "Agricola",
+   "GenBank",
+   "Medline",
+   "GBIF",
+   "WHO health reports",
+   "Human Disease - OBO",
+   "Core Software Ontology",
+   "Cell System Ontology",
+   "CIA Factbook",
+   "Evidence Code",
+   "WordNet",
+   "Ontology for Geographic Information ISO 19115:2003",
+   "General Ontology for Linguistic Description"
+ )
> evaluations <- c(
+   2,3,4,3,1,1,3,
+   2.5,4,4,3,2,5,3,
+   5,4,1,2,5,5,5,
+   5,4,4,4,5,5,5,
+   5,5,3,5,5,5,5,
+   5,5,3,4,4,5,5,
+   5,4,2,4,5,5,5,
+   3,4,2,4,4,5,5,
+   3,4,4,4,4,4,1,
+   0.5,2,4,3,1,5,5,
+   2,2,5,1,1,5,3,
+   2,3,5,1,1,5,5,
+   5,2,2,1,2,5,5,
+   2.5,2,1,3,3,1,3,
+   5,4,3,3,2,1,3,
+   5,4,2,3,2,5,3,
+   5,4,2,1,2,5,1,
+   4,4,2,1,4,5,5,
+   5,4,2,2,3,5,3,
+   4,4,3,3,2,5,5,
+   4,4,1,1,2,5,1,
+   5,4,2,1,4,5,3,
+   5,5,5,1,5,5,3,
+   2,2,5,1,1,5,1,
+   5,3,2,1,2,5,5,
+   5,3,2,1,2,5,5,
+   1,2,4,3,2,5,1,
+   2,2,5,3,2,5,5,
+   1,2,5,3,2,5,5,
+   2,2,5,1,1,5,5,
+   2,2,5,1,2,5,5,
+   1,2,4,1,1,5,5,
+   5,5,4,1,4,5,5,
+   5,3,1,1,3,5,5,
+   5,5,2,1,5,5,5,
+   5,2,2,1,2,5,5,
+   5,3,2,1,2,5,5,
+   1.5,2,2,1,1,1,3,
+   5,3,2,1,3,5,5,
+   5,3,2,1,3,5,3
+ )
> r<-7*(0:39)
> s<-sort(c(r+1,r+2,r+3,r+5))
> selected<-evaluations[s]
> x <- matrix(evaluations, 40, 7, byrow=T)
> factanal(x, 1)

Call:
factanal(x = x, factors = 1)

Uniquenesses:
[1] 0.431 0.269 0.836 0.940 0.284 0.958 0.984

Loadings:
     Factor1
[1,]  0.755 
[2,]  0.855 
[3,] -0.405 
[4,]  0.244 
[5,]  0.846 
[6,]  0.206 
[7,]  0.128 

               Factor1
SS loadings      2.299
Proportion Var   0.328

Test of the hypothesis that 1 factor is sufficient.
The chi square statistic is 34.61 on 14 degrees of freedom.
The p-value is 0.00168 
> factanal(x, 2)

Call:
factanal(x = x, factors = 2)

Uniquenesses:
[1] 0.005 0.251 0.606 0.618 0.268 0.955 0.977

Loadings:
     Factor1 Factor2
[1,]  0.992         
[2,]  0.594   0.629 
[3,] -0.619   0.102 
[4,] -0.148   0.600 
[5,]  0.537   0.666 
[6,]  0.202         
[7,]  0.150         

               Factor1 Factor2
SS loadings      2.094   1.226
Proportion Var   0.299   0.175
Cumulative Var   0.299   0.474

Test of the hypothesis that 2 factors are sufficient.
The chi square statistic is 12.27 on 8 degrees of freedom.
The p-value is 0.139 
> factanal(x, 3)

Call:
factanal(x = x, factors = 3)

Uniquenesses:
[1] 0.005 0.258 0.552 0.566 0.262 0.005 0.919

Loadings:
     Factor1 Factor2 Factor3
[1,]  0.932   0.218   0.281 
[2,]  0.480   0.688   0.195 
[3,] -0.667                 
[4,] -0.184   0.609  -0.171 
[5,]  0.401   0.719   0.246 
[6,]                  0.995 
[7,]                  0.274 

               Factor1 Factor2 Factor3
SS loadings       1.75   1.409   1.274
Proportion Var    0.25   0.201   0.182
Cumulative Var    0.25   0.451   0.633

Test of the hypothesis that 3 factors are sufficient.
The chi square statistic is 4.64 on 3 degrees of freedom.
The p-value is 0.2 
> x <- matrix(selected, 40, 4, byrow=T)
> factanal(x, 1)

Call:
factanal(x = x, factors = 1)

Uniquenesses:
[1] 0.390 0.276 0.807 0.327

Loadings:
     Factor1
[1,]  0.781 
[2,]  0.851 
[3,] -0.439 
[4,]  0.820 

               Factor1
SS loadings       2.20
Proportion Var    0.55

Test of the hypothesis that 1 factor is sufficient.
The chi square statistic is 12.95 on 2 degrees of freedom.
The p-value is 0.00154 
> cov(x)
          [,1]       [,2]       [,3]       [,4]
[1,]  2.481410  1.0916667 -1.2948718  1.2891026
[2,]  1.091667  1.1275641 -0.3589744  1.0660256
[3,] -1.294872 -0.3589744  1.8461538 -0.5641026
[4,]  1.289103  1.0660256 -0.5641026  1.8660256
> cor(x)
           [,1]       [,2]       [,3]       [,4]
[1,]  1.0000000  0.6526347 -0.6049835  0.5990725
[2,]  0.6526347  1.0000000 -0.2488049  0.7349169
[3,] -0.6049835 -0.2488049  1.0000000 -0.3039242
[4,]  0.5990725  0.7349169 -0.3039242  1.0000000
> prcomp(x)
Standard deviations:
[1] 2.2085742 1.2262670 0.7933596 0.5569594

Rotation:
            PC1        PC2        PC3        PC4
[1,]  0.6630407 -0.1399413 -0.6618071 -0.3206320
[2,]  0.3737721  0.3831627 -0.1121577  0.8371985
[3,] -0.4194978  0.7440946 -0.4726879 -0.2165894
[4,]  0.4946646  0.5290803  0.5709625 -0.3865007
>