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_linear_regression.wasp
Title produced by softwareLinear Regression Graphical Model Validation
Date of computationSat, 24 Nov 2012 19:32:20 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2012/Nov/24/t1353803628b4rtux1d4xa2tqh.htm/, Retrieved Sun, 28 Apr 2024 19:57:20 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=192545, Retrieved Sun, 28 Apr 2024 19:57:20 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Bivariate Data Series] [Bivariate dataset] [2008-01-05 23:51:08] [74be16979710d4c4e7c6647856088456]
- RMPD  [Blocked Bootstrap Plot - Central Tendency] [Colombia Coffee] [2008-01-07 10:26:26] [74be16979710d4c4e7c6647856088456]
-  M D    [Blocked Bootstrap Plot - Central Tendency] [Mini-tutorial - B...] [2010-11-15 22:35:54] [6f0e7a2d1a07390e3505a2db8288f975]
-   PD      [Blocked Bootstrap Plot - Central Tendency] [Mini-tutorial - B...] [2010-11-15 23:41:51] [6f0e7a2d1a07390e3505a2db8288f975]
- RMPD        [Linear Regression Graphical Model Validation] [Workshop 6 - Regr...] [2010-11-16 19:26:47] [8b017ffbf7b0eded54d8efebfb3e4cfa]
-    D          [Linear Regression Graphical Model Validation] [Simple Linear Reg...] [2010-11-28 10:16:25] [8b017ffbf7b0eded54d8efebfb3e4cfa]
-    D            [Linear Regression Graphical Model Validation] [Simple Linear Reg...] [2010-11-28 13:21:55] [8b017ffbf7b0eded54d8efebfb3e4cfa]
- R  D                [Linear Regression Graphical Model Validation] [Linear Regression] [2012-11-25 00:32:20] [4289cf790da1cc09a0cb8798de13fde3] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
1579
2146
2462
3695
4831
5134
6250
5760
6249
2917
1741
2359
1511
2059
2635
2867
4403
5720
4502
5749
5627
2846
1762
2429
1169
2154
2249
2687
4359
5382
4459
6398
4596
3024
1887
2070
1351
2218
2461
3028
4784
4975
4607
6249
4809
3157
1910
2228
1594
2467
2222
3607
4685
4962
5770
5480
5000
3228
1993
2288
1580
2111
2192
3601
4665
4876
5813
5589
5331
3075
2002
2306
1507
1992
2487
3490
4647
5594
5611
5788
6204
3013
1931
2549
1504
2090
2702
2939
4500
6208
6415
5657
5964
3163
1997
2422
1.376
2.202
2.683
3.303
5.202
5.231
4.880
7.998
4.977
3.531
2.025
2.205
1.442
2.238
2.179
3.218
5.139
4.990
4.914
6.084
5.672
3.548
1.793
2.086
Dataseries Y:
Download CSV
Histogram 
9.769
9.321
9.939
9.336
10.195
9.464
10.010
10.213
9.563
9.890
9.305
9.391
9.928
8.686
9.843
9.627
10.074
9.503
10.119
10.000
9.313
9.866
9.172
9.241
9.659
8.904
9.755
9.080
9.435
8.971
10.063
9.793
9.454
9.759
8.820
9.403
9.676
8.642
9.402
9.610
9.294
9.448
10.319
9.548
9.801
9.596
8.923
9.746
9.829
9.125
9.782
9.441
9.162
9.915
10.444
10.209
9.985
9.842
9.429
10.132
9.849
9.172
10.313
9.819
9.955
10.048
10.082
10.541
10.208
10.233
9.439
9.963
10.158
9.225
10.474
9.757
10.490
10.281
10.444
10.640
10.695
10.786
9.832
9.747
10.411
9.511
10.402
9.701
10.540
10.112
10.915
11.183
10.384
10.834
9.886
10.216
10.943
9.867
10.203
10.837
10.573
10.647
11.502
10.656
10.866
10.835
9.945
10.331
10.718
9.462
10.579
10.633
10.346
10.757
11.207
11.013
11.015
10.765
10.042
10.661

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time7 seconds
R Server'Sir Maurice George Kendall' @ kendall.wessa.net

\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 & 7 seconds \tabularnewline
R Server & 'Sir Maurice George Kendall' @ kendall.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=192545&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]7 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Maurice George Kendall' @ kendall.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=192545&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=192545&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 time7 seconds
R Server'Sir Maurice George Kendall' @ kendall.wessa.net






Simple Linear Regression
StatisticsEstimateS.D.T-STAT (H0: coeff=0)P-value (two-sided)
constant term10.1171616006650.0929290753095364108.8697113036570
slope-4.97955681710431e-052.63604557014289e-05-1.889025316369780.0613430141474736

\begin{tabular}{lllllllll}
\hline
Simple Linear Regression \tabularnewline
Statistics & Estimate & S.D. & T-STAT (H0: coeff=0) & P-value (two-sided) \tabularnewline
constant term & 10.117161600665 & 0.0929290753095364 & 108.869711303657 & 0 \tabularnewline
slope & -4.97955681710431e-05 & 2.63604557014289e-05 & -1.88902531636978 & 0.0613430141474736 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=192545&T=1

