Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1.3067 1.2894 1.2770 1.2208 1.2565 1.3406 1.3569 1.3686 1.4272 1.4614 1.4914 1.4816 1.4562 1.4268 1.4088 1.4016 1.3650 1.3190 1.3050 1.2785 1.3239 1.3449 1.2732 1.3322 1.4369 1.4975 1.5770 1.5553 1.5557 1.5750 1.5527 1.4748 1.4718 1.4570 1.4684 1.4227 1.3896 1.3622 1.3716 1.3419 1.3511 1.3516 1.3242 1.3074 1.2999 1.3213 1.2881 1.2611 1.2727 1.2811 1.2684 1.2650 1.2770 1.2271 1.2020 1.1938 1.2103 1.1856 1.1786 1.2015 1.2256 1.2292 1.2037 1.2165 1.2694 1.2938 1.3201 1.3014 1.3119 1.3408 1.2991 1.2490 1.2218 1.2176 1.2266 1.2138 1.2007 1.1985 1.2262 1.2646 1.2613 1.2286 1.1702 1.1692 1.1222 1.1139 1.1372 1.1663 1.1582 1.0848 1.0807 1.0773 1.0622 1.0183 1.0014 0.9811 0.9808 0.9778 0.9922 0.9554 0.9170 0.8858 0.8758 0.8700 0.8833 0.8924 0.8883 0.9059 0.9111 0.9005 0.8607 0.8532 0.8742 0.8920 0.9095 0.9217 0.9383
Data Y:
8.7000 8.9000 8.9000 8.1000 8.0000 8.3000 8.5000 8.7000 8.6000 8.3000 7.9000 7.9000 8.1000 8.3000 8.1000 7.4000 7.3000 7.7000 8.0000 8.0000 7.7000 6.9000 6.6000 6.9000 7.5000 7.9000 7.7000 6.5000 6.1000 6.4000 6.8000 7.1000 7.3000 7.2000 7.0000 7.0000 7.0000 7.3000 7.5000 7.2000 7.7000 8.0000 7.9000 8.0000 8.0000 7.9000 7.9000 8.0000 8.1000 8.1000 8.2000 8.0000 8.3000 8.5000 8.6000 8.7000 8.7000 8.5000 8.4000 8.5000 8.7000 8.7000 8.6000 7.9000 8.1000 8.2000 8.5000 8.6000 8.5000 8.3000 8.2000 8.7000 9.3000 9.3000 8.8000 7.4000 7.2000 7.5000 8.3000 8.8000 8.9000 8.6000 8.4000 8.4000 8.4000 8.4000 8.3000 7.6000 7.6000 7.9000 8.0000 8.2000 8.3000 8.2000 8.1000 8.0000 7.8000 7.6000 7.5000 6.8000 6.9000 7.1000 7.3000 7.4000 7.6000 7.6000 7.5000 7.5000 6.8000 6.4000 6.2000 6.0000 6.3000 6.3000 6.1000 6.1000 6.3000
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
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation