Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
3 0 5 7 4 3 4 0 5 14 0 0 0 4 7 52 42 47 42 77 7 0 8 36 191 11 32 20 87 66 54 54 74 41 38 41 25 30 29 53 31 63 55 36
Data Y:
5 0 59 24 21 64 58 41 162 32 20 22 16 75 50 116 447 455 447 802 34 22 63 549 1396 603 387 199 495 608 883 344 903 412 908 766 779 904 440 648 360 259 448 255
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