Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
235.1 37 280.7 30 264.6 47 240.7 35 201.4 30 240.8 43 241.1 82 223.8 40 206.1 47 174.7 19 203.3 52 220.5 136 299.5 80 347.4 42 338.3 54 327.7 66 351.6 81 396.6 63 438.8 137 395.6 72 363.5 107 378.8 58 357 36 369 52 464.8 79 479.1 77 431.3 54 366.5 84 326.3 48 355.1 96 331.6 83 261.3 66 249 61 205.5 53 235.6 30 240.9 74 264.9 69 253.8 59 232.3 42 193.8 65 177 70 213.2 100 207.2 63 180.6 105 188.6 82 175.4 81 199 75 179.6 102 225.8 121 234 98 200.2 76 183.6 77 178.2 63 203.2 37 208.5 35 191.8 23 172.8 40 148 29 159.4 37 154.5 51 213.2 20 196.4 28 182.8 13 176.4 22 153.6 25 173.2 13 171 16 151.2 13 161.9 16 157.2 17 201.7 9 236.4 17 356.1 25 398.3 14 403.7 8 384.6 7 365.8 10 368.1 7 367.9 10 347 3 343.3 292.9 311.5 300.9 366.9 356.9 329.7 316.2 269 289.3 266.2 253.6 233.8 228.4 253.6 260.1 306.6 309.2 309.5 271 279.9 317.9 298.4 246.7 227.3 209.1 259.9 266 320.6 308.5 282.2 262.7 263.5 313.1 284.3 252.6 250.3 246.5 312.7 333.2 446.4 511.6 515.5 506.4 483.2 522.3 509.8 460.7 405.8 375 378.5 406.8 467.8 469.8 429.8 355.8 332.7 378 360.5 334.7 319.5 323.1 363.6 352.1 411.9 388.6 416.4 360.7 338 417.2 388.4 371.1 331.5 353.7 396.7 447 533.5 565.4 542.3 488.7 467.1 531.3 496.1 444 403.4 386.3 394.1 404.1 462.1 448.1 432.3 386.3 395.2 421.9 382.9 384.2 345.5 323.4 372.6 376 462.7 487 444.2 399.3 394.9 455.4 414 375.5 347 339.4 385.8 378.8 451.8 446.1 422.5 383.1 352.8 445.3 367.5 355.1 326.2 319.8 331.8 340.9 394.1 417.2 369.9 349.2 321.4 405.7 342.9 316.5 284.2 270.9 288.8 278.8 324.4 310.9 299 273 279.3 359.2 305 282.1 250.3 246.5 257.9 266.5 315.9 318.4 295.4 266.4 245.8 362.8 324.9 294.2 289.5 295.2 290.3 272 307.4 328.7 292.9 249.1 230.4 361.5 321.7 277.2 260.7 251 257.6 241.8 287.5 292.3 274.7 254.2 230 339 318.2 287 295.8 284 271 262.7 340.6 379.4 373.3 355.2 338.4 466.9 451 422 429.2 425.9 460.7 463.6 541.4 544.2 517.5 469.4 439.4 549 533 506.1 484 457 481.5 469.5 544.7 541.2 521.5 469.7 434.4 542.6 517.3 485.7 465.8 447 426.6 411.6 467.5 484.5 451.2 417.4 379.9 484.7 455 420.8 416.5 376.3 405.6 405.8 500.8 514 475.5 430.1 414.4 538 526 488.5 520.2 504.4 568.5 610.6 818 830.9 835.9 782 762.3 856.9 820.9 769.6 752.2 724.4 723.1 719.5 817.4 803.3 752.5 689 630.4 765.5 757.7 732.2 702.6 683.3 709.5 702.2 784.8 810.9 755.6 656.8 615.1 745.3 694.1 675.7 643.7 622.1 634.6 588 689.7 673.9 647.9 568.8 545.7 632.6 643.8 593.1 579.7 546 562.9 572.5
Names of X columns:
Werkloosheid Militairen
Column number of first sample
Column number of second sample
Confidence
Alternative
two.sided
two.sided
less
greater
Are observations paired?
unpaired
unpaired
paired
Null Hypothesis
Chart options
Title:
R Code
par1 <- as.numeric(par1) #column number of first sample par2 <- as.numeric(par2) #column number of second sample par3 <- as.numeric(par3) #confidence (= 1 - alpha) if (par5 == 'unpaired') paired <- FALSE else paired <- TRUE par6 <- as.numeric(par6) #H0 z <- t(y) if (par1 == par2) stop('Please, select two different column numbers') if (par1 < 1) stop('Please, select a column number greater than zero for the first sample') if (par2 < 1) stop('Please, select a column number greater than zero for the second sample') if (par1 > length(z[1,])) stop('The column number for the first sample should be smaller') if (par2 > length(z[1,])) stop('The column number for the second sample should be smaller') if (par3 <= 0) stop('The confidence level should be larger than zero') if (par3 >= 1) stop('The confidence level should be smaller than zero') (r.t <- t.test(z[,par1],z[,par2],var.equal=TRUE,alternative=par4,paired=paired,mu=par6,conf.level=par3)) (v.t <- var.test(z[,par1],z[,par2],conf.level=par3)) (r.w <- t.test(z[,par1],z[,par2],var.equal=FALSE,alternative=par4,paired=paired,mu=par6,conf.level=par3)) (w.t <- wilcox.test(z[,par1],z[,par2],alternative=par4,paired=paired,mu=par6,conf.level=par3)) (ks.t <- ks.test(z[,par1],z[,par2],alternative=par4)) m1 <- mean(z[,par1],na.rm=T) m2 <- mean(z[,par2],na.rm=T) mdiff <- m1 - m2 newsam1 <- z[!is.na(z[,par1]),par1] newsam2 <- z[,par2]+mdiff newsam2 <- newsam2[!is.na(newsam2)] (ks1.t <- ks.test(newsam1,newsam2,alternative=par4)) mydf <- data.frame(cbind(z[,par1],z[,par2])) colnames(mydf) <- c('Variable 1','Variable 2') bitmap(file='test1.png') boxplot(mydf, notch=TRUE, ylab='value',main=main) dev.off() bitmap(file='test2.png') qqnorm(z[,par1],main='Normal QQplot - Variable 1') qqline(z[,par1]) dev.off() bitmap(file='test3.png') qqnorm(z[,par2],main='Normal QQplot - Variable 2') qqline(z[,par2]) dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,paste('Two Sample t-test (',par5,')',sep=''),2,TRUE) a<-table.row.end(a) if(!paired){ a<-table.row.start(a) a<-table.element(a,'Mean of Sample 1',header=TRUE) a<-table.element(a,r.t$estimate[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Mean of Sample 2',header=TRUE) a<-table.element(a,r.t$estimate[[2]]) a<-table.row.end(a) } else { a<-table.row.start(a) a<-table.element(a,'Difference: Mean1 - Mean2',header=TRUE) a<-table.element(a,r.t$estimate) a<-table.row.end(a) } a<-table.row.start(a) a<-table.element(a,'t-stat',header=TRUE) a<-table.element(a,r.t$statistic[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'df',header=TRUE) a<-table.element(a,r.t$parameter[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value',header=TRUE) a<-table.element(a,r.t$p.value) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'H0 value',header=TRUE) a<-table.element(a,r.t$null.value[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Alternative',header=TRUE) a<-table.element(a,r.t$alternative) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'CI Level',header=TRUE) a<-table.element(a,attr(r.t$conf.int,'conf.level')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'CI',header=TRUE) a<-table.element(a,paste('[',r.t$conf.int[1],',',r.t$conf.int[2],']',sep='')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'F-test to compare two variances',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'F-stat',header=TRUE) a<-table.element(a,v.t$statistic[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'df',header=TRUE) a<-table.element(a,v.t$parameter[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value',header=TRUE) a<-table.element(a,v.t$p.value) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'H0 value',header=TRUE) a<-table.element(a,v.t$null.value[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Alternative',header=TRUE) a<-table.element(a,v.t$alternative) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'CI Level',header=TRUE) a<-table.element(a,attr(v.t$conf.int,'conf.level')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'CI',header=TRUE) a<-table.element(a,paste('[',v.t$conf.int[1],',',v.t$conf.int[2],']',sep='')) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab') a<-table.start() a<-table.row.start(a) a<-table.element(a,paste('Welch Two Sample t-test (',par5,')',sep=''),2,TRUE) a<-table.row.end(a) if(!paired){ a<-table.row.start(a) a<-table.element(a,'Mean of Sample 1',header=TRUE) a<-table.element(a,r.w$estimate[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Mean of Sample 2',header=TRUE) a<-table.element(a,r.w$estimate[[2]]) a<-table.row.end(a) } else { a<-table.row.start(a) a<-table.element(a,'Difference: Mean1 - Mean2',header=TRUE) a<-table.element(a,r.w$estimate) a<-table.row.end(a) } a<-table.row.start(a) a<-table.element(a,'t-stat',header=TRUE) a<-table.element(a,r.w$statistic[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'df',header=TRUE) a<-table.element(a,r.w$parameter[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value',header=TRUE) a<-table.element(a,r.w$p.value) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'H0 value',header=TRUE) a<-table.element(a,r.w$null.value[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Alternative',header=TRUE) a<-table.element(a,r.w$alternative) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'CI Level',header=TRUE) a<-table.element(a,attr(r.w$conf.int,'conf.level')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'CI',header=TRUE) a<-table.element(a,paste('[',r.w$conf.int[1],',',r.w$conf.int[2],']',sep='')) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable1.tab') a<-table.start() a<-table.row.start(a) a<-table.element(a,paste('Wicoxon rank sum test with continuity correction (',par5,')',sep=''),2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'W',header=TRUE) a<-table.element(a,w.t$statistic[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value',header=TRUE) a<-table.element(a,w.t$p.value) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'H0 value',header=TRUE) a<-table.element(a,w.t$null.value[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Alternative',header=TRUE) a<-table.element(a,w.t$alternative) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Kolmogorov-Smirnov Test to compare <i>Distributions</i> of two Samples',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'KS Statistic',header=TRUE) a<-table.element(a,ks.t$statistic[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value',header=TRUE) a<-table.element(a,ks.t$p.value) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Kolmogorov-Smirnov Test to compare <i>Distributional Shape</i> of two Samples',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'KS Statistic',header=TRUE) a<-table.element(a,ks1.t$statistic[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value',header=TRUE) a<-table.element(a,ks1.t$p.value) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable2.tab')
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