Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
112.52 112.39 112.24 112.10 109.85 111.89 111.88 111.48 110.98 110.42 107.90 109.46 109.11 109.26 109.99 110.17 110.28 109.13 110.15 109.39 108.45 108.23 107.44 104.86 106.23 105.85 104.95 104.46 104.66 103.05 104.16 104.08 104.20 103.68 103.69 101.29 103.03 102.90 102.68 102.98 103.47 101.72 102.82 102.74 102.38 101.81 101.88 99.60 100.93 100.85 100.93 101.10 101.10 99.31 100.33 99.99 99.82 99.65 99.06 96.92 98.20 98.54 98.71 98.20 98.29 96.67 97.69 97.78 97.44 96.92 96.84 95.05 96.33 96.33 96.16 96.50 96.33 94.71 95.82 95.47 95.82 95.99 95.73 93.77 94.71 94.62 94.79 94.88 94.79 93.43 94.37 94.62 94.45 94.37 94.20 92.66 93.51 93.60 93.60 93.77 93.60 92.41 93.60 93.34 92.92 92.07 91.89 90.27 91.72 91.98 91.81 91.98 91.30 89.93 90.87 90.53 90.27 90.10 89.68 87.89
Data Y:
108.06 108.22 107.96 107.16 107.01 107.64 108.35 108.45 108.08 106.51 105.52 105.82 105.07 105.88 106.32 106.49 105.96 105.82 105.94 106.25 105.81 105.44 104.20 103.32 103.35 103.92 103.76 103.60 102.89 102.77 103.55 104.09 104.05 103.47 102.16 101.50 101.75 102.10 102.12 102.30 101.81 101.32 101.68 102.02 102.06 101.58 100.78 100.28 100.06 100.46 100.78 100.78 99.90 99.58 99.90 100.22 100.30 100.14 99.34 98.53 98.05 98.85 99.26 99.10 98.29 98.13 98.45 99.18 99.01 98.69 97.89 97.33 96.84 97.41 97.81 97.97 97.16 96.92 97.00 97.49 97.57 97.49 96.60 95.88 95.31 95.48 95.96 96.04 95.07 94.91 94.99 95.31 95.48 94.83 93.87 93.38 92.34 92.50 92.74 92.66 91.70 91.45 91.53 91.86 91.61 90.89 89.85 89.04 87.83 88.24 88.32 88.00 87.19 86.87 87.11 87.11 86.87 86.55 85.58 85.10
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