Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1280 1024 1120 1024 1280 1280 1280 1024 1280 1280 1280 1280 1280 1688 1440 1600 1280 1280 1280 1176 1280 1503 1440 1366 1280 1024 1280 2560 1280 1024 1280 1280 1440 1280 1440 1024 1440 1143 1280 1440 1280 1366 1024 1408 1366 1176 1920 1257 1280 1280 1440 1680 1440 1024 1140 1280 1280 1280 1280 1280 1440 1280 1152 1280 1280 1440 1280 1280 1440 1280 1280 1440 1280 1280 1600 1024 1366 1280 1280 1440 1366 1280 1024 1280 1440 1280 1280 1408 1280 1600 1600 1680 1440 1440 917 1280 1760 1280 1280 1280 1024 1366 1440 1280 1280 1920 1024 1024 1600 1117 1440 983 1024 1024 1280 1440 1280 1280 1280 1440 1280 1024 1024 1152 1280 1024 1366 1680 1680 1280 1366 1024 1440 1024 1280 1280 1280 1024 1280
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
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation