Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
-50.74706514 -48.93785253 -47.15852374 -45.17543128 -46.94276659 -47.29735115 -42.25443623 -45.636011 -40.5706558 -34.94664303 -32.38276206 -32.7877051 -26.14063339 -22.09709979 -21.64764034 -16.98759233 -17.90735002 -15.42376615 -11.53425643 -13.32206639 -9.014273186 -4.253762759 -0.425986315 2.843751224 5.304276145 3.352359662 4.868775794 2.749804351 0.226872673 0.617183142 -0.16183631 0.436478416 2.548934314 6.596089819 10.3692672 9.720903334 12.32020087 12.14404607 12.94839943 13.76345412 16.31332119 16.81965467 15.68801627 17.19533256 18.78454744 24.92140923 27.36016736 27.69815372 30.48854798 30.43069115 30.29093572 28.97275052 27.2951215 25.22718201 26.17030799 25.74716525 29.830356 38.0397072 56.54216423 42.59868849
Data Y:
-7.384479833 -348.812296 656.3387963 168.4268193 163.9341261 -219.1238071 131.7351533 68.87952107 -499.4813227 88.460018 -3.676899151 -250.2346886 174.1014888 -231.4249098 544.4571966 124.3024887 245.2234366 98.23681742 372.6587033 36.0872455 -337.1878898 -72.72145156 -121.8431801 -137.6524521 19.966092 -447.5428034 769.4438158 243.4453361 -164.2602175 454.3943467 -69.44942508 405.8522526 -358.8682021 -83.29905889 50.36917358 -94.12503992 -307.9726466 -171.4972374 581.7501227 244.4978459 74.07234434 218.1732414 41.97296758 230.9245802 -379.8087485 -230.0475763 -14.66583292 -134.2125562 -277.3250263 -389.3945325 682.2209027 -261.6970045 148.0023313 158.047103 -258.1716876 153.2005099 -667.3263955 -162.5924918 -58.96570282 -588.4092144
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