Home » date » 2011 » Nov » 22 »

*The author of this computation has been verified* R Software Module: /rwasp_edauni.wasp (opens new window with default values) Title produced by software: Univariate Explorative Data Analysis Date of computation: Tue, 22 Nov 2011 12:11:26 -0500   Cite this page as follows: Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2011/Nov/22/t132198190770mibkusdazjum3.htm/, Retrieved Tue, 18 Jun 2013 21:43:30 +0000   Original text written by user:   IsPrivate? No (this computation is public)   User-defined keywords:   System-generated keywords (parent): t1289984313fiq9ol6jc8ggb97 (pk = 96536) Estimated Impact 18   Dataseries X: » Textfile « » CSV « » Stem and Leaf « » Histogram « » Kernel Density « » Harrell-Davis Quantiles « » Central Tendency « » Variability « 23 20 20 21 24 22 23 20 25 23 27 27 22 24 25 22 28 28 27 25 16 28 21 24 27 14 14 27 20 21 22 21 12 20 24 19 28 23 27 22 27 26 22 21 19 24 19 26 22 28 21 23 28 10 24 21 21 24 24 25 25 23 21 16 17 25 24 23 25 23 28 26 22 19 26 18 18 25 27 12 15 21 23 22 21 24 27 22 28 26 10 19 22 21 24 25 21 20 21 24 23 18 24 24 19 20 18 20 27 23 26 23 17 21 25 23 27 24 20 27 21 24 21 15 25 25 22 24 21 22 23 22 20 23 25 23 22 25 26 22 24 24 25 20 26 21 26 21 22 16 26 28 18 25 23 21 20 25 22 21 16 18   Output produced by software: Summary of computational transaction Raw Input view raw input (R code) Raw Output view raw output of R engine Computing time 4 seconds R Server 'George Udny Yule' @ yule.wessa.net Descriptive Statistics # observations 162 minimum 10 Q1 21 median 23 mean 22.3950617283951 Q3 25 maximum 28   Charts produced by software: http://www.freestatistics.org/blog/date/2011/Nov/22/t132198190770mibkusdazjum3/1tn4g1321981882.png (opens in new window) http://www.freestatistics.org/blog/date/2011/Nov/22/t132198190770mibkusdazjum3/1tn4g1321981882.ps (opens in new window) Click here to open pdf file. http://www.freestatistics.org/blog/date/2011/Nov/22/t132198190770mibkusdazjum3/226pw1321981882.png (opens in new window) http://www.freestatistics.org/blog/date/2011/Nov/22/t132198190770mibkusdazjum3/226pw1321981882.ps (opens in new window) Click here to open pdf file. http://www.freestatistics.org/blog/date/2011/Nov/22/t132198190770mibkusdazjum3/3hu7e1321981882.png (opens in new window) http://www.freestatistics.org/blog/date/2011/Nov/22/t132198190770mibkusdazjum3/3hu7e1321981882.ps (opens in new window) Click here to open pdf file. http://www.freestatistics.org/blog/date/2011/Nov/22/t132198190770mibkusdazjum3/4jfzx1321981882.png (opens in new window) http://www.freestatistics.org/blog/date/2011/Nov/22/t132198190770mibkusdazjum3/4jfzx1321981882.ps (opens in new window) Click here to open pdf file. http://www.freestatistics.org/blog/date/2011/Nov/22/t132198190770mibkusdazjum3/58w0z1321981882.png (opens in new window) http://www.freestatistics.org/blog/date/2011/Nov/22/t132198190770mibkusdazjum3/58w0z1321981882.ps (opens in new window) Click here to open pdf file. http://www.freestatistics.org/blog/date/2011/Nov/22/t132198190770mibkusdazjum3/65w931321981882.png (opens in new window) http://www.freestatistics.org/blog/date/2011/Nov/22/t132198190770mibkusdazjum3/65w931321981882.ps (opens in new window) Click here to open pdf file. http://www.freestatistics.org/blog/date/2011/Nov/22/t132198190770mibkusdazjum3/7fg0r1321981882.png (opens in new window) http://www.freestatistics.org/blog/date/2011/Nov/22/t132198190770mibkusdazjum3/7fg0r1321981882.ps (opens in new window) Click here to open pdf file.   Parameters (Session): par1 = 0 ; par2 = 36 ;   Parameters (R input): par1 = 0 ; par2 = 36 ;   R code (references can be found in the software module): par1 <- as.numeric(par1) par2 <- as.numeric(par2) x <- as.ts(x) library(lattice) bitmap(file='pic1.png') plot(x,type='l',main='Run Sequence Plot',xlab='time or index',ylab='value') grid() dev.off() bitmap(file='pic2.png') hist(x) grid() dev.off() bitmap(file='pic3.png') if (par1 > 0) { densityplot(~x,col='black',main=paste('Density Plot bw = ',par1),bw=par1) } else { densityplot(~x,col='black',main='Density Plot') } dev.off() bitmap(file='pic4.png') qqnorm(x) qqline(x) grid() dev.off() if (par2 > 0) { bitmap(file='lagplot1.png') dum <- cbind(lag(x,k=1),x) dum dum1 <- dum[2:length(x),] dum1 z <- as.data.frame(dum1) z plot(z,main='Lag plot (k=1), lowess, and regression line') lines(lowess(z)) abline(lm(z)) dev.off() if (par2 > 1) { bitmap(file='lagplotpar2.png') dum <- cbind(lag(x,k=par2),x) dum dum1 <- dum[(par2+1):length(x),] dum1 z <- as.data.frame(dum1) z mylagtitle <- 'Lag plot (k=' mylagtitle <- paste(mylagtitle,par2,sep='') mylagtitle <- paste(mylagtitle,'), and lowess',sep='') plot(z,main=mylagtitle) lines(lowess(z)) dev.off() } bitmap(file='pic5.png') acf(x,lag.max=par2,main='Autocorrelation Function') grid() dev.off() } summary(x) load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Descriptive Statistics',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'# observations',header=TRUE) a<-table.element(a,length(x)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'minimum',header=TRUE) a<-table.element(a,min(x)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Q1',header=TRUE) a<-table.element(a,quantile(x,0.25)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'median',header=TRUE) a<-table.element(a,median(x)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'mean',header=TRUE) a<-table.element(a,mean(x)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Q3',header=TRUE) a<-table.element(a,quantile(x,0.75)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'maximum',header=TRUE) a<-table.element(a,max(x)) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')