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:59:00 -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/t1256752786izlte0fn7eiywhv.htm/, Retrieved Sun, 05 May 2024 22:06:46 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=51670, Retrieved Sun, 05 May 2024 22:06:46 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsWorkshop 4 Part 2.3
Estimated Impact97

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] [shw-ws4p2] [2009-10-28 17:59:00] [5b5bced41faf164488f2c271c918b21f] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
6395841
4822416
10252804
7387524
7441984
5541316
7273809
7027801
4272489
6974881
6446521
5262436
7354944
5354596
9560464
7166329
7912969
7118224
8637721
6848689
4977361
6155361
5861241
5798464
6553600
4410000
10989225
7845601
5774409
9144576
6285049
8880400
4888521
6105841
6728836
6012304
4981824
5631129
9778129
7851204
6974881
7767369
6859161
7873636
4809249
5396329
6395841
5817744
5116644
4639716
10432900
5267025
7371225
7469289
5368489
7452900
3659569
5712100
6170256
3841600
Dataseries Y:
Download CSV
Histogram 
226622650401
225249906025
221312852721
212752485001
206773916176
207595051876
267130821409
275826636864
273502850625
268930402225
259324359121
262387768644
269531258896
267298306081
260031664489
259210302129
250857734449
257019594841
324128678329
336068321796
334074752064
319749535296
299585454336
307789724944
316209405625
314557209316
308393630224
295499872801
288006102244
294547169284
352277860900
373031442169
375294687769
373717032976
353034423889
354565466116
349121448225
347367605641
341556087184
328443610000
322006311936
323792864784
385311940225
395495085456
394675445824
374687221689
354505923216
356577373881
352133054464
348184965184
336166880401
329711382025
328071200625
328262535364
383863267489
391636904481
384295847056
345303140625
320064010564
310554311076

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 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=51670&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=51670&T=0

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






Model: Y[t] = c + b X[t] + e[t]
c351005341553.397
b-5618.25486015437

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=51670&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]
c351005341553.397
b-5618.25486015437






Descriptive Statistics about e[t]
# observations60
minimum-112277764128.746
Q1-40836169650.213
median10610561299.7231
mean9.5367431640625e-06
Q337568485756.1802
maximum88725837626.689

\begin{tabular}{lllllllll}
\hline
Descriptive Statistics about e[t] \tabularnewline
# observations & 60 \tabularnewline
minimum & -112277764128.746 \tabularnewline
Q1 & -40836169650.213 \tabularnewline
median & 10610561299.7231 \tabularnewline
mean & 9.5367431640625e-06 \tabularnewline
Q3 & 37568485756.1802 \tabularnewline
maximum & 88725837626.689 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=51670&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]-112277764128.746[/C][/ROW]
[ROW][C]Q1[/C][C]-40836169650.213[/C][/ROW]
[ROW][C]median[/C][C]10610561299.7231[/C][/ROW]
[ROW][C]mean[/C][C]9.5367431640625e-06[/C][/ROW]
[ROW][C]Q3[/C][C]37568485756.1802[/C][/ROW]
[ROW][C]maximum[/C][C]88725837626.689[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=51670&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=51670&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-112277764128.746
Q1-40836169650.213
median10610561299.7231
mean9.5367431640625e-06
Q337568485756.1802
maximum88725837626.689


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