Mercurial > repos > devteam > linear_regression

changeset 0:5cbb241afdc5 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Imported from capsule None
author devteam
date Tue, 01 Apr 2014 09:12:59 -0400
parents
children
files linear_regression.py linear_regression.xml test-data/regr_inp.tabular test-data/regr_out.pdf test-data/regr_out.tabular tool_dependencies.xml
diffstat 6 files changed, 1844 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/linear_regression.py	Tue Apr 01 09:12:59 2014 -0400
@@ -0,0 +1,146 @@
+#!/usr/bin/env python
+
+import sys
+from rpy import *
+import numpy
+
+def stop_err(msg):
+    sys.stderr.write(msg)
+    sys.exit()
+
+infile = sys.argv[1]
+y_col = int(sys.argv[2])-1
+x_cols = sys.argv[3].split(',')
+outfile = sys.argv[4]
+outfile2 = sys.argv[5]
+
+print "Predictor columns: %s; Response column: %d" % ( x_cols, y_col+1 )
+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." )
+
+y_vals = []
+x_vals = []
+
+for k, col in enumerate(x_cols):
+    x_cols[k] = int(col)-1
+    x_vals.append([])
+
+NA = 'NA'
+for ind, line in enumerate( file( infile )):
+    if line and not line.startswith( '#' ):
+        try:
+            fields = line.split("\t")
+            try:
+                yval = float(fields[y_col])
+            except:
+                yval = r('NA')
+            y_vals.append(yval)
+            for k, col in enumerate(x_cols):
+                try:
+                    xval = float(fields[col])
+                except:
+                    xval = r('NA')
+                x_vals[k].append(xval)
+        except:
+            pass
+
+x_vals1 = numpy.asarray(x_vals).transpose()
+
+dat = r.list(x=array(x_vals1), y=y_vals)
+
+set_default_mode(NO_CONVERSION)
+try:
+    linear_model = r.lm(r("y ~ x"), data = r.na_exclude(dat))
+except RException, rex:
+    stop_err("Error performing linear regression on the input data.\nEither the response column or one of the predictor columns contain only non-numeric or invalid values.")
+set_default_mode(BASIC_CONVERSION)
+
+coeffs = linear_model.as_py()['coefficients']
+yintercept = coeffs['(Intercept)']
+summary = r.summary(linear_model)
+
+co = summary.get('coefficients', 'NA')
+"""
+if len(co) != len(x_vals)+1:
+    stop_err("Stopped performing linear regression on the input data, since one of the predictor columns contains only non-numeric or invalid values.")
+"""
+
+try:
+    yintercept = r.round(float(yintercept), digits=10)
+    pvaly = r.round(float(co[0][3]), digits=10)
+except:
+    pass
+
+print >> fout, "Y-intercept\t%s" % (yintercept)
+print >> fout, "p-value (Y-intercept)\t%s" % (pvaly)
+
+if len(x_vals) == 1:    #Simple linear  regression case with 1 predictor variable
+    try:
+        slope = r.round(float(coeffs['x']), digits=10)
+    except:
+        slope = 'NA'
+    try:
+        pval = r.round(float(co[1][3]), digits=10)
+    except:
+        pval = 'NA'
+    print >> fout, "Slope (c%d)\t%s" % ( x_cols[0]+1, slope )
+    print >> fout, "p-value (c%d)\t%s" % ( x_cols[0]+1, pval )
+else:    #Multiple regression case with >1 predictors
+    ind = 1
+    while ind < len(coeffs.keys()):
+        try:
+            slope = r.round(float(coeffs['x'+str(ind)]), digits=10)
+        except:
+            slope = 'NA'
+        print >> fout, "Slope (c%d)\t%s" % ( x_cols[ind-1]+1, slope )
+        try:
+            pval = r.round(float(co[ind][3]), digits=10)
+        except:
+            pval = 'NA'
+        print >> fout, "p-value (c%d)\t%s" % ( x_cols[ind-1]+1, pval )
+        ind += 1
+
+rsq = summary.get('r.squared','NA')
+adjrsq = summary.get('adj.r.squared','NA')
+fstat = summary.get('fstatistic','NA')
+sigma = summary.get('sigma','NA')
+
+try:
+    rsq = r.round(float(rsq), digits=5)
+    adjrsq = r.round(float(adjrsq), digits=5)
+    fval = r.round(fstat['value'], digits=5)
+    fstat['value'] = str(fval)
+    sigma = r.round(float(sigma), digits=10)
+except:
+    pass
+
+print >> fout, "R-squared\t%s" % (rsq)
+print >> fout, "Adjusted R-squared\t%s" % (adjrsq)
+print >> fout, "F-statistic\t%s" % (fstat)
+print >> fout, "Sigma\t%s" % (sigma)
+
+r.pdf( outfile2, 8, 8 )
+if len(x_vals) == 1:    #Simple linear  regression case with 1 predictor variable
+    sub_title =  "Slope = %s; Y-int = %s" % ( slope, yintercept )
+    try:
+        r.plot(x=x_vals[0], y=y_vals, xlab="X", ylab="Y", sub=sub_title, main="Scatterplot with regression")
+        r.abline(a=yintercept, b=slope, col="red")
+    except:
+        pass
+else:
+    r.pairs(dat, main="Scatterplot Matrix", col="blue")
+try:
+    r.plot(linear_model)
+except:
+    pass
+r.dev_off()
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/linear_regression.xml	Tue Apr 01 09:12:59 2014 -0400
@@ -0,0 +1,76 @@
+<tool id="LinearRegression1" name="Perform Linear Regression" version="1.0.1">
+  <description> </description>
+  <requirements>
+    <requirement type="package" version="2.11.0">R</requirement>
+    <requirement type="package" version="1.7.1">numpy</requirement>
+    <requirement type="package" version="1.0.3">rpy</requirement>
+  </requirements>
+  <command interpreter="python">
+    linear_regression.py 
+      $input1
+      $response_col
+      $predictor_cols
+      $out_file1
+      $out_file2
+      1>/dev/null
+  </command>
+  <inputs>
+    <param format="tabular" name="input1" type="data" label="Select data" help="Dataset missing? See TIP below."/>
+    <param name="response_col" label="Response column (Y)" type="data_column" data_ref="input1" numerical="True"/>
+    <param name="predictor_cols" label="Predictor columns (X)" type="data_column" data_ref="input1" numerical="True" multiple="true" >
+        <validator type="no_options" message="Please select at least one column."/>
+    </param>
+  </inputs>
+  <outputs>
+    <data format="input" name="out_file1" metadata_source="input1" />
+    <data format="pdf" name="out_file2" />
+  </outputs>
+  <requirements>
+    <requirement type="python-module">rpy</requirement>
+  </requirements>
+  <tests>
+    <test>
+        <param name="input1" value="regr_inp.tabular"/>
+        <param name="response_col" value="3"/>
+        <param name="predictor_cols" value="1,2"/>
+        <output name="out_file1" file="regr_out.tabular"/>
+        <output name="out_file2" file="regr_out.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 uses the 'lm' function from R statistical package to perform linear regression on the input data. It outputs two files, one containing the summary statistics of the performed regression, and the other containing diagnostic plots to check whether model assumptions are satisfied.   
+
+*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 predictor and response 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:
+
+  - sigma: the square root of the estimated variance of the random error (standard error of the residiuals)
+  - R-squared: the fraction of variance explained by the model
+  - Adjusted R-squared: the above R-squared statistic adjusted, penalizing for the number of the predictors (p)
+  - p-value: p-value for the t-test of the null hypothesis that the corresponding slope is equal to zero against the two-sided alternative.
+
+
+  </help>
+</tool>
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/regr_inp.tabular	Tue Apr 01 09:12:59 2014 -0400
@@ -0,0 +1,10 @@
+7	33	42
+4	41	33
+16	7	75
+3	49	28
+21	5	91
+8	31	55
+7	35	52
+5	30	16
+15	10	69
+20	10	94
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/regr_out.pdf	Tue Apr 01 09:12:59 2014 -0400
@@ -0,0 +1,1590 @@
+%PDF-1.4
+%ρ\r
+1 0 obj
+<<
+/CreationDate (D:20110308105934)
+/ModDate (D:20110308105934)
+/Title (R Graphics Output)
+/Producer (R 2.11.0)
+/Creator (R)
+>>
+endobj
+2 0 obj
+<<
+/Type /Catalog
+/Pages 3 0 R
+>>
+endobj
+5 0 obj
+<<
+/Type /Page
+/Parent 3 0 R
+/Contents 6 0 R
+/Resources 4 0 R
+>>
+endobj
+6 0 obj
+<<
+/Length 7 0 R
+>>
+stream
+1 J 1 j q
+Q q 38.02 358.66 166.66 160.32 re W n
+Q q
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+42.77 363.41 m
+199.92 363.41 l
+199.92 514.22 l
+42.77 514.22 l
+42.77 363.41 l
+S
+Q q 42.77 363.41 157.15 150.82 re W n
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 16.00 0.00 -0.00 16.00 117.34 436.34 Tm (y) Tj
+ET
+Q q 209.42 363.41 157.15 150.82 re W n
+Q q 204.67 358.66 166.66 160.32 re W n
+Q q
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+209.42 363.41 m
+366.58 363.41 l
+366.58 514.22 l
+209.42 514.22 l
+209.42 363.41 l
+S
+231.41 514.22 m 352.67 514.22 l S
+231.41 514.22 m 231.41 518.98 l S
+271.83 514.22 m 271.83 518.98 l S
+312.25 514.22 m 312.25 518.98 l S
+352.67 514.22 m 352.67 518.98 l S
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 8.00 0.00 -0.00 8.00 229.19 525.63 Tm (5) Tj
+ET
+BT
+/F2 1 Tf 8.00 0.00 -0.00 8.00 267.38 525.63 Tm (10) Tj
+ET
+BT
+/F2 1 Tf 8.00 0.00 -0.00 8.00 307.80 525.63 Tm (15) Tj
+ET
+BT
+/F2 1 Tf 8.00 0.00 -0.00 8.00 348.22 525.63 Tm (20) Tj
+ET
+Q q 209.42 363.41 157.15 150.82 re W n
+0.000 0.000 1.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 245.63 413.83 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 221.37 397.72 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 318.38 472.91 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 213.29 388.76 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 358.80 501.55 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 253.71 437.10 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 245.63 431.73 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 229.46 367.28 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 310.30 462.17 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 350.72 506.93 Tm (l) Tj 0 Tr
+ET
+Q q 376.08 363.41 157.15 150.82 re W n
+Q q 371.33 358.66 166.66 160.32 re W n
+Q q
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+376.08 363.41 m
+533.23 363.41 l
+533.23 514.22 l
+376.08 514.22 l
+376.08 363.41 l
+S
+533.23 376.16 m 533.23 483.57 l S
+533.23 376.16 m 537.98 376.16 l S
+533.23 411.96 m 537.98 411.96 l S
+533.23 447.77 m 537.98 447.77 l S
+533.23 483.57 m 537.98 483.57 l S
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 0.00 8.00 -8.00 0.00 550.34 371.71 Tm (20) Tj
+ET
+BT
+/F2 1 Tf 0.00 8.00 -8.00 0.00 550.34 407.51 Tm (40) Tj
+ET
+BT
+/F2 1 Tf 0.00 8.00 -8.00 0.00 550.34 443.32 Tm (60) Tj
+ET
+BT
+/F2 1 Tf 0.00 8.00 -8.00 0.00 550.34 479.13 Tm (80) Tj
+ET
+Q q 376.08 363.41 157.15 150.82 re W n
+0.000 0.000 1.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 472.54 413.83 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 499.00 397.72 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 386.56 472.91 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 525.46 388.76 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 379.95 501.55 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 465.93 437.10 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 479.16 431.73 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 462.62 367.28 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 396.48 462.17 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 396.48 506.93 Tm (l) Tj 0 Tr
+ET
+Q q 42.77 203.09 157.15 150.82 re W n
+Q q 38.02 198.34 166.66 160.32 re W n
+Q q
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+42.77 203.09 m
+199.92 203.09 l
+199.92 353.90 l
+42.77 353.90 l
+42.77 203.09 l
+S
+42.77 224.19 m 42.77 340.56 l S
+42.77 224.19 m 38.02 224.19 l S
+42.77 262.98 m 38.02 262.98 l S
+42.77 301.77 m 38.02 301.77 l S
+42.77 340.56 m 38.02 340.56 l S
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 0.00 8.00 -8.00 0.00 31.36 221.97 Tm (5) Tj
+ET
+BT
+/F2 1 Tf 0.00 8.00 -8.00 0.00 31.36 258.53 Tm (10) Tj
+ET
+BT
+/F2 1 Tf 0.00 8.00 -8.00 0.00 31.36 297.32 Tm (15) Tj
+ET
+BT
+/F2 1 Tf 0.00 8.00 -8.00 0.00 31.36 336.11 Tm (20) Tj
+ET
+Q q 42.77 203.09 157.15 150.82 re W n
+0.000 0.000 1.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 95.14 237.99 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 78.35 214.72 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 156.70 307.82 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 69.02 206.96 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 186.55 346.61 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 119.39 245.75 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 113.79 237.99 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 46.63 222.48 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 145.51 300.06 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 192.14 338.85 Tm (l) Tj 0 Tr
+ET
+Q q 209.42 203.09 157.15 150.82 re W n
+Q q 204.67 198.34 166.66 160.32 re W n
+Q q
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+209.42 203.09 m
+366.58 203.09 l
+366.58 353.90 l
+209.42 353.90 l
+209.42 203.09 l
+S
+Q q 209.42 203.09 157.15 150.82 re W n
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 16.00 0.00 -0.00 16.00 277.33 272.87 Tm (x.1) Tj
+ET
+Q q 376.08 203.09 157.15 150.82 re W n
+Q q 371.33 198.34 166.66 160.32 re W n
+Q q
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+376.08 203.09 m
+533.23 203.09 l
+533.23 353.90 l
+376.08 353.90 l
+376.08 203.09 l
+S
+Q q 376.08 203.09 157.15 150.82 re W n
+0.000 0.000 1.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 472.54 237.99 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 499.00 214.72 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 386.56 307.82 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 525.46 206.96 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 379.95 346.61 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 465.93 245.75 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 479.16 237.99 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 462.62 222.48 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 396.48 300.06 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 396.48 338.85 Tm (l) Tj 0 Tr
+ET
+Q q 42.77 42.77 157.15 150.82 re W n
+Q q 38.02 38.02 166.66 160.32 re W n
+Q q
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+42.77 42.77 m
+199.92 42.77 l
+199.92 193.58 l
+42.77 193.58 l
+42.77 42.77 l
+S
+56.05 42.77 m 167.98 42.77 l S
+56.05 42.77 m 56.05 38.02 l S
+93.36 42.77 m 93.36 38.02 l S
+130.67 42.77 m 130.67 38.02 l S
+167.98 42.77 m 167.98 38.02 l S
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 8.00 0.00 -0.00 8.00 51.60 25.66 Tm (20) Tj
+ET
+BT
+/F2 1 Tf 8.00 0.00 -0.00 8.00 88.91 25.66 Tm (40) Tj
+ET
+BT
+/F2 1 Tf 8.00 0.00 -0.00 8.00 126.22 25.66 Tm (60) Tj
+ET
+BT
+/F2 1 Tf 8.00 0.00 -0.00 8.00 163.53 25.66 Tm (80) Tj
+ET
+Q q 42.77 42.77 157.15 150.82 re W n
+0.000 0.000 1.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 95.14 135.51 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 78.35 160.90 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 156.70 52.99 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 69.02 186.29 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 186.55 46.64 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 119.39 129.16 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 113.79 141.85 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 46.63 125.98 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 145.51 62.51 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 192.14 62.51 Tm (l) Tj 0 Tr
+ET
+Q q 209.42 42.77 157.15 150.82 re W n
+Q q 204.67 38.02 166.66 160.32 re W n
+Q q
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+209.42 42.77 m
+366.58 42.77 l
+366.58 193.58 l
+209.42 193.58 l
+209.42 42.77 l
+S
+Q q 209.42 42.77 157.15 150.82 re W n
+0.000 0.000 1.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 245.63 135.51 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 221.37 160.90 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 318.38 52.99 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 213.29 186.29 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 358.80 46.64 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 253.71 129.16 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 245.63 141.85 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 229.46 125.98 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 310.30 62.51 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 4.94 0 0 4.94 350.72 62.51 Tm (l) Tj 0 Tr
+ET
+Q q 376.08 42.77 157.15 150.82 re W n
+Q q 371.33 38.02 166.66 160.32 re W n
+Q q
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+376.08 42.77 m
+533.23 42.77 l
+533.23 193.58 l
+376.08 193.58 l
+376.08 42.77 l
+S
+398.44 42.77 m 530.72 42.77 l S
+398.44 42.77 m 398.44 38.02 l S
+431.51 42.77 m 431.51 38.02 l S
+464.58 42.77 m 464.58 38.02 l S
+497.65 42.77 m 497.65 38.02 l S
+530.72 42.77 m 530.72 38.02 l S
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 8.00 0.00 -0.00 8.00 393.99 25.66 Tm (10) Tj
+ET
+BT
+/F2 1 Tf 8.00 0.00 -0.00 8.00 427.06 25.66 Tm (20) Tj
+ET
+BT
+/F2 1 Tf 8.00 0.00 -0.00 8.00 460.13 25.66 Tm (30) Tj
+ET
+BT
+/F2 1 Tf 8.00 0.00 -0.00 8.00 493.20 25.66 Tm (40) Tj
+ET
+BT
+/F2 1 Tf 8.00 0.00 -0.00 8.00 526.27 25.66 Tm (50) Tj
+ET
+533.23 64.22 m 533.23 191.17 l S
+533.23 64.22 m 537.98 64.22 l S
+533.23 95.96 m 537.98 95.96 l S
+533.23 127.70 m 537.98 127.70 l S
+533.23 159.43 m 537.98 159.43 l S
+533.23 191.17 m 537.98 191.17 l S
+BT
+/F2 1 Tf 0.00 8.00 -8.00 0.00 550.34 59.77 Tm (10) Tj
+ET
+BT
+/F2 1 Tf 0.00 8.00 -8.00 0.00 550.34 91.51 Tm (20) Tj
+ET
+BT
+/F2 1 Tf 0.00 8.00 -8.00 0.00 550.34 123.25 Tm (30) Tj
+ET
+BT
+/F2 1 Tf 0.00 8.00 -8.00 0.00 550.34 154.99 Tm (40) Tj
+ET
+BT
+/F2 1 Tf 0.00 8.00 -8.00 0.00 550.34 186.72 Tm (50) Tj
+ET
+Q q 376.08 42.77 157.15 150.82 re W n
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 16.00 0.00 -0.00 16.00 443.98 112.55 Tm (x.2) Tj
+ET
+Q q
+BT
+0.000 0.000 0.000 rg
+/F3 1 Tf 14.00 0.00 -0.00 14.00 228.88 549.39 Tm (Scatterplot Matrix) Tj
+ET
+Q
+endstream
+endobj
+7 0 obj
+9369
+endobj
+8 0 obj
+<<
+/Type /Page
+/Parent 3 0 R
+/Contents 9 0 R
+/Resources 4 0 R
+>>
+endobj
+9 0 obj
+<<
+/Length 10 0 R
+>>
+stream
+1 J 1 j q
+Q q 59.04 73.44 486.72 443.52 re W n
+Q q
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+91.71 73.44 m 481.63 73.44 l S
+91.71 73.44 m 91.71 66.24 l S
+156.70 73.44 m 156.70 66.24 l S
+221.68 73.44 m 221.68 66.24 l S
+286.67 73.44 m 286.67 66.24 l S
+351.66 73.44 m 351.66 66.24 l S
+416.64 73.44 m 416.64 66.24 l S
+481.63 73.44 m 481.63 66.24 l S
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 85.04 47.52 Tm (30) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 150.02 47.52 Tm (40) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 215.01 47.52 Tm (50) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 280.00 47.52 Tm (60) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 344.98 47.52 Tm (70) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 409.97 47.52 Tm (80) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 474.96 47.52 Tm (90) Tj
+ET
+59.04 147.85 m 59.04 487.37 l S
+59.04 147.85 m 51.84 147.85 l S
+59.04 232.73 m 51.84 232.73 l S
+59.04 317.61 m 51.84 317.61 l S
+59.04 402.49 m 51.84 402.49 l S
+59.04 487.37 m 51.84 487.37 l S
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 137.67 Tm (-10) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 225.89 Tm (-5) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 314.27 Tm (0) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 399.15 Tm (5) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 480.69 Tm (10) Tj
+ET
+59.04 73.44 m
+545.76 73.44 l
+545.76 516.96 l
+59.04 516.96 l
+59.04 73.44 l
+S
+Q q
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 268.54 18.72 Tm [(Fitted v) 25 (alues)] TJ
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 12.96 268.86 Tm (Residuals) Tj
+ET
+Q q 59.04 73.44 486.72 443.52 re W n
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 166.43 315.79 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 101.40 332.89 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 349.78 397.05 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 111.06 222.76 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 524.77 211.55 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 192.03 469.62 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 178.19 454.85 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 74.10 115.59 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 330.06 346.71 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 516.81 283.29 Tm (l) Tj 0 Tr
+ET
+1.000 0.000 0.000 RG
+77.07 141.77 m
+104.36 215.63 l
+114.02 240.45 l
+169.39 389.46 l
+181.15 422.49 l
+194.99 463.49 l
+333.02 381.33 l
+352.74 378.55 l
+519.77 253.25 l
+527.73 247.63 l
+S
+Q q
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 99.28 4.32 Tm [(function \(f) 30 (or) -25 (m) 10 (ula, data, subset, w) 10 (eights) 15 (, na.action, method = "qr", model = ...)] TJ
+ET
+Q q
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 251.72 523.44 Tm (Residuals vs Fitted) Tj
+ET
+Q q
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 9.00 0.00 -0.00 9.00 80.67 113.16 Tm (8) Tj
+ET
+BT
+/F2 1 Tf 9.00 0.00 -0.00 9.00 186.39 470.06 Tm (6) Tj
+ET
+BT
+/F2 1 Tf 9.00 0.00 -0.00 9.00 172.54 455.29 Tm (7) Tj
+ET
+Q q 59.04 73.44 486.72 443.52 re W n
+0.745 0.745 0.745 RG
+0.75 w
+[ 0.00 4.00] 0 d
+1 J
+1 j
+10.00 M
+59.04 317.61 m 545.76 317.61 l S
+Q
+endstream
+endobj
+10 0 obj
+3098
+endobj
+11 0 obj
+<<
+/Type /Page
+/Parent 3 0 R
+/Contents 12 0 R
+/Resources 4 0 R
+>>
+endobj
+12 0 obj
+<<
+/Length 13 0 R
+>>
+stream
+1 J 1 j q
+Q q 59.04 73.44 486.72 443.52 re W n
+Q q 59.04 73.44 486.72 443.52 re W n
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 281.58 323.65 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 317.30 340.84 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 394.93 412.99 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 203.95 210.25 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 153.67 203.79 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 524.77 469.29 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 445.20 457.77 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 74.10 87.27 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 354.14 356.10 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 244.74 284.34 Tm (l) Tj 0 Tr
+ET
+Q q
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+83.86 73.44 m 520.94 73.44 l S
+83.86 73.44 m 83.86 66.24 l S
+156.71 73.44 m 156.71 66.24 l S
+229.55 73.44 m 229.55 66.24 l S
+302.40 73.44 m 302.40 66.24 l S
+375.25 73.44 m 375.25 66.24 l S
+448.09 73.44 m 448.09 66.24 l S
+520.94 73.44 m 520.94 66.24 l S
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 72.02 47.52 Tm (-1.5) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 144.86 47.52 Tm (-1.0) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 217.71 47.52 Tm (-0.5) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 294.06 47.52 Tm (0.0) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 366.91 47.52 Tm (0.5) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 439.75 47.52 Tm (1.0) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 512.60 47.52 Tm (1.5) Tj
+ET
+59.04 100.87 m 59.04 437.81 l S
+59.04 100.87 m 51.84 100.87 l S
+59.04 213.18 m 51.84 213.18 l S
+59.04 325.50 m 51.84 325.50 l S
+59.04 437.81 m 51.84 437.81 l S
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 94.03 Tm (-2) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 206.34 Tm (-1) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 322.16 Tm (0) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 434.48 Tm (1) Tj
+ET
+59.04 73.44 m
+545.76 73.44 l
+545.76 516.96 l
+59.04 516.96 l
+59.04 73.44 l
+S
+Q q
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 245.71 18.72 Tm (Theoretical Quantiles) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 12.96 234.26 Tm [(Standardiz) 15 (ed residuals)] TJ
+ET
+Q q 59.04 73.44 486.72 443.52 re W n
+0.498 0.498 0.498 RG
+0.75 w
+[ 0.00 4.00] 0 d
+1 J
+1 j
+10.00 M
+59.04 105.88 m 545.76 526.86 l S
+Q q
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 99.28 4.32 Tm [(function \(f) 30 (or) -25 (m) 10 (ula, data, subset, w) 10 (eights) 15 (, na.action, method = "qr", model = ...)] TJ
+ET
+Q q
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 268.41 523.44 Tm [(Nor) -25 (mal Q-Q)] TJ
+ET
+Q q
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 9.00 0.00 -0.00 9.00 80.67 87.71 Tm (8) Tj
+ET
+BT
+/F2 1 Tf 9.00 0.00 -0.00 9.00 519.13 469.73 Tm (6) Tj
+ET
+BT
+/F2 1 Tf 9.00 0.00 -0.00 9.00 439.56 458.21 Tm (7) Tj
+ET
+Q
+endstream
+endobj
+13 0 obj
+2849
+endobj
+14 0 obj
+<<
+/Type /Page
+/Parent 3 0 R
+/Contents 15 0 R
+/Resources 4 0 R
+>>
+endobj
+15 0 obj
+<<
+/Length 16 0 R
+>>
+stream
+1 J 1 j q
+Q q 59.04 73.44 486.72 443.52 re W n
+Q q
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+91.71 73.44 m 481.63 73.44 l S
+91.71 73.44 m 91.71 66.24 l S
+156.70 73.44 m 156.70 66.24 l S
+221.68 73.44 m 221.68 66.24 l S
+286.67 73.44 m 286.67 66.24 l S
+351.66 73.44 m 351.66 66.24 l S
+416.64 73.44 m 416.64 66.24 l S
+481.63 73.44 m 481.63 66.24 l S
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 85.04 47.52 Tm (30) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 150.02 47.52 Tm (40) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 215.01 47.52 Tm (50) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 280.00 47.52 Tm (60) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 344.98 47.52 Tm (70) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 409.97 47.52 Tm (80) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 474.96 47.52 Tm (90) Tj
+ET
+59.04 89.87 m 59.04 515.16 l S
+59.04 89.87 m 51.84 89.87 l S
+59.04 231.63 m 51.84 231.63 l S
+59.04 373.39 m 51.84 373.39 l S
+59.04 515.16 m 51.84 515.16 l S
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 81.53 Tm (0.0) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 223.29 Tm (0.5) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 365.05 Tm (1.0) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 506.82 Tm (1.5) Tj
+ET
+59.04 73.44 m
+545.76 73.44 l
+545.76 516.96 l
+59.04 516.96 l
+59.04 73.44 l
+S
+Q q
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 268.54 18.72 Tm [(Fitted v) 25 (alues)] TJ
+ET
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+12.17 229.30 m
+11.30 230.85 l
+15.84 232.40 l
+4.26 234.47 l
+4.26 360.55 l
+S
+15.84 236.28 m
+6.77 236.28 l
+S
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 236.83 Tm (S) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 244.84 Tm (t) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 248.17 Tm (a) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 254.84 Tm (n) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 261.52 Tm (d) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 268.19 Tm (a) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 274.86 Tm (r) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 278.86 Tm (d) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 285.53 Tm (i) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 288.19 Tm (z) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 294.19 Tm (e) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 300.86 Tm (d) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 307.54 Tm ( ) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 310.87 Tm (r) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 314.87 Tm (e) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 321.54 Tm (s) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 327.54 Tm (i) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 330.20 Tm (d) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 336.88 Tm (u) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 343.55 Tm (a) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 350.22 Tm (l) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 15.61 352.88 Tm (s) Tj
+ET
+15.84 359.44 m
+6.77 359.44 l
+S
+Q q 59.04 73.44 486.72 443.52 re W n
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 166.43 110.35 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 101.40 200.59 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 349.78 341.20 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 111.06 371.22 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 524.77 379.25 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 192.03 410.95 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 178.19 397.96 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 74.10 497.94 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 330.06 241.42 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 516.81 253.39 Tm (l) Tj 0 Tr
+ET
+1.000 0.000 0.000 RG
+77.07 461.38 m
+104.36 370.36 l
+114.02 345.81 l
+169.39 376.18 l
+181.15 388.28 l
+194.99 403.24 l
+333.02 304.25 l
+352.74 305.93 l
+519.77 321.89 l
+527.73 323.06 l
+S
+Q q
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 99.28 4.32 Tm [(function \(f) 30 (or) -25 (m) 10 (ula, data, subset, w) 10 (eights) 15 (, na.action, method = "qr", model = ...)] TJ
+ET
+Q q
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 261.21 523.44 Tm (Scale-Location) Tj
+ET
+Q q
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 9.00 0.00 -0.00 9.00 80.67 498.38 Tm (8) Tj
+ET
+BT
+/F2 1 Tf 9.00 0.00 -0.00 9.00 186.39 411.40 Tm (6) Tj
+ET
+BT
+/F2 1 Tf 9.00 0.00 -0.00 9.00 172.54 398.40 Tm (7) Tj
+ET
+Q
+endstream
+endobj
+16 0 obj
+4330
+endobj
+17 0 obj
+<<
+/Type /Page
+/Parent 3 0 R
+/Contents 18 0 R
+/Resources 4 0 R
+>>
+endobj
+18 0 obj
+<<
+/Length 19 0 R
+>>
+stream
+1 J 1 j q
+Q q 59.04 73.44 486.72 443.52 re W n
+Q q
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+77.07 73.44 m 461.50 73.44 l S
+77.07 73.44 m 77.07 66.24 l S
+173.17 73.44 m 173.17 66.24 l S
+269.28 73.44 m 269.28 66.24 l S
+365.39 73.44 m 365.39 66.24 l S
+461.50 73.44 m 461.50 66.24 l S
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 68.73 47.52 Tm (0.0) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 164.83 47.52 Tm (0.1) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 260.94 47.52 Tm (0.2) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 357.05 47.52 Tm (0.3) Tj
+ET
+BT
+/F2 1 Tf 12.00 0.00 -0.00 12.00 453.16 47.52 Tm (0.4) Tj
+ET
+59.04 128.38 m 59.04 440.64 l S
+59.04 128.38 m 51.84 128.38 l S
+59.04 232.47 m 51.84 232.47 l S
+59.04 336.55 m 51.84 336.55 l S
+59.04 440.64 m 51.84 440.64 l S
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 121.54 Tm (-2) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 225.63 Tm (-1) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 333.22 Tm (0) Tj
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 41.76 437.30 Tm (1) Tj
+ET
+59.04 73.44 m
+545.76 73.44 l
+545.76 516.96 l
+59.04 516.96 l
+59.04 73.44 l
+S
+Q q
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 277.78 18.72 Tm [(Le) 30 (v) 25 (er) 10 (age)] TJ
+ET
+BT
+/F2 1 Tf 0.00 12.00 -12.00 0.00 12.96 234.26 Tm [(Standardiz) 15 (ed residuals)] TJ
+ET
+Q q 59.04 73.44 486.72 443.52 re W n
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 200.58 334.65 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 279.49 350.59 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 403.96 417.45 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 524.77 229.56 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 461.00 223.57 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 186.13 469.62 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 216.88 458.94 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 490.02 115.59 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 341.39 364.73 Tm (l) Tj 0 Tr
+ET
+BT
+/F1 1 Tf 1 Tr 7.48 0 0 7.48 520.06 298.23 Tm (l) Tj 0 Tr
+ET
+1.000 0.000 0.000 RG
+189.10 433.95 m
+203.55 424.65 l
+219.85 414.34 l
+282.45 379.59 l
+344.35 374.86 l
+406.92 335.40 l
+463.96 275.83 l
+492.98 246.97 l
+523.03 225.37 l
+527.73 222.96 l
+S
+0.745 0.745 0.745 RG
+0.75 w
+[ 0.00 4.00] 0 d
+59.04 336.55 m 545.76 336.55 l S
+77.07 73.44 m 77.07 516.96 l S
+Q q
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 99.28 4.32 Tm [(function \(f) 30 (or) -25 (m) 10 (ula, data, subset, w) 10 (eights) 15 (, na.action, method = "qr", model = ...)] TJ
+ET
+Q q 59.04 73.44 486.72 443.52 re W n
+1.000 0.000 0.000 RG
+0.75 w
+[ 3.00 5.00] 0 d
+1 J
+1 j
+10.00 M
+289.33 576.00 m
+290.46 575.18 l
+295.10 571.89 l
+299.74 568.69 l
+304.38 565.59 l
+309.02 562.57 l
+313.67 559.63 l
+318.31 556.76 l
+322.95 553.97 l
+327.59 551.24 l
+332.23 548.59 l
+336.88 545.99 l
+341.52 543.46 l
+346.16 540.98 l
+350.80 538.56 l
+355.44 536.19 l
+360.09 533.87 l
+364.73 531.60 l
+369.37 529.38 l
+374.01 527.20 l
+378.65 525.07 l
+383.29 522.97 l
+387.94 520.92 l
+392.58 518.90 l
+397.22 516.92 l
+401.86 514.98 l
+406.50 513.07 l
+411.15 511.20 l
+415.79 509.35 l
+420.43 507.54 l
+425.07 505.75 l
+429.71 504.00 l
+434.36 502.27 l
+439.00 500.57 l
+443.64 498.90 l
+448.28 497.25 l
+452.92 495.62 l
+457.56 494.02 l
+462.21 492.44 l
+466.85 490.89 l
+471.49 489.35 l
+476.13 487.84 l
+480.77 486.34 l
+485.42 484.87 l
+490.06 483.41 l
+494.70 481.97 l
+499.34 480.55 l
+503.98 479.15 l
+508.63 477.76 l
+513.27 476.39 l
+517.91 475.04 l
+522.55 473.70 l
+527.19 472.37 l
+531.83 471.06 l
+536.48 469.77 l
+541.12 468.48 l
+545.76 467.22 l
+S
+197.66 0.00 m
+202.26 7.16 l
+206.90 14.00 l
+211.55 20.50 l
+216.19 26.69 l
+220.83 32.60 l
+225.47 38.24 l
+230.11 43.64 l
+234.75 48.81 l
+239.40 53.78 l
+244.04 58.55 l
+248.68 63.13 l
+253.32 67.55 l
+257.96 71.81 l
+262.61 75.92 l
+267.25 79.90 l
+271.89 83.74 l
+276.53 87.45 l
+281.17 91.05 l
+285.82 94.54 l
+290.46 97.93 l
+295.10 101.22 l
+299.74 104.41 l
+304.38 107.52 l
+309.02 110.54 l
+313.67 113.48 l
+318.31 116.35 l
+322.95 119.14 l
+327.59 121.86 l
+332.23 124.52 l
+336.88 127.12 l
+341.52 129.65 l
+346.16 132.13 l
+350.80 134.55 l
+355.44 136.92 l
+360.09 139.24 l
+364.73 141.51 l
+369.37 143.73 l
+374.01 145.91 l
+378.65 148.04 l
+383.29 150.14 l
+387.94 152.19 l
+392.58 154.21 l
+397.22 156.18 l
+401.86 158.13 l
+406.50 160.04 l
+411.15 161.91 l
+415.79 163.76 l
+420.43 165.57 l
+425.07 167.35 l
+429.71 169.11 l
+434.36 170.83 l
+439.00 172.53 l
+443.64 174.21 l
+448.28 175.86 l
+452.92 177.48 l
+457.56 179.08 l
+462.21 180.66 l
+466.85 182.22 l
+471.49 183.76 l
+476.13 185.27 l
+480.77 186.77 l
+485.42 188.24 l
+490.06 189.70 l
+494.70 191.14 l
+499.34 192.56 l
+503.98 193.96 l
+508.63 195.35 l
+513.27 196.72 l
+517.91 198.07 l
+522.55 199.41 l
+527.19 200.73 l
+531.83 202.04 l
+536.48 203.34 l
+541.12 204.62 l
+545.76 205.89 l
+S
+424.78 576.00 m
+425.07 575.84 l
+429.71 573.36 l
+434.36 570.92 l
+439.00 568.51 l
+443.64 566.14 l
+448.28 563.81 l
+452.92 561.51 l
+457.56 559.25 l
+462.21 557.02 l
+466.85 554.81 l
+471.49 552.64 l
+476.13 550.50 l
+480.77 548.39 l
+485.42 546.30 l
+490.06 544.24 l
+494.70 542.21 l
+499.34 540.20 l
+503.98 538.21 l
+508.63 536.25 l
+513.27 534.31 l
+517.91 532.40 l
+522.55 530.51 l
+527.19 528.63 l
+531.83 526.78 l
+536.48 524.95 l
+541.12 523.13 l
+545.76 521.34 l
+S
+291.37 0.00 m
+295.10 3.74 l
+299.74 8.26 l
+304.38 12.65 l
+309.02 16.92 l
+313.67 21.08 l
+318.31 25.13 l
+322.95 29.08 l
+327.59 32.94 l
+332.23 36.70 l
+336.88 40.36 l
+341.52 43.95 l
+346.16 47.45 l
+350.80 50.88 l
+355.44 54.23 l
+360.09 57.50 l
+364.73 60.71 l
+369.37 63.86 l
+374.01 66.94 l
+378.65 69.96 l
+383.29 72.92 l
+387.94 75.82 l
+392.58 78.67 l
+397.22 81.47 l
+401.86 84.22 l
+406.50 86.92 l
+411.15 89.57 l
+415.79 92.18 l
+420.43 94.75 l
+425.07 97.27 l
+429.71 99.75 l
+434.36 102.19 l
+439.00 104.60 l
+443.64 106.96 l
+448.28 109.30 l
+452.92 111.59 l
+457.56 113.86 l
+462.21 116.09 l
+466.85 118.29 l
+471.49 120.46 l
+476.13 122.61 l
+480.77 124.72 l
+485.42 126.81 l
+490.06 128.87 l
+494.70 130.90 l
+499.34 132.91 l
+503.98 134.89 l
+508.63 136.85 l
+513.27 138.79 l
+517.91 140.71 l
+522.55 142.60 l
+527.19 144.48 l
+531.83 146.33 l
+536.48 148.16 l
+541.12 149.98 l
+545.76 151.77 l
+S
+69.84 87.84 m 91.44 87.84 l S
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 102.24 83.53 Tm (Cook's distance) Tj
+ET
+Q q
+0.000 0.000 0.000 RG
+0.75 w
+[] 0 d
+1 J
+1 j
+10.00 M
+545.76 151.77 m 545.76 516.96 l S
+545.76 151.77 m 545.76 151.77 l S
+545.76 205.89 m 545.76 205.89 l S
+545.76 467.22 m 545.76 467.22 l S
+BT
+1.000 0.000 0.000 rg
+/F2 1 Tf 9.00 0.00 -0.00 9.00 549.36 148.54 Tm (1) Tj
+ET
+BT
+/F2 1 Tf 9.00 0.00 -0.00 9.00 549.36 202.66 Tm (0.5) Tj
+ET
+BT
+/F2 1 Tf 9.00 0.00 -0.00 9.00 549.36 463.98 Tm (0.5) Tj
+ET
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 12.00 0.00 -0.00 12.00 242.10 523.44 Tm [(Residuals vs Le) 30 (v) 25 (er) 10 (age)] TJ
+ET
+Q q
+BT
+0.000 0.000 0.000 rg
+/F2 1 Tf 9.00 0.00 -0.00 9.00 496.58 113.16 Tm (8) Tj
+ET
+BT
+/F2 1 Tf 9.00 0.00 -0.00 9.00 519.13 227.13 Tm (4) Tj
+ET
+BT
+/F2 1 Tf 9.00 0.00 -0.00 9.00 467.56 221.14 Tm (5) Tj
+ET
+Q
+endstream
+endobj
+19 0 obj
+6890
+endobj
+3 0 obj
+<<
+/Type /Pages
+/Kids [
+5 0 R
+8 0 R
+11 0 R
+14 0 R
+17 0 R
+]
+/Count 5
+/MediaBox [0 0 576 576]
+>>
+endobj
+4 0 obj
+<<
+/ProcSet [/PDF /Text]
+/Font << /F1 21 0 R /F2 22 0 R /F3 23 0 R >>
+/ExtGState << >>
+>>
+endobj
+20 0 obj
+<<
+/Type /Encoding
+/BaseEncoding /WinAnsiEncoding
+/Differences [ 45/minus 96/quoteleft
+144/dotlessi /grave /acute /circumflex /tilde /macron /breve /dotaccent
+/dieresis /.notdef /ring /cedilla /.notdef /hungarumlaut /ogonek /caron /space]
+>>
+endobj
+21 0 obj
+<<
+/Type /Font
+/Subtype /Type1
+/Name /F1
+/BaseFont /ZapfDingbats
+>>
+endobj
+22 0 obj <<
+/Type /Font
+/Subtype /Type1
+/Name /F2
+/BaseFont /Helvetica
+/Encoding 20 0 R
+>> endobj
+23 0 obj <<
+/Type /Font
+/Subtype /Type1
+/Name /F3
+/BaseFont /Helvetica-Bold
+/Encoding 20 0 R
+>> endobj
+xref
+0 24
+0000000000 65535 f 
+0000000021 00000 n 
+0000000164 00000 n 
+0000027531 00000 n 
+0000027641 00000 n 
+0000000213 00000 n 
+0000000293 00000 n 
+0000009715 00000 n 
+0000009735 00000 n 
+0000009815 00000 n 
+0000012967 00000 n 
+0000012988 00000 n 
+0000013070 00000 n 
+0000015974 00000 n 
+0000015995 00000 n 
+0000016077 00000 n 
+0000020462 00000 n 
+0000020483 00000 n 
+0000020565 00000 n 
+0000027510 00000 n 
+0000027746 00000 n 
+0000028004 00000 n 
+0000028088 00000 n 
+0000028186 00000 n 
+trailer
+<<
+/Size 24
+/Info 1 0 R
+/Root 2 0 R
+>>
+startxref
+28289
+%%EOF
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/regr_out.tabular	Tue Apr 01 09:12:59 2014 -0400
@@ -0,0 +1,10 @@
+Y-intercept	-28.1169126094
+p-value (Y-intercept)	0.2942632048
+Slope (c1)	5.7471341245
+p-value (c1)	0.001732877
+Slope (c2)	0.9042745374
+p-value (c2)	0.1146675468
+R-squared	0.93881
+Adjusted R-squared	0.92132
+F-statistic	{'dendf': 7.0, 'value': '53.69463', 'numdf': 2.0}
+Sigma	7.4345007113
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_dependencies.xml	Tue Apr 01 09:12:59 2014 -0400
@@ -0,0 +1,12 @@
+<?xml version="1.0"?>
+<tool_dependency>
+  <package name="numpy" version="1.7.1">
+      <repository changeset_revision="55a7a5e9d63f" name="package_numpy_1_7" owner="devteam" prior_installation_required="False" toolshed="http://testtoolshed.g2.bx.psu.edu" />
+    </package>
+    <package name="rpy" version="1.0.3">
+      <repository changeset_revision="c0eb80864491" name="package_rpy_1_0_3" owner="devteam" prior_installation_required="False" toolshed="http://testtoolshed.g2.bx.psu.edu" />
+    </package>
+    <package name="R" version="2.11.0">
+      <repository changeset_revision="2c0a13200a73" name="package_r_2_11_0" owner="devteam" prior_installation_required="False" toolshed="http://testtoolshed.g2.bx.psu.edu" />
+    </package>
+</tool_dependency>