Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
0.18428321652041 0.00710294820826641 0.00116535686783743 0.108729652178155 0.0226027559903587 0.124926004055343 0.101432440208012 0.175027945010452 0.062637379076165 -0.0399529666703403 0.12319625970099 0.0520964104695329 0.0741331631496151 0.0786957318355064 0.174609618553996 0.244943379543757 0.059233510470993 -0.190572109424889 0.168110836840757 0.0619037066351169 0.0028916512338603 -0.0455102815583758 -0.116673237656353 0.0585961274231594 -0.0259904869472595 -0.38609466229135 -0.00630486571274935 0.108071358409122 -0.0628380702948525 -0.108938212257777 -0.222915326468656 0.114551148378627 -0.05663840395505 0.21268911402888 -0.163736199827149 -0.0970470642759202 -0.344081807315471 0.120216028668961 0.00553481619705711 -0.0695365440812045 -0.137723372960688 -0.0331822099960538 -0.0519876598333341 0.0428565265578385 -0.150658573328081 0.142080265356049 0.171890801091025 0.117819622815556 0.19812576837063 -0.168129122408109 0.0334278027906954 0.0341629287679525 -0.0448059132193436 0.172818741897772 0.000264903859908037 0.0406610626265162 -0.0818108725970905 -0.0340018005629163 -0.238100003197897 0.00397269541951085 0.0189627413480861 -0.113053693027687 -0.212327952679913 0.17133339855926 0.104615602654068 0.0591425553751558 -0.0210562617949824 0.151200479401401 -0.273735841794059 -0.108631289474458 -0.0378621921715031 0.103860617232709 0.0983542185986925 0.0103329441280026 0.110219929923621 0.35782815360866 0.13351098127867 0.0663861640050136 -0.215302598317922 -0.131173201874798 0.109061224968608 0.0580779813126851 -0.0621925064159906 -0.0739719635176505 -0.316637561783232 0.00528577565867649 0.179146880263971 0.0460731097071163 0.0331282476898352 -0.256960651581475 -0.18522873356683 -0.191442896365065 -0.0727792767500551 -0.0636102018938309 -0.036909363748502 0.104790703021071 0.0762554967286991 0.0469421885705287 -0.0103024654292262 -0.0813499403171658 0.201634739282139 -0.0131210674600976 0.143276895835489 0.0125717196543428 0.0823586484894057 0.118164122692306 -0.353133656579894 -0.29898269307436 0.162450395485489 -0.104803390059382 0.127237688427226 -0.185866116614663 -0.0738365697199846 -0.0692042094692913 -0.0554190036354777 0.0865118390707002 -0.218077961751833 -0.0840119461651836 -0.115182302562131 0.0385281570884781 0.103971130963065 0.0390919775238502 -0.144634021606588 0.00387258757976661 -0.192081871803229 -0.0778488445947311 0.193471878953273 -0.0937084898096633 0.149521593287202 -0.0124618768302845 -0.0322009156170574 0.264754711345494 -0.180559758939988 -0.370275351860741 0.161459008505465 -0.310527635010566 -0.114678483337918 0.214716249229012 0.162412274833786 0.115856388849552 -0.119759590941546 0.0572554618915419 -0.173654888485154 -0.00449361735518098 -0.073785909159818 0.22996418528914 -0.351961151089054 -0.0190265195732244 -0.208142786031066 -0.0218964300189377 0.166784478379786 0.169470419732816 -0.0568361897331848 0.105711924964955 0.159760607311472 0.118586739034524 0.0163259889618812 -0.0910949736837612 -0.166423405357875 -0.0486179136523914 -0.0652975372246972 -0.211296693244894 -0.0188080034595025 0.0728513645580383 -0.0216314623818401 -0.0393605905500316 0.0732184496268977 0.142606854007517 0.230264901552793 -0.17765370911218 -0.0648366049447725 0.263595177749921 0.0359022816230272 -0.0412107660544421 -0.0412107660544421 0.0760788626708624 0.0803742085384491 0.18959900560191 0.245764462044252 -0.154866420165325 0.206987155672928 0.190054437409483 -0.0335814517548916 0.19256840584008 0.171711336771808 0.128443224537509 0.12910559425248 0.127449669965053 0.0167873002500829 0.167835094375898 -0.027788271388026 -0.307982401412077 -0.0833466616670683 -0.104435232710914 -0.327740095639977 0.0242646599257948 0.001750336994704 -0.127413156580059 -0.0425325470784345 -0.148628836062361 0.0837688546564615 -0.0687395648357738 -0.0839019857359938 0.262753608830748 -0.0732092377998709 0.0721122234881066 0.180483702518402 0.0206222957163605 -0.441809713262175 -0.0885037131586943 0.175924471951188 0.127640567064511 0.0580821483701045 0.179862227747864 -0.014792911557665 -0.0983958147575908 -0.0792947305261701 0.138910890677312 0.16245577581659 0.0265702132032087 -0.106524504438134 -0.254099743989774 -0.11083619484143 -0.341286273923522 -0.284333592785057 0.0984697956529047 -0.0800113908073176 0.0291538503295974 0.152834215132548 -0.0191376904535237 -0.0347911453237658 0.0990354144948576 0.00685885057672529 0.0435164077429748 -0.0437837970360358 0.109898369289943 0.0937454153718923 -0.0806082064769133 -0.144510621026447 0.135886463474767 0.0126767074010823 0.0803158030782113 -0.0383048899704326 -0.0216314623818401 -0.168359226669814 -0.102967113990753 -0.103360472601493 0.112360425417033 0.0676322818196781 -0.0619347965284489 0.11687284394656 0.272457795440399 -0.0892970766494202 0.125500313013645 -0.0854622488862674 -0.0388471191074806 -0.199703244467343 0.234958355545567 -0.0847219378049103 -0.187094292102355 -0.0393605905500316 0.0784389477571962 0.131338314864342 -0.0383048899704326 -0.165572795659583 0.189178128705204 0.179702233658296 -0.245405493033845 0.0889430188286013 0.0990755489441131 0.109898829599236 0.00998583316884095 -0.34807552629503 -0.0101848133811379 0.0823199909017991 0.149317636174818 0.0776051496313301 0.194778341324341
Chart options
Title:
Y-axis minimum
Y-axis maximum
R Code
geomean <- function(x) { return(exp(mean(log(x)))) } harmean <- function(x) { return(1/mean(1/x)) } quamean <- function(x) { return(sqrt(mean(x*x))) } winmean <- function(x) { x <-sort(x[!is.na(x)]) n<-length(x) denom <- 3 nodenom <- n/denom if (nodenom>40) denom <- n/40 sqrtn = sqrt(n) roundnodenom = floor(nodenom) win <- array(NA,dim=c(roundnodenom,2)) for (j in 1:roundnodenom) { win[j,1] <- (j*x[j+1]+sum(x[(j+1):(n-j)])+j*x[n-j])/n win[j,2] <- sd(c(rep(x[j+1],j),x[(j+1):(n-j)],rep(x[n-j],j)))/sqrtn } return(win) } trimean <- function(x) { x <-sort(x[!is.na(x)]) n<-length(x) denom <- 3 nodenom <- n/denom if (nodenom>40) denom <- n/40 sqrtn = sqrt(n) roundnodenom = floor(nodenom) tri <- array(NA,dim=c(roundnodenom,2)) for (j in 1:roundnodenom) { tri[j,1] <- mean(x,trim=j/n) tri[j,2] <- sd(x[(j+1):(n-j)]) / sqrt(n-j*2) } return(tri) } midrange <- function(x) { return((max(x)+min(x))/2) } q1 <- function(data,n,p,i,f) { np <- n*p; i <<- floor(np) f <<- np - i qvalue <- (1-f)*data[i] + f*data[i+1] } q2 <- function(data,n,p,i,f) { np <- (n+1)*p i <<- floor(np) f <<- np - i qvalue <- (1-f)*data[i] + f*data[i+1] } q3 <- function(data,n,p,i,f) { np <- n*p i <<- floor(np) f <<- np - i if (f==0) { qvalue <- data[i] } else { qvalue <- data[i+1] } } q4 <- function(data,n,p,i,f) { np <- n*p i <<- floor(np) f <<- np - i if (f==0) { qvalue <- (data[i]+data[i+1])/2 } else { qvalue <- data[i+1] } } q5 <- function(data,n,p,i,f) { np <- (n-1)*p i <<- floor(np) f <<- np - i if (f==0) { qvalue <- data[i+1] } else { qvalue <- data[i+1] + f*(data[i+2]-data[i+1]) } } q6 <- function(data,n,p,i,f) { np <- n*p+0.5 i <<- floor(np) f <<- np - i qvalue <- data[i] } q7 <- function(data,n,p,i,f) { np <- (n+1)*p i <<- floor(np) f <<- np - i if (f==0) { qvalue <- data[i] } else { qvalue <- f*data[i] + (1-f)*data[i+1] } } q8 <- function(data,n,p,i,f) { np <- (n+1)*p i <<- floor(np) f <<- np - i if (f==0) { qvalue <- data[i] } else { if (f == 0.5) { qvalue <- (data[i]+data[i+1])/2 } else { if (f < 0.5) { qvalue <- data[i] } else { qvalue <- data[i+1] } } } } midmean <- function(x,def) { x <-sort(x[!is.na(x)]) n<-length(x) if (def==1) { qvalue1 <- q1(x,n,0.25,i,f) qvalue3 <- q1(x,n,0.75,i,f) } if (def==2) { qvalue1 <- q2(x,n,0.25,i,f) qvalue3 <- q2(x,n,0.75,i,f) } if (def==3) { qvalue1 <- q3(x,n,0.25,i,f) qvalue3 <- q3(x,n,0.75,i,f) } if (def==4) { qvalue1 <- q4(x,n,0.25,i,f) qvalue3 <- q4(x,n,0.75,i,f) } if (def==5) { qvalue1 <- q5(x,n,0.25,i,f) qvalue3 <- q5(x,n,0.75,i,f) } if (def==6) { qvalue1 <- q6(x,n,0.25,i,f) qvalue3 <- q6(x,n,0.75,i,f) } if (def==7) { qvalue1 <- q7(x,n,0.25,i,f) qvalue3 <- q7(x,n,0.75,i,f) } if (def==8) { qvalue1 <- q8(x,n,0.25,i,f) qvalue3 <- q8(x,n,0.75,i,f) } midm <- 0 myn <- 0 roundno4 <- round(n/4) round3no4 <- round(3*n/4) for (i in 1:n) { if ((x[i]>=qvalue1) & (x[i]<=qvalue3)){ midm = midm + x[i] myn = myn + 1 } } midm = midm / myn return(midm) } (arm <- mean(x)) sqrtn <- sqrt(length(x)) (armse <- sd(x) / sqrtn) (armose <- arm / armse) (geo <- geomean(x)) (har <- harmean(x)) (qua <- quamean(x)) (win <- winmean(x)) (tri <- trimean(x)) (midr <- midrange(x)) midm <- array(NA,dim=8) for (j in 1:8) midm[j] <- midmean(x,j) midm bitmap(file='test1.png') lb <- win[,1] - 2*win[,2] ub <- win[,1] + 2*win[,2] if ((ylimmin == '') | (ylimmax == '')) plot(win[,1],type='b',main=main, xlab='j', pch=19, ylab='Winsorized Mean(j/n)', ylim=c(min(lb),max(ub))) else plot(win[,1],type='l',main=main, xlab='j', pch=19, ylab='Winsorized Mean(j/n)', ylim=c(ylimmin,ylimmax)) lines(ub,lty=3) lines(lb,lty=3) grid() dev.off() bitmap(file='test2.png') lb <- tri[,1] - 2*tri[,2] ub <- tri[,1] + 2*tri[,2] if ((ylimmin == '') | (ylimmax == '')) plot(tri[,1],type='b',main=main, xlab='j', pch=19, ylab='Trimmed Mean(j/n)', ylim=c(min(lb),max(ub))) else plot(tri[,1],type='l',main=main, xlab='j', pch=19, ylab='Trimmed Mean(j/n)', ylim=c(ylimmin,ylimmax)) lines(ub,lty=3) lines(lb,lty=3) grid() dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Central Tendency - Ungrouped Data',4,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Measure',header=TRUE) a<-table.element(a,'Value',header=TRUE) a<-table.element(a,'S.E.',header=TRUE) a<-table.element(a,'Value/S.E.',header=TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Arithmetic Mean',header=TRUE) a<-table.element(a,signif(arm,6)) a<-table.element(a, signif(armse,6)) a<-table.element(a,signif(armose,6)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'Geometric Mean',header=TRUE) a<-table.element(a,signif(geo,6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'Harmonic Mean',header=TRUE) a<-table.element(a,signif(har,6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'Quadratic Mean',header=TRUE) a<-table.element(a,signif(qua,6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) for (j in 1:length(win[,1])) { a<-table.row.start(a) mylabel <- paste('Winsorized Mean (',j) mylabel <- paste(mylabel,'/') mylabel <- paste(mylabel,length(win[,1])) mylabel <- paste(mylabel,')') a<-table.element(a, mylabel,header=TRUE) a<-table.element(a,signif(win[j,1],6)) a<-table.element(a,signif(win[j,2],6)) a<-table.element(a,signif(win[j,1]/win[j,2],6)) a<-table.row.end(a) } for (j in 1:length(tri[,1])) { a<-table.row.start(a) mylabel <- paste('Trimmed Mean (',j) mylabel <- paste(mylabel,'/') mylabel <- paste(mylabel,length(tri[,1])) mylabel <- paste(mylabel,')') a<-table.element(a, mylabel,header=TRUE) a<-table.element(a,signif(tri[j,1],6)) a<-table.element(a,signif(tri[j,2],6)) a<-table.element(a,signif(tri[j,1]/tri[j,2],6)) a<-table.row.end(a) } a<-table.row.start(a) a<-table.element(a, 'Median',header=TRUE) a<-table.element(a,signif(median(x),6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'Midrange',header=TRUE) a<-table.element(a,signif(midr,6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'Weighted Average at Xnp',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[1],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'Weighted Average at X(n+1)p',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[2],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'Empirical Distribution Function',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[3],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'Empirical Distribution Function - Averaging',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[4],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'Empirical Distribution Function - Interpolation',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[5],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'Closest Observation',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[6],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'True Basic - Statistics Graphics Toolkit',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[7],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'MS Excel (old versions)',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[8],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Number of observations',header=TRUE) a<-table.element(a,signif(length(x),6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.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