Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*The author of this computation has been verified*

R Software Module

rwasp_edauni.wasp

Title produced by software

Univariate Explorative Data Analysis

Date of computation

Sun, 26 Oct 2008 05:52:50 -0600

Cite this page as follows

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2008/Oct/26/t12250220287jk2g2wgp8mlbdk.htm/, Retrieved Sat, 18 May 2024 01:10:57 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=18858, Retrieved Sat, 18 May 2024 01:10:57 +0000

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords

Estimated Impact

177

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)

F     [Tukey lambda PPCC Plot] [Tukey lambda tot ...] [2008-10-24 09:02:28] [e1a46c1dcfccb0cb690f79a1a409b517]
F RMPD  [Univariate Explorative Data Analysis] [Univariate Explor...] [2008-10-24 12:02:31] [e1a46c1dcfccb0cb690f79a1a409b517]
-   PD    [Univariate Explorative Data Analysis] [UEDA - Vlaams gew...] [2008-10-26 09:55:38] [46c5a5fbda57fdfa1d4ef48658f82a0c]
F    D      [Univariate Explorative Data Analysis] [4plot prijs kledi...] [2008-10-26 10:30:49] [46c5a5fbda57fdfa1d4ef48658f82a0c]
F    D          [Univariate Explorative Data Analysis] [taak 2/1] [2008-10-26 11:52:50] [dbfa7caa6871c163dec68ca05d48bb00] [Current]
F    D            [Univariate Explorative Data Analysis] [taak 2/2] [2008-10-26 14:28:17] [29647dffafb5b58c12a48dbf6cba2b57] 

Feedback Forum

2008-11-02 14:50:17 [Evelyn Ongena] [reply] 
Als we naar de lag plots kijken, is er volgens mij geen autocorrelatie. Afgaande op het histogram en de density plot is de verdeling redelijk normaal. Bovendien dient de R code veranderd te worden door invoeging van x<-x-gevonden gemiddelde om de 4de voorwaarde voor het model na te gaan.
2008-11-03 08:07:01 [Thomas Baken] [reply] 
Volgens mij kunnen we hier inderdaad niet spreken over autocorrelatie. En om al te spreken van totaal geen normaalverdeling... In het histogram vind ik dat er toch enige normaalverdeling terug te vinden is, echter in het density plot bemerken we rechts een kleine afwijking. De determinerende component tenslotte is niet constant. Dit kunnen we zien door te kijken naar het Run Sequence Plot welke een licht stijgend verloop kent en niet vlak is. 

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

0,914835165
1,126410835
1,058324496
1,02034588
1,15625
1,239454094
0,960614793
0,925138632
1,071734475
1,393602225
1,065887354
1,060063224
1,037344398
1,098792536
1,187203791
1,157407407
1,1375
1,332889481
0,912488605
0,975728155
1,224116931
1,477941176
0,998962656
1,021208908
1,075555556
1,1
1,273324573
1,173333333
1,189189189
1,455639098
0,995884774
1,030818278
1,201982652
1,375886525
1,102505695
1,082681564
0,975903614
1,155581948
1,295605859
1,056521739
1,204207921
1,331053352
0,9749499
1,081111111
1,17087846
1,343922652
1,244923858
1,10882016
1,063736264
1,208489388
1,315217391
1,12037037
1,299328859
1,359550562
1,047619048
1,187730061
1,134818288
1,384835479
1,150831354
1,070562293
0,968095713

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'George Udny Yule' @ 72.249.76.132

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 2 seconds \tabularnewline
R Server & 'George Udny Yule' @ 72.249.76.132 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=18858&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]2 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'George Udny Yule' @ 72.249.76.132[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=18858&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=18858&T=0

As an alternative you can also use a QR Code:

The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'George Udny Yule' @ 72.249.76.132

Descriptive Statistics
# observations	61
minimum	0.912488605
Q1	1.056521739
median	1.126410835
mean	1.14533954259016
Q3	1.224116931
maximum	1.477941176

\begin{tabular}{lllllllll}
\hline
Descriptive Statistics \tabularnewline
# observations & 61 \tabularnewline
minimum & 0.912488605 \tabularnewline
Q1 & 1.056521739 \tabularnewline
median & 1.126410835 \tabularnewline
mean & 1.14533954259016 \tabularnewline
Q3 & 1.224116931 \tabularnewline
maximum & 1.477941176 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=18858&T=1

[TABLE]
[ROW][C]Descriptive Statistics[/C][/ROW]
[ROW][C]# observations[/C][C]61[/C][/ROW]
[ROW][C]minimum[/C][C]0.912488605[/C][/ROW]
[ROW][C]Q1[/C][C]1.056521739[/C][/ROW]
[ROW][C]median[/C][C]1.126410835[/C][/ROW]
[ROW][C]mean[/C][C]1.14533954259016[/C][/ROW]
[ROW][C]Q3[/C][C]1.224116931[/C][/ROW]
[ROW][C]maximum[/C][C]1.477941176[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=18858&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=18858&T=1

As an alternative you can also use a QR Code:

The GUIDs for individual cells are displayed in the table below:

Descriptive Statistics
# observations	61
minimum	0.912488605
Q1	1.056521739
median	1.126410835
mean	1.14533954259016
Q3	1.224116931
maximum	1.477941176

Figure 1

PNG link

Postscript link

PDF link

Figure 2

PNG link

Postscript link

PDF link

Figure 3

PNG link

Postscript link

PDF link

Figure 4

PNG link

Postscript link

PDF link

Figure 5

PNG link

Postscript link

PDF link

Figure 6

PNG link

Postscript link

PDF link

Figure 7

PNG link

Postscript link

PDF link

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')

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code