[TABLE]
[ROW][C]Simple Linear Regression[/C][/ROW]
[ROW][C]Statistics[/C][C]Estimate[/C][C]S.D.[/C][C]T-STAT (H0: coeff=0)[/C][C]P-value (two-sided)[/C][/ROW]
[ROW][C]constant term[/C][C]10.117161600665[/C][C]0.0929290753095364[/C][C]108.869711303657[/C][C]0[/C][/ROW]
[ROW][C]slope[/C][C]-4.97955681710431e-05[/C][C]2.63604557014289e-05[/C][C]-1.88902531636978[/C][C]0.0613430141474736[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=192545&T=1

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

Simple Linear Regression
StatisticsEstimateS.D.T-STAT (H0: coeff=0)P-value (two-sided)
constant term10.1171616006650.0929290753095364108.8697113036570
slope-4.97955681710431e-052.63604557014289e-05-1.889025316369780.0613430141474736


Figures (Output of Computation)

Figure 1
PNG link  
QR Code
Postscript link  
QR Code
PDF link  
QR Code


Figure 2
PNG link  
QR Code
Postscript link  
QR Code
PDF link  
QR Code


Figure 3
PNG link  
QR Code
Postscript link  
QR Code
PDF link  
QR Code


Figure 4
PNG link  
QR Code
Postscript link  
QR Code
PDF link  
QR Code


Figure 5
PNG link  
QR Code
Postscript link  
QR Code
PDF link  
QR Code


Figure 6
PNG link  
QR Code
Postscript link  
QR Code
PDF link  
QR Code


Figure 7
PNG link  
QR Code
Postscript link  
QR Code
PDF link  
QR Code


Figure 8
PNG link  
QR Code
Postscript link  
QR Code
PDF link  
QR Code


Figure 9
PNG link  
QR Code
Postscript link  
QR Code
PDF link  
QR Code


Input Parameters & R Code

Parameters (Session):
par1 = 0 ;
Parameters (R input):
par1 = 0 ;
R code (references can be found in the software module):
par1 <- as.numeric(par1)
library(lattice)
z <- as.data.frame(cbind(x,y))
m <- lm(y~x)
summary(m)
bitmap(file='test1.png')
plot(z,main='Scatterplot, lowess, and regression line')
lines(lowess(z),col='red')
abline(m)
grid()
dev.off()
bitmap(file='test2.png')
m2 <- lm(m$fitted.values ~ x)
summary(m2)
z2 <- as.data.frame(cbind(x,m$fitted.values))
names(z2) <- list('x','Fitted')
plot(z2,main='Scatterplot, lowess, and regression line')
lines(lowess(z2),col='red')
abline(m2)
grid()
dev.off()
bitmap(file='test3.png')
m3 <- lm(m$residuals ~ x)
summary(m3)
z3 <- as.data.frame(cbind(x,m$residuals))
names(z3) <- list('x','Residuals')
plot(z3,main='Scatterplot, lowess, and regression line')
lines(lowess(z3),col='red')
abline(m3)
grid()
dev.off()
bitmap(file='test4.png')
m4 <- lm(m$fitted.values ~ m$residuals)
summary(m4)
z4 <- as.data.frame(cbind(m$residuals,m$fitted.values))
names(z4) <- list('Residuals','Fitted')
plot(z4,main='Scatterplot, lowess, and regression line')
lines(lowess(z4),col='red')
abline(m4)
grid()
dev.off()
bitmap(file='test5.png')
myr <- as.ts(m$residuals)
z5 <- as.data.frame(cbind(lag(myr,1),myr))
names(z5) <- list('Lagged Residuals','Residuals')
plot(z5,main='Lag plot')
m5 <- lm(z5)
summary(m5)
abline(m5)
grid()
dev.off()
bitmap(file='test6.png')
hist(m$residuals,main='Residual Histogram',xlab='Residuals')
dev.off()
bitmap(file='test7.png')
if (par1 > 0)
{
densityplot(~m$residuals,col='black',main=paste('Density Plot   bw = ',par1),bw=par1)
} else {
densityplot(~m$residuals,col='black',main='Density Plot')
}
dev.off()
bitmap(file='test8.png')
acf(m$residuals,main='Residual Autocorrelation Function')
dev.off()
bitmap(file='test9.png')
qqnorm(x)
qqline(x)
grid()
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Simple Linear Regression',5,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Statistics',1,TRUE)
a<-table.element(a,'Estimate',1,TRUE)
a<-table.element(a,'S.D.',1,TRUE)
a<-table.element(a,'T-STAT (H0: coeff=0)',1,TRUE)
a<-table.element(a,'P-value (two-sided)',1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'constant term',header=TRUE)
a<-table.element(a,m$coefficients[[1]])
sd <- sqrt(vcov(m)[1,1])
a<-table.element(a,sd)
tstat <- m$coefficients[[1]]/sd
a<-table.element(a,tstat)
pval <- 2*(1-pt(abs(tstat),length(x)-2))
a<-table.element(a,pval)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'slope',header=TRUE)
a<-table.element(a,m$coefficients[[2]])
sd <- sqrt(vcov(m)[2,2])
a<-table.element(a,sd)
tstat <- m$coefficients[[2]]/sd
a<-table.element(a,tstat)
pval <- 2*(1-pt(abs(tstat),length(x)-2))
a<-table.element(a,pval)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')