Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
19.6000 16.0000 17.7000 19.8000 17.0000 17.4000 18.9000 15.7000 15.2000 15.8000 16.0000 16.1000 16.2000 12.5000 14.8000 15.4000 13.6000 14.2000 15.0000 14.1000 13.7000 14.4000 15.6000 19.7000 20.4000 16.1000 20.1000 20.6000 19.3000 20.5000 19.2000 19.0000 18.7000 16.5000 19.0000 20.5000 18.4000 16.2000 18.1000 19.3000 18.3000 17.2000 19.6000 17.2000 17.4000 16.0000 18.5000 18.4000 18.2000 14.9000 16.3000 18.3000 18.0000 15.9000 19.6000 16.6000 16.2000 16.6000 17.5000 16.2000 17.5000 13.8000 14.9000 17.2000 15.6000 16.2000 17.4000 15.1000 14.5000 15.1000 15.5000 15.9000 15.9000 12.3000 14.4000 16.0000 13.9000 14.7000 16.2000 13.8000 13.2000 13.5000 13.5000 15.0000 14.5000 10.5000 13.7000 13.9000 13.4000 14.0000 14.3000 13.3000 13.2000 12.6000 13.7000 15.6000 14.4000 11.0000 13.7000 13.8000 14.3000 14.0000 14.6000 13.1000 13.2000 11.6000 13.3000 14.4000 13.3000 11.3000 13.2000 14.1000 14.0000 12.9000 15.2000 13.6000 13.7000
Data Y:
18.9000 16.6000 17.2000 19.2000 17.1000 17.7000 18.7000 15.9000 16.0000 16.8000 16.0000 16.8000 16.3000 13.6000 14.3000 15.5000 13.9000 14.3000 15.8000 14.5000 15.1000 15.8000 17.2000 20.4000 21.3000 18.2000 20.2000 21.1000 19.7000 21.5000 20.2000 19.0000 20.2000 18.0000 19.5000 20.3000 18.0000 16.4000 17.8000 18.5000 18.2000 16.7000 19.1000 16.8000 17.5000 16.2000 17.9000 17.7000 17.2000 15.7000 15.2000 17.7000 17.4000 15.9000 19.7000 16.7000 16.9000 18.0000 17.6000 15.2000 16.5000 14.7000 14.1000 16.9000 15.2000 15.4000 16.8000 14.8000 14.1000 15.0000 14.8000 15.0000 15.1000 12.8000 13.0000 15.7000 12.8000 13.9000 15.4000 13.2000 12.7000 13.5000 12.8000 13.9000 13.3000 10.7000 12.3000 12.9000 12.5000 13.0000 13.9000 13.1000 13.1000 13.0000 12.8000 14.2000 13.0000 11.2000 12.1000 12.9000 13.2000 13.2000 13.5000 12.4000 12.4000 11.6000 12.6000 13.1000 12.3000 11.4000 11.8000 13.4000 13.6000 12.9000 14.5000 13.3000 13.5000
Chart options
Title:
Label y-axis:
Label x-axis:
R Code
bitmap(file='test1.png') 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) par(mar=c(0,4,1,1)) barplot(histx$counts, axes=FALSE, ylim=c(0, maxcounts), space=0) par(mar=c(4,0,1,1)) barplot(histy$counts, axes=FALSE, xlim=c(0, maxcounts), space=0, horiz=TRUE) dev.off() lx = length(x) makebiased = (lx-1)/lx varx = var(x)*makebiased vary = var(y)*makebiased corxy <- cor.test(x,y,method='pearson') cxy <- as.matrix(corxy$estimate)[1,1] load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Pearson Product Moment Correlation - Ungrouped Data',3,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Statistic',1,TRUE) a<-table.element(a,'Variable X',1,TRUE) a<-table.element(a,'Variable Y',1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,hyperlink('http://www.xycoon.com/arithmetic_mean.htm','Mean',''),header=TRUE) a<-table.element(a,mean(x)) a<-table.element(a,mean(y)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,hyperlink('http://www.xycoon.com/biased.htm','Biased Variance',''),header=TRUE) a<-table.element(a,varx) a<-table.element(a,vary) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,hyperlink('http://www.xycoon.com/biased1.htm','Biased Standard Deviation',''),header=TRUE) a<-table.element(a,sqrt(varx)) a<-table.element(a,sqrt(vary)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,hyperlink('http://www.xycoon.com/covariance.htm','Covariance',''),header=TRUE) a<-table.element(a,cov(x,y),2) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,hyperlink('http://www.xycoon.com/pearson_correlation.htm','Correlation',''),header=TRUE) a<-table.element(a,cxy,2) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,hyperlink('http://www.xycoon.com/coeff_of_determination.htm','Determination',''),header=TRUE) a<-table.element(a,cxy*cxy,2) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,hyperlink('http://www.xycoon.com/ttest_statistic.htm','T-Test',''),header=TRUE) a<-table.element(a,as.matrix(corxy$statistic)[1,1],2) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value (2 sided)',header=TRUE) a<-table.element(a,(p2 <- as.matrix(corxy$p.value)[1,1]),2) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value (1 sided)',header=TRUE) a<-table.element(a,p2/2,2) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Degrees of Freedom',header=TRUE) a<-table.element(a,lx-2,2) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Number of Observations',header=TRUE) a<-table.element(a,lx,2) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation