Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_meanplot.wasp
Title produced by softwareMean Plot
Date of computationThu, 25 Nov 2010 22:06:58 +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/Nov/25/t1290722738pplf666zbwdni8c.htm/, Retrieved Thu, 28 Mar 2024 11:48:26 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=101592, Retrieved Thu, 28 Mar 2024 11:48:26 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact116
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Central Tendency] [Arabica Price in ...] [2008-01-19 12:03:37] [74be16979710d4c4e7c6647856088456]
- RMPD  [Kernel Density Estimation] [] [2010-11-24 12:06:47] [c895532cb7349383dee5125244983cc8]
- RMP     [Mean Plot] [] [2010-11-24 12:10:12] [c895532cb7349383dee5125244983cc8]
-    D      [Mean Plot] [] [2010-11-25 17:36:24] [c895532cb7349383dee5125244983cc8]
-               [Mean Plot] [mini tutorial 5] [2010-11-25 22:06:58] [87bb5e10c18d96bd329dff2d857096c8] [Current]
Feedback Forum

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




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

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



Parameters (Session):
par1 = 12 ;
Parameters (R input):
par1 = 12 ; par2 = ; par3 = ; par4 = ; par5 = ; par6 = ; par7 = ; par8 = ; par9 = ; par10 = ; par11 = ; par12 = ; par13 = ; par14 = ; par15 = ; par16 = ; par17 = ; par18 = ; par19 = ; par20 = ;
R code (references can be found in the software module):
par1 <- as.numeric(par1)
(n <- length(x))
(np <- floor(n / par1))
arr <- array(NA,dim=c(par1,np+1))
darr <- array(NA,dim=c(par1,np+1))
ari <- array(0,dim=par1)
dx <- diff(x)
j <- 0
for (i in 1:n)
{
j = j + 1
ari[j] = ari[j] + 1
arr[j,ari[j]] <- x[i]
darr[j,ari[j]] <- dx[i]
if (j == par1) j = 0
}
ari
arr
darr
arr.mean <- array(NA,dim=par1)
arr.median <- array(NA,dim=par1)
arr.midrange <- array(NA,dim=par1)
for (j in 1:par1)
{
arr.mean[j] <- mean(arr[j,],na.rm=TRUE)
arr.median[j] <- median(arr[j,],na.rm=TRUE)
arr.midrange[j] <- (quantile(arr[j,],0.75,na.rm=TRUE) + quantile(arr[j,],0.25,na.rm=TRUE)) / 2
}
overall.mean <- mean(x)
overall.median <- median(x)
overall.midrange <- (quantile(x,0.75) + quantile(x,0.25)) / 2
bitmap(file='plot1.png')
plot(arr.mean,type='b',ylab='mean',main='Mean Plot',xlab='Periodic Index')
mtext(paste('#blocks = ',np))
abline(overall.mean,0)
dev.off()
bitmap(file='plot2.png')
plot(arr.median,type='b',ylab='median',main='Median Plot',xlab='Periodic Index')
mtext(paste('#blocks = ',np))
abline(overall.median,0)
dev.off()
bitmap(file='plot3.png')
plot(arr.midrange,type='b',ylab='midrange',main='Midrange Plot',xlab='Periodic Index')
mtext(paste('#blocks = ',np))
abline(overall.midrange,0)
dev.off()
bitmap(file='plot4.png')
z <- data.frame(t(arr))
names(z) <- c(1:par1)
(boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Periodic Subseries'))
dev.off()
bitmap(file='plot4b.png')
z <- data.frame(t(darr))
names(z) <- c(1:par1)
(boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Differenced Periodic Subseries'))
dev.off()
bitmap(file='plot5.png')
z <- data.frame(arr)
names(z) <- c(1:np)
(boxplot(z,notch=TRUE,col='grey',xlab='Block Index',ylab='Value',main='Notched Box Plots - Sequential Blocks'))
dev.off()
bitmap(file='plot6.png')
z <- data.frame(cbind(arr.mean,arr.median,arr.midrange))
names(z) <- list('mean','median','midrange')
(boxplot(z,notch=TRUE,col='grey',ylab='Overall Central Tendency',main='Notched Box Plots'))
dev.off()