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 computationTue, 08 Dec 2009 13:18:56 -0700
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/Dec/08/t1260303578r1z2jz7anpwaici.htm/, Retrieved Sat, 27 Apr 2024 22:45:39 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=64839, Retrieved Sat, 27 Apr 2024 22:45:39 +0000
QR Codes:

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

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 4] [2009-10-29 12:50:48] [1646a2766cb8c4a6f9d3b2fffef409b3]
- RM D    [Bivariate Explorative Data Analysis] [] [2009-12-08 20:18:56] [d41d8cd98f00b204e9800998ecf8427e] [Current]
-    D      [Bivariate Explorative Data Analysis] [] [2009-12-08 20:28:08] [74be16979710d4c4e7c6647856088456]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
12,59443229
12,58959203
12,57954817
12,55931771
12,54901675
12,57020869
12,66432942
12,69284735
12,6946864
12,70162827
12,6681701
12,66835592
12,65631033
12,65164269
12,64236113
12,62132305
12,61259053
12,62489811
12,71732618
12,72990728
12,72466404
12,68828581
12,65691575
12,65869166
12,65465471
12,64324312
12,62136596
12,60816206
12,60626722
12,61695969
12,70613556
12,72213354
12,69959747
12,6414816
12,59783091
12,58249678
12,58293936
12,559114
12,52987568
12,51563607
12,4874134
12,48101137
12,60083401
12,62345974
12,56721659
12,52970205
12,4928443
12,49557681
12,50016289
12,47809017
12,45114206
12,44176908
12,4041158
12,43217581
12,54413291
12,56114581
12,50978554
12,47256683
12,45138074
12,46988202
12,48333689
12,48480605
12,48095819
12,47828061
12,45358485
12,48454115
12,59085874
Dataseries Y:
Download CSV
Histogram 
12,49654956
12,49637379
12,48664324
12,46405424
12,44840687
12,44850104
12,53273622
12,56147528
12,56584977
12,55335058
12,53064218
12,53507059
12,53223477
12,53214456
12,52466091
12,50673264
12,49552814
12,48745869
12,56813953
12,58171752
12,58553595
12,57136138
12,54780602
12,55198589
12,54322322
12,54400096
12,53118055
12,52531152
12,52216723
12,51102104
12,57695321
12,58013652
12,58553595
12,53719488
12,50553259
12,49072889
12,50259304
12,48437083
12,44979499
12,44253905
12,4133185
12,36892252
12,46286758
12,47224887
12,44769638
12,43178489
12,40255184
12,41756863
12,4233631
12,40917266
12,37615292
12,37617401
12,32631488
12,33199579
12,42091786
12,42252587
12,4170181
12,40941341
12,41415535
12,45201433
12,48491199
12,5000362
12,51739556
12,52074834
12,49661313
12,51353161
12,58693764

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'Gwilym Jenkins' @ 72.249.127.135

\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 & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=64839&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]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=64839&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=64839&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'Gwilym Jenkins' @ 72.249.127.135






Model: Y[t] = c + b X[t] + e[t]
c4.61920095570182
b0.625572230955777

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=64839&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]
c4.61920095570182
b0.625572230955777






Descriptive Statistics about e[t]
# observations67
minimum-0.064429122797953
Q1-0.0180898515098264
median-0.00197372882437226
mean-9.01141917616217e-19
Q30.0103120593755818
maximum0.0954815446082772

\begin{tabular}{lllllllll}
\hline
Descriptive Statistics about e[t] \tabularnewline
# observations & 67 \tabularnewline
minimum & -0.064429122797953 \tabularnewline
Q1 & -0.0180898515098264 \tabularnewline
median & -0.00197372882437226 \tabularnewline
mean & -9.01141917616217e-19 \tabularnewline
Q3 & 0.0103120593755818 \tabularnewline
maximum & 0.0954815446082772 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=64839&T=2

[TABLE]
[ROW][C]Descriptive Statistics about e[t][/C][/ROW]
[ROW][C]# observations[/C][C]67[/C][/ROW]
[ROW][C]minimum[/C][C]-0.064429122797953[/C][/ROW]
[ROW][C]Q1[/C][C]-0.0180898515098264[/C][/ROW]
[ROW][C]median[/C][C]-0.00197372882437226[/C][/ROW]
[ROW][C]mean[/C][C]-9.01141917616217e-19[/C][/ROW]
[ROW][C]Q3[/C][C]0.0103120593755818[/C][/ROW]
[ROW][C]maximum[/C][C]0.0954815446082772[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=64839&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=64839&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]
# observations67
minimum-0.064429122797953
Q1-0.0180898515098264
median-0.00197372882437226
mean-9.01141917616217e-19
Q30.0103120593755818
maximum0.0954815446082772


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