FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_edabi.wasp
Title produced by softwareBivariate Explorative Data Analysis
Date of computationWed, 28 Oct 2009 11:49:47 -0600
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2009/Oct/28/t125675226609lanlsbd6q2tpd.htm/, Retrieved Mon, 29 Apr 2024 08:39:54 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=51656, Retrieved Mon, 29 Apr 2024 08:39:54 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsWSH 5 e(t) en e'(t)
Estimated Impact145

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Bivariate Explorative Data Analysis] [Workshop 5] [2009-10-28 17:49:47] [e7a989b306049c061a54f626f1127c12] [Current]
-    D    [Bivariate Explorative Data Analysis] [SHW_W5.4] [2009-10-30 15:14:38] [f966872135bb25240f339c0c372beeec]
- RMPD      [Kendall tau Rank Correlation] [Kendall tau Rank ...] [2009-11-08 23:21:57] [8733f8ed033058987ec00f5e71b74854]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
0.614147829
0.634181616
0.379802399
0.193130101
0.886550717
0.648704427
0.383092091
0.327602232
0.570959498
0.813290589
0.408848021
0.419912206
0.584456134
0.496715428
0.489025402
0.596799895
-0.064369874
0.15780073
0.598936712
0.376808342
0.467880973
-0.177655342
-0.055315805
0.122471358
-0.160869014
-0.424133366
-0.497477961
-0.510974598
-0.54993153
-0.271033725
-0.28209791
-0.381113969
-0.441046395
-0.56432764
-0.745362262
-1.242030336
-1.36420094
-1.570780324
-1.539724585
-1.589619001
-1.00508352
-0.67159533
-0.386075908
-0.043829284
-0.047118976
0.407471409
-0.132427231
0.276711305
0.659067736
0.705756928
0.936559264
1.607429175
2.013628457
1.124439241
1.085355608
0.41427453
0.267585339
-1.221223776
-1.253390156
-2.476586934
Dataseries Y:
Download CSV
Histogram 
26.35404199
8.085095727
-7.521619976
-6.126795599
-21.64657344
-11.37621019
-1.961731053
-11.93468212
-15.31268572
-24.17929035
-29.24423181
-21.72069533
-33.86259266
-33.99518648
-51.38119969
-49.91755214
-41.33813154
-38.65224142
-36.96271575
-47.0097887
-34.27848879
-35.69413875
-48.75297994
-41.07768723
-43.08415676
-47.12075538
-54.03110663
-25.58119969
3.422928203
13.26445627
11.2409198
29.99070365
25.55281112
17.60068563
-3.895678825
-0.996972731
4.217137145
-15.2026407
-0.486867659
11.41017527
2.460514505
-18.79715587
15.10296722
26.25872825
49.79883933
51.09964089
46.8928021
5.755588042
4.912273719
-17.36702379
-3.184153774
20.06082372
28.13124702
50.92188047
72.90955685
84.1704935
77.24979101
69.189779
60.70777059
56.93327945

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 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 & 4 seconds \tabularnewline
R Server & 'George Udny Yule' @ 72.249.76.132 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=51656&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]4 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=51656&T=0

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

As an alternative you can also use a QR Code:  
QR Code


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

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'George Udny Yule' @ 72.249.76.132






Model: Y[t] = c + b X[t] + e[t]
c-2.66668199412352e-10
b-4.76090275302302

\begin{tabular}{lllllllll}
\hline
Model: Y[t] = c + b X[t] + e[t] \tabularnewline
c & -2.66668199412352e-10 \tabularnewline
b & -4.76090275302302 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=51656&T=1

