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, 03 Nov 2009 15:53:47 -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/03/t1257288895gelbq1uiton6b7h.htm/, Retrieved Wed, 01 May 2024 22:41:46 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=53465, Retrieved Wed, 01 May 2024 22:41:46 +0000
QR Codes:

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

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] [] [2009-11-03 22:53:47] [54e293c1fb7c46e2abc5c1dda68d8adb] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
258881,9875
255059,9883
247697,9873
244778,9879
240900,9888
239932,9876
270246,9882
283892,9885
282347,986
273569,9859
254755,9873
254353,9863
255842,9886
254489,9882
251994,9864
246338,9879
244018,9884
245952,9859
279805,9879
283110,9877
281096,9855
275963,9864
270693,9865
271900,9855
274411,9877
272432,9881
268360,9865
268585,9875
264767,9875
269973,9865
304743,989
309364,9879
308346,9857
298426,9868
289230,9865
291974,9869
294911,9875
293487,9873
290554,9853
284735,9877
281817,9866
287853,986
316262,9885
325411,9873
326010,9856
328281,9854
317479,9858
317538,9865
313736,987
312275,9872
309390,9854
302949,9874
300315,9868
304034,9857
333475,9874
337697,9871
335931,9864
323930,9852
313926,9859
314484,9864
313217,9867
309663,9868
302962,9864
298988,9856
298422,9873
301630,986
329764,987
335082,9882
327615,9853
309118,9855
295915,9872
291412,9859
291458,4609
284594,461
276391,4593
272482,4605
264897,461
263206,4597
296722,4612
Dataseries Y:
Download CSV
Histogram 
208069,33
204924,3308
200204,3298
198002,3304
195100,3314
191309,3301
213682,3307
225694,331
230493,3285
229175,3283
216142,3298
216627,3287
220120,3312
220029,3308
218358,3289
214826,3304
210619,331
209587,3283
236955,3304
241995,3303
241792,3279
242535,3288
238459,329
240251,3279
244666,3303
244490,3306
241486,3289
240455,33
236003,33
236911,3289
264493,3316
270263,3304
269559,3281
266951,3293
258027,3289
262727,3294
267327,33
267280,3298
264691,3277
258777,3302
254758,3291
254782,3285
277181,3311
285265,3298
286516,328
282956,3278
276601,3282
277829,329
277042,3294
277017,3297
274951,3278
270064,3299
267054,3293
264907,3281
287173,3299
291100,3296
292214,3288
288100,3275
281391,3283
282570,3289
280104,3292
280322,3292
276750,3289
275130,328
274266,3298
271225,3284
289716,3295
290640,3307
292214,3276
278420,3278
269740,3297
265775,3283
268948,3304
264090,3305
255112,3287
253267,33
245971,3306
235286,3291
258470,3308

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

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

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

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

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






Descriptive Statistics about e[t]
# observations79
minimum-22994.3654293957
Q1-4538.09274498655
median1308.28316866395
mean-2.4891200212096e-13
Q36178.15265032349
maximum16123.7157378378

\begin{tabular}{lllllllll}
\hline
Descriptive Statistics about e[t] \tabularnewline
# observations & 79 \tabularnewline
minimum & -22994.3654293957 \tabularnewline
Q1 & -4538.09274498655 \tabularnewline
median & 1308.28316866395 \tabularnewline
mean & -2.4891200212096e-13 \tabularnewline
Q3 & 6178.15265032349 \tabularnewline
maximum & 16123.7157378378 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=53465&T=2

[TABLE]
[ROW][C]Descriptive Statistics about e[t][/C][/ROW]
[ROW][C]# observations[/C][C]79[/C][/ROW]
[ROW][C]minimum[/C][C]-22994.3654293957[/C][/ROW]
[ROW][C]Q1[/C][C]-4538.09274498655[/C][/ROW]
[ROW][C]median[/C][C]1308.28316866395[/C][/ROW]
[ROW][C]mean[/C][C]-2.4891200212096e-13[/C][/ROW]
[ROW][C]Q3[/C][C]6178.15265032349[/C][/ROW]
[ROW][C]maximum[/C][C]16123.7157378378[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=53465&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=53465&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]
# observations79
minimum-22994.3654293957
Q1-4538.09274498655
median1308.28316866395
mean-2.4891200212096e-13
Q36178.15265032349
maximum16123.7157378378


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