Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
30 27 18 28 27 28 21 22 20 23 18 18 25 25 25 24 13 17 4 16 21 20 22 0 18 29 15 22 22 17 27 15 26 18 27 17 19 13 16 2 26 22 23 18 18 19 22 14 23 20 13 16 7 17 19 23 19 16 20 25 17 12 24 14 18 19 16 19 4 20 24 17 22 19 22 23 17 24 20 25 20 21 22 25 28 29 20 20 19 19 26 10 17 30 22 23 16 18 25 18 24 23 24 15 20 26 23 23 22 15 22 10 20 23 27 23 25 20 24 23 22 21 25 27 23 23 19 15 20 16 25 25 19 16 19 19 23 21 19 20 3 23 15 24 24 24 28 23 25 25 20 25 23 20 16 23 11 23 29 16 23 20 4 24 16 3 23 20 19 24 27
Data Y:
16.6 14.85 11.75 18.45 19.9 18.45 15 11.35 18.1 13.4 13.9 15.25 16.1 17.35 13.15 12.15 18.2 13.6 14.1 14.9 16.25 15.65 14.6 19.2 13.2 15.65 7.65 15.2 11.85 11.4 19.9 15.15 16.85 12.6 12.35 16.65 13.95 15.7 15.35 15.1 17.75 16.65 19.1 13.35 18.4 16.15 18.4 15.6 16.35 17.65 11.7 14.35 14.75 9.9 16.85 15.6 17.1 19.1 7.6 14.75 13.6 11.9 16.35 14.35 17.75 19.3 17.1 19.05 18.55 19.1 13.35 17.6 16.1 11.95 7.7 14.6 12.6 18.9 11.6 14.6 13.85 15.9 10.95 15.1 15.95 14.6 17.6 15.35 12.9 12.6 10.35 15.4 9.6 14.85 13.6 12.65 11.9 16.6 11.2 15.85 11.15 15.6 13.1 12.4 14.9 11.2 14.6 14.75 7.85 7.85 10.95 9.95 14.9 13.4 16.85 10.95 12.2 15.2 8.1 4.35 12.7 18.1 17.85 17.1 16.1 14.7 10.6 12.6 16.2 13.6 14.1 14.5 14.75 14.8 12.45 12.65 17.35 8.6 16.1 17.75 15.25 17.65 13.6 18.25 16 18.25 18.95 16.1 15.4 15.4 13.35 19.1 19.25 12.75 9.85 15.25 12.4 18.15 12.35 15.6 18.4 12.85 9.5 4.5 13.6 11.7 17.75 14.05 13.35 11.85 13.2
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