Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
96 75 70 134 72 8 169 1 88 98 106 122 57 139 87 176 114 121 103 135 123 99 74 103 158 116 102 132 62 150 143 50 141 48 141 83 112 79 33 149 126 85 84 68 50 101 20 101 150 118 99 8 88 21 30 97 163 132 161 89 160 139 104 103 66 163 93 85 150 143 107 22 85 91 131 140 156 81 137 102 72 161 30 120 49 71 76 85 146 165 89 168 48 149 75 107 116 165 155 165 121 156 86 13 113 112 133 169 30 121 82 148 12 146 23 84 163 4 81 18 118 76 55 62 16 98 137 50 152 163 142 77 59 94 128 63 127 59 118 110 45 96 128 41 146 147 121 185 0 4 0 0 0 0 85 157 0 0 7 12 0 37 0 62
Data Y:
186099 113854 99776 106194 100792 47552 250931 6853 115466 110896 169351 94853 72591 101345 113713 165354 164263 135213 111669 134163 140303 150773 111848 102509 96785 116136 158376 153990 64057 230054 184531 114198 198299 33750 189723 100826 188355 104470 58391 164808 134097 80238 133252 54518 121850 79367 56968 106314 191889 104864 160792 15049 191179 25109 45824 129711 210012 194679 197680 81180 197765 214738 96252 124527 153242 145707 113963 134904 114268 94333 102204 23824 111563 91313 89770 100125 165278 181712 80906 75881 83963 175721 68580 136323 55792 25157 100922 118845 170492 81716 115750 105590 92795 82390 135599 127667 163073 211381 189944 226168 117495 195894 80684 19630 88634 139292 128602 135848 178377 106330 178303 116938 5841 106020 24610 74151 232241 6622 127097 13155 160501 91502 24469 88229 13983 80716 157384 122975 191469 231257 258287 122531 61394 86480 195791 18284 147581 72558 147341 114651 100187 130332 134218 10901 145758 75767 134969 169216 0 7953 0 0 0 0 105406 174586 0 0 4245 21509 7670 15673 0 75882
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