Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
221 219 214 210 207 206 217 231 234 233 228 226 227 225 219 215 210 206 215 228 229 222 215 212 211 208 205 201 198 198 210 224 226 222 216 215 215 214 211 207 203 200 209 223 225 216 206 203 203 201 197 192 187 184 194 203 197 191 182 175 163 155 151 156 154 153 167 177 171 169 160 151 139 130 126 130 127 122 129 135 142 156 157 165 170 169 162 148 143 146 175 181 178 166 161 164 173 174 167 156 148 150 174 181 183 178 176 184 193 192 182 163 157 167 205 219 214 198 183 184 192 196 194 185 181 184 206 210 208 197 189 190 191 190 187 184 183 184 203 208 205 195 189 188 190 190 190 193 185 173 176 170 163 170 171 173 171 162 152 142 136 146 179 191 181 170 161 168 180 182 176 164 154 160 189 196 186 171 169 181 198 202 196 183 173 175 198 203 197 191 182 172 158 147 143 146 147 152 177 184 174 162 157 155 159 158 156 157 156 158 173 179 172 169 168 172 180 182 182 181 178 178 196 199 192 187 184 184 188 183 176 168 163 166 189 195 192 189 187 187 190 187 179 168 160 161 177 182 176
Data Y:
191 189 184 179 175 171 179 191 195 195 193 193 195 193 187 181 176 169 174 185 186 182 178 178 179 178 174 171 168 167 175 187 191 188 185 185 187 188 186 183 179 176 183 198 203 198 192 191 194 194 192 188 182 175 178 181 171 164 159 160 163 159 148 139 129 124 136 146 143 141 135 134 135 134 136 142 142 135 140 146 155 170 167 166 160 156 156 160 156 150 157 158 167 189 197 199 193 188 186 190 186 181 190 189 192 201 200 206 208 202 190 171 163 167 195 208 208 197 189 192 199 202 200 191 190 180 194 196 199 200 199 205 207 211 210 208 201 186 177 168 173 181 185 186 189 186 181 182 176 165 176 174 168 165 162 170 179 178 169 160 151 159 191 195 184 162 152 162 188 202 209 204 193 191 202 204 206 211 214 224 224 222 219 218 213 213 229 225 220 212 204 204 202 195 186 175 170 171 196 202 200 191 186 186 193 193 188 185 182 180 194 204 216 233 241 243 241 233 228 225 219 217 235 237 238 235 234 239 248 248 247 246 240 233 242 239 238 238 238 240 249 251 253 251 246 247 260 260 259
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,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,'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') 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