Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
210.5719011 211.8997081 210.0843833 216.9714191 212.5325759 213.3266004 211.6320939 213.8495876 209.7836603 207.8165281 209.5503598 197.7604812 218.7340218 222.8200937 222.6256364 227.2831783 230.9728651 223.1882391 226.6908408 227.3001901 222.3150821 223.6556123 227.8932506 219.077975 237.9620698 239.3620698 240.1949375 247.6320939 250.0126723 241.7479499 247.8635157 248.1117083 241.5191789 244.6548894 252.0347448 237.994072 260.2839013 263.301395 261.6975884 263.8558533 268.7959015 249.7504089 270.4517565 270.1824554 266.7364317 273.7251544 269.6086247 263.3813488 283.9826964 285.5836603 285.0489138 291.8846243 289.7983114 285.6105526 292.9853472 298.2038049 291.9685765 288.9378776 284.5719503 273.9959999
Data Y:
-489.655367 -447.6900293 -537.3677823 -481.4102379 -503.7985479 -581.5953765 -494.8534187 -446.578747 -524.8000885 -579.7086071 -653.407384 -847.3550966 -531.7339354 -455.9861327 -562.2051181 -515.2549594 -509.8456254 -612.3288157 -470.5160829 -457.8351585 -551.5732196 -562.0428617 -516.3417287 -798.8804249 -494.0911623 -494.2911623 -549.4996808 -494.2534187 -504.6475737 -606.2817381 -521.1265498 -486.5573154 -630.0568175 -546.8751679 -544.5016292 -915.6313089 -528.6226531 -480.2479814 -598.5478913 -565.6654263 -566.6899392 -844.5683272 -486.7818283 -499.9872656 -576.1595813 -552.459989 -618.924919 -862.4963291 -504.7098302 -523.2000885 -624.5719965 -541.8903469 -612.415585 -655.7054357 -555.9580406 -526.4736273 -618.4004061 -605.1949688 -683.0163103 -952.4118256
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