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 computationSat, 24 Oct 2009 10:40:11 -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/24/t1256402515kgxc7j4dsnxqrbi.htm/, Retrieved Fri, 03 May 2024 08:58:43 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=50193, Retrieved Fri, 03 May 2024 08:58:43 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact123

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Bivariate Data Series] [Bivariate dataset] [2008-01-05 23:51:08] [74be16979710d4c4e7c6647856088456]
- RMPD    [Bivariate Explorative Data Analysis] [bivariate EDA] [2009-10-24 16:40:11] [ea241b681aafed79da4b5b99fad98471] [Current]
-    D      [Bivariate Explorative Data Analysis] [X(t) als Y(t) / X(t)] [2009-10-25 10:46:15] [cd6314e7e707a6546bd4604c9d1f2b69]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
204114.4389
203580.4908
201576.0985
197317.2185
194799.2051
196998.8823
215441.3247
221696.6115
222368.1301
221900.245
215672.5449
216139.7041
214473.9749
213934.5821
212137.4531
208025.5812
205976.9049
206547.5142
225316.2784
228274.2272
228037.5622
222188.0905
216121.5549
216752.0565
215397.0408
214186.1294
210457.2046
208426.6774
207907.6117
207968.2299
224892.3142
227158.0543
225018.9953
213297.9099
205354.752
202281.0116
203479.5816
199224.6943
192988.2827
190899.6783
185498.1276
181005.8495
201586.988
204936.9585
196612.3054
191317.4718
185065.4519
186666.933
187633.9199
184236.0356
178813.432
177949.8953
170368.6315
173353.441
191792.6166
193684.1212
188005.2515
183732.5781
182168.8472
187350.4301
191735.6282
193338.5613
194675.7909
194753.1063
190056.8317
194636.5887
212914.9627
216553.5048
210720.0043
Dataseries Y:
Download CSV
Histogram 
562325
560854
555332
543599
536662
542722
593530
610763
612613
611324
594167
595454
590865
589379
584428
573100
567456
569028
620735
628884
628232
612117
595404
597141
593408
590072
579799
574205
572775
572942
619567
625809
619916
587625
565742
557274
560576
548854
531673
525919
511038
498662
555362
564591
541657
527070
509846
514258
516922
507561
492622
490243
469357
477580
528379
533590
517945
506174
501866
516141
528222
532638
536322
536535
523597
536214
586570
596594
580523

Tables (Output of Computation)





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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=50193&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 time8 seconds
R Server'George Udny Yule' @ 72.249.76.132






Model: Y[t] = c + b X[t] + e[t]
c0.000129070089478634
b2.75494963942146

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

[TABLE]
[ROW][C]Model: Y[t] = c + b X[t] + e[t][/C][/ROW]
[ROW][C]c[/C][C]0.000129070089478634[/C][/ROW]
[ROW][C]b[/C][C]2.75494963942146[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=50193&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=50193&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]
c0.000129070089478634
b2.75494963942146






Descriptive Statistics about e[t]
# observations69
minimum-0.000141184117480032
Q1-7.31128829036564e-05
median-1.24737024314189e-05
mean8.05074510385669e-22
Q36.42119308685378e-05
maximum0.000153131023167299

\begin{tabular}{lllllllll}
\hline
Descriptive Statistics about e[t] \tabularnewline
# observations & 69 \tabularnewline
minimum & -0.000141184117480032 \tabularnewline
Q1 & -7.31128829036564e-05 \tabularnewline
median & -1.24737024314189e-05 \tabularnewline
mean & 8.05074510385669e-22 \tabularnewline
Q3 & 6.42119308685378e-05 \tabularnewline
maximum & 0.000153131023167299 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=50193&T=2

[TABLE]
[ROW][C]Descriptive Statistics about e[t][/C][/ROW]
[ROW][C]# observations[/C][C]69[/C][/ROW]
[ROW][C]minimum[/C][C]-0.000141184117480032[/C][/ROW]
[ROW][C]Q1[/C][C]-7.31128829036564e-05[/C][/ROW]
[ROW][C]median[/C][C]-1.24737024314189e-05[/C][/ROW]
[ROW][C]mean[/C][C]8.05074510385669e-22[/C][/ROW]
[ROW][C]Q3[/C][C]6.42119308685378e-05[/C][/ROW]
[ROW][C]maximum[/C][C]0.000153131023167299[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=50193&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=50193&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]
# observations69
minimum-0.000141184117480032
Q1-7.31128829036564e-05
median-1.24737024314189e-05
mean8.05074510385669e-22
Q36.42119308685378e-05
maximum0.000153131023167299


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