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