[TABLE]
[ROW][C]Model: Y[t] = c + b X[t] + e[t][/C][/ROW]
[ROW][C]c[/C][C]-2.66668199412352e-10[/C][/ROW]
[ROW][C]b[/C][C]-4.76090275302302[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=51656&T=1

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

As an alternative you can also use a QR Code:  
QR Code


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

Model: Y[t] = c + b X[t] + e[t]
c-2.66668199412352e-10
b-4.76090275302302






Descriptive Statistics about e[t]
# observations60
minimum-56.3995508238265
Q1-28.7804385056452
median-5.46036982896015
mean4.41799374861773e-15
Q324.1022894159272
maximum86.142814250651

\begin{tabular}{lllllllll}
\hline
Descriptive Statistics about e[t] \tabularnewline
# observations & 60 \tabularnewline
minimum & -56.3995508238265 \tabularnewline
Q1 & -28.7804385056452 \tabularnewline
median & -5.46036982896015 \tabularnewline
mean & 4.41799374861773e-15 \tabularnewline
Q3 & 24.1022894159272 \tabularnewline
maximum & 86.142814250651 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=51656&T=2

[TABLE]
[ROW][C]Descriptive Statistics about e[t][/C][/ROW]
[ROW][C]# observations[/C][C]60[/C][/ROW]
[ROW][C]minimum[/C][C]-56.3995508238265[/C][/ROW]
[ROW][C]Q1[/C][C]-28.7804385056452[/C][/ROW]
[ROW][C]median[/C][C]-5.46036982896015[/C][/ROW]
[ROW][C]mean[/C][C]4.41799374861773e-15[/C][/ROW]
[ROW][C]Q3[/C][C]24.1022894159272[/C][/ROW]
[ROW][C]maximum[/C][C]86.142814250651[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=51656&T=2

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

As an alternative you can also use a QR Code:  
QR Code


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

Descriptive Statistics about e[t]
# observations60
minimum-56.3995508238265
Q1-28.7804385056452
median-5.46036982896015
mean4.41799374861773e-15
Q324.1022894159272
maximum86.142814250651


Figures (Output of Computation)

Input Parameters & R Code

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)
y <- as.ts(y)
mylm <- lm(y~x)
cbind(mylm$resid)
library(lattice)
bitmap(file='pic1.png')
plot(y,type='l',main='Run Sequence Plot of Y[t]',xlab='time or index',ylab='value')
grid()
dev.off()
bitmap(file='pic1a.png')
plot(x,type='l',main='Run Sequence Plot of X[t]',xlab='time or index',ylab='value')
grid()
dev.off()
bitmap(file='pic1b.png')
plot(x,y,main='Scatter Plot',xlab='X[t]',ylab='Y[t]')
grid()
dev.off()
bitmap(file='pic1c.png')
plot(mylm$resid,type='l',main='Run Sequence Plot of e[t]',xlab='time or index',ylab='value')
grid()
dev.off()
bitmap(file='pic2.png')
hist(mylm$resid,main='Histogram of e[t]')
dev.off()
bitmap(file='pic3.png')
if (par1 > 0)
{
densityplot(~mylm$resid,col='black',main=paste('Density Plot of e[t]   bw = ',par1),bw=par1)
} else {
densityplot(~mylm$resid,col='black',main='Density Plot of e[t]')
}
dev.off()
bitmap(file='pic4.png')
qqnorm(mylm$resid,main='QQ plot of e[t]')
qqline(mylm$resid)
grid()
dev.off()
if (par2 > 0)
{
bitmap(file='pic5.png')
acf(mylm$resid,lag.max=par2,main='Residual Autocorrelation Function')
grid()
dev.off()
}
summary(x)
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Model: Y[t] = c + b X[t] + e[t]',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'c',1,TRUE)
a<-table.element(a,mylm$coeff[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'b',1,TRUE)
a<-table.element(a,mylm$coeff[[2]])
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,'Descriptive Statistics about e[t]',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'# observations',header=TRUE)
a<-table.element(a,length(mylm$resid))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'minimum',header=TRUE)
a<-table.element(a,min(mylm$resid))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Q1',header=TRUE)
a<-table.element(a,quantile(mylm$resid,0.25))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'median',header=TRUE)
a<-table.element(a,median(mylm$resid))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'mean',header=TRUE)
a<-table.element(a,mean(mylm$resid))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Q3',header=TRUE)
a<-table.element(a,quantile(mylm$resid,0.75))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'maximum',header=TRUE)
a<-table.element(a,max(mylm$resid))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')