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)
factanal(x, 2)
factanal(x, 3)
x <- matrix(selected, 40, 4, byrow=T)
factanal(x, 1)
cov(x)
cor(x)
prcomp(x)