Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
-0.317563243 0.793263584 0.245684805 -0.027819099 0.072180901 0.414346267 0.893263584 0.366767488 0.245684805 0.666767488 1.004090411 1.204090411 1.283007729 0.356511632 0.330015536 0.956511632 1.01975968 0.998676998 0.314346267 0.651098219 0.893263584 0.572180901 0.788421142 0.688421142 0.183007729 -0.39648056 -0.307307387 -0.22839007 -0.586224705 -0.612720801 -0.844630311 -0.907878359 -1.57653982 -1.787366647 -1.26029958 -1.297051532 -1.434374455 -0.760870551 -0.229532013 0.075881401 0.539129449 0.512633353 1.143971891 0.822889208 1.569897016 2.375310429 2.960212132 2.460212132 2.307219939 1.522889208 -0.077110792 -0.324118599 -0.713291772 -0.807878359 -2.023547628 -1.881382262 -2.907878359 -3.518134214 -3.733232512 -2.669984464
Data Y:
33.97260143 16.58856422 10.35739026 3.126967938 -1.873032062 -5.949839504 -10.41143578 -12.18101346 -16.64260974 -16.18101346 11.20452701 19.20452701 19.74293073 12.97335305 6.20377538 -2.026646945 -10.64185811 -12.10345439 -15.9498395 -17.33462834 -21.41143578 -19.87303206 4.050912124 10.05091212 11.74293073 6.434197706 0.818234915 -5.643361364 -5.720168806 -8.489746481 -8.56730555 -9.952094388 -12.64486462 -9.260827411 8.279079566 12.6638684 15.27832794 9.508750263 4.815980031 1.123961426 -0.491249737 -1.260827411 -2.953597644 -4.415193923 -6.954349272 -6.646367876 9.970346542 13.97034654 17.43119119 11.58480608 7.584806077 3.123961426 -0.260075783 -2.952094388 -6.105709271 -8.182516713 -12.95209439 -13.79772788 1.818986543 6.20377538
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