Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
13768040.14 17487530.67 16198106.13 17535166.38 16571771.60 16198892.67 16554237.93 19554176.37 15903762.33 18003781.65 18329610.38 16260733.42 14851949.20 18174068.44 18406552.23 18466459.42 16016524.60 17428458.32 17167191.42 19629987.60 17183629.01 18344657.85 19301440.71 18147463.68 16192909.22 18374420.60 20515191.95 18957217.20 16471529.53 18746813.27 19009453.59 19211178.55 20547653.75 19325754.03 20605542.58 20056915.06 16141449.72 20359793.22 19711553.27 15638580.70 14384486.00 13855616.12 14308336.46 15290621.44 14423755.53 13779681.49 15686348.94 14733828.17 12522497.94 16189383.57 16059123.25 16007123.26 15806842.33 15159951.13 15692144.17 18908869.11 16969881.42 16997477.78 19858875.65 17681170.13
Data Y:
14731798.37 16471559.62 15213975.95 17637387.4 17972385.83 16896235.55 16697955.94 19691579.52 15930700.75 17444615.98 17699369.88 15189796.81 15672722.75 17180794.3 17664893.45 17862884.98 16162288.88 17463628.82 16772112.17 19106861.48 16721314.25 18161267.85 18509941.2 17802737.97 16409869.75 17967742.04 20286602.27 19537280.81 18021889.62 20194317.23 19049596.62 20244720.94 21473302.24 19673603.19 21053177.29 20159479.84 18203628.31 21289464.94 20432335.71 17180395.07 15816786.32 15071819.75 14521120.61 15668789.39 14346884.11 13881008.13 15465943.69 14238232.92 13557713.21 16127590.29 16793894.2 16014007.43 16867867.15 16014583.21 15878594.85 18664899.14 17962530.06 17332692.2 19542066.35 17203555.19
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