Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
-1756.606034 4515.24258 6350.757752 187.6662664 7179.987129 -1766.376102 3341.770864 -6153.4031 20418.81063 5825.391929 -10079.64907 3756.072201 -1726.742767 7013.734425 10008.75491 -3889.520808 2853.072056 4867.000198 9517.403357 -19081.5299 2090.060331 7667.210345 18234.87677 4200.880538 1264.162256 2584.899596 1795.964144 6713.258676 -0.068849117 7108.528605 5387.158619 -6468.095299 2883.307014 -4432.345942 -1782.780676 6472.577676 1013.79339 2177.704119 2867.612346 -8285.22375 1348.880462 3408.921922 -8620.730165 2859.197602 4446.042697 9696.016672 1717.131294 -3150.034481 -11543.52083 884.4316352 1621.208212 -2631.418695 1761.402548 -2649.643223 -2776.023462 -7699.097349 -533.0422272 -11156.45315 595.2518785 -619.038755 -4126.313963 -1418.221331 -4510.069121
Data Y:
9.251916481 22.28951742 -37.15509397 34.22187354 -175.779264 104.9707065 -61.23944789 -396.7006521 -396.4104855 -259.9810144 127.9713638 -36.01972185 204.1305246 44.23781198 -159.6936733 203.1240259 121.7477886 -16.20225966 -243.8846 391.6693185 -164.2219687 293.7998866 150.1830112 -191.0192529 251.8229631 -153.3772568 8.18066887 276.1462458 7.505698792 290.2551301 154.6986628 -344.4860962 -268.7252068 308.7804581 53.41379783 309.7151277 86.25318818 -34.86211092 196.5575758 226.4231797 -187.9949744 -58.58222587 209.148432 145.9303467 129.4478425 -33.08964353 -386.0328906 144.2386434 117.0703219 129.9208501 -152.2074543 -55.09481442 144.0193516 -16.89369247 144.6275477 -160.1503109 329.6155122 138.2292641 9.036902416 50.52401291 313.8066156 236.8949084 110.0260622
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