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 computationMon, 02 Nov 2009 09:30:30 -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/Nov/02/t125717950768u28y31cal70pl.htm/, Retrieved Sat, 04 May 2024 00:22:21 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=52765, Retrieved Sat, 04 May 2024 00:22:21 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Univariate Data Series] [WS2] [2009-10-12 17:24:11] [4f76e114ed5e444b1133aad392380aad]
- RMPD  [Bivariate Explorative Data Analysis] [] [2009-10-28 18:11:46] [4f76e114ed5e444b1133aad392380aad]
-  M D      [Bivariate Explorative Data Analysis] [ln op reeksen toe...] [2009-11-02 16:30:30] [ea241b681aafed79da4b5b99fad98471] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
11.57467886
11.60253762
11.66544052
11.72235136
11.65106
11.67780411
11.74772878
11.78474489
11.6997867
11.72634994
11.77385023
11.80157956
11.71983394
11.73940744
11.77681933
11.7982247
11.70729087
11.71213691
11.76222787
11.7871879
11.695472
11.71994776
11.77614321
11.79451904
11.72595403
11.73551686
11.80294979
11.82449739
11.75594223
11.76714124
11.82912939
11.84736707
11.76821089
11.78755269
11.84519603
11.87147365
11.80520694
11.82858256
11.88355065
11.91355892
11.84076051
11.87240248
11.91861721
11.94929475
11.88371631
11.91010299
11.96081129
11.98970608
11.93670189
11.97041993
12.0034545
12.04337135
11.99701073
12.0183264
12.04707691
12.07071673
12.02193344
12.04496449
12.07909118
12.11206961
12.0565164
12.08872951
12.11315111
12.14467585
12.0831307
12.11873212
12.14403233
12.17895835
12.11743813
12.15222873
12.17460778
12.20856454
12.14495228
12.17423104
12.1915984
12.23144795
12.16395293
12.189435
12.20499206
12.24221169
12.18638511
12.21018418
12.2343831
12.2666455
12.21011447
12.22432014
12.24185476
12.25918678
12.18031885
12.2918833
Dataseries Y:
Download CSV
Histogram 
4.608165695
4.62399194
4.625952725
4.636668853
4.641502115
4.665324109
4.679349584
4.699570861
4.715816706
4.744062185
4.753590191
4.773223771
4.790819533
4.836281907
4.852811209
4.868303386
4.873669439
4.894850261
4.899331225
4.90897164
4.914124394
4.935193099
4.934473933
4.93878119
4.932313327
4.948050419
4.95088529
4.954417614
4.955827058
4.973279508
4.97397131
4.980176087
4.989071116
5.007296393
5.009968405
5.011301739
5.016617366
5.029129877
5.031091303
5.036952602
5.042779747
5.05560866
5.065123479
5.073297055
5.076423035
5.094976443
5.098035484
5.102302483
5.098035484
5.108971195
5.109575242
5.11679549
5.120983351
5.139321635
5.14107859
5.147494477
5.14633101
5.158480421
5.159055299
5.157905213
5.158480421
5.170483995
5.173887288
5.182906515
5.188502501
5.200153118
5.202907182
5.208939555
5.21384821
5.227358278
5.23324537
5.242804657
5.245443877
5.256974403
5.260615499
5.26631057
5.268888556
5.286244785
5.294811227
5.30529298
5.313205979
5.32933142
5.333201768
5.346154696
5.352331955
5.372032404
5.381738975
5.37296091
5.351384087
5.359412466

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=52765&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=52765&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=52765&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]
c-6.67687199305052
b0.981718152841215

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=52765&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-6.67687199305052
b0.981718152841215






Descriptive Statistics about e[t]
# observations90
minimum-0.194504278044386
Q1-0.0351358846960069
median0.00642475872363349
mean2.98385237476992e-18
Q30.0415675771100746
maximum0.112876851873847

\begin{tabular}{lllllllll}
\hline
Descriptive Statistics about e[t] \tabularnewline
# observations & 90 \tabularnewline
minimum & -0.194504278044386 \tabularnewline
Q1 & -0.0351358846960069 \tabularnewline
median & 0.00642475872363349 \tabularnewline
mean & 2.98385237476992e-18 \tabularnewline
Q3 & 0.0415675771100746 \tabularnewline
maximum & 0.112876851873847 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=52765&T=2

[TABLE]
[ROW][C]Descriptive Statistics about e[t][/C][/ROW]
[ROW][C]# observations[/C][C]90[/C][/ROW]
[ROW][C]minimum[/C][C]-0.194504278044386[/C][/ROW]
[ROW][C]Q1[/C][C]-0.0351358846960069[/C][/ROW]
[ROW][C]median[/C][C]0.00642475872363349[/C][/ROW]
[ROW][C]mean[/C][C]2.98385237476992e-18[/C][/ROW]
[ROW][C]Q3[/C][C]0.0415675771100746[/C][/ROW]
[ROW][C]maximum[/C][C]0.112876851873847[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=52765&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=52765&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]
# observations90
minimum-0.194504278044386
Q1-0.0351358846960069
median0.00642475872363349
mean2.98385237476992e-18
Q30.0415675771100746
maximum0.112876851873847


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