Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1440 1920 1280 1280 1280 1366 1440 953 1280 1920 1366 1366 1280 1280 1280 1024 1366 1600 1024 1920 1280 1280 1024 1024 1311 1280 1280 1440 1229 1280 1600 1280 1366 1152 1280 1280 1280 1280 1280 1280 1366 1680 1280 1366 1280 1280 1366 1024 1680 1366 1280 1280 1366 1280 1024 1280 1280 1366 1440 1366 1280 1440 1920 1366 1280 1280 1440 1280 1680 1280 1366 1280 1280 1280 1366 1920 1920 1229 1280 1920 1280 1024 1280 1024 1280 1280 1280 1680 1360 1024 1920 1366 1366 1280 1600 1024 1366 1280 1024 1280 1366 1440 1200 1440 1680 1280 1400 1366 1280 1440 1280 1440 1440 1920 1920 1366 1366 1517 1680 1280 1280 1280 1280 1366 320 1920 1280 1440 1366 1280 1366 1366 1280 1536 1280 1280 1680 1280 1440 1280 1360 1920 1680 1366 1280 1280 1440 1280 1280 1680 320 1366 1280 1366 1024 1280 1366 1280 1366 1280 1117 1280
Data Y:
900 1080 800 1024 1024 768 900 762 800 1080 768 768 800 800 800 600 768 900 600 1080 800 800 600 768 737 800 800 900 737 800 900 800 768 720 800 800 800 960 800 800 768 1050 800 768 1024 1024 768 768 1050 768 1024 1024 768 800 768 800 1024 768 900 768 800 900 1080 768 800 800 900 800 1050 800 768 800 800 800 768 1080 1080 691 800 1080 800 768 800 768 1024 1024 800 1050 768 768 1080 768 768 800 900 768 768 800 600 800 768 900 1600 900 1050 800 1050 768 800 900 800 900 900 1080 1080 768 768 948 945 1024 800 1024 1024 768 480 1080 800 900 768 800 768 768 800 864 800 800 1050 800 900 800 768 1080 1050 768 800 800 900 800 1024 1050 480 768 800 768 768 1024 768 1024 768 800 698 800
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