Free Statistics

of Irreproducible Research!

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, 17 Nov 2011 14:33:39 -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/Nov/17/t13215584740wmthgyzs75pmrx.htm/, Retrieved Fri, 29 Mar 2024 14:04:18 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=145252, Retrieved Fri, 29 Mar 2024 14:04:18 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact52
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Linear Regression Graphical Model Validation] [] [2011-11-17 19:33:39] [1118fb1265e4c78f2f623b6bb1012fba] [Current]
Feedback Forum

Post a new message
Dataseries X:
9256885
10333752
10720612
11907834
9690590
10609448
11977254
13244600
11240162
11576825
11807926
13588068
11691588
12562273
12438117
12928814
11548838
12776845
11997156
13330088
11435231
12572825
11345956
13470822
10775234
12519177
11934407
13812730
11072269
12875557
12839205
14102759
10644606
13551116
12790714
13436225
10543361
13588113
13910100
13688184
9460957
12983955
12072813
13410011
11595556
12711511
13192150
11446149
11595556
12188025
13224051
13065877
9609338
12607393
12501842
12689869
10650116
12591209
13880292
13514738
9758680
11408706
11300823
11987461
9644061
11129888
11516053
10674049
8831349
10896072
10658420
11588748
9034816
10147881
11040127
11968843
9204166
9828988
10599002
10383324
9344956
10236895
9834064
10900904
9573817
10054837
9099351
10830753
8890756
10240704
9407973
10746225
8685878
10120696
9299008
8928589
7499163
8814962
9215177
10158324
8373334
9485738
9333946
10910330
8816168
8923734
Dataseries Y:
420535353
488878453
498041239
564072495
457511553
534305471
580509081
648770437
585105127
630821361
635048475
684750205
624579278
708146235
679748791
727949555
711704976
823818767
751357794
852585004
769520139
832010662
718771367
872040833
743812598
818570610
886246705
976236126
819407290
1026133596
953889829
1060970548
831565599
1082691438
1037403716
1100450194
877855090
1190802579
1145401598
1124835873
828092945
1113564364
1079819494
1220707994
1027966674
1255047042
1229675334
1039095435
1027966674
1243925959
1291053814
1369959521
1032152788
1312286401
1353955493
1423890029
1198277525
1459407892
1612018595
1528714149
1169666679
1403102002
1420269219
1510080957
1275608999
1471568598
1483636755
1482944107
1274746307
1636335358
1618764870
1787598977
1436519436
1673395701
1798072680
1944544177
1530417042
1726170719
1781938866
1819651590
2023539182
2310956246
2296138765
2489571272
2342282823
2622730494
2443617370
2758409400
2534556428
2871825467
2811245838
2962815256
2623492171
3127575470
2817693911
2647347460
2229602701
2625109728
2803549524
2950848868
2646481372
2928423819
2993317766
3259160572
2879505487
2938057843




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'George Udny Yule' @ yule.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 & 'George Udny Yule' @ yule.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=145252&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]'George Udny Yule' @ yule.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=145252&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=145252&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 Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'George Udny Yule' @ yule.wessa.net







Simple Linear Regression
StatisticsEstimateS.D.T-STAT (H0: coeff=0)P-value (two-sided)
constant term4519092466.02605459873337.3474449.826819906742682.22044604925031e-16
slope-270.24655896521340.5519258069166-6.66421023386071.30206267989763e-09

\begin{tabular}{lllllllll}
\hline
Simple Linear Regression \tabularnewline
Statistics & Estimate & S.D. & T-STAT (H0: coeff=0) & P-value (two-sided) \tabularnewline
constant term & 4519092466.02605 & 459873337.347444 & 9.82681990674268 & 2.22044604925031e-16 \tabularnewline
slope & -270.246558965213 & 40.5519258069166 & -6.6642102338607 & 1.30206267989763e-09 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=145252&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]4519092466.02605[/C][C]459873337.347444[/C][C]9.82681990674268[/C][C]2.22044604925031e-16[/C][/ROW]
[ROW][C]slope[/C][C]-270.246558965213[/C][C]40.5519258069166[/C][C]-6.6642102338607[/C][C]1.30206267989763e-09[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=145252&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=145252&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
StatisticsEstimateS.D.T-STAT (H0: coeff=0)P-value (two-sided)
constant term4519092466.02605459873337.3474449.826819906742682.22044604925031e-16
slope-270.24655896521340.5519258069166-6.66421023386071.30206267989763e-09



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