Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1280 1024 1024 768 1120 700 1024 768 1280 800 1280 1024 1280 800 1024 768 1280 800 1280 1024 1280 800 1280 800 1280 1024 1688 949 1440 900 1600 1200 1280 800 1280 800 1280 768 1176 735 1280 800 1503 845 1440 900 1366 768 1280 768 1024 768 1280 800 2560 1440 1280 768 1024 768 1280 1024 1280 800 1440 900 1280 800 1440 900 1024 768 1440 900 1143 857 1280 800 1440 900 1280 800 1366 768 1024 768 1408 880 1366 768 1176 735 1920 1200 1257 785 1280 800 1280 800 1440 900 1680 1050 1440 900 1024 768 1140 641 1280 1024 1280 800 1280 800 1280 800 1280 800 1440 900 1280 800 1152 864 1280 1024 1280 800 1440 900 1280 800 1280 1024 1440 900 1280 800 1280 800 1440 900 1280 800 1280 1024 1600 900 1024 768 1366 768 1280 800 1280 800 1440 900 1366 768 1280 800 1024 768 1280 800 1440 900 1280 800 1280 800 1408 880 1280 800 1600 900 1600 900 1680 1050 1440 900 1440 900 917 550 1280 800 1760 990 1280 800 1280 800 1280 800 1024 768 1366 768 1440 900 1280 800 1280 1024 1920 1080 1024 768 1024 768 1600 900 1117 698 1440 900 983 737 1024 768 1024 640 1280 800 1440 900 1280 800 1280 800 1280 800 1440 900 1280 800 1024 768 1024 768 1152 864 1280 768 1024 768 1366 768 1680 1050 1680 1050 1280 800 1366 768 1024 768 1440 900 1024 768 1280 800 1280 800 1280 800 1024 768 1280 800
Data Y:
1024 768 700 768 800 1024 800 768 800 1024 800 800 1024 949 900 1200 800 800 768 735 800 845 900 768 768 768 800 1440 768 768 1024 800 900 800 900 768 900 857 800 900 800 768 768 880 768 735 1200 785 800 800 900 1050 900 768 641 1024 800 800 800 800 900 800 864 1024 800 900 800 1024 900 800 800 900 800 1024 900 768 768 800 800 900 768 800 768 800 900 800 800 880 800 900 900 1050 900 900 550 800 990 800 800 800 768 768 900 800 1024 1080 768 768 900 698 900 737 768 640 800 900 800 800 800 900 800 768 768 864 768 768 768 1050 1050 800 768 768 900 768 800 800 800 768 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