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 computationThu, 22 Dec 2011 15:28:15 -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/2011/Dec/22/t1324585970nph3cz84q2a65c1.htm/, Retrieved Mon, 29 Apr 2024 16:09:39 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=159953, Retrieved Mon, 29 Apr 2024 16:09:39 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Linear Regression Graphical Model Validation] [Colombia Coffee -...] [2008-02-26 10:22:06] [74be16979710d4c4e7c6647856088456]
- RM D  [Linear Regression Graphical Model Validation] [Regression model 1] [2011-11-15 15:47:03] [74b1e5a3104ff0b2404b2865a63336ad]
-   PD    [Linear Regression Graphical Model Validation] [9-plot] [2011-11-15 20:19:09] [74b1e5a3104ff0b2404b2865a63336ad]
-   PD        [Linear Regression Graphical Model Validation] [Enkelvoudige Line...] [2011-12-22 20:28:15] [f9bdb25068ab2a4592adc645515299ca] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
140824
110459
105079
112098
43929
76173
187326
22807
144408
66485
79089
81625
68788
103297
69446
114948
167949
125081
125818
136588
112431
103037
82317
118906
83515
104581
103129
83243
37110
113344
139165
86652
112302
69652
119442
69867
101629
70168
31081
103925
92622
79011
93487
64520
93473
114360
33032
96125
151911
89256
95676
5950
149695
32551
31701
100087
169707
150491
120192
95893
151715
176225
59900
104767
114799
72128
143592
89626
131072
126817
81351
22618
88977
92059
81897
108146
126372
249771
71154
71571
55918
160141
38692
102812
56622
15986
123534
108535
93879
144551
56750
127654
65594
59938
146975
165904
169265
183500
165986
184923
140358
149959
57224
43750
48029
104978
100046
101047
197426
160902
147172
109432
1168
83248
25162
45724
110529
855
101382
14116
89506
135356
116066
144244
8773
102153
117440
104128
134238
134047
279488
79756
66089
102070
146760
154771
165933
64593
92280
67150
128692
124089
125386
37238
140015
150047
154451
156349
0
6023
0
0
0
0
84601
68946
0
0
1644
6179
3926
52789
0
100350
Dataseries Y:
Download CSV
Histogram 
279055
209884
233939
222117
179751
70849
568125
33186
227332
258874
351915
260484
204003
368577
269455
395936
335567
423110
182016
267365
279428
508849
206722
200004
257139
270815
296850
307100
184160
393860
327660
252512
373013
115602
430118
273950
428077
251349
115658
388812
343783
207021
214344
182398
157164
459455
78800
217932
368086
215843
244765
24188
399093
65029
101097
300488
369627
367127
374193
270099
391871
315924
291391
295075
276201
267432
215924
256641
260919
182961
256967
73566
272362
220707
228835
371391
398210
220401
229333
217623
200046
483074
145943
295224
80953
180759
179344
415550
369093
180679
299505
292260
199481
282361
329281
234577
297995
305984
416463
414359
297080
318283
222281
43287
223456
258249
299566
321797
174736
169579
354041
303273
23668
196743
61857
207339
431443
21054
252805
31961
360401
251240
187003
180842
38214
278173
358276
211775
445926
348017
441946
210700
126320
316128
466139
162279
412099
173802
292443
283913
243609
387072
246963
173260
346748
176654
264767
314070
1
14688
98
455
0
0
284420
410509
0
203
7199
46660
17547
121550
969
242258

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=159953&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'AstonUniversity' @ aston.wessa.net






Simple Linear Regression
StatisticsEstimateS.D.T-STAT (H0: coeff=0)P-value (two-sided)
constant term84571.05046732914380.8279635965.880819288111542.26738112729663e-08
slope1.682931165543520.1328601628025312.66693589744320

\begin{tabular}{lllllllll}
\hline
Simple Linear Regression \tabularnewline
Statistics & Estimate & S.D. & T-STAT (H0: coeff=0) & P-value (two-sided) \tabularnewline
constant term & 84571.050467329 & 14380.827963596 & 5.88081928811154 & 2.26738112729663e-08 \tabularnewline
slope & 1.68293116554352 & 0.13286016280253 & 12.6669358974432 & 0 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=159953&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]84571.050467329[/C][C]14380.827963596[/C][C]5.88081928811154[/C][C]2.26738112729663e-08[/C][/ROW]
[ROW][C]slope[/C][C]1.68293116554352[/C][C]0.13286016280253[/C][C]12.6669358974432[/C][C]0[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=159953&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=159953&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 term84571.05046732914380.8279635965.880819288111542.26738112729663e-08
slope1.682931165543520.1328601628025312.66693589744320


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 = 1 ; par2 = 2 ; par3 = Pearson Chi-Squared ;
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')