Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*The author of this computation has been verified*

R Software Module

rwasp_linear_regression.wasp

Title produced by software

Linear Regression Graphical Model Validation

Date of computation

Tue, 16 Nov 2010 20:19:12 +0000

Cite this page as follows

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2010/Nov/16/t1289938701y9we50yzbh93a40.htm/, Retrieved Thu, 02 May 2024 08:09:16 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=96397, Retrieved Thu, 02 May 2024 08:09:16 +0000

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords

Estimated Impact

123

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]
-  MPD  [Bivariate Data Series] [Schermbreedte en ...] [2010-11-11 17:54:48] [6bc4f9343b7ea3ef5a59412d1f72bb2b]
- RMPD      [Linear Regression Graphical Model Validation] [] [2010-11-16 20:19:12] [6b31f806e9ccc1f74a26091056f791cb] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

Dataseries Y:

Download CSV

Histogram

61.10986116
54.12031648
54.95749446
51.99594378
56.24546553
64.77452146
54.33699821
53.72119053
48.83342072
55.93379725
64.65918024
61.60053636
63.05659263
67.15166067
60.02845855
69.52301201
65.01551552
60.87993666
70.72557348
73.09285805
73.29160052
72.66843666
68.83636849
76.67414394
70.41408315
77.9401298
73.20981757
77.41384989
79.4630442
75.42057899
82.68120243
94.03299383
84.0698862
85.67449188
81.89949935
92.9720738
88.85764153
80.49481348
91.71796544
82.26348598
82.94116291
90.99841989
91.31788491
96.95308859
90.04825345
100.9221214
99.89956639
93.01828481
95.97742386
103.5617168
110.7469177
108.3696381
107.0012411
110.7524482
98.56007743
100.5808127
101.9883776
101.8791562
108.9640668
104.5917605
119.5561823
119.3482511
105.3661433
121.4417904
117.747784
116.6637445
119.6190875
115.2413454
106.7405737
113.7508379
124.9325395
124.9165153
134.0975866
132.0328883
124.4484784
115.0077499
122.6887566
125.0399372
119.2917291
131.6979609
132.3570185
136.9068086
147.3595766
135.5406186
137.5528819
128.6595377
135.701906
128.428847
134.1136873
133.3455986
139.8282064
152.4304094
148.7181493
142.8992652
152.7563365
139.7011827
146.5597201
146.5053758
145.5923464
140.0328201
154.0235601
154.9190928
158.101761
154.2821531
163.6434212
163.7318741
161.4377475
156.7145855
153.561049
153.8774688
169.2235043
156.3795439
165.490775
154.0411596
154.9816657
169.586015
160.6091489
160.6353333
163.6734875
166.1373986
175.3861465
181.8040035
172.5480095
175.5379656
177.4873292
178.8761947
178.375176
176.7772914
181.8905564
180.3690141
180.2465823
189.4586351
190.2788307
169.6704531
186.9818197
179.4104796
182.1682369
188.2755288
176.4747358
189.2416564
192.4300356
203.944427
203.3403548
188.1282826
196.354948
197.9087125
201.9367353
193.445003
198.8122889
196.4723325
204.5398687
205.0704692
204.3924263
207.9987893
200.9225201
214.2768736
208.0597487
208.5093111
201.22446
203.3434561
212.5706014
210.5406441
210.887768
209.6129575
209.313722
215.7808384
226.4811578
220.9966323
205.9677281
210.4273103
216.9371659
226.5364144
226.6246165
227.5501626
230.5976031
218.4966592
225.660519
219.7215854
224.2551636
223.4610044
241.6074588
220.6620768
231.2874518
239.7293117
235.6812993
239.2962581
230.904901
236.5510765
248.5657307
226.3445746
242.5403835
242.0843806
245.795533
244.5513355
244.957443
242.9284176
246.9975039
243.5141125
242.3741637
249.0753006
252.3467369
248.3848335
256.0723248
256.1781746
250.3571524
258.7296195
256.3099685
256.0184129
245.31261
248.7418796
263.5634464
268.8057802
259.432041
280.6897114
263.5020284
262.910557
254.8534255
262.914486
262.4803197
265.5528724
268.3350219
274.2269561
285.0879902
287.4846281
276.9574892
270.2328838
273.343679
282.3265576
271.3683696
276.7731473
284.984373
285.2336093
286.4188091
283.9944259
285.9733102
284.4558291
281.7885504
287.3979095
290.8197966
281.4066935
293.8918414
295.4459854
296.4138982
299.309
299.1331287
297.5508534
299.870089
298.1574167
293.7846644
289.1280973
302.7972545
311.6582427
306.8750536
300.9725842
309.210378
314.0951678
302.0571523
307.1139977
306.9729935
307.7220354
306.3991609
310.2296239
327.1883072
323.0723765
319.3126762
320.7490193
316.7516682
314.156943
318.148116
313.1329884
324.4617348
323.299863
323.5194599
325.766307
315.8600699
314.4862665
318.7478783
330.3323151
324.4215348
319.5462879
333.9837493
329.2936641
332.9469371
332.2435291
331.1389228
329.5469119
340.7498356
335.7682144
336.8801173
333.1690483
343.3672947
339.7428617
347.8462662
338.0019601
339.0321535
349.7254835
339.9088441
341.8686896
348.8204737
351.4017628
354.3571353
354.7671573
354.6799216
346.3211461
352.8674728
353.7544304
352.3393992
360.1987502
354.4125537
353.6133866
366.9320672
364.0824985
363.1392851
363.203798
358.7129274
357.8204599
366.2509726
365.6025863
370.8290372
366.5541817
374.8778657
374.9878806
375.9899812
377.0300208
376.3060598
381.8212896
380.7436544
370.7164768
368.0023819
376.9856375
387.8136742
383.166558
386.1462067
388.3909748
383.3338225
379.9889392
383.332495
394.2459788
391.6334594
390.6626585
385.1222833
397.8269152
394.4597945
395.180818
393.6250717
396.190475
395.3262294
392.0564838
402.4440354
393.7467618
394.3378194
410.4003906
398.2528185
407.9129263
408.5478406
403.4111597
398.4169705
406.8649612
406.8042436
405.1187658
417.866504
416.6827187
412.3299343
409.3072168
413.2115411
416.0710076
408.932371
404.8629156
413.4135646
406.9434141
426.4176986
420.0718548
430.8056774
419.9182
431.2393615
429.5679042
420.0486899
428.0971546
430.9662678
425.1897311
434.9097032
443.4206608
435.6845397
445.0445079
433.9489262
423.4262595
432.0529621
446.8513244
432.7851253
433.2442466
447.2221488
446.0185703
456.7722076
441.2170262
444.4915125
439.5461461
443.1393853
444.2187665
443.1337644
447.4276883
453.5096647
458.1323521
453.2685384
461.58333
454.3458605
460.8233618
456.7880646
451.2040437
450.7924225
450.3340228

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	8 seconds
R Server	'RServer@AstonUniversity' @ vre.aston.ac.uk

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=96397&T=0

