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 13:47:06 -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/t1256759283r3wwv84ua7hrn7d.htm/, Retrieved Mon, 06 May 2024 03:02:40 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=51787, Retrieved Mon, 06 May 2024 03:02:40 +0000
QR Codes:

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

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] [Workshop4 Part2 b...] [2009-10-28 19:47:06] [5ed0eef5d4509bbfdac0ae6d87f3b4bf] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
12,28411643
12,27179952
12,25231194
12,22609584
12,20725694
12,22423668
12,39452537
12,46961698
12,40214917
12,40680698
12,3612585
12,37673061
12,3683743
12,34737254
12,33351554
12,30898116
12,2883068
12,29888881
12,4546734
12,49066871
12,42369699
12,39302513
12,34254658
12,35365178
12,29809105
12,29825521
12,26721478
12,27732762
12,26028617
12,26177695
12,38859418
12,3910407
12,34835315
12,24910819
12,19562829
12,17868089
12,17860381
12,13181235
12,09013416
12,05877841
12,02755584
12,03203468
12,21750702
12,21794193
12,11459312
12,06379454
12,0224142
12,05115417
12,04179334
12,00873085
11,99482061
11,98050115
11,92504174
11,97107157
12,13634736
12,13923281
12,06870533
12,03979371
12,02471263
12,08974667
12,13571456
12,15069735
12,16297202
12,14900482
12,13352883
12,1823382
12,31730986
12,33226026
12,27431991
Dataseries Y:
Download CSV
Histogram 
12,46443275
12,46321948
12,45411952
12,43120622
12,42071214
12,42742259
12,47848625
12,4852789
12,51740655
12,51517736
12,49711407
12,49107112
12,48346572
12,48659792
12,48137592
12,46452542
12,45813013
12,45868567
12,51166585
12,52266282
12,54708456
12,51603222
12,50508034
12,49958892
12,52976355
12,51446502
12,50066578
12,48145185
12,48226016
12,49089438
12,5485305
12,56867548
12,57564456
12,54697079
12,52043803
12,50254844
12,51587084
12,50723964
12,47015466
12,4770459
12,44037768
12,38280779
12,4624111
12,49395095
12,48096959
12,4609629
12,42451717
12,42232853
12,44262983
12,43657307
12,38628364
12,38718097
12,34057294
12,34374981
12,4292122
12,44515254
12,43391713
12,40055331
12,39453365
12,40486203
12,41530387
12,41845045
12,41956578
12,42514793
12,393904
12,40068091
12,48097718
12,49928316
12,48350739

Tables (Output of Computation)





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

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






Model: Y[t] = c + b X[t] + e[t]
c9.11790896447152
b0.274002359799523

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=51787&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]
c9.11790896447152
b0.274002359799523






Descriptive Statistics about e[t]
# observations69
minimum-0.0552174691474005
Q1-0.0234069895553007
median-0.00324978433310424
mean-9.48983796920198e-20
Q30.0252410741502926
maximum0.0742576927906003

\begin{tabular}{lllllllll}
\hline
Descriptive Statistics about e[t] \tabularnewline
# observations & 69 \tabularnewline
minimum & -0.0552174691474005 \tabularnewline
Q1 & -0.0234069895553007 \tabularnewline
median & -0.00324978433310424 \tabularnewline
mean & -9.48983796920198e-20 \tabularnewline
Q3 & 0.0252410741502926 \tabularnewline
maximum & 0.0742576927906003 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=51787&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.0552174691474005[/C][/ROW]
[ROW][C]Q1[/C][C]-0.0234069895553007[/C][/ROW]
[ROW][C]median[/C][C]-0.00324978433310424[/C][/ROW]
[ROW][C]mean[/C][C]-9.48983796920198e-20[/C][/ROW]
[ROW][C]Q3[/C][C]0.0252410741502926[/C][/ROW]
[ROW][C]maximum[/C][C]0.0742576927906003[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=51787&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=51787&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.0552174691474005
Q1-0.0234069895553007
median-0.00324978433310424
mean-9.48983796920198e-20
Q30.0252410741502926
maximum0.0742576927906003


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