Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
25.64 27.97 27.62 23.31 29.07 29.58 28.63 29.92 32.68 31.54 32.43 26.54 25.85 27.60 25.71 25.38 28.57 27.64 25.36 25.90 26.29 21.74 19.20 19.32 19.82 20.36 24.31 25.97 25.61 24.67 25.59 26.09 28.37 27.34 24.46 27.46 30.23 32.33 29.87 24.87 25.48 27.28 28.24 29.58 26.95 29.08 28.76 29.59 30.70 30.52 32.67 33.19 37.13 35.54 37.75 41.84 42.94 49.14 44.61 40.22 44.23 45.85 53.38 53.26 51.80 55.30 57.81 63.96 63.77 59.15 56.12 57.42 63.52 61.71 63.01 68.18 72.03 69.75 74.41 74.33 64.24 60.03 59.44 62.50 55.04 58.34 61.92 67.65 67.68 70.30 75.26 71.44 76.36 81.71 92.60 90.60 92.23 94.09 102.79 109.65 124.05 132.69 135.81 116.07 101.42 75.73 55.48 43.80 45.29 44.01 47.48 51.07 57.84 69.04 65.61 72.87 68.41 73.25 77.43 75.28 77.33 74.31 79.70 85.47 77.98 75.69 75.20 77.21 77.85 83.53 85.99
Data Y:
104.8 105.2 105.6 105.8 106.1 106.5 106.71 106.68 107.41 107.15 107.5 107.22 107.11 107.57 107.81 108.75 109.43 109.62 109.54 109.53 109.84 109.67 109.79 109.56 110.22 110.4 110.69 110.72 110.89 110.58 110.94 110.91 111.22 111.09 111 111.06 111.55 112.32 112.64 112.36 112.04 112.37 112.59 112.89 113.22 112.85 113.06 112.99 113.32 113.74 113.91 114.52 114.96 114.91 115.3 115.44 115.52 116.08 115.94 115.56 115.88 116.66 117.41 117.68 117.85 118.21 118.92 119.03 119.17 118.95 118.92 118.9 118.92 119.44 119.40 119.98 120.43 120.41 120.82 120.97 120.63 120.38 120.68 120.84 120.90 121.56 121.57 122.12 121.97 121.96 122.48 122.33 122.44 123.08 124.23 124.58 125.08 125.98 126.90 127.19 128.33 129.04 129.72 128.92 129.13 128.90 128.13 127.85 127.98 128.42 127.68 127.95 127.85 127.61 127.53 127.92 127.59 127.65 127.98 128.19 128.77 129.31 129.80 130.24 130.76 130.75 130.81 130.89 131.30 131.49 131.65
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