Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
4.000 50.000 102.909 40.500 27.500 31.500 60.500 47.500 2.694.420 1.675.000 81.564 490.531 44.500 935.321 1.283.190 1.453.640 1.453.640 169.769 491.182 212.400 1.066.550 5.123.210 107.589 50.000 2.515.000 134.926 77.799 638.394 2.539.000 330.314 1.251.320 846.000 20.000 1.100.000 648.000 638.500 195.000 667.000 692.864 1.105.000 226.641 1.708.000 594.500 679.000 667.047 932.159 239.000 26.304 934.000 390.453 428.000 140.000 1.296.520 85.500 538.000 6.305.000 399.050 1.091.950 248.000 182.200 6.556.000 771.000 4.450.000 1.120.820 1.934.000 4.000.000 707.001 1.031.210 1.298.750 1.945.000 163.000 233.500 204.000 1.396.000 475.882 359.900 612.264 92.300 1.019.000 97.250 366.000 259.538 473.000 300.000 1.797.000 1.631.000 981.000 1.480.630 4.146.000 60.250 167.987 198.604 1.025.130 222.000 879.166 1.900.000 1.460.000 914.874 2.885.000 1.235.770 4.207.960 711.250 2.347.000 1.811.350 202.532 872.074 125.520 785.000 82.500 5.200.000 2.325.000 420.385 66.000 2.885.000 1.153.000 1.499.340 431.492 350.000 3.424.470 1.413.670 5.013.090 91.260 567.000 6.950.420 278.900 899.000 1.038.240 2.146.550 381.157 2.202.810 2.862.000 214.733 720.482 402.500 849.330 596.833 210.000 1.083.750 187.940 1.453.000 1.212.700 800.430 318.032 2.100.000 399.217 40.000 679.000 4.000.000 1.810.000 2.065.890 2.615.000 352.000 146.000 176.000 3.447.010 115.000 1.026.600 708.000 226.054 1.222.000 434.724 254.832 868.964 875.000 1.313.000 1.215.000 2.187.000 1.766.000 1.158.000 13.089.000 1.754.420 981.469 69.400 956.550 3.034.000 59.000 116.807 683.000 967.000 782.400 862.500 926.247 1.595.000 1.153.470 430.000 2.025.420 6.774.000 2.467.400 1.328.000 36.600 93.000 279.200 544.847 178.500 2.382.290 840.000 342.000 853.000 113.000 65.000 160.797 434.000
Data Y:
715 955.646 1.375.649 500.000 4.706.424 1.043.000 893.367 1.015.236 105.568.000 678.367.000 2.761.149 54.839.000 3.117.144 324.610.000 1.088.642.000 182.023.000 182.023.000 5.829.596 48.692.339 10.254.799 225.438.000 1.668.100.000 47.133.850 7.833.000 374.800.000 13.641.221 38.806.307 62.463.000 714.151.000 25.526.452 305.288.000 27.322.000 1.600.000 1.437.977.000 27.601.000 827.633.000 22.011.000 174.407.000 118.105.000 449.764.000 24.153.310 263.312.000 214.452.000 635.487.000 66.179.000 143.000.000 80.139.887 6.652.015 200.700.000 59.913.000 75.543.000 14.456.000 655.968.000 3.583.600 84.020.000 5.317.858.000 52.956.000 336.389.000 48.889.000 37.514.000 6.273.516.000 128.819.000 1.401.307.000 555.869.000 912.508.000 11.459.000.000 83.814.000 307.002.000 67.885.000 233.388.000 11.525.000 47.915.000 172.660.729 429.686.000 116.171.000 21.856.000 147.470.000 10.200.000 199.795.000 29.471.835 155.002.000 101.186.000 80.944.000 41.129.205 531.872.000 657.238.000 147.144.000 407.657.000 1.392.842.000 6.437.036 28.041.000 3.154.230 663.656.000 26.921.000 552.456.000 1.361.597.000 666.931.000 350.201.000 374.145.000 248.815.000 3.123.897.000 251.355.000 1.762.400.000 811.130.000 25.136.000 161.930.000 5.827.041 296.557.000 9.548.849 2.719.546.000 3.566.500.000 107.946.000 34.700.000 2.967.930.000 458.963.000 252.832.000 224.211.000 46.177.000 3.350.727.000 219.372.000 397.754.000 38.055.000 97.511.941 3.081.775.000 64.208.000 1.021.898.000 1.069.178.000 543.058.000 49.200.000 2.696.895.000 5.127.242.000 77.755.000 120.473.000 207.004.000 312.636.000 111.597.000 51.454.000 1.404.567.000 87.131.376 474.828.000 680.041.000 246.314.000 57.389.000 2.414.300.000 33.726.952 2.833.345 115.427.000 9.078.000.000 1.165.781.000 1.465.720.000 1.963.300.000 84.080.000 16.598.200 51.748.825 1.872.529.000 18.208.000 1.022.672.000 275.386.000 62.155.250 915.349.000 143.729.000 51.700.000 190.304.000 245.958.000 344.567.000 1.352.591.000 3.369.407.000 330.947.000 241.821.000 53.095.000.000 525.930.000 420.795.000 30.594.342 742.838.000 3.219.184.000 9.765.000 3.026.000 92.673.000 2.405.711.000 106.442.000 183.224.468 327.579.000 2.293.231.000 616.318.000 56.078.000 737.423.000 12.119.000.000 2.811.901.000 780.290.000 5.482.439 34.469.000 5.851.600 292.693.000 27.777.678 3.349.512.000 1.421.529.000 144.415.000 325.562.000 42.566.440 200.000 4.477.670 13.980.000
Chart options
Title:
Label y-axis:
Label x-axis:
R Code
library(psychometric) x <- x[!is.na(y)] y <- y[!is.na(y)] y <- y[!is.na(x)] x <- x[!is.na(x)] 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, sub=main) 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', na.rm = T) 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,'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,'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,'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,'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,'Correlation',header=TRUE) a<-table.element(a,cxy,2) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Determination',header=TRUE) a<-table.element(a,cxy*cxy,2) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'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,'95% CI of Correlation',header=TRUE) a<-table.element(a,paste('[',CIr(r=cxy, n = lx, level = .95)[1],', ', CIr(r=cxy, n = lx, level = .95)[2],']',sep=''),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') library(moments) library(nortest) jarque.x <- jarque.test(x) jarque.y <- jarque.test(y) if(lx>7) { ad.x <- ad.test(x) ad.y <- ad.test(y) } a<-table.start() a<-table.row.start(a) a<-table.element(a,'Normality Tests',1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,paste('<pre>',RC.texteval('jarque.x'),'</pre>',sep='')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,paste('<pre>',RC.texteval('jarque.y'),'</pre>',sep='')) a<-table.row.end(a) if(lx>7) { a<-table.row.start(a) a<-table.element(a,paste('<pre>',RC.texteval('ad.x'),'</pre>',sep='')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,paste('<pre>',RC.texteval('ad.y'),'</pre>',sep='')) a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable1.tab') library(car) bitmap(file='test2.png') qqPlot(x,main='QQplot of variable x') dev.off() bitmap(file='test3.png') qqPlot(y,main='QQplot of variable y') dev.off()
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