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 computationThu, 10 Dec 2009 12:12:18 -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/10/t12604727839jr8pgf80kfggcw.htm/, Retrieved Fri, 19 Apr 2024 03:02:00 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=65742, Retrieved Fri, 19 Apr 2024 03:02:00 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Bivariate Kernel Density Estimation] [] [2009-12-10 17:31:58] [023d83ebdf42a2acf423907b4076e8a1]
- RMP     [Bivariate Explorative Data Analysis] [] [2009-12-10 19:12:18] [9f6463b67b1eb7bae5c03a796abf0348] [Current]
- RMPD      [Harrell-Davis Quantiles] [Harrell Davis 2] [2009-12-11 21:04:05] [830e13ac5e5ac1e5b21c6af0c149b21d]
- RMPD      [Harrell-Davis Quantiles] [Harrell Davis 2] [2009-12-11 21:04:05] [74be16979710d4c4e7c6647856088456]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
94.204
97.217
95.118
93.688
93.140
91.516
90.957
90.372
89.749
85.813
86.026
83.933
83.602
83.384
76.369
60.808
48.071
42.604
41.402
62.121
79.739
79.006
74.472
75.956
75.041
74.873
72.922
70.472
71.423
71.363
73.297
72.081
70.488
65.544
64.450
61.698
61.352
61.072
63.722
61.987
53.802
47.818
50.998
58.438
60.143
61.854
70.987
70.389
72.175
70.243
69.616
69.443
70.833
71.059
72.218
72.647
73.299
73.756
75.557
78.172
75.624
76.959
74.994
76.841
78.043
75.187
73.387
70.798
68.722
68.396
68.466
67.675
65.248
62.974
59.801
57.894
58.592
59.249
59.554
59.753
60.877
60.532
58.452
56.955
56.437
55.588
56.702
57.062
57.826
58.755
60.250
61.142
60.690
58.495
56.020
55.814
56.489
56.587
55.714
55.611
56.093
55.929
54.181
54.810
56.189
57.427
59.432
58.951
Dataseries Y:
Download CSV
Histogram 
99.026
101.851
99.958
97.875
97.927
95.149
94.551
93.999
93.297
89.901
89.742
87.096
86.863
86.718
80.020
63.483
51.289
44.071
43.654
66.115
84.518
83.395
78.307
80.049
78.346
78.317
75.918
73.739
74.530
74.179
76.974
75.408
73.336
69.210
67.286
64.606
64.159
64.423
66.411
64.270
56.521
50.599
54.751
62.227
63.932
65.391
75.744
74.590
76.035
74.427
73.354
73.081
75.309
75.463
75.910
76.151
76.882
78.632
80.137
82.490
79.896
81.303
79.344
81.355
82.328
79.669
77.249
75.101
72.520
72.438
72.653
71.429
69.189
66.451
63.354
61.379
61.880
62.274
62.429
63.905
63.917
64.295
61.930
60.440
59.353
58.695
60.569
60.386
60.938
61.795
63.304
64.270
63.492
61.333
59.341
58.412
58.725
59.277
58.562
57.858
58.790
58.243
57.044
57.339
59.429
60.575
61.950
61.712

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=65742&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]
c0.934888350894207
b1.03794723808285

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=65742&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.934888350894207
b1.03794723808285






Descriptive Statistics about e[t]
# observations108
minimum-1.08459248217577
Q1-0.396885529995194
median0.0240943138495157
mean1.10397882322436e-17
Q30.395795066496338
maximum1.14227515706736

\begin{tabular}{lllllllll}
\hline
Descriptive Statistics about e[t] \tabularnewline
# observations & 108 \tabularnewline
minimum & -1.08459248217577 \tabularnewline
Q1 & -0.396885529995194 \tabularnewline
median & 0.0240943138495157 \tabularnewline
mean & 1.10397882322436e-17 \tabularnewline
Q3 & 0.395795066496338 \tabularnewline
maximum & 1.14227515706736 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=65742&T=2

[TABLE]
[ROW][C]Descriptive Statistics about e[t][/C][/ROW]
[ROW][C]# observations[/C][C]108[/C][/ROW]
[ROW][C]minimum[/C][C]-1.08459248217577[/C][/ROW]
[ROW][C]Q1[/C][C]-0.396885529995194[/C][/ROW]
[ROW][C]median[/C][C]0.0240943138495157[/C][/ROW]
[ROW][C]mean[/C][C]1.10397882322436e-17[/C][/ROW]
[ROW][C]Q3[/C][C]0.395795066496338[/C][/ROW]
[ROW][C]maximum[/C][C]1.14227515706736[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=65742&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=65742&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]
# observations108
minimum-1.08459248217577
Q1-0.396885529995194
median0.0240943138495157
mean1.10397882322436e-17
Q30.395795066496338
maximum1.14227515706736


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