R version 3.3.2 (2016-10-31) -- "Sincere Pumpkin Patch"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

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

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

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

> 
> RC.capture <- function (expression, collapse = NULL) {
+     resultConn <- textConnection('RC.resultText', open = 'w', local=TRUE)
+     sink(resultConn)
+     on.exit(function() {
+         sink()
+         close(resultConn)
+     })
+     expression
+     on.exit(NULL)
+     sink()
+     close(resultConn)
+     return(paste(c(RC.resultText, ''), collapse = collapse, sep = ''))
+ }
> RC.texteval <- function (sourceText, collapse = NULL, echo = TRUE) {
+     sourceConn <- textConnection(sourceText, open = 'r')
+     on.exit(close(sourceConn))
+     result <- RC.capture(source(file = sourceConn, local = FALSE, echo = echo, print.eval = TRUE), collapse = collapse)
+     on.exit(NULL)
+     close(sourceConn)
+     res <- ''
+     for(i in 1:length(result)) {
+       if (result[i]!='') res <- paste(res,result[i],'
+ ',sep='')
+     }
+     return(res)
+ }
> y <- c(4,3,3,1,0,0,2,1,0,1)
> x <- c(0,1,1,3,9,5,2,4,2,2)
> ylab = 'Number of Pets'
> xlab = 'Number of Siblings'
> main = 'Scatterplot and histograms'
> #'GNU S' R Code compiled by R2WASP v. 1.2.327 (Thu, 24 Sep 2015 13:03:42 +0100)
> #Author: root
> #To cite this work: Wessa P., (2015), Pearson Correlation (v1.0.10) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_correlation.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #
> library(psychometric)
Loading required package: multilevel
Loading required package: nlme
Loading required package: MASS
> x <- x[!is.na(y)]
> y <- y[!is.na(y)]
> y <- y[!is.na(x)]
> x <- x[!is.na(x)]
> postscript(file="/var/wessaorg/rcomp/tmp/1ii2m1489518526.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=0.13888888888889,height=0.13888888888889) 
> histx <- hist(x, plot=FALSE)
> histy <- hist(y, plot=FALSE)
> maxcounts <- max(c(histx$counts, histx$counts))
> xrange <- c(min(x),max(x))
> yrange <- c(min(y),max(y))
> nf <- layout(matrix(c(2,0,1,3),2,2,byrow=TRUE), c(3,1), c(1,3), TRUE)
> par(mar=c(4,4,1,1))
> plot(x, y, xlim=xrange, ylim=yrange, xlab=xlab, ylab=ylab, sub=main)
Error in plot.new() : figure margins too large
Calls: plot -> plot.default -> plot.new
Execution halted