Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
110.089 99.966 94.195 119.588 125.946 118.906 108.717 118.357 115.422 130.937 111.802 120.357 121.12 108.752 98.973 116.721 123.857 114.056 98.309 110.251 108.538 122.526 118.394 116.691 116.014 111.784 95.164 121.028 129.89 116.102 102.055 113.562 113.071 126.486 119.472 117.141 121.925 112.688 90.974 126.398 130.401 116.873 111.917 113.919 117.931 139.332 127.781 118.103 136.984 119.566 107.238 139.389 139.798 138.074 129.739 125.098 129.341 149.083 137.721 130.126 140.499 118.165 115.932 144.662 136.159 147.339 138.807 135.275 135.84 161.702 132.606 147.419 154.865 131.545 126.212 150.316 154.524 154.28 131.059 147.168 141.8 162.022 143.924 151.406 159.601 143.513 132.302 152.021 171.573 161.591 134.057 172.247 173.384 173.706 188.178 165.932 179.795
Data Y:
34.679 32.566 37.111 35.977 41.811 43.465 36.466 39.188 35.318 39.76 35.102 38.829 36.304 35.199 33.761 30.579 36.467 35.855 30.077 35.45 32.956 33.903 35.319 37.387 34.77 33.383 32.644 34.524 40.451 36.605 35.849 38.199 34.092 36.435 36.733 36.452 34.213 35.314 32.192 35.523 37.877 33.927 35.5 36.819 35.191 42.832 40.462 38.249 47.931 39.777 39.834 41.529 42.861 43.481 39.493 41.667 42.179 45.282 42.322 44.709 51.212 41.708 43.986 44.914 42.518 49.019 47.452 46.669 42.104 52.736 40.953 46.101 47.181 41.855 43.147 44.226 50.592 50.436 44.366 49.226 44.548 50.587 47.689 51.718 49.893 52.021 48.384 48.765 58.453 56.863 47.538 66.204 55.972 56.435 63.573 59.091 61.658
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