Repository 'principal_component_analysis'
hg clone https://toolshed.g2.bx.psu.edu/repos/devteam/principal_component_analysis

Changeset 0:f568051cdf2e (2014-05-19)
Commit message:
Imported from capsule None
added:
pca.py
pca.xml
test-data/1.tabular
test-data/iris.tabular
test-data/pca_out1.tabular
test-data/pca_out2.pdf
test-data/pca_out3.tabular
test-data/pca_out4.pdf
test-data/pca_out5.tabular
test-data/pca_out6.pdf
tool_dependencies.xml
b
diff -r 000000000000 -r f568051cdf2e pca.py
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/pca.py Mon May 19 12:34:43 2014 -0400
[
@@ -0,0 +1,129 @@
+#!/usr/bin/env python
+
+import sys, string
+from rpy import *
+import numpy
+
+def stop_err(msg):
+    sys.stderr.write(msg)
+    sys.exit()
+
+infile = sys.argv[1]
+x_cols = sys.argv[2].split(',')
+method = sys.argv[3]
+outfile = sys.argv[4]
+outfile2 = sys.argv[5]
+
+if method == 'svd':
+    scale = center = "FALSE"
+    if sys.argv[6] == 'both':
+        scale = center = "TRUE"
+    elif sys.argv[6] == 'center':
+        center = "TRUE"
+    elif sys.argv[6] == 'scale':
+        scale = "TRUE"
+    
+fout = open(outfile,'w')
+elems = []
+for i, line in enumerate( file ( infile )):
+    line = line.rstrip('\r\n')
+    if len( line )>0 and not line.startswith( '#' ):
+        elems = line.split( '\t' )
+        break 
+    if i == 30:
+        break # Hopefully we'll never get here...
+
+if len( elems )<1:
+    stop_err( "The data in your input dataset is either missing or not formatted properly." )
+
+x_vals = []
+
+for k,col in enumerate(x_cols):
+    x_cols[k] = int(col)-1
+    x_vals.append([])
+
+NA = 'NA'
+skipped = 0
+for ind,line in enumerate( file( infile )):
+    if line and not line.startswith( '#' ):
+        try:
+            fields = line.strip().split("\t")
+            valid_line = True
+            for k,col in enumerate(x_cols):
+                try:
+                    xval = float(fields[col])
+                except:
+                    skipped += 1 
+                    valid_line = False
+                    break
+            if valid_line:
+                for k,col in enumerate(x_cols):
+                    xval = float(fields[col])
+                    x_vals[k].append(xval)
+        except:
+            skipped += 1
+
+x_vals1 = numpy.asarray(x_vals).transpose()
+dat= r.list(array(x_vals1))
+
+set_default_mode(NO_CONVERSION)
+try:
+    if method == "cor":
+        pc = r.princomp(r.na_exclude(dat), cor = r("TRUE"))
+    elif method == "cov":
+        pc = r.princomp(r.na_exclude(dat), cor = r("FALSE"))
+    elif method=="svd":
+        pc = r.prcomp(r.na_exclude(dat), center = r(center), scale = r(scale))
+except RException, rex:
+    stop_err("Encountered error while performing PCA on the input data: %s" %(rex))
+
+set_default_mode(BASIC_CONVERSION)
+summary = r.summary(pc, loadings="TRUE")
+ncomps = len(summary['sdev'])
+
+if type(summary['sdev']) == type({}):
+    comps_unsorted = summary['sdev'].keys()
+    comps=[]
+    sd = summary['sdev'].values()
+    for i in range(ncomps):
+        sd[i] = summary['sdev'].values()[comps_unsorted.index('Comp.%s' %(i+1))]
+        comps.append('Comp.%s' %(i+1))
+elif type(summary['sdev']) == type([]):
+    comps=[]
+    for i in range(ncomps):
+        comps.append('Comp.%s' %(i+1))
+        sd = summary['sdev']
+
+print >>fout, "#Component\t%s" %("\t".join(["%s" % el for el in range(1,ncomps+1)]))
+print >>fout, "#Std. deviation\t%s" %("\t".join(["%.4g" % el for el in sd]))
+total_var = 0
+vars = []
+for s in sd:
+    var = s*s
+    total_var += var
+    vars.append(var)
+for i,var in enumerate(vars):
+    vars[i] = vars[i]/total_var
+       
+print >>fout, "#Proportion of variance explained\t%s" %("\t".join(["%.4g" % el for el in vars]))
+
+print >>fout, "#Loadings\t%s" %("\t".join(["%s" % el for el in range(1,ncomps+1)]))
+xcolnames = ["c%d" %(el+1) for el in x_cols]
+if 'loadings' in summary: #in case of princomp
+    loadings = 'loadings'
+elif 'rotation' in summary: #in case of prcomp
+    loadings = 'rotation'
+for i,val in enumerate(summary[loadings]):
+    print >>fout, "%s\t%s" %(xcolnames[i], "\t".join(["%.4g" % el for el in val]))
+
+print >>fout, "#Scores\t%s" %("\t".join(["%s" % el for el in range(1,ncomps+1)]))
+if 'scores' in summary: #in case of princomp
+    scores = 'scores'
+elif 'x' in summary: #in case of prcomp
+    scores = 'x'
+for obs,sc in enumerate(summary[scores]):
+    print >>fout, "%s\t%s" %(obs+1, "\t".join(["%.4g" % el for el in sc]))
+
+r.pdf( outfile2, 8, 8 )
+r.biplot(pc)
+r.dev_off()
\ No newline at end of file
b
diff -r 000000000000 -r f568051cdf2e pca.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/pca.xml Mon May 19 12:34:43 2014 -0400
b
@@ -0,0 +1,102 @@
+<tool id="pca1" name="Principal Component Analysis" version="1.0.2">
+  <description> </description>
+  <requirements>
+    <requirement type="package" version="1.0.3">rpy</requirement>
+    <requirement type="package" version="2.11.0">R</requirement>
+    <requirement type="package" version="1.7.1">numpy</requirement>
+  </requirements>
+  <command interpreter="python">
+    pca.py 
+      $input1
+      $var_cols
+      $methodChoice.method
+      $out_file1
+      $out_file2
+      #if $methodChoice.method == "svd":
+      $methodChoice.scale
+      #end if
+  </command>
+  <inputs>
+    <param format="tabular" name="input1" type="data" label="Select data" help="Dataset missing? See TIP below."/>
+    <param name="var_cols" label="Select columns containing input variables " type="data_column" data_ref="input1" numerical="True" multiple="true" >
+        <validator type="no_options" message="Please select at least one column."/>
+    </param>
+    <conditional name="methodChoice">
+        <param name="method" type="select" label="Method" help="The correlation matrix can only be used if there are no constant variables">
+            <option value="cor" selected="true">Eigenvectors of Correlation (princomp)</option>
+            <option value="cov">Eigenvectors of Covariance (princomp)</option>
+            <option value="svd">Singular Value Decomposition (prcomp)</option>
+        </param>
+        <when value="cor" />
+        <when value="cov" />
+        <when value="svd">
+            <param name="scale" type="select" label="Centering and Scaling" help="Can be used to center and/or scale variables">
+                <option value="none" selected="true">None</option>
+                <option value="center">Center only</option>
+                <option value="scale">Scale only</option>
+                <option value="both">Center and Scale</option>
+            </param>        
+        </when>
+    </conditional>
+  </inputs>
+  <outputs>
+    <data format="input" name="out_file1" metadata_source="input1" />
+    <data format="pdf" name="out_file2" />
+  </outputs>
+  <tests>
+    <test>
+        <param name="input1" value="iris.tabular"/>
+        <param name="var_cols" value="1,2,3,4"/>
+        <param name="method" value="cor"/>
+        <output name="out_file1" file="pca_out1.tabular"/>
+        <output name="out_file2" file="pca_out2.pdf"/>
+    </test>
+    <test>
+        <param name="input1" value="iris.tabular"/>
+        <param name="var_cols" value="1,2,3,4"/>
+        <param name="method" value="cov"/>
+        <output name="out_file1" file="pca_out3.tabular"/>
+        <output name="out_file2" file="pca_out4.pdf"/>
+    </test>
+    <test>
+        <param name="input1" value="iris.tabular"/>
+        <param name="var_cols" value="1,2,3,4"/>
+        <param name="method" value="svd"/>
+        <param name="scale" value="both"/>
+        <output name="out_file1" file="pca_out5.tabular"/>
+        <output name="out_file2" file="pca_out6.pdf"/>
+    </test>
+  </tests>
+  <help>
+
+
+.. class:: infomark
+
+**TIP:** If your data is not TAB delimited, use *Edit Datasets-&gt;Convert characters*
+
+-----
+
+.. class:: infomark
+
+**What it does**
+
+This tool performs Principal Component Analysis on the given numeric input data using functions from R statistical package - 'princomp' function (for Eigenvector based solution) and 'prcomp' function (for Singular value decomposition based solution). It outputs two files, one containing the summary statistics of PCA, and the other containing biplots of the observations and principal components.   
+
+*R Development Core Team (2009). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL http://www.R-project.org.*
+
+-----
+
+.. class:: warningmark
+
+**Note**
+
+- This tool currently treats all variables as continuous numeric variables. Running the tool on categorical variables might result in incorrect results. Rows containing non-numeric (or missing) data in any of the chosen columns will be skipped from the analysis.
+
+- The summary statistics in the output are described below:
+
+  - Std. deviation: Standard deviations of the principal components
+  - Loadings: a list of eigen-vectors/variable loadings
+  - Scores: Scores of the input data on the principal components
+
+  </help>
+</tool>
b
diff -r 000000000000 -r f568051cdf2e test-data/1.tabular
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/1.tabular Mon May 19 12:34:43 2014 -0400
b
@@ -0,0 +1,6 @@
+chr22 1000 NM_17
+chr22 2000 NM_18
+chr10 2200 NM_10
+chr10 hap test
+chr10 1200 NM_11
+chr22 1600 NM_19
\ No newline at end of file
b
diff -r 000000000000 -r f568051cdf2e test-data/iris.tabular
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/iris.tabular Mon May 19 12:34:43 2014 -0400
b
@@ -0,0 +1,151 @@
+5.1 3.5 1.4 0.2 Iris-setosa
+4.9 3.0 1.4 0.2 Iris-setosa
+4.7 3.2 1.3 0.2 Iris-setosa
+4.6 3.1 1.5 0.2 Iris-setosa
+5.0 3.6 1.4 0.2 Iris-setosa
+5.4 3.9 1.7 0.4 Iris-setosa
+4.6 3.4 1.4 0.3 Iris-setosa
+5.0 3.4 1.5 0.2 Iris-setosa
+4.4 2.9 1.4 0.2 Iris-setosa
+4.9 3.1 1.5 0.1 Iris-setosa
+5.4 3.7 1.5 0.2 Iris-setosa
+4.8 3.4 1.6 0.2 Iris-setosa
+4.8 3.0 1.4 0.1 Iris-setosa
+4.3 3.0 1.1 0.1 Iris-setosa
+5.8 4.0 1.2 0.2 Iris-setosa
+5.7 4.4 1.5 0.4 Iris-setosa
+5.4 3.9 1.3 0.4 Iris-setosa
+5.1 3.5 1.4 0.3 Iris-setosa
+5.7 3.8 1.7 0.3 Iris-setosa
+5.1 3.8 1.5 0.3 Iris-setosa
+5.4 3.4 1.7 0.2 Iris-setosa
+5.1 3.7 1.5 0.4 Iris-setosa
+4.6 3.6 1.0 0.2 Iris-setosa
+5.1 3.3 1.7 0.5 Iris-setosa
+4.8 3.4 1.9 0.2 Iris-setosa
+5.0 3.0 1.6 0.2 Iris-setosa
+5.0 3.4 1.6 0.4 Iris-setosa
+5.2 3.5 1.5 0.2 Iris-setosa
+5.2 3.4 1.4 0.2 Iris-setosa
+4.7 3.2 1.6 0.2 Iris-setosa
+4.8 3.1 1.6 0.2 Iris-setosa
+5.4 3.4 1.5 0.4 Iris-setosa
+5.2 4.1 1.5 0.1 Iris-setosa
+5.5 4.2 1.4 0.2 Iris-setosa
+4.9 3.1 1.5 0.1 Iris-setosa
+5.0 3.2 1.2 0.2 Iris-setosa
+5.5 3.5 1.3 0.2 Iris-setosa
+4.9 3.1 1.5 0.1 Iris-setosa
+4.4 3.0 1.3 0.2 Iris-setosa
+5.1 3.4 1.5 0.2 Iris-setosa
+5.0 3.5 1.3 0.3 Iris-setosa
+4.5 2.3 1.3 0.3 Iris-setosa
+4.4 3.2 1.3 0.2 Iris-setosa
+5.0 3.5 1.6 0.6 Iris-setosa
+5.1 3.8 1.9 0.4 Iris-setosa
+4.8 3.0 1.4 0.3 Iris-setosa
+5.1 3.8 1.6 0.2 Iris-setosa
+4.6 3.2 1.4 0.2 Iris-setosa
+5.3 3.7 1.5 0.2 Iris-setosa
+5.0 3.3 1.4 0.2 Iris-setosa
+7.0 3.2 4.7 1.4 Iris-versicolor
+6.4 3.2 4.5 1.5 Iris-versicolor
+6.9 3.1 4.9 1.5 Iris-versicolor
+5.5 2.3 4.0 1.3 Iris-versicolor
+6.5 2.8 4.6 1.5 Iris-versicolor
+5.7 2.8 4.5 1.3 Iris-versicolor
+6.3 3.3 4.7 1.6 Iris-versicolor
+4.9 2.4 3.3 1.0 Iris-versicolor
+6.6 2.9 4.6 1.3 Iris-versicolor
+5.2 2.7 3.9 1.4 Iris-versicolor
+5.0 2.0 3.5 1.0 Iris-versicolor
+5.9 3.0 4.2 1.5 Iris-versicolor
+6.0 2.2 4.0 1.0 Iris-versicolor
+6.1 2.9 4.7 1.4 Iris-versicolor
+5.6 2.9 3.6 1.3 Iris-versicolor
+6.7 3.1 4.4 1.4 Iris-versicolor
+5.6 3.0 4.5 1.5 Iris-versicolor
+5.8 2.7 4.1 1.0 Iris-versicolor
+6.2 2.2 4.5 1.5 Iris-versicolor
+5.6 2.5 3.9 1.1 Iris-versicolor
+5.9 3.2 4.8 1.8 Iris-versicolor
+6.1 2.8 4.0 1.3 Iris-versicolor
+6.3 2.5 4.9 1.5 Iris-versicolor
+6.1 2.8 4.7 1.2 Iris-versicolor
+6.4 2.9 4.3 1.3 Iris-versicolor
+6.6 3.0 4.4 1.4 Iris-versicolor
+6.8 2.8 4.8 1.4 Iris-versicolor
+6.7 3.0 5.0 1.7 Iris-versicolor
+6.0 2.9 4.5 1.5 Iris-versicolor
+5.7 2.6 3.5 1.0 Iris-versicolor
+5.5 2.4 3.8 1.1 Iris-versicolor
+5.5 2.4 3.7 1.0 Iris-versicolor
+5.8 2.7 3.9 1.2 Iris-versicolor
+6.0 2.7 5.1 1.6 Iris-versicolor
+5.4 3.0 4.5 1.5 Iris-versicolor
+6.0 3.4 4.5 1.6 Iris-versicolor
+6.7 3.1 4.7 1.5 Iris-versicolor
+6.3 2.3 4.4 1.3 Iris-versicolor
+5.6 3.0 4.1 1.3 Iris-versicolor
+5.5 2.5 4.0 1.3 Iris-versicolor
+5.5 2.6 4.4 1.2 Iris-versicolor
+6.1 3.0 4.6 1.4 Iris-versicolor
+5.8 2.6 4.0 1.2 Iris-versicolor
+5.0 2.3 3.3 1.0 Iris-versicolor
+5.6 2.7 4.2 1.3 Iris-versicolor
+5.7 3.0 4.2 1.2 Iris-versicolor
+5.7 2.9 4.2 1.3 Iris-versicolor
+6.2 2.9 4.3 1.3 Iris-versicolor
+5.1 2.5 3.0 1.1 Iris-versicolor
+5.7 2.8 4.1 1.3 Iris-versicolor
+6.3 3.3 6.0 2.5 Iris-virginica
+5.8 2.7 5.1 1.9 Iris-virginica
+7.1 3.0 5.9 2.1 Iris-virginica
+6.3 2.9 5.6 1.8 Iris-virginica
+6.5 3.0 5.8 2.2 Iris-virginica
+7.6 3.0 6.6 2.1 Iris-virginica
+4.9 2.5 4.5 1.7 Iris-virginica
+7.3 2.9 6.3 1.8 Iris-virginica
+6.7 2.5 5.8 1.8 Iris-virginica
+7.2 3.6 6.1 2.5 Iris-virginica
+6.5 3.2 5.1 2.0 Iris-virginica
+6.4 2.7 5.3 1.9 Iris-virginica
+6.8 3.0 5.5 2.1 Iris-virginica
+5.7 2.5 5.0 2.0 Iris-virginica
+5.8 2.8 5.1 2.4 Iris-virginica
+6.4 3.2 5.3 2.3 Iris-virginica
+6.5 3.0 5.5 1.8 Iris-virginica
+7.7 3.8 6.7 2.2 Iris-virginica
+7.7 2.6 6.9 2.3 Iris-virginica
+6.0 2.2 5.0 1.5 Iris-virginica
+6.9 3.2 5.7 2.3 Iris-virginica
+5.6 2.8 4.9 2.0 Iris-virginica
+7.7 2.8 6.7 2.0 Iris-virginica
+6.3 2.7 4.9 1.8 Iris-virginica
+6.7 3.3 5.7 2.1 Iris-virginica
+7.2 3.2 6.0 1.8 Iris-virginica
+6.2 2.8 4.8 1.8 Iris-virginica
+6.1 3.0 4.9 1.8 Iris-virginica
+6.4 2.8 5.6 2.1 Iris-virginica
+7.2 3.0 5.8 1.6 Iris-virginica
+7.4 2.8 6.1 1.9 Iris-virginica
+7.9 3.8 6.4 2.0 Iris-virginica
+6.4 2.8 5.6 2.2 Iris-virginica
+6.3 2.8 5.1 1.5 Iris-virginica
+6.1 2.6 5.6 1.4 Iris-virginica
+7.7 3.0 6.1 2.3 Iris-virginica
+6.3 3.4 5.6 2.4 Iris-virginica
+6.4 3.1 5.5 1.8 Iris-virginica
+6.0 3.0 4.8 1.8 Iris-virginica
+6.9 3.1 5.4 2.1 Iris-virginica
+6.7 3.1 5.6 2.4 Iris-virginica
+6.9 3.1 5.1 2.3 Iris-virginica
+5.8 2.7 5.1 1.9 Iris-virginica
+6.8 3.2 5.9 2.3 Iris-virginica
+6.7 3.3 5.7 2.5 Iris-virginica
+6.7 3.0 5.2 2.3 Iris-virginica
+6.3 2.5 5.0 1.9 Iris-virginica
+6.5 3.0 5.2 2.0 Iris-virginica
+6.2 3.4 5.4 2.3 Iris-virginica
+5.9 3.0 5.1 1.8 Iris-virginica
+
b
diff -r 000000000000 -r f568051cdf2e test-data/pca_out1.tabular
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/pca_out1.tabular Mon May 19 12:34:43 2014 -0400
b
@@ -0,0 +1,159 @@
+#Component 1 2 3 4
+#Std. deviation 1.706 0.9598 0.3839 0.1436
+#Proportion of variance explained 0.7277 0.2303 0.03684 0.005152
+#Loadings 1 2 3 4
+c1 0.5224 -0.3723 0.721 0.262
+c2 -0.2634 -0.9256 -0.242 -0.1241
+c3 0.5813 -0.02109 -0.1409 -0.8012
+c4 0.5656 -0.06542 -0.6338 0.5235
+#Scores 1 2 3 4
+1 -2.265 -0.5057 0.1219 0.02307
+2 -2.086 0.6554 0.2273 0.1032
+3 -2.368 0.3185 -0.05148 0.02783
+4 -2.304 0.5754 -0.09886 -0.06631
+5 -2.389 -0.6748 -0.02143 -0.0374
+6 -2.071 -1.519 -0.03068 0.004399
+7 -2.446 -0.07456 -0.3422 -0.0381
+8 -2.234 -0.2476 0.08257 -0.02551
+9 -2.342 1.095 -0.1536 -0.02679
+10 -2.189 0.4486 0.2466 -0.03991
+11 -2.163 -1.071 0.264 0.0153
+12 -2.327 -0.1586 -0.1002 -0.1346
+13 -2.224 0.7091 0.2232 0.002631
+14 -2.64 0.9383 -0.1896 -0.01942
+15 -2.192 -1.89 0.4695 0.1928
+16 -2.251 -2.722 -0.0326 0.04713
+17 -2.203 -1.514 0.001363 0.1866
+18 -2.19 -0.5143 0.03862 0.09191
+19 -1.894 -1.431 0.3707 0.05953
+20 -2.34 -1.158 -0.1374 -0.03983
+21 -1.915 -0.4305 0.416 0.01036
+22 -2.205 -0.9525 -0.1647 0.05773
+23 -2.774 -0.4895 -0.3388 0.01785
+24 -1.82 -0.1068 -0.04006 0.1503
+25 -2.228 -0.1622 -0.1242 -0.2712
+26 -1.957 0.6079 0.2986 0.04384
+27 -2.052 -0.266 -0.09209 0.0666
+28 -2.168 -0.552 0.2013 0.009261
+29 -2.14 -0.3366 0.2653 0.08354
+30 -2.269 0.3149 -0.07552 -0.1088
+31 -2.145 0.4839 0.06786 -0.04838
+32 -1.832 -0.4453 0.2654 0.2391
+33 -2.608 -1.828 -0.05142 -0.2319
+34 -2.438 -2.185 0.07935 -0.05102
+35 -2.189 0.4486 0.2466 -0.03991
+36 -2.211 0.1843 0.2186 0.1686
+37 -2.044 -0.685 0.4794 0.1956
+38 -2.189 0.4486 0.2466 -0.03991
+39 -2.436 0.8822 -0.2016 -0.009961
+40 -2.171 -0.2927 0.1699 0.00624
+41 -2.287 -0.468 -0.04074 0.1057
+42 -1.872 2.328 0.1945 0.2917
+43 -2.558 0.4538 -0.3136 -0.06741
+44 -1.964 -0.4974 -0.3148 0.1755
+45 -2.133 -1.171 -0.2528 -0.1532
+46 -2.075 0.6919 0.05656 0.1403
+47 -2.381 -1.151 -0.0621 -0.1542
+48 -2.398 0.3624 -0.1469 -0.04948
+49 -2.227 -1.025 0.1766 -0.01644
+50 -2.206 -0.03224 0.1466 0.04878
+51 1.104 -0.8631 0.6856 0.03498
+52 0.7325 -0.5986 0.09407 0.004454
+53 1.242 -0.6148 0.5548 0.009673
+54 0.3973 1.758 0.01857 0.0674
+55 1.073 0.2118 0.3974 0.1055
+56 0.3845 0.5911 -0.1268 -0.2405
+57 0.7487 -0.7787 -0.1487 -0.0783
+58 -0.4979 1.849 -0.2556 -0.03939
+59 0.9262 -0.03033 0.5955 -0.0291
+60 0.004968 1.029 -0.5429 -0.02835
+61 -0.1247 2.658 0.03981 0.01614
+62 0.4387 0.05888 -0.2067 0.03985
+63 0.5516 1.773 0.7614 0.04835
+64 0.7172 0.1854 0.0673 -0.1646
+65 -0.03726 0.4328 -0.1981 0.109
+66 0.8759 -0.51 0.5035 0.1051
+67 0.348 0.1906 -0.4928 -0.1921
+68 0.1534 0.7907 0.2986 -0.2043
+69 1.215 1.633 0.4794 0.2282
+70 0.1569 1.303 0.1686 -0.05041
+71 0.7383 -0.4025 -0.6168 -0.08445
+72 0.4724 0.4166 0.2627 0.1142
+73 1.228 0.9409 0.3667 -0.008448
+74 0.6294 0.4168 0.29 -0.2735
+75 0.7005 0.06349 0.4448 0.04408
+76 0.8735 -0.2507 0.4721 0.1021
+77 1.254 0.08262 0.7268 0.04083
+78 1.358 -0.3288 0.2615 0.06701
+79 0.6621 0.2243 -0.08737 -0.03635
+80 -0.04728 1.057 0.3153 0.06601
+81 0.1215 1.564 0.1452 -0.007875
+82 0.01412 1.573 0.2366 -0.03115
+83 0.236 0.7759 0.148 0.02446
+84 1.057 0.6369 -0.1068 -0.1834
+85 0.2214 0.2808 -0.6676 -0.2556
+86 0.4318 -0.8551 -0.4507 -0.1111
+87 1.049 -0.5222 0.3961 0.0373
+88 1.036 1.392 0.6854 0.1391
+89 0.06707 0.2126 -0.2941 -0.1475
+90 0.2754 1.33 -0.09344 0.009948
+91 0.2723 1.119 -0.09817 -0.2698
+92 0.6232 -0.02754 0.0193 -0.1477
+93 0.33 0.9889 0.196 0.007627
+94 -0.3736 2.018 -0.1122 0.02108
+95 0.2829 0.854 -0.1341 -0.1069
+96 0.08905 0.1749 -0.1314 -0.2301
+97 0.2244 0.3805 -0.1588 -0.1326
+98 0.5739 0.1537 0.27 -0.01941
+99 -0.457 1.539 -0.1961 0.2009
+100 0.2522 0.5959 -0.09475 -0.0583
+101 1.848 -0.8717 -1.003 -0.05107
+102 1.153 0.7013 -0.5315 -0.04041
+103 2.206 -0.5545 0.2055 0.0593
+104 1.439 0.05001 -0.1634 -0.2358
+105 1.868 -0.2912 -0.394 -0.01678
+106 2.754 -0.7884 0.5862 -0.1009
+107 0.3584 1.56 -0.991 -0.133
+108 2.303 -0.4095 0.6542 -0.2372
+109 2.002 0.7239 0.3941 -0.08499
+110 2.268 -1.921 -0.3925 0.1029
+111 1.366 -0.6939 -0.2833 0.107
+112 1.599 0.4282 -0.0233 0.05894
+113 1.884 -0.4143 -0.02455 0.1463
+114 1.253 1.167 -0.5821 0.09968
+115 1.464 0.4441 -1.004 0.275
+116 1.592 -0.677 -0.6367 0.1906
+117 1.471 -0.2532 -0.03666 -0.1554
+118 2.437 -2.557 0.1342 -0.2757
+119 3.309 0.002361 0.7069 0.04676
+120 1.254 1.718 0.2646 -0.06307
+121 2.04 -0.9074 -0.2319 0.1671
+122 0.9739 0.5712 -0.8295 0.02732
+123 2.898 -0.3978 0.8609 -0.1261
+124 1.329 0.4868 0.004707 0.1406
+125 1.704 -1.014 -0.296 -0.06274
+126 1.958 -1.003 0.4228 -0.2185
+127 1.172 0.3189 -0.1307 0.1257
+128 1.02 -0.06554 -0.338 -0.009069
+129 1.786 0.1933 -0.27 0.03121
+130 1.865 -0.5554 0.7175 -0.2076
+131 2.435 -0.2467 0.7302 -0.01679
+132 2.316 -2.626 0.4996 -0.2132
+133 1.86 0.1847 -0.3533 0.1
+134 1.111 0.296 0.1827 -0.1857
+135 1.197 0.8172 0.1632 -0.4884
+136 2.801 -0.8447 0.547 0.2963
+137 1.58 -1.072 -0.9434 0.03361
+138 1.347 -0.4223 -0.18 -0.2159
+139 0.9234 -0.01923 -0.4174 0.004744
+140 1.854 -0.6724 0.01482 0.1949
+141 2.016 -0.6104 -0.4259 0.2468
+142 1.903 -0.686 -0.1278 0.4692
+143 1.153 0.7013 -0.5315 -0.04041
+144 2.043 -0.8647 -0.3353 0.04428
+145 2.002 -1.049 -0.6293 0.2126
+146 1.871 -0.3828 -0.2545 0.3889
+147 1.558 0.9053 0.02538 0.2213
+148 1.521 -0.2668 -0.1793 0.1189
+149 1.376 -1.016 -0.9314 0.02415
+150 0.9593 0.02228 -0.5288 -0.1637
b
diff -r 000000000000 -r f568051cdf2e test-data/pca_out2.pdf
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/pca_out2.pdf Mon May 19 12:34:43 2014 -0400
[
b'@@ -0,0 +1,705 @@\n+%PDF-1.4\n+%\x81\xe2\x81\xe3\x81\xcf\x81\xd3\\r\n+1 0 obj\n+<<\n+/CreationDate (D:20110308103314)\n+/ModDate (D:20110308103314)\n+/Title (R Graphics Output)\n+/Producer (R 2.11.0)\n+/Creator (R)\n+>>\n+endobj\n+2 0 obj\n+<<\n+/Type /Catalog\n+/Pages 3 0 R\n+>>\n+endobj\n+5 0 obj\n+<<\n+/Type /Page\n+/Parent 3 0 R\n+/Contents 6 0 R\n+/Resources 4 0 R\n+>>\n+endobj\n+6 0 obj\n+<<\n+/Length 7 0 R\n+>>\n+stream\r\n+1 J 1 j q\n+Q q\n+0.000 0.000 0.000 RG\n+0.75 w\n+[] 0 d\n+1 J\n+1 j\n+10.00 M\n+125.41 73.44 m 484.30 73.44 l S\n+125.41 73.44 m 125.41 66.24 l S\n+215.13 73.44 m 215.13 66.24 l S\n+304.85 73.44 m 304.85 66.24 l S\n+394.58 73.44 m 394.58 66.24 l S\n+484.30 73.44 m 484.30 66.24 l S\n+BT\n+0.000 0.000 0.000 rg\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 113.57 47.52 Tm (-0.2) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 203.29 47.52 Tm (-0.1) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 296.51 47.52 Tm (0.0) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 386.24 47.52 Tm (0.1) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 475.96 47.52 Tm (0.2) Tj\n+ET\n+80.64 118.21 m 80.64 477.10 l S\n+80.64 118.21 m 73.44 118.21 l S\n+80.64 207.93 m 73.44 207.93 l S\n+80.64 297.65 m 73.44 297.65 l S\n+80.64 387.38 m 73.44 387.38 l S\n+80.64 477.10 m 73.44 477.10 l S\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 63.36 106.37 Tm (-0.2) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 63.36 196.09 Tm (-0.1) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 63.36 289.31 Tm (0.0) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 63.36 379.04 Tm (0.1) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 63.36 468.76 Tm (0.2) Tj\n+ET\n+80.64 73.44 m\n+524.16 73.44 l\n+524.16 516.96 l\n+80.64 516.96 l\n+80.64 73.44 l\n+S\n+Q q\n+BT\n+0.000 0.000 0.000 rg\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 281.60 18.72 Tm [(Comp) 35 (.1)] TJ\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 34.56 274.40 Tm [(Comp) 35 (.2)] TJ\n+ET\n+Q q 80.64 73.44 443.52 443.52 re W n\n+BT\n+0.000 0.000 0.000 rg\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 204.28 254.84 Tm (1) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 211.93 343.46 Tm (2) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 199.84 317.86 Tm (3) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 202.58 337.35 Tm (4) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 198.95 242.14 Tm (5) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 212.61 177.65 Tm (6) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 196.50 287.84 Tm (7) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 205.60 274.65 Tm (8) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 200.96 377.14 Tm (9) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 204.20 327.79 Tm (10) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 205.29 211.72 Tm (11) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 198.25 281.33 Tm (12) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 202.68 347.67 Tm (13) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 184.84 365.05 Tm (14) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 204.05 149.30 Tm (15) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 201.51 85.76 Tm (16) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 203.60 177.90 Tm (17) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 204.14 254.30 Tm (18) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 216.85 184.32 Tm (19) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 197.71 205.16 Tm (20) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 215.97 260.58 Tm (21) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 203.52 220.74 Tm (22) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 179.06 256.19 Tm (23) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 220.02 285.29 Tm (24) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 202.51 281.17 Tm (25) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 214.15 339.95 Tm (26) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 210.07 273.13 Tm (27) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 205.08 251.42 Tm (28) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 206.28 267.86 Tm (29) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 200.76 317.58 Tm (30) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 206.10 330.49 Tm (31) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 219.52 259.56 Tm (32) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 186.1'..b'.00 -0.00 12.00 346.26 355.92 Tm (135) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 415.11 229.07 Tm (136) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 362.70 211.69 Tm (137) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 352.69 261.32 Tm (138) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 334.50 292.08 Tm (139) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 374.43 242.23 Tm (140) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 381.42 246.85 Tm (141) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 376.56 241.07 Tm (142) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 344.36 347.08 Tm (143) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 382.58 227.44 Tm (144) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 380.80 213.52 Tm (145) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 375.16 264.33 Tm (146) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 361.77 362.54 Tm (147) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 360.15 273.19 Tm (148) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 353.95 215.98 Tm (149) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 336.04 295.25 Tm (150) Tj\n+ET\n+Q q\n+Q q\n+1.000 0.000 0.000 RG\n+0.75 w\n+[] 0 d\n+1 J\n+1 j\n+10.00 M\n+137.82 516.96 m 471.89 516.96 l S\n+137.82 516.96 m 137.82 524.16 l S\n+221.33 516.96 m 221.33 524.16 l S\n+304.85 516.96 m 304.85 524.16 l S\n+388.37 516.96 m 388.37 524.16 l S\n+471.89 516.96 m 471.89 524.16 l S\n+BT\n+0.000 0.000 0.000 rg\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 127.64 534.24 Tm (-10) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 214.49 534.24 Tm (-5) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 301.52 534.24 Tm (0) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 385.04 534.24 Tm (5) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 465.22 534.24 Tm (10) Tj\n+ET\n+524.16 130.62 m 524.16 464.69 l S\n+524.16 130.62 m 531.36 130.62 l S\n+524.16 214.13 m 531.36 214.13 l S\n+524.16 297.65 m 531.36 297.65 l S\n+524.16 381.17 m 531.36 381.17 l S\n+524.16 464.69 m 531.36 464.69 l S\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 550.08 120.44 Tm (-10) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 550.08 207.29 Tm (-5) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 550.08 294.32 Tm (0) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 550.08 377.84 Tm (5) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 550.08 458.02 Tm (10) Tj\n+ET\n+0.000 0.000 0.000 RG\n+80.64 73.44 m\n+524.16 73.44 l\n+524.16 516.96 l\n+80.64 516.96 l\n+80.64 73.44 l\n+S\n+Q q 80.64 73.44 443.52 443.52 re W n\n+BT\n+1.000 0.000 0.000 rg\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 473.26 220.33 Tm [(V) 70 (ar 1)] TJ\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 199.01 111.70 Tm [(V) 70 (ar 2)] TJ\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 493.81 289.32 Tm [(V) 70 (ar 3)] TJ\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 488.35 280.59 Tm [(V) 70 (ar 4)] TJ\n+ET\n+1.000 0.000 0.000 RG\n+0.75 w\n+[] 0 d\n+1 J\n+1 j\n+10.00 M\n+304.85 297.65 m 450.72 239.17 l S\n+443.59 238.15 m\n+450.72 239.17 l\n+446.27 244.83 l\n+S\n+304.85 297.65 m 231.32 152.26 l S\n+230.92 159.45 m\n+231.32 152.26 l\n+237.34 156.20 l\n+S\n+304.85 297.65 m 467.16 294.34 l S\n+460.85 290.87 m\n+467.16 294.34 l\n+461.00 298.07 l\n+S\n+304.85 297.65 m 462.79 287.38 l S\n+456.33 284.19 m\n+462.79 287.38 l\n+456.80 291.38 l\n+S\n+Q\n+endstream\n+endobj\n+7 0 obj\n+12856\n+endobj\n+3 0 obj\n+<<\n+/Type /Pages\n+/Kids [\n+5 0 R\n+]\n+/Count 1\n+/MediaBox [0 0 576 576]\n+>>\n+endobj\n+4 0 obj\n+<<\n+/ProcSet [/PDF /Text]\n+/Font <</F2 9 0 R >>\n+/ExtGState << >>\n+>>\n+endobj\n+8 0 obj\n+<<\n+/Type /Encoding\n+/BaseEncoding /WinAnsiEncoding\n+/Differences [ 45/minus 96/quoteleft\n+144/dotlessi /grave /acute /circumflex /tilde /macron /breve /dotaccent\n+/dieresis /.notdef /ring /cedilla /.notdef /hungarumlaut /ogonek /caron /space]\n+>>\n+endobj\n+9 0 obj <<\n+/Type /Font\n+/Subtype /Type1\n+/Name /F2\n+/BaseFont /Helvetica\n+/Encoding 8 0 R\n+>> endobj\n+xref\n+0 10\n+0000000000 65535 f \n+0000000021 00000 n \n+0000000164 00000 n \n+0000013223 00000 n \n+0000013306 00000 n \n+0000000213 00000 n \n+0000000293 00000 n \n+0000013202 00000 n \n+0000013387 00000 n \n+0000013644 00000 n \n+trailer\n+<<\n+/Size 10\n+/Info 1 0 R\n+/Root 2 0 R\n+>>\n+startxref\n+13740\n+%%EOF\n'
b
diff -r 000000000000 -r f568051cdf2e test-data/pca_out3.tabular
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/pca_out3.tabular Mon May 19 12:34:43 2014 -0400
b
@@ -0,0 +1,159 @@
+#Component 1 2 3 4
+#Std. deviation 2.049 0.4905 0.2793 0.1534
+#Proportion of variance explained 0.9246 0.05302 0.01719 0.005183
+#Loadings 1 2 3 4
+c1 0.3616 -0.6565 0.581 0.3173
+c2 -0.08227 -0.7297 -0.5964 -0.3241
+c3 0.8566 0.1758 -0.07252 -0.4797
+c4 0.3588 0.07471 -0.5491 0.7511
+#Scores 1 2 3 4
+1 -2.684 -0.3266 0.02151 0.001006
+2 -2.715 0.1696 0.2035 0.0996
+3 -2.89 0.1373 -0.02471 0.0193
+4 -2.746 0.3111 -0.03767 -0.07596
+5 -2.729 -0.3339 -0.09623 -0.06313
+6 -2.28 -0.7478 -0.1743 -0.02715
+7 -2.821 0.0821 -0.2643 -0.0501
+8 -2.626 -0.1704 0.0158 -0.04628
+9 -2.888 0.5708 -0.02734 -0.02662
+10 -2.674 0.1067 0.1915 -0.05589
+11 -2.507 -0.6519 0.06927 -0.01661
+12 -2.613 -0.02152 -0.1077 -0.1577
+13 -2.787 0.2277 0.2003 -0.007235
+14 -3.225 0.5033 -0.06841 -0.02195
+15 -2.644 -1.186 0.1445 0.157
+16 -2.384 -1.345 -0.2837 0.001926
+17 -2.623 -0.8181 -0.1453 0.1647
+18 -2.648 -0.3191 -0.03339 0.07612
+19 -2.199 -0.8792 0.1145 0.02533
+20 -2.587 -0.5205 -0.2196 -0.06908
+21 -2.311 -0.3979 0.2337 -0.01532
+22 -2.543 -0.44 -0.2148 0.03844
+23 -3.216 -0.1416 -0.2996 0.001857
+24 -2.303 -0.1055 -0.04568 0.1472
+25 -2.356 0.03121 -0.1294 -0.3016
+26 -2.508 0.1391 0.2471 0.03538
+27 -2.469 -0.1379 -0.1013 0.05597
+28 -2.562 -0.3747 0.07236 -0.01524
+29 -2.64 -0.3193 0.1393 0.06514
+30 -2.633 0.1901 -0.04647 -0.1246
+31 -2.588 0.1974 0.07128 -0.06048
+32 -2.41 -0.4181 0.1384 0.2308
+33 -2.648 -0.82 -0.2306 -0.2848
+34 -2.597 -1.1 -0.1636 -0.09896
+35 -2.674 0.1067 0.1915 -0.05589
+36 -2.867 -0.07719 0.1568 0.1625
+37 -2.625 -0.6068 0.2612 0.1759
+38 -2.674 0.1067 0.1915 -0.05589
+39 -2.982 0.4803 -0.07972 -0.01105
+40 -2.59 -0.2361 0.0739 -0.01456
+41 -2.77 -0.2711 -0.08424 0.09236
+42 -2.852 0.9329 0.341 0.3227
+43 -2.998 0.3343 -0.199 -0.07587
+44 -2.406 -0.1959 -0.2707 0.1738
+45 -2.209 -0.4427 -0.3035 -0.1859
+46 -2.716 0.2427 0.09052 0.143
+47 -2.538 -0.5104 -0.1719 -0.1922
+48 -2.84 0.2206 -0.09006 -0.06039
+49 -2.543 -0.5863 0.01118 -0.04833
+50 -2.704 -0.115 0.0827 0.0341
+51 1.285 -0.6854 0.4061 0.01929
+52 0.9324 -0.3192 0.01713 -6.758e-06
+53 1.464 -0.5042 0.3383 -0.0008576
+54 0.181 0.8256 0.1771 0.09578
+55 1.087 -0.07539 0.3065 0.1134
+56 0.6404 0.4173 -0.04119 -0.2427
+57 1.095 -0.2839 -0.17 -0.08497
+58 -0.7515 1.001 -0.01567 -0.01651
+59 1.043 -0.229 0.4148 -0.03752
+60 -0.01019 0.7206 -0.2834 -0.005946
+61 -0.5111 1.262 0.2665 0.04891
+62 0.5111 0.1023 -0.1323 0.0501
+63 0.2623 0.5479 0.6919 0.06148
+64 0.984 0.1244 0.06216 -0.169
+65 -0.1749 0.2518 -0.09366 0.1249
+66 0.9276 -0.4682 0.3132 0.1004
+67 0.6596 0.352 -0.3284 -0.189
+68 0.2345 0.3319 0.2703 -0.212
+69 0.9424 0.5418 0.4973 0.2606
+70 0.04325 0.5815 0.233 -0.03956
+71 1.116 0.08421 -0.4598 -0.07721
+72 0.3568 0.06682 0.2275 0.1241
+73 1.296 0.3276 0.3475 0.003246
+74 0.9205 0.1824 0.2316 -0.2868
+75 0.714 -0.1504 0.3204 0.04294
+76 0.8996 -0.3296 0.3148 0.1011
+77 1.331 -0.2447 0.5212 0.03751
+78 1.557 -0.2674 0.1646 0.07035
+79 0.8125 0.1623 -0.03634 -0.02968
+80 -0.3073 0.3651 0.3153 0.07653
+81 -0.07034 0.7025 0.2418 0.009095
+82 -0.1919 0.6775 0.3039 -0.01805
+83 0.135 0.3117 0.175 0.03418
+84 1.379 0.4212 -0.01548 -0.1776
+85 0.5873 0.4833 -0.4446 -0.2524
+86 0.8072 -0.1951 -0.3895 -0.1166
+87 1.22 -0.408 0.2366 0.03164
+88 0.8129 0.3707 0.6129 0.1577
+89 0.2452 0.2667 -0.1896 -0.1473
+90 0.1645 0.6797 0.0578 0.03097
+91 0.463 0.6695 0.02405 -0.2684
+92 0.8902 0.03381 0.009768 -0.1534
+93 0.2289 0.4023 0.2274 0.01862
+94 -0.7071 1.008 0.1021 0.04762
+95 0.3555 0.5032 -0.01789 -0.09807
+96 0.3311 0.2112 -0.08381 -0.2387
+97 0.3752 0.2916 -0.07907 -0.1312
+98 0.6417 -0.01907 0.2042 -0.02051
+99 -0.9085 0.7516 0.007737 0.2336
+100 0.2978 0.347 -0.01218 -0.05078
+101 2.532 0.01184 -0.7585 -0.0326
+102 1.414 0.5749 -0.2964 -0.0157
+103 2.616 -0.3419 0.1121 0.06596
+104 1.971 0.1811 -0.1065 -0.2369
+105 2.35 0.04188 -0.2841 -0.001313
+106 3.397 -0.5472 0.3519 -0.1112
+107 0.5194 1.191 -0.5467 -0.0988
+108 2.932 -0.3524 0.4237 -0.2554
+109 2.32 0.2455 0.3499 -0.07626
+110 2.918 -0.7804 -0.4217 0.1077
+111 1.662 -0.242 -0.2428 0.1194
+112 1.802 0.2162 0.0377 0.07871
+113 2.165 -0.2153 -0.03315 0.1627
+114 1.345 0.7764 -0.2829 0.1405
+115 1.585 0.5393 -0.6306 0.3275
+116 1.905 -0.1188 -0.4801 0.2171
+117 1.949 -0.04073 -0.04273 -0.1578
+118 3.489 -1.172 -0.1293 -0.3116
+119 3.795 -0.2533 0.517 0.05645
+120 1.298 0.761 0.3449 -0.04267
+121 2.428 -0.3768 -0.2186 0.1839
+122 1.198 0.6056 -0.5126 0.0595
+123 3.499 -0.4568 0.5769 -0.1378
+124 1.388 0.204 0.06351 0.1638
+125 2.276 -0.3334 -0.2847 -0.06223
+126 2.614 -0.5584 0.2084 -0.2404
+127 1.258 0.1791 -0.04698 0.1476
+128 1.291 0.1164 -0.2316 0.003084
+129 2.123 0.2109 -0.1535 0.05261
+130 2.388 -0.4625 0.452 -0.2299
+131 2.841 -0.3727 0.501 -0.02022
+132 3.232 -1.371 0.1184 -0.2545
+133 2.159 0.2183 -0.2084 0.1277
+134 1.443 0.1438 0.1541 -0.1899
+135 1.78 0.5015 0.1758 -0.5035
+136 3.077 -0.6858 0.3364 0.3106
+137 2.145 -0.1389 -0.7342 0.05177
+138 1.905 -0.04805 -0.1605 -0.222
+139 1.169 0.1645 -0.2825 0.01933
+140 2.108 -0.3715 -0.02744 0.21
+141 2.314 -0.1826 -0.3229 0.2759
+142 1.922 -0.4093 -0.1155 0.5041
+143 1.414 0.5749 -0.2964 -0.0157
+144 2.563 -0.276 -0.2913 0.05618
+145 2.419 -0.3035 -0.5043 0.2382
+146 1.944 -0.1874 -0.1793 0.4251
+147 1.526 0.375 0.1206 0.2557
+148 1.764 -0.07852 -0.1308 0.1363
+149 1.902 -0.1159 -0.7229 0.04087
+150 1.39 0.2829 -0.3623 -0.1563
b
diff -r 000000000000 -r f568051cdf2e test-data/pca_out4.pdf
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/pca_out4.pdf Mon May 19 12:34:43 2014 -0400
[
b'@@ -0,0 +1,705 @@\n+%PDF-1.4\n+%\x81\xe2\x81\xe3\x81\xcf\x81\xd3\\r\n+1 0 obj\n+<<\n+/CreationDate (D:20110308103327)\n+/ModDate (D:20110308103327)\n+/Title (R Graphics Output)\n+/Producer (R 2.11.0)\n+/Creator (R)\n+>>\n+endobj\n+2 0 obj\n+<<\n+/Type /Catalog\n+/Pages 3 0 R\n+>>\n+endobj\n+5 0 obj\n+<<\n+/Type /Page\n+/Parent 3 0 R\n+/Contents 6 0 R\n+/Resources 4 0 R\n+>>\n+endobj\n+6 0 obj\n+<<\n+/Length 7 0 R\n+>>\n+stream\r\n+1 J 1 j q\n+Q q\n+0.000 0.000 0.000 RG\n+0.75 w\n+[] 0 d\n+1 J\n+1 j\n+10.00 M\n+123.42 73.44 m 498.23 73.44 l S\n+123.42 73.44 m 123.42 66.24 l S\n+217.12 73.44 m 217.12 66.24 l S\n+310.82 73.44 m 310.82 66.24 l S\n+404.53 73.44 m 404.53 66.24 l S\n+498.23 73.44 m 498.23 66.24 l S\n+BT\n+0.000 0.000 0.000 rg\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 111.57 47.52 Tm (-0.2) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 205.28 47.52 Tm (-0.1) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 302.48 47.52 Tm (0.0) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 396.19 47.52 Tm (0.1) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 489.89 47.52 Tm (0.2) Tj\n+ET\n+80.64 116.22 m 80.64 491.03 l S\n+80.64 116.22 m 73.44 116.22 l S\n+80.64 209.92 m 73.44 209.92 l S\n+80.64 303.62 m 73.44 303.62 l S\n+80.64 397.33 m 73.44 397.33 l S\n+80.64 491.03 m 73.44 491.03 l S\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 63.36 104.37 Tm (-0.2) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 63.36 198.08 Tm (-0.1) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 63.36 295.28 Tm (0.0) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 63.36 388.99 Tm (0.1) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 63.36 482.69 Tm (0.2) Tj\n+ET\n+80.64 73.44 m\n+524.16 73.44 l\n+524.16 516.96 l\n+80.64 516.96 l\n+80.64 73.44 l\n+S\n+Q q\n+BT\n+0.000 0.000 0.000 rg\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 281.60 18.72 Tm [(Comp) 35 (.1)] TJ\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 34.56 274.40 Tm [(Comp) 35 (.2)] TJ\n+ET\n+Q q 80.64 73.44 443.52 443.52 re W n\n+BT\n+0.000 0.000 0.000 rg\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 207.24 248.47 Tm (1) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 206.08 325.85 Tm (2) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 199.56 320.94 Tm (3) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 204.92 347.93 Tm (4) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 205.58 247.53 Tm (5) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 222.34 182.89 Tm (6) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 202.14 312.30 Tm (7) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 209.40 272.94 Tm (8) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 199.63 388.55 Tm (9) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 204.29 316.16 Tm (10) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 210.54 197.73 Tm (11) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 206.56 296.05 Tm (12) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 200.05 335.04 Tm (13) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 183.70 377.90 Tm (14) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 205.42 114.51 Tm (15) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 215.12 89.78 Tm (16) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 206.21 171.81 Tm (17) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 205.25 249.75 Tm (18) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 222.02 162.39 Tm (19) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 207.52 218.34 Tm (20) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 217.86 237.35 Tm (21) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 209.17 230.78 Tm (22) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 184.05 277.43 Tm (23) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 218.14 282.95 Tm (24) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 216.16 304.39 Tm (25) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 210.49 321.21 Tm (26) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 211.94 277.90 Tm (27) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 208.45 241.08 Tm (28) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 205.56 249.72 Tm (29) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 205.82 329.17 Tm (30) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 207.48 330.31 Tm (31) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 214.14 234.31 Tm (32) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 205.2'..b'-0.00 12.00 367.28 377.73 Tm (135) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 415.72 192.56 Tm (136) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 380.93 277.86 Tm (137) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 371.96 292.03 Tm (138) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 344.47 325.18 Tm (139) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 379.53 241.58 Tm (140) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 387.25 270.93 Tm (141) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 372.61 235.57 Tm (142) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 353.63 389.19 Tm (143) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 396.55 256.36 Tm (144) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 391.17 252.18 Tm (145) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 373.42 270.29 Tm (146) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 357.80 357.90 Tm (147) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 366.70 287.27 Tm (148) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 371.84 281.45 Tm (149) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 352.72 343.64 Tm (150) Tj\n+ET\n+Q q\n+Q q\n+1.000 0.000 0.000 RG\n+0.75 w\n+[] 0 d\n+1 J\n+1 j\n+10.00 M\n+127.58 516.96 m 494.07 516.96 l S\n+127.58 516.96 m 127.58 524.16 l S\n+219.20 516.96 m 219.20 524.16 l S\n+310.82 516.96 m 310.82 524.16 l S\n+402.45 516.96 m 402.45 524.16 l S\n+494.07 516.96 m 494.07 524.16 l S\n+BT\n+0.000 0.000 0.000 rg\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 117.40 534.24 Tm (-20) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 209.03 534.24 Tm (-10) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 307.49 534.24 Tm (0) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 395.78 534.24 Tm (10) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 487.40 534.24 Tm (20) Tj\n+ET\n+524.16 120.38 m 524.16 486.87 l S\n+524.16 120.38 m 531.36 120.38 l S\n+524.16 212.00 m 531.36 212.00 l S\n+524.16 303.62 m 531.36 303.62 l S\n+524.16 395.25 m 531.36 395.25 l S\n+524.16 486.87 m 531.36 486.87 l S\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 550.08 110.20 Tm (-20) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 550.08 201.83 Tm (-10) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 550.08 300.29 Tm (0) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 550.08 388.58 Tm (10) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 550.08 480.20 Tm (20) Tj\n+ET\n+0.000 0.000 0.000 RG\n+80.64 73.44 m\n+524.16 73.44 l\n+524.16 516.96 l\n+80.64 516.96 l\n+80.64 73.44 l\n+S\n+Q q 80.64 73.44 443.52 443.52 re W n\n+BT\n+1.000 0.000 0.000 rg\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 380.03 263.27 Tm [(V) 70 (ar 1)] TJ\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 277.99 259.24 Tm [(V) 70 (ar 2)] TJ\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 493.81 309.11 Tm [(V) 70 (ar 3)] TJ\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 379.40 303.52 Tm [(V) 70 (ar 4)] TJ\n+ET\n+1.000 0.000 0.000 RG\n+0.75 w\n+[] 0 d\n+1 J\n+1 j\n+10.00 M\n+310.82 303.62 m 377.32 274.71 l S\n+370.17 273.90 m\n+377.32 274.71 l\n+373.04 280.50 l\n+S\n+310.82 303.62 m 295.70 271.49 l S\n+295.09 278.67 m\n+295.70 271.49 l\n+301.61 275.60 l\n+S\n+310.82 303.62 m 468.35 311.36 l S\n+462.30 307.46 m\n+468.35 311.36 l\n+461.95 314.65 l\n+S\n+310.82 303.62 m 376.82 306.91 l S\n+370.77 303.01 m\n+376.82 306.91 l\n+370.41 310.20 l\n+S\n+Q\n+endstream\n+endobj\n+7 0 obj\n+12861\n+endobj\n+3 0 obj\n+<<\n+/Type /Pages\n+/Kids [\n+5 0 R\n+]\n+/Count 1\n+/MediaBox [0 0 576 576]\n+>>\n+endobj\n+4 0 obj\n+<<\n+/ProcSet [/PDF /Text]\n+/Font <</F2 9 0 R >>\n+/ExtGState << >>\n+>>\n+endobj\n+8 0 obj\n+<<\n+/Type /Encoding\n+/BaseEncoding /WinAnsiEncoding\n+/Differences [ 45/minus 96/quoteleft\n+144/dotlessi /grave /acute /circumflex /tilde /macron /breve /dotaccent\n+/dieresis /.notdef /ring /cedilla /.notdef /hungarumlaut /ogonek /caron /space]\n+>>\n+endobj\n+9 0 obj <<\n+/Type /Font\n+/Subtype /Type1\n+/Name /F2\n+/BaseFont /Helvetica\n+/Encoding 8 0 R\n+>> endobj\n+xref\n+0 10\n+0000000000 65535 f \n+0000000021 00000 n \n+0000000164 00000 n \n+0000013228 00000 n \n+0000013311 00000 n \n+0000000213 00000 n \n+0000000293 00000 n \n+0000013207 00000 n \n+0000013392 00000 n \n+0000013649 00000 n \n+trailer\n+<<\n+/Size 10\n+/Info 1 0 R\n+/Root 2 0 R\n+>>\n+startxref\n+13745\n+%%EOF\n'
b
diff -r 000000000000 -r f568051cdf2e test-data/pca_out5.tabular
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/pca_out5.tabular Mon May 19 12:34:43 2014 -0400
b
@@ -0,0 +1,159 @@
+#Component 1 2 3 4
+#Std. deviation 1.706 0.9598 0.3839 0.1436
+#Proportion of variance explained 0.7277 0.2303 0.03684 0.005152
+#Loadings 1 2 3 4
+c1 0.5224 -0.3723 0.721 0.262
+c2 -0.2634 -0.9256 -0.242 -0.1241
+c3 0.5813 -0.02109 -0.1409 -0.8012
+c4 0.5656 -0.06542 -0.6338 0.5235
+#Scores 1 2 3 4
+1 -2.257 -0.504 0.1215 0.023
+2 -2.079 0.6532 0.2265 0.1029
+3 -2.36 0.3174 -0.05131 0.02773
+4 -2.297 0.5734 -0.09853 -0.06609
+5 -2.381 -0.6725 -0.02136 -0.03727
+6 -2.064 -1.513 -0.03058 0.004384
+7 -2.438 -0.07431 -0.3411 -0.03797
+8 -2.226 -0.2468 0.0823 -0.02542
+9 -2.334 1.091 -0.153 -0.0267
+10 -2.181 0.4471 0.2457 -0.03977
+11 -2.156 -1.067 0.2631 0.01525
+12 -2.32 -0.1581 -0.09983 -0.1341
+13 -2.217 0.7068 0.2225 0.002622
+14 -2.631 0.9351 -0.1889 -0.01936
+15 -2.185 -1.884 0.4679 0.1921
+16 -2.244 -2.713 -0.03249 0.04697
+17 -2.195 -1.509 0.001359 0.186
+18 -2.183 -0.5126 0.03849 0.0916
+19 -1.888 -1.426 0.3695 0.05933
+20 -2.332 -1.154 -0.137 -0.0397
+21 -1.908 -0.429 0.4146 0.01033
+22 -2.197 -0.9493 -0.1642 0.05754
+23 -2.765 -0.4879 -0.3377 0.01779
+24 -1.814 -0.1064 -0.03993 0.1498
+25 -2.221 -0.1616 -0.1238 -0.2703
+26 -1.95 0.6059 0.2976 0.04369
+27 -2.045 -0.2651 -0.09179 0.06638
+28 -2.161 -0.5502 0.2006 0.00923
+29 -2.133 -0.3355 0.2644 0.08326
+30 -2.261 0.3138 -0.07526 -0.1085
+31 -2.137 0.4823 0.06763 -0.04822
+32 -1.826 -0.4438 0.2645 0.2383
+33 -2.599 -1.822 -0.05125 -0.2311
+34 -2.43 -2.178 0.07908 -0.05085
+35 -2.181 0.4471 0.2457 -0.03977
+36 -2.204 0.1837 0.2179 0.1681
+37 -2.038 -0.6827 0.4778 0.195
+38 -2.181 0.4471 0.2457 -0.03977
+39 -2.428 0.8792 -0.2009 -0.009928
+40 -2.163 -0.2917 0.1694 0.006219
+41 -2.279 -0.4664 -0.0406 0.1054
+42 -1.865 2.32 0.1939 0.2907
+43 -2.549 0.4523 -0.3125 -0.06719
+44 -1.958 -0.4957 -0.3137 0.175
+45 -2.126 -1.168 -0.2519 -0.1527
+46 -2.068 0.6896 0.05637 0.1398
+47 -2.373 -1.147 -0.06189 -0.1537
+48 -2.39 0.3612 -0.1464 -0.04931
+49 -2.219 -1.022 0.1761 -0.01639
+50 -2.199 -0.03213 0.1461 0.04862
+51 1.1 -0.8602 0.6833 0.03486
+52 0.73 -0.5966 0.09375 0.004439
+53 1.238 -0.6128 0.553 0.009641
+54 0.396 1.752 0.01851 0.06717
+55 1.069 0.2111 0.3961 0.1052
+56 0.3832 0.5891 -0.1264 -0.2397
+57 0.7462 -0.7761 -0.1482 -0.07804
+58 -0.4962 1.843 -0.2547 -0.03926
+59 0.9231 -0.03023 0.5935 -0.02901
+60 0.004951 1.026 -0.5411 -0.02825
+61 -0.1243 2.649 0.03968 0.01609
+62 0.4373 0.05868 -0.206 0.03972
+63 0.5498 1.767 0.7588 0.04819
+64 0.7148 0.1848 0.06708 -0.164
+65 -0.03713 0.4314 -0.1974 0.1087
+66 0.873 -0.5083 0.5018 0.1048
+67 0.3468 0.19 -0.4912 -0.1914
+68 0.1529 0.7881 0.2976 -0.2036
+69 1.211 1.628 0.4778 0.2275
+70 0.1564 1.299 0.168 -0.05024
+71 0.7358 -0.4011 -0.6147 -0.08417
+72 0.4708 0.4152 0.2618 0.1139
+73 1.224 0.9378 0.3655 -0.00842
+74 0.6273 0.4154 0.289 -0.2726
+75 0.6981 0.06328 0.4433 0.04393
+76 0.8706 -0.2499 0.4706 0.1018
+77 1.25 0.08234 0.7244 0.04069
+78 1.354 -0.3277 0.2606 0.06679
+79 0.6599 0.2236 -0.08708 -0.03623
+80 -0.04712 1.054 0.3143 0.06579
+81 0.1211 1.558 0.1448 -0.007849
+82 0.01407 1.568 0.2358 -0.03105
+83 0.2352 0.7733 0.1475 0.02438
+84 1.053 0.6348 -0.1064 -0.1828
+85 0.2207 0.2799 -0.6653 -0.2547
+86 0.4303 -0.8523 -0.4492 -0.1108
+87 1.046 -0.5205 0.3948 0.03717
+88 1.032 1.388 0.6831 0.1387
+89 0.06684 0.2119 -0.2931 -0.147
+90 0.2745 1.325 -0.09313 0.009914
+91 0.2714 1.116 -0.09784 -0.2689
+92 0.6211 -0.02745 0.01924 -0.1472
+93 0.3289 0.9856 0.1953 0.007601
+94 -0.3724 2.011 -0.1118 0.02101
+95 0.282 0.8511 -0.1337 -0.1065
+96 0.08876 0.1743 -0.131 -0.2294
+97 0.2236 0.3792 -0.1582 -0.1321
+98 0.572 0.1532 0.2691 -0.01934
+99 -0.4555 1.534 -0.1955 0.2002
+100 0.2514 0.5939 -0.09443 -0.0581
+101 1.842 -0.8688 -0.9994 -0.0509
+102 1.149 0.699 -0.5297 -0.04028
+103 2.199 -0.5526 0.2048 0.0591
+104 1.434 0.04984 -0.1628 -0.235
+105 1.862 -0.2902 -0.3927 -0.01673
+106 2.745 -0.7858 0.5843 -0.1005
+107 0.3572 1.555 -0.9877 -0.1325
+108 2.295 -0.4081 0.652 -0.2364
+109 1.995 0.7214 0.3928 -0.08471
+110 2.26 -1.915 -0.3912 0.1026
+111 1.361 -0.6916 -0.2823 0.1067
+112 1.594 0.4268 -0.02323 0.05875
+113 1.878 -0.4129 -0.02447 0.1458
+114 1.249 1.163 -0.5802 0.09935
+115 1.459 0.4427 -1.001 0.2741
+116 1.586 -0.6748 -0.6345 0.19
+117 1.466 -0.2523 -0.03654 -0.1549
+118 2.429 -2.548 0.1338 -0.2747
+119 3.298 0.002353 0.7046 0.0466
+120 1.25 1.712 0.2637 -0.06286
+121 2.034 -0.9044 -0.2311 0.1666
+122 0.9707 0.5693 -0.8267 0.02723
+123 2.888 -0.3965 0.8581 -0.1257
+124 1.325 0.4851 0.004692 0.1401
+125 1.699 -1.011 -0.295 -0.06253
+126 1.951 -1 0.4214 -0.2177
+127 1.168 0.3178 -0.1302 0.1253
+128 1.016 -0.06532 -0.3369 -0.009038
+129 1.78 0.1926 -0.2691 0.0311
+130 1.859 -0.5535 0.7151 -0.2069
+131 2.427 -0.2458 0.7278 -0.01674
+132 2.308 -2.617 0.498 -0.2124
+133 1.854 0.1841 -0.3522 0.09971
+134 1.108 0.295 0.182 -0.1851
+135 1.193 0.8144 0.1627 -0.4868
+136 2.792 -0.8419 0.5452 0.2953
+137 1.575 -1.069 -0.9402 0.0335
+138 1.343 -0.4208 -0.1794 -0.2152
+139 0.9203 -0.01917 -0.416 0.004728
+140 1.847 -0.6702 0.01477 0.1942
+141 2.009 -0.6084 -0.4245 0.2459
+142 1.897 -0.6837 -0.1274 0.4676
+143 1.149 0.699 -0.5297 -0.04028
+144 2.036 -0.8618 -0.3341 0.04413
+145 1.995 -1.045 -0.6272 0.2119
+146 1.864 -0.3815 -0.2537 0.3876
+147 1.553 0.9023 0.0253 0.2206
+148 1.516 -0.2659 -0.1787 0.1185
+149 1.372 -1.013 -0.9283 0.02407
+150 0.9561 0.02221 -0.527 -0.1631
b
diff -r 000000000000 -r f568051cdf2e test-data/pca_out6.pdf
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/pca_out6.pdf Mon May 19 12:34:43 2014 -0400
[
b'@@ -0,0 +1,705 @@\n+%PDF-1.4\n+%\x81\xe2\x81\xe3\x81\xcf\x81\xd3\\r\n+1 0 obj\n+<<\n+/CreationDate (D:20110308103341)\n+/ModDate (D:20110308103341)\n+/Title (R Graphics Output)\n+/Producer (R 2.11.0)\n+/Creator (R)\n+>>\n+endobj\n+2 0 obj\n+<<\n+/Type /Catalog\n+/Pages 3 0 R\n+>>\n+endobj\n+5 0 obj\n+<<\n+/Type /Page\n+/Parent 3 0 R\n+/Contents 6 0 R\n+/Resources 4 0 R\n+>>\n+endobj\n+6 0 obj\n+<<\n+/Length 7 0 R\n+>>\n+stream\r\n+1 J 1 j q\n+Q q\n+0.000 0.000 0.000 RG\n+0.75 w\n+[] 0 d\n+1 J\n+1 j\n+10.00 M\n+124.81 73.44 m 484.90 73.44 l S\n+124.81 73.44 m 124.81 66.24 l S\n+214.83 73.44 m 214.83 66.24 l S\n+304.85 73.44 m 304.85 66.24 l S\n+394.88 73.44 m 394.88 66.24 l S\n+484.90 73.44 m 484.90 66.24 l S\n+BT\n+0.000 0.000 0.000 rg\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 112.96 47.52 Tm (-0.2) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 202.99 47.52 Tm (-0.1) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 296.51 47.52 Tm (0.0) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 386.54 47.52 Tm (0.1) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 476.56 47.52 Tm (0.2) Tj\n+ET\n+80.64 117.61 m 80.64 477.70 l S\n+80.64 117.61 m 73.44 117.61 l S\n+80.64 207.63 m 73.44 207.63 l S\n+80.64 297.65 m 73.44 297.65 l S\n+80.64 387.68 m 73.44 387.68 l S\n+80.64 477.70 m 73.44 477.70 l S\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 63.36 105.76 Tm (-0.2) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 63.36 195.79 Tm (-0.1) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 63.36 289.31 Tm (0.0) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 63.36 379.34 Tm (0.1) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 63.36 469.36 Tm (0.2) Tj\n+ET\n+80.64 73.44 m\n+524.16 73.44 l\n+524.16 516.96 l\n+80.64 516.96 l\n+80.64 73.44 l\n+S\n+Q q\n+BT\n+0.000 0.000 0.000 rg\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 290.73 18.72 Tm (PC1) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 34.56 283.53 Tm (PC2) Tj\n+ET\n+Q q 80.64 73.44 443.52 443.52 re W n\n+BT\n+0.000 0.000 0.000 rg\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 204.28 254.84 Tm (1) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 211.93 343.46 Tm (2) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 199.84 317.86 Tm (3) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 202.58 337.35 Tm (4) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 198.95 242.14 Tm (5) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 212.61 177.65 Tm (6) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 196.50 287.84 Tm (7) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 205.60 274.65 Tm (8) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 200.96 377.14 Tm (9) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 204.20 327.79 Tm (10) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 205.29 211.72 Tm (11) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 198.25 281.33 Tm (12) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 202.68 347.67 Tm (13) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 184.84 365.05 Tm (14) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 204.05 149.30 Tm (15) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 201.51 85.76 Tm (16) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 203.60 177.90 Tm (17) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 204.14 254.30 Tm (18) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 216.85 184.32 Tm (19) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 197.71 205.16 Tm (20) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 215.97 260.58 Tm (21) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 203.52 220.74 Tm (22) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 179.06 256.19 Tm (23) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 220.02 285.29 Tm (24) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 202.51 281.17 Tm (25) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 214.15 339.95 Tm (26) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 210.07 273.13 Tm (27) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 205.08 251.42 Tm (28) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 206.28 267.86 Tm (29) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 200.76 317.58 Tm (30) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 206.10 330.49 Tm (31) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 219.52 259.56 Tm (32) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 186.19 153.99 Tm (33) Tj\n+E'..b'.00 -0.00 12.00 346.26 355.92 Tm (135) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 415.11 229.07 Tm (136) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 362.70 211.69 Tm (137) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 352.69 261.32 Tm (138) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 334.50 292.08 Tm (139) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 374.43 242.23 Tm (140) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 381.42 246.85 Tm (141) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 376.56 241.07 Tm (142) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 344.36 347.08 Tm (143) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 382.58 227.44 Tm (144) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 380.80 213.52 Tm (145) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 375.16 264.33 Tm (146) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 361.77 362.54 Tm (147) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 360.15 273.19 Tm (148) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 353.95 215.98 Tm (149) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 336.04 295.25 Tm (150) Tj\n+ET\n+Q q\n+Q q\n+1.000 0.000 0.000 RG\n+0.75 w\n+[] 0 d\n+1 J\n+1 j\n+10.00 M\n+137.82 516.96 m 471.89 516.96 l S\n+137.82 516.96 m 137.82 524.16 l S\n+221.33 516.96 m 221.33 524.16 l S\n+304.85 516.96 m 304.85 524.16 l S\n+388.37 516.96 m 388.37 524.16 l S\n+471.89 516.96 m 471.89 524.16 l S\n+BT\n+0.000 0.000 0.000 rg\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 127.64 534.24 Tm (-10) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 214.49 534.24 Tm (-5) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 301.52 534.24 Tm (0) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 385.04 534.24 Tm (5) Tj\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 465.22 534.24 Tm (10) Tj\n+ET\n+524.16 130.62 m 524.16 464.69 l S\n+524.16 130.62 m 531.36 130.62 l S\n+524.16 214.13 m 531.36 214.13 l S\n+524.16 297.65 m 531.36 297.65 l S\n+524.16 381.17 m 531.36 381.17 l S\n+524.16 464.69 m 531.36 464.69 l S\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 550.08 120.44 Tm (-10) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 550.08 207.29 Tm (-5) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 550.08 294.32 Tm (0) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 550.08 377.84 Tm (5) Tj\n+ET\n+BT\n+/F2 1 Tf 0.00 12.00 -12.00 0.00 550.08 458.02 Tm (10) Tj\n+ET\n+0.000 0.000 0.000 RG\n+80.64 73.44 m\n+524.16 73.44 l\n+524.16 516.96 l\n+80.64 516.96 l\n+80.64 73.44 l\n+S\n+Q q 80.64 73.44 443.52 443.52 re W n\n+BT\n+1.000 0.000 0.000 rg\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 473.26 220.33 Tm [(V) 70 (ar 1)] TJ\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 199.01 111.70 Tm [(V) 70 (ar 2)] TJ\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 493.81 289.32 Tm [(V) 70 (ar 3)] TJ\n+ET\n+BT\n+/F2 1 Tf 12.00 0.00 -0.00 12.00 488.35 280.59 Tm [(V) 70 (ar 4)] TJ\n+ET\n+1.000 0.000 0.000 RG\n+0.75 w\n+[] 0 d\n+1 J\n+1 j\n+10.00 M\n+304.85 297.65 m 450.72 239.17 l S\n+443.59 238.15 m\n+450.72 239.17 l\n+446.27 244.83 l\n+S\n+304.85 297.65 m 231.32 152.26 l S\n+230.92 159.45 m\n+231.32 152.26 l\n+237.34 156.20 l\n+S\n+304.85 297.65 m 467.16 294.34 l S\n+460.85 290.87 m\n+467.16 294.34 l\n+461.00 298.07 l\n+S\n+304.85 297.65 m 462.79 287.38 l S\n+456.33 284.19 m\n+462.79 287.38 l\n+456.80 291.38 l\n+S\n+Q\n+endstream\n+endobj\n+7 0 obj\n+12834\n+endobj\n+3 0 obj\n+<<\n+/Type /Pages\n+/Kids [\n+5 0 R\n+]\n+/Count 1\n+/MediaBox [0 0 576 576]\n+>>\n+endobj\n+4 0 obj\n+<<\n+/ProcSet [/PDF /Text]\n+/Font <</F2 9 0 R >>\n+/ExtGState << >>\n+>>\n+endobj\n+8 0 obj\n+<<\n+/Type /Encoding\n+/BaseEncoding /WinAnsiEncoding\n+/Differences [ 45/minus 96/quoteleft\n+144/dotlessi /grave /acute /circumflex /tilde /macron /breve /dotaccent\n+/dieresis /.notdef /ring /cedilla /.notdef /hungarumlaut /ogonek /caron /space]\n+>>\n+endobj\n+9 0 obj <<\n+/Type /Font\n+/Subtype /Type1\n+/Name /F2\n+/BaseFont /Helvetica\n+/Encoding 8 0 R\n+>> endobj\n+xref\n+0 10\n+0000000000 65535 f \n+0000000021 00000 n \n+0000000164 00000 n \n+0000013201 00000 n \n+0000013284 00000 n \n+0000000213 00000 n \n+0000000293 00000 n \n+0000013180 00000 n \n+0000013365 00000 n \n+0000013622 00000 n \n+trailer\n+<<\n+/Size 10\n+/Info 1 0 R\n+/Root 2 0 R\n+>>\n+startxref\n+13718\n+%%EOF\n'
b
diff -r 000000000000 -r f568051cdf2e tool_dependencies.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_dependencies.xml Mon May 19 12:34:43 2014 -0400
b
@@ -0,0 +1,12 @@
+<?xml version="1.0"?>
+<tool_dependency>
+  <package name="numpy" version="1.7.1">
+      <repository changeset_revision="0c288abd2a1e" name="package_numpy_1_7" owner="devteam" toolshed="http://toolshed.g2.bx.psu.edu" />
+    </package>
+    <package name="rpy" version="1.0.3">
+      <repository changeset_revision="82170c94ca7c" name="package_rpy_1_0_3" owner="devteam" toolshed="http://toolshed.g2.bx.psu.edu" />
+    </package>
+    <package name="R" version="2.11.0">
+      <repository changeset_revision="5824d2b3bc8b" name="package_r_2_11_0" owner="devteam" toolshed="http://toolshed.g2.bx.psu.edu" />
+    </package>
+</tool_dependency>