Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
85 99 98 99 97 99 99 99 99 99 99 99 98 95 98 98 99 96 88 94 88 98 98 96 98 98 98 96 94 100 98 99 97 97 99 99 98 98 99 98 98 99 98 98 98 99 98 99 85 95 98 98 98 99 99 98 99 99 92 99 98 99 98 98 98 99 99 98 98 98 98 98 98 98 99 97 98 96 99 99 98 98 98 99 99 99 99 99 99 99 99 99 99 99 99 98 98 99 99 98 98 88 98 98 98 84 99 96 99 99 100 85 60 96 99 94 98 99 96 97 97 97 98 98 99 99 97 96 98 99 99 99 98 99 98 98 98 98 98 98 99 99 92 99 99 96 98 98 96 98 96 98 99 98 98 98 98 97 98 96 97 85 92 70 98 99 99 96 99 99 99 98 99 98 99 98 99 99 98 97 98 91 99 99 99 99 97 98 98 98 96 88 95 98 96 92 98 98 97 98 97 96 98 98 97 98 98 98 98 98 98 98 99 97 98 98 99 98 98 98 99 99 98 98 99 99 98 82 98 98 98 99 99 98 98 98 98 98 93 99 98 99 98 98 98 80 99 98 98 97 98 100 98 98 99 98 100 98 96 86 97 98 94 88 97 96 99 98 98 98 98 97 99 99 96 97 98 99 98 99 98 99 98 99 99 99 98 98 99 96 97 98 99 98 99 99 99 99 99 98 98 99 99 98 99 98 97 99 98 98 98 97 99 98 99 66 99 98 99 99 99 98 99 99 98 98 99 99 99 98 100 98 97 98 99 93 98 99 98 99 96 99 98 99 93 99 98 99 99 97 99 98 99 99 98 96 98 99 98 99 99 98 99 96 97 99 98 96 99 90 98 97 98 97 99 98 99 99 98 94 98 99 98 97 98 96 98 99 85 96 98 96 92 98 99 98 100 97 100 98 96 99 99 98 98 98 100 98 100 97 96 90 98 97 97
Data Y:
28 0 0 0 4 0 0 0 0 0 0 0 0 12 10 0 0 4 20 10 18 0 0 12 0 0 0 0 0 0 0 0 0 0 0 0 0 22 0 0 0 0 0 0 0 0 0 0 32 12 0 0 0 0 0 0 0 0 18 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 20 0 0 0 0 0 28 30 10 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 12 0 16 16 24 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 14 0 0 0 0 0 0 0 0 26 40 0 0 0 12 0 0 0 0 0 0 0 0 0 0 12 0 12 24 0 0 0 0 0 0 0 0 0 19 12 28 0 0 0 0 0 31 0 0 0 0 0 0 0 12 0 0 12 0 0 0 0 0 0 20 0 0 0 0 2 27 0 0 0 0 0 0 12 24 0 0 28 35 6 12 24 7 0 0 1 12 20 19 25 12 9 24 0 0 18 19 18 0 12 8 12 21 14 4 9 0 0 0 20 0 0 0 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 26 12 12 0 0 22 0 0 5 0 0 0 0 14 2 12 0 0 0 0 32 3 28 20 0 35 29 0 26 26 0 0 0 29 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 12 0 12 0 0 32 0 0 0 0 12 0 22 4 12 12 12 0 24 0 0 26 34 0 0 7 0 0 0 0 0 0 0 0 0 12 0 0 0 3 12 0 12 28 21 0 15 12
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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation