R version 3.0.2 (2013-09-25) -- "Frisbee Sailing" Copyright (C) 2013 The R Foundation for Statistical Computing Platform: i686-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > > RC.capture <- function (expression, collapse = NULL) { + resultConn <- textConnection('RC.resultText', open = 'w', local=TRUE) + sink(resultConn) + on.exit(function() { + sink() + close(resultConn) + }) + expression + on.exit(NULL) + sink() + close(resultConn) + return(paste(c(RC.resultText, ''), collapse = collapse, sep = '')) + } > RC.texteval <- function (sourceText, collapse = NULL, echo = TRUE) { + sourceConn <- textConnection(sourceText, open = 'r') + on.exit(close(sourceConn)) + result <- RC.capture(source(file = sourceConn, local = FALSE, echo = echo, print.eval = TRUE), collapse = collapse) + on.exit(NULL) + close(sourceConn) + res <- '' + for(i in 1:length(result)) { + if (result[i]!='') res <- paste(res,result[i],' + ',sep='') + } + return(res) + } > x <- c(0.00550118759611825 + ,0.100183846454512 + ,0.407011937175347 + ,-0.16737394013665 + ,0.0874322043971727 + ,0.411921594926703 + ,-0.16303809760263 + ,-0.347349364549994 + ,-0.0296049307932174 + ,0.323806680916916 + ,-0.239140348405293 + ,-0.0990421431826477 + ,0.0467509732678099 + ,0.163478509066741 + ,-0.350096927982722 + ,0.32436525705737 + ,0.311383604987517 + ,-0.0960371667588798 + ,0.124548267674515 + ,-0.511624147113553 + ,0.314584415989742 + ,0.287739894396906 + ,0.282204368241592 + ,-0.25961081368395 + ,0.139880423671516 + ,0.0468251561703096 + ,0.141297672568196 + ,0.131779095699905 + ,0.188261365013337 + ,0.34246868819404 + ,0.195820162553774 + ,-0.179557156234377 + ,0.232001180091104 + ,0.0900659140956214 + ,-0.1338451790818 + ,-0.00791389868766919 + ,-0.0425554884595751 + ,0.482693233899655 + ,0.282396103736813 + ,0.156155004245922 + ,-0.114835337644616 + ,0.00659340969595179 + ,-0.224997452667842 + ,-0.127774776710956 + ,0.0620755029470993 + ,-0.0998584806098916 + ,0.267864642498592 + ,0.389744472117138 + ,0.0761067268738662 + ,-0.135349114066045 + ,-0.199236081602781 + ,0.202107333568291 + ,0.47797312151875 + ,-0.330634318524646 + ,-0.00413610310225685 + ,0.14024473328936 + ,-0.223728014463631 + ,0.223086730239358 + ,0.073422994370519 + ,0.225451904074953 + ,-0.0398059712998328 + ,-0.0698047210357811 + ,-0.259657758635374 + ,0.167444170403979 + ,-0.0135427590292624 + ,0.192334314572137) > #'GNU S' R Code compiled by R2WASP v. 1.2.327 () > #Author: root > #To cite this work: Wessa P., (2013), Skewness and Kurtosis Test (v1.0.3) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_skewness_kurtosis.wasp/ > #Source of accompanying publication: > # > library(moments) > agostino <- agostino.test(x) > anscombe <- anscombe.test(x) > geary <- geary(x) > jarque <- jarque.test(x) > > #Note: the /var/fisher/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/fisher/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Skewness and Kurtosis Test',1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,paste('
',RC.texteval('agostino'),'
',sep='')) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,paste('
',RC.texteval('anscombe'),'
',sep='')) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,paste('
',RC.texteval('jarque'),'
',sep='')) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,paste('
',RC.texteval('geary'),'
',sep='')) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/fisher/rcomp/tmp/11w8n1386839573.tab") > > > > proc.time() user system elapsed 0.753 0.237 0.972