Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_edauni.wasp
Title produced by softwareUnivariate Explorative Data Analysis
Date of computationSun, 23 Nov 2008 05:17:49 -0700
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2008/Nov/23/t1227442707iflwyaj7sx6za42.htm/, Retrieved Mon, 29 Apr 2024 14:21:14 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=25229, Retrieved Mon, 29 Apr 2024 14:21:14 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsVan Dooren Leen
Estimated Impact218
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Multiple Regression] [Q1 The Seatbeltlaw] [2007-11-14 19:27:43] [8cd6641b921d30ebe00b648d1481bba0]
F    D  [Multiple Regression] [The Seatbelt Law ...] [2008-11-15 12:00:57] [93834488277b53a4510bfd06084ae13b]
- RMP       [Univariate Explorative Data Analysis] [Seatbelt Law Q2] [2008-11-23 12:17:49] [d175f84d503eb4f2a43145d5e67795b5] [Current]
Feedback Forum

Post a new message
Dataseries X:
1687	0
1508	0
1507	0
1385	0
1632	0
1511	0
1559	0
1630	0
1579	0
1653	0
2152	0
2148	0
1752	0
1765	0
1717	0
1558	0
1575	0
1520	0
1805	0
1800	0
1719	0
2008	0
2242	0
2478	0
2030	0
1655	0
1693	0
1623	0
1805	0
1746	0
1795	0
1926	0
1619	0
1992	0
2233	0
2192	0
2080	0
1768	0
1835	0
1569	0
1976	0
1853	0
1965	0
1689	0
1778	0
1976	0
2397	0
2654	0
2097	0
1963	0
1677	0
1941	0
2003	0
1813	0
2012	0
1912	0
2084	0
2080	0
2118	0
2150	0
1608	0
1503	0
1548	0
1382	0
1731	0
1798	0
1779	0
1887	0
2004	0
2077	0
2092	0
2051	0
1577	0
1356	0
1652	0
1382	0
1519	0
1421	0
1442	0
1543	0
1656	0
1561	0
1905	0
2199	0
1473	0
1655	0
1407	0
1395	0
1530	0
1309	0
1526	0
1327	0
1627	0
1748	0
1958	0
2274	0
1648	0
1401	0
1411	0
1403	0
1394	0
1520	0
1528	0
1643	0
1515	0
1685	0
2000	0
2215	0
1956	0
1462	0
1563	0
1459	0
1446	0
1622	0
1657	0
1638	0
1643	0
1683	0
2050	0
2262	0
1813	0
1445	0
1762	0
1461	0
1556	0
1431	0
1427	0
1554	0
1645	0
1653	0
2016	0
2207	0
1665	0
1361	0
1506	0
1360	0
1453	0
1522	0
1460	0
1552	0
1548	0
1827	0
1737	0
1941	0
1474	0
1458	0
1542	0
1404	0
1522	0
1385	0
1641	0
1510	0
1681	0
1938	0
1868	0
1726	0
1456	0
1445	0
1456	0
1365	0
1487	0
1558	0
1488	0
1684	0
1594	0
1850	0
1998	0
2079	0
1494	0
1057	1
1218	1
1168	1
1236	1
1076	1
1174	1
1139	1
1427	1
1487	1
1483	1
1513	1
1357	1
1165	1
1282	1
1110	1
1297	1
1185	1
1222	1
1284	1
1444	1
1575	1
1737	1
1763	1




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135

\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 & 1 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=25229&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]1 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=25229&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=25229&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 time1 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135



Parameters (Session):
par1 = 0 ; par2 = 0 ;
Parameters (R input):
par1 = 0 ; par2 = 0 ;
R code (references can be found in the software module):
par1 <- as.numeric(par1)
par2 <- as.numeric(par2)
x <- as.ts(x)
library(lattice)
bitmap(file='pic1.png')
plot(x,type='l',main='Run Sequence Plot',xlab='time or index',ylab='value')
grid()
dev.off()
bitmap(file='pic2.png')
hist(x)
grid()
dev.off()
bitmap(file='pic3.png')
if (par1 > 0)
{
densityplot(~x,col='black',main=paste('Density Plot bw = ',par1),bw=par1)
} else {
densityplot(~x,col='black',main='Density Plot')
}
dev.off()
bitmap(file='pic4.png')
qqnorm(x)
qqline(x)
grid()
dev.off()
if (par2 > 0)
{
bitmap(file='lagplot1.png')
dum <- cbind(lag(x,k=1),x)
dum
dum1 <- dum[2:length(x),]
dum1
z <- as.data.frame(dum1)
z
plot(z,main='Lag plot (k=1), lowess, and regression line')
lines(lowess(z))
abline(lm(z))
dev.off()
if (par2 > 1) {
bitmap(file='lagplotpar2.png')
dum <- cbind(lag(x,k=par2),x)
dum
dum1 <- dum[(par2+1):length(x),]
dum1
z <- as.data.frame(dum1)
z
mylagtitle <- 'Lag plot (k='
mylagtitle <- paste(mylagtitle,par2,sep='')
mylagtitle <- paste(mylagtitle,'), and lowess',sep='')
plot(z,main=mylagtitle)
lines(lowess(z))
dev.off()
}
bitmap(file='pic5.png')
acf(x,lag.max=par2,main='Autocorrelation Function')
grid()
dev.off()
}
summary(x)
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Descriptive Statistics',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'# observations',header=TRUE)
a<-table.element(a,length(x))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'minimum',header=TRUE)
a<-table.element(a,min(x))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Q1',header=TRUE)
a<-table.element(a,quantile(x,0.25))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'median',header=TRUE)
a<-table.element(a,median(x))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'mean',header=TRUE)
a<-table.element(a,mean(x))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Q3',header=TRUE)
a<-table.element(a,quantile(x,0.75))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'maximum',header=TRUE)
a<-table.element(a,max(x))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')