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 computationTue, 21 Dec 2010 14:30:35 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2010/Dec/21/t1292942150bpeg3gpxpbwqx6k.htm/, Retrieved Wed, 08 May 2024 21:08:02 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=113621, Retrieved Wed, 08 May 2024 21:08:02 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsSpaargeld Huishoudens NL vs Uitgegeven Leningen Banken NL
Estimated Impact128

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]
-  M D  [Linear Regression Graphical Model Validation] [ws6vb] [2010-11-13 11:31:05] [c7506ced21a6c0dca45d37c8a93c80e0]
-    D    [Linear Regression Graphical Model Validation] [tt] [2010-11-13 12:51:28] [4a7069087cf9e0eda253aeed7d8c30d6]
-    D      [Linear Regression Graphical Model Validation] [W6tutorial] [2010-11-13 14:33:02] [c7506ced21a6c0dca45d37c8a93c80e0]
-    D          [Linear Regression Graphical Model Validation] [Spaargeld NL vs L...] [2010-12-21 14:30:35] [a8abc7260f3c847aeac0a796e7895a2e] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
143827
145191
146832
148577
149873
151847
153252
154292
155657
156523
156416
156693
160312
160438
160882
161668
164391
168556
169738
170387
171294
172202
172651
172770
178366
180014
181067
182586
184957
186417
188599
189490
190264
191221
191110
190674
195438
196393
197172
198760
200945
203845
204613
205487
206100
206315
206291
207801
211653
211325
211893
212056
214696
217455
218884
219816
219984
219062
218550
218179
222218
222196
223393
223292
226236
228831
228745
229140
229270
229359
230006
228810
232677
232961
234629
235660
240024
243554
244368
244356
245126
246321
246797
246735
251083
251786
252732
255051
259022
261698
263891
265247
262228
263429
264305
266371
273248
275472
278146
279506
283991
286794
288703
289285
288869
286942
285833
284095
289229
289389
290793
291454
294733
293853
294056
293982
293075
292391
Dataseries Y:
Download CSV
Histogram 
829461
837669
854793
850092
848783
846150
828543
830389
848989
841106
854616
832714
839290
840572
869186
856979
872126
868281
862455
881177
886924
886842
916407
890606
900409
920169
922871
920004
945772
937507
941691
958256
963509
970266
972853
982168
999892
1002099
1017611
1029782
1047956
1047689
1060054
1067078
1072366
1081823
1087601
1089905
1116316
1111355
1124250
1140597
1151683
1137532
967532
972994
999207
1007982
1015892
994850
987503
986743
1020674
1024067
1040444
1019081
1027828
1021010
1025563
1044756
1062545
1070425
1100087
1093596
1109143
1113855
1129275
1131996
1144103
1167830
1153194
1175008
1175805
1173456
1187498
1202958
1206229
1249533
1279743
1283496
1282942
1284739
1337169
1314087
1306144
1200391
1265445
1259329
1219342
1227626
1232874
1241046
1244172
1237838
1212801
1234237
1224699
1237432
1248847
1256543
1252434
1265176
1314670
1299329
1216744
1225275
1193478
1207226

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

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






Simple Linear Regression
StatisticsEstimateS.D.T-STAT (H0: coeff=0)P-value (two-sided)
constant term371942.19331974423555.432517587015.79008124949710
slope3.130181538589950.10425239987236530.0250310057340

\begin{tabular}{lllllllll}
\hline
Simple Linear Regression \tabularnewline
Statistics & Estimate & S.D. & T-STAT (H0: coeff=0) & P-value (two-sided) \tabularnewline
constant term & 371942.193319744 & 23555.4325175870 & 15.7900812494971 & 0 \tabularnewline
slope & 3.13018153858995 & 0.104252399872365 & 30.025031005734 & 0 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=113621&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]371942.193319744[/C][C]23555.4325175870[/C][C]15.7900812494971[/C][C]0[/C][/ROW]
[ROW][C]slope[/C][C]3.13018153858995[/C][C]0.104252399872365[/C][C]30.025031005734[/C][C]0[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=113621&T=1

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


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