Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
5 4 5 6 6 6 7 8 7 8 7 8 8 9 9 8 9 9 10 11 12 13 13 13 14 14 15 15 16 16 17 18 19 20 22 20 22 25 24 25 28 26 27 26 25 27 28 30 31 32 34 34 33 32 34 36 37 40 38 38 36 40 40 42 44 45 47 49 47 49 52 50 50 57 58 58 58 61 61 64 68 40 34 46 36 34 45 55 50 56 72 76 78 77 90 88 97 93 84 67 72 75 71 75 90 78 73 62 65 61 58 33 39 56 79 82 79 73 87 85 83 82 83 92 95 97 87 84 84 89 103 106 109 106 105 115 120 124 121 131 139 133 119 123 120 128 134 126 115 106 99 100 99 99 100 100 108 109 115 114 108 113 118 122 118 121 118 121 121 112 119 116 110 111 106 108
Data Y:
8.7 10.6 9.0 9.2 8.3 7.6 9.0 8.2 9.0 8.7 9.1 8.0 7.7 9.9 8.4 9.0 8.8 8.3 8.8 9.8 8.1 8.9 7.4 8.9 9.8 8.5 9.7 7.7 8.8 9.6 9.6 7.8 9.4 9.4 8.8 10.1 9.0 8.0 7.8 10.0 8.9 8.9 8.6 9.0 9.3 8.9 7.0 9.3 8.4 9.0 8.9 9.4 8.3 8.8 7.5 7.5 8.2 8.1 8.2 8.4 9.3 9.2 8.3 8.6 9.2 9.5 9.2 9.4 8.8 8.4 9.3 9.1 8.9 9.3 8.9 8.8 8.6 9.7 10.4 9.7 10.0 9.8 9.3 9.5 8.4 9.8 8.6 9.7 10.5 8.6 9.3 9.1 9.1 9.9 9.3 9.7 8.9 10.0 9.1 9.7 9.1 10.3 9.7 9.6 9.8 9.8 9.8 8.7 9.1 8.9 10.1 9.5 10.3 9.3 10.2 10.1 10.3 9.8 9.8 9.3 9.9 9.3 9.2 8.4 10.1 9.6 10.6 10.0 10.4 8.6 8.4 9.6 9.1 10.0 10.0 9.3 9.6 9.6 9.9 9.3 9.7 10.0 9.9 10.4 9.8 9.4 9.0 9.4 9.7 10.5 10.5 9.9 8.9 9.4 9.2 10.6 11.3 11.2 10.0 10.6 10.1 11.1 10.9 9.1 10.8 10.6 11.1 11.2 10.7 11.2 11.1 10.7 11.0 11.4 11.5 10.9
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