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, 05 Nov 2009 08:12:25 -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/05/t1257434047h2873ns9ggocw05.htm/, Retrieved Thu, 02 May 2024 17:00:11 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=54130, Retrieved Thu, 02 May 2024 17:00:11 +0000
QR Codes:

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

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-04 19:36:16] [3425351e86519d261a643e224a0c8ee1]
-    D    [Bivariate Explorative Data Analysis] [] [2009-11-05 15:12:25] [17416e80e7873ecccac25c455c5f767e] [Current]
- RMPD      [Partial Correlation] [] [2009-11-05 15:35:37] [3425351e86519d261a643e224a0c8ee1]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
3278,397949
3153,10596
3161,180194
3130,213256
3118,064786
3081,668613
3000,397031
2956,091623
2971,220093
3357,120296
3416,494173
3409,204939
3290,467184
3164,480194
3179,278664
3149,671725
3127,894021
3053,761674
3009,816266
3002,797031
3002,777031
3351,371061
3395,956469
3351,201061
3113,895551
2951,152389
2884,609276
2914,296215
2825,955399
2700,008409
2656,033001
2544,764481
2455,843665
2871,210807
2945,173154
2774,860756
2664,282236
2539,325247
2569,262185
2591,929889
2526,026777
2415,688257
2393,710553
2239,106625
2306,409737
2677,461471
2721,226879
2602,859124
2513,748308
2483,441369
2586,930655
2675,951471
2713,377644
2734,965348
2742,604583
2646,204532
2734,905348
3112,706317
3187,138664
3068,300909
Dataseries Y:
Download CSV
Histogram 
-19,05358195
-43,46849122
-49,82643774
-42,99465168
-40,97875866
-28,48902609
-39,65161444
-52,80393537
-37,91982839
15,96695292
35,10338081
33,56132732
10,08847154
-12,72643774
-5,142330764
7,389455289
1,563294829
-4,657240038
8,990439042
13,14838556
9,648385555
65,72489943
69,07722035
40,22489943
2,479187856
-16,04598885
-15,82447023
-28,85625629
-29,26089813
-36,3758074
-29,02812832
-42,85893062
-46,76357246
18,79142279
18,91195766
12,74472747
14,11392517
19,99901589
19,46722984
34,3933903
42,81490892
50,18410662
61,25794616
76,97482294
109,3533043
162,9559787
117,7082996
61,83544379
-31,96919806
-84,137412
-102,2486632
-87,24402135
-88,33375391
-76,90759345
-65,46553997
-56,11223529
-16,90759345
21,13713437
50,25766924
26,18481345

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24

\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 & 4 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=54130&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]4 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=54130&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=54130&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 time4 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24






Model: Y[t] = c + b X[t] + e[t]
c-1.49548939446923
b0.000518073906807889

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

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

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






Descriptive Statistics about e[t]
# observations60
minimum-102.093395076608
Q1-38.400376062162
median1.89987658244984
mean5.96860523967753e-16
Q327.8906713021142
maximum163.064345169861

\begin{tabular}{lllllllll}
\hline
Descriptive Statistics about e[t] \tabularnewline
# observations & 60 \tabularnewline
minimum & -102.093395076608 \tabularnewline
Q1 & -38.400376062162 \tabularnewline
median & 1.89987658244984 \tabularnewline
mean & 5.96860523967753e-16 \tabularnewline
Q3 & 27.8906713021142 \tabularnewline
maximum & 163.064345169861 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=54130&T=2

[TABLE]
[ROW][C]Descriptive Statistics about e[t][/C][/ROW]
[ROW][C]# observations[/C][C]60[/C][/ROW]
[ROW][C]minimum[/C][C]-102.093395076608[/C][/ROW]
[ROW][C]Q1[/C][C]-38.400376062162[/C][/ROW]
[ROW][C]median[/C][C]1.89987658244984[/C][/ROW]
[ROW][C]mean[/C][C]5.96860523967753e-16[/C][/ROW]
[ROW][C]Q3[/C][C]27.8906713021142[/C][/ROW]
[ROW][C]maximum[/C][C]163.064345169861[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=54130&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=54130&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]
# observations60
minimum-102.093395076608
Q1-38.400376062162
median1.89987658244984
mean5.96860523967753e-16
Q327.8906713021142
maximum163.064345169861


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