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