As an alternative you can also use a QR Code:

The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	8 seconds
R Server	'RServer@AstonUniversity' @ vre.aston.ac.uk

Simple Linear Regression
Statistics	Estimate	S.D.	T-STAT (H0: coeff=0)	P-value (two-sided)
constant term	-44.4760358159129	0.66426852711803	-66.9549045306646	0
slope	0.099841303648979	0.000205980926142623	484.711402743417	0

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=96397&T=1

As an alternative you can also use a QR Code:

The GUIDs for individual cells are displayed in the table below:

Simple Linear Regression
Statistics	Estimate	S.D.	T-STAT (H0: coeff=0)	P-value (two-sided)
constant term	-44.4760358159129	0.66426852711803	-66.9549045306646	0
slope	0.099841303648979	0.000205980926142623	484.711402743417	0

Figure 1

PNG link

Postscript link

PDF link

Figure 2

PNG link

Postscript link

PDF link

Figure 3

PNG link

Postscript link

PDF link

Figure 4

PNG link

Postscript link

PDF link

Figure 5

PNG link

Postscript link

PDF link

Figure 6

PNG link

Postscript link

PDF link

Figure 7

PNG link

Postscript link

PDF link

Figure 8

PNG link

Postscript link

PDF link

Figure 9

PNG link

Postscript link

PDF link

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

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code