Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
14.69 20.15 11.78 21.47 17.54 21.47 18.55 14.33 8.78 17.66 19.64 13.25 16.17 9.66 20.45 24.21 19.65 21.33 14.47 16.22 16.66 19.87 14.55 19.41 20.63 9.12 8.36 19.36 25.78 14.66 27.55 22.63 30.33 22.22 11.99 9.47 10.33 15.36 14.99 17.31 18.94 16.54 18.21 11.74 17.41 9.99 21.14 19.70 14.87 19.65 27.65 13.45 20.00 13.47 18.54 20.00 11.47 28.65 24.14 19.87 17.65 13.69 15.47 11.63 10.24 13.54 11.00 19.74 21.25 29.54 11.41 10.87 8.47 12.69 10.14 18.63 22.54 17.65 16.41 17.65 19.87 21.45 22.65 20.69 23.74 20.96 18.35 17.84 20.54 21.63 27.61 19.47 23.52 25.64 19.33 20.64 18.98 15.67 23.74 22.89 24.66 29.40 19.87 17.55 15.94 21.74 23.51 24.84 18.47 21.01 16.74 27.45 19.54 17.74 21.65 24.74 26.89 23.22 20.16 21.88
Data Y:
11.54 25.47 34.78 29.21 43.9 33.47 55.10 48.13 33.25 25.33 26.98 88.74 36.54 11.78 29.11 32.47 36.99 27.55 18.55 14.33 8.78 40.55 22.63 20.47 60.03 50.14 30.63 36.99 27.55 18.55 14.33 21.63 36.12 23.69 54.12 45.16 69.41 41.96 13.25 26.98 88.74 36.54 11.78 29.11 32.47 55.17 17.66 19.64 44.47 50.17 21.63 36.12 23.69 33.25 25.33 26.98 14.33 8.78 17.66 19.64 13.25 16.17 9.66 20.45 24.21 19.65 21.25 33.64 25.44 60.74 71.33 55.21 66.23 19.65 23.52 25.64 19.33 20.64 18.98 15.67 20.00 11.47 28.65 24.14 19.87 17.65 30.54 28.57 29.65 31.22 32.00 19.64 41.88 45.98 34.12 30.88 25.96 20.33 26.98 27.12 21.66 27.96 28.63 29.10 34.25 39.77 41.23 14.10 19.54 8.78 24.55 34.87 31.47 21.65 19.77 9.23 7.55 22.44 16.89 47.55
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