changeset 0:cf431604ec3e draft default tip

Imported from capsule None
author devteam
date Tue, 01 Apr 2014 10:52:17 -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 10:52:17 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 10:52:17 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 10:52:17 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 10:52:17 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 10:52:17 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 10:52:17 2014 -0400
@@ -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" prior_installation_required="False" 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" prior_installation_required="False" toolshed="http://toolshed.g2.bx.psu.edu" />
+    </package>
+    <package name="R" version="2.11.0">
+      <repository changeset_revision="6e1b17857732" name="package_r_2_11_0" owner="devteam" prior_installation_required="False" toolshed="http://toolshed.g2.bx.psu.edu" />
+    </package>
+</tool_dependency>