Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
24 47 36 47 55 66 57 52 67 54 46 49 28 50 52 45 64 63 39 42 54 36 49 60 33 46 48 52 35 40 63 32 45 51 48 38 59 53 50 30 61 61 47 51 44 54 53 58 49 53 48 40 55 55 52 50 37 47 51 45 57 28 53 60 38 44 40 48 48 62 58 47 24 41 47 37 54 53 55 57 58 58 35 43 44 45 61 53 41 38 41 55 51 54 39 48 54 59 47 44 55 48 59 48 51 65 46 44 38 58 50 19 54 60 50 50 48 32 69 55 52 53 41 53 56 55 13 52 46 56 27 53 45 52 31 38 41 44 49 51 40 49 51 43 56 38 42 55 55 44 49 51 42 57 43 48 59 55 50 30 54 49 48 59 43 49 56 49 50 49 29 44 57 57 49 54 60 58 62 43 54 49 64 55 54 64 45 59 50 57 65 60 55 46 61 43 55 47 38 59 48 51 55 71 54 60 42 49 52 61 52 52 54 38 64 46 59
Data Y:
52 62 40 48 40 77 61 44 66 49 40 53 59 58 53 58 49 51 58 48 53 39 47 51 52 57 37 46 48 33 28 45 56 45 53 47 39 62 59 55 50 45 59 39 41 53 44 61 41 54 36 47 54 54 57 47 32 66 41 54 54 40 60 53 25 58 50 56 48 15 43 52 49 46 56 62 48 50 35 50 67 48 42 48 58 31 60 43 50 57 30 44 65 49 70 48 54 52 61 61 66 53 56 44 50 73 52 53 53 48 45 44 50 31 50 43 64 56 46 41 58 46 52 37 50 46 51 54 47 35 57 44 49 45 57 58 63 46 52 53 57 27 59 52 51 45 53 38 60 51 45 39 52 39 48 38 57 47 62 57 45 55 43 53 51 28 41 49 50 42 40 41 42 46 49 52 54 50 61 52 49 39 61 41 63 62 46 47 65 44 37 66 73 38 49 64 56 42 36 45 45 47 59 48 53 50 49 65 47 60 46 50 55 48 43 66 65 52 38 44 61 41 38 64 59 54 36 34 48 55 44 46 53 20 49 35 62 41 36 54 59 54 54 52 39 50 59 53 43 55 56 45 42 42 53 51 38 51 61 57 55 62 23 52 48 55 43 44 70 31 45 49 66 49 46 50 48 58 62 55 49
Chart options
Title:
Label y-axis:
Label x-axis:
R Code
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,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') 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') qq.plot(x,main='QQplot of variable x') dev.off() bitmap(file='test3.png') qq.plot(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