Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_smp.wasp
Title produced by softwareStandard Deviation-Mean Plot
Date of computationWed, 04 Dec 2013 13:22:48 -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/2013/Dec/04/t1386181551qv6kkhnlan7tgnq.htm/, Retrieved Thu, 28 Mar 2024 20:23:36 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=230744, Retrieved Thu, 28 Mar 2024 20:23:36 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact73
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Univariate Data Series] [data set] [2008-12-01 19:54:57] [b98453cac15ba1066b407e146608df68]
- RMP   [Standard Deviation-Mean Plot] [Unemployment] [2010-11-29 10:34:47] [b98453cac15ba1066b407e146608df68]
- R PD      [Standard Deviation-Mean Plot] [WS10] [2013-12-04 18:22:48] [f1e366d257cd544a6a94e0c7cf247a26] [Current]
Feedback Forum

Post a new message
Dataseries X:
1356,876334
1267,205322
1296,765977
1088,768988
974,9813884
901,3596576
697,4384802
619,518781
446,3076104
279,7073872
388,5496968
290,3616669
278,7900756
433,3678289
436,029714
642,0796877
800,8666721
1005,060868
986,1196671
1320,620732
1293,652614
1315,564345
1438,509653
1311,68492
1281,851742
1229,184547
1157,483378
1049,930474
720,8534348
756,8141819
587,8054359
440,0632504
476,8424862
376,9041092
408,9909215
335,6074609
510,7385144
784,0262709
868,9048671
930,5891366
1232,387738
1433,973937
1441,341759
1560,894767
1570,121112
1539,736335
1509,516021
1264,261299
1280,1511
1104,198071
867,0732077
832,9038416
551,1656185
578,9382082
548,7690841
526,7505509
543,7872669
387,4285031
583,0257696
858,5103593
835,2939603
1234,576103
1307,198241
1550,368145
1513,93417
1607,433912
1768,513076
1778,655313
1668,394241
1383,577733
1271,618447
1140,579064
998,7028247
804,5968421
743,1134557
661,133487
635,0605076
334,5931783
502,1494165
450,8321675
700,3681247
799,3587889
1068,293463
1229,161035
1369,983001
1586,670588
1649,664042
1871,680788
1816,035687
1919,694672
1817,473725
1707,303794
1606,456459
1446,737592
1005,045182
913,67958
694,7422417
647,2193014
603,2233095
561,3616255
477,9286605
648,8516322
676,4550378
1016,256499
1139,927185
1269,321991
1546,025319
1753,304467
1958,22963
2094,065774
1977,440699
2027,739055
1842,633768
1852,639287
1706,985147
1471,71997
1159,935346
888,3265125
747,4889418
573,6363
696,6030203
736,0360216
625,9859074
668,0129992
878,1341576
925,7355157
1384,068603
1504,264376
1677,525248
1935,129054
1995,385689
2158,635275
2278,786793
2004,468115
2055,40499
1967,431468
1770,334016
1699,341376
1476,810831
970,7915825
997,474492
652,0843291
708,677725
656,5517255
770,8461221
901,3881145
832,0016083
1028,128263
1296,071512
1785,233557
1804,294853
2165,319782
2274,425758
2444,192672
2254,896853
2475,21364
2096,602738
2180,683121
2069,506046
1776,663609
1612,449679
1322,927331
914,8881288
749,3597793
985,6467572
880,7232311
819,4581121
760,5542696
985,0475722
1295,758114
1683,545809
1787,214405
1972,790239
2134,412385
2263,189757
2432,266646
2369,774565
2360,715385
2432,992026
2437,706015
1921,308588
1960,743986
1481,140265
1436,91711
1299,84893
985,736701
805,8439046
736,9963963
898,3310247
1065,843227
1113,559476
1211,918524
1462,155718
1801,38398
2045,731198
2405,387519
2376,281401
2542,194732
2535,172687
2518,76254
2467,712934
2483,547793
2111,5245
1987,469819
1900,284914
1545,174509
1241,063215
1244,17204
1001,510714
1004,787213
1064,602774
1250,752996
1244,505238
1565,513983
1904,45973
2200,393469
2489,302736
2382,022655
2740,656964
2692,201467
2707,985974
2793,734301
2909,368555
2611,616903
2367,927913
2036,392713
1794,090508
1658,074271
1450,206111
1159,117531
1144,297116
1188,107289
1082,864685
1152,30001
1370,470201
1560,414285
1878,82661
2213,171169
2306,794427
2878,274067
3043,563683
3163,511145
3218,690493
3177,601319
2941,518453
2922,823356
2325,88093
2414,009426
1874,014382
1866,949762
1330,453111
1413,735476
988,2221384
1046,533458
1398,652585
1444,958119
1637,387943
1926,364707
2088,421239
2250,892557
2803,393628
2920,006864
3028,392467
3284,085805
3116,539112
3120,244711
2946,135624
2740,67997
2810,140837
2233,823039
2114,011418
2076,542437
1571,481405
1567,958655
1217,444043
1136,643661
1481,287919
1528,477337
1797,720002
1819,469046
2395,358987
2707,314708
2921,060577
3246,433366
3051,477548
3444,481052
3432,580153
3421,04804
3381,830601
3169,679136
2872,11508
2601,76563
2143,011809
2159,598897
1645,12941
1622,54215
1350,064447
1372,486802
1475,078831
1356,08523
1974,230192
1835,493872
2495,225532
2813,660064
3199,803837
3128,428211
3614,588977
3741,816099
3711,62971
3501,116488
3620,888016
3348,935195




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=230744&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 time3 seconds
R Server'Herman Ole Andreas Wold' @ wold.wessa.net







Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
1800.653440791667400.8355816842221077.1689468
2938.528898116667411.2024252884751159.7195774
3735.19428515355.597109572321946.2442811
41220.54097975359.5302869929071059.3825976
5721.891798408333268.05408418852892.7225969
61421.678533775280.116692791589943.3613527
7743.946940916667260.21840037403894.5678567
81559.2020925325.5222185109021006.015092
9919.551439133333420.8587251257371275.3758065
101525.07003585833544.5106716537911520.429474
111265.45965556667576.4385644220521532.6493676
121436.51312025833606.1879680223041626.7024639
131669.334394575666.4105646427431704.3675179
141347.45523350833567.9870265170091431.3233417
152012.95107651667485.2685380611361452.6584428
161243.18234438333399.3217703557081223.7475897
172228.11040175347.327287062611080.039014
181430.60173291667393.827597809141198.882755
192431.94791333333404.7915371720741251.294284
201615.40362508333580.02118381321795.409382
212474.39596133333704.2308758316461888.237382
222068.10929253333795.0307601481982295.8636666
232220.97040933333727.2224130156521983.60105
242603.89239458333769.3573288639151963.193133
252095.78233525745.3437193263642031.766154
263082.15134941667665.1063016461511906.322227

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 800.653440791667 & 400.835581684222 & 1077.1689468 \tabularnewline
2 & 938.528898116667 & 411.202425288475 & 1159.7195774 \tabularnewline
3 & 735.19428515 & 355.597109572321 & 946.2442811 \tabularnewline
4 & 1220.54097975 & 359.530286992907 & 1059.3825976 \tabularnewline
5 & 721.891798408333 & 268.05408418852 & 892.7225969 \tabularnewline
6 & 1421.678533775 & 280.116692791589 & 943.3613527 \tabularnewline
7 & 743.946940916667 & 260.21840037403 & 894.5678567 \tabularnewline
8 & 1559.2020925 & 325.522218510902 & 1006.015092 \tabularnewline
9 & 919.551439133333 & 420.858725125737 & 1275.3758065 \tabularnewline
10 & 1525.07003585833 & 544.510671653791 & 1520.429474 \tabularnewline
11 & 1265.45965556667 & 576.438564422052 & 1532.6493676 \tabularnewline
12 & 1436.51312025833 & 606.187968022304 & 1626.7024639 \tabularnewline
13 & 1669.334394575 & 666.410564642743 & 1704.3675179 \tabularnewline
14 & 1347.45523350833 & 567.987026517009 & 1431.3233417 \tabularnewline
15 & 2012.95107651667 & 485.268538061136 & 1452.6584428 \tabularnewline
16 & 1243.18234438333 & 399.321770355708 & 1223.7475897 \tabularnewline
17 & 2228.11040175 & 347.32728706261 & 1080.039014 \tabularnewline
18 & 1430.60173291667 & 393.82759780914 & 1198.882755 \tabularnewline
19 & 2431.94791333333 & 404.791537172074 & 1251.294284 \tabularnewline
20 & 1615.40362508333 & 580.0211838132 & 1795.409382 \tabularnewline
21 & 2474.39596133333 & 704.230875831646 & 1888.237382 \tabularnewline
22 & 2068.10929253333 & 795.030760148198 & 2295.8636666 \tabularnewline
23 & 2220.97040933333 & 727.222413015652 & 1983.60105 \tabularnewline
24 & 2603.89239458333 & 769.357328863915 & 1963.193133 \tabularnewline
25 & 2095.78233525 & 745.343719326364 & 2031.766154 \tabularnewline
26 & 3082.15134941667 & 665.106301646151 & 1906.322227 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=230744&T=1

[TABLE]
[ROW][C]Standard Deviation-Mean Plot[/C][/ROW]
[ROW][C]Section[/C][C]Mean[/C][C]Standard Deviation[/C][C]Range[/C][/ROW]
[ROW][C]1[/C][C]800.653440791667[/C][C]400.835581684222[/C][C]1077.1689468[/C][/ROW]
[ROW][C]2[/C][C]938.528898116667[/C][C]411.202425288475[/C][C]1159.7195774[/C][/ROW]
[ROW][C]3[/C][C]735.19428515[/C][C]355.597109572321[/C][C]946.2442811[/C][/ROW]
[ROW][C]4[/C][C]1220.54097975[/C][C]359.530286992907[/C][C]1059.3825976[/C][/ROW]
[ROW][C]5[/C][C]721.891798408333[/C][C]268.05408418852[/C][C]892.7225969[/C][/ROW]
[ROW][C]6[/C][C]1421.678533775[/C][C]280.116692791589[/C][C]943.3613527[/C][/ROW]
[ROW][C]7[/C][C]743.946940916667[/C][C]260.21840037403[/C][C]894.5678567[/C][/ROW]
[ROW][C]8[/C][C]1559.2020925[/C][C]325.522218510902[/C][C]1006.015092[/C][/ROW]
[ROW][C]9[/C][C]919.551439133333[/C][C]420.858725125737[/C][C]1275.3758065[/C][/ROW]
[ROW][C]10[/C][C]1525.07003585833[/C][C]544.510671653791[/C][C]1520.429474[/C][/ROW]
[ROW][C]11[/C][C]1265.45965556667[/C][C]576.438564422052[/C][C]1532.6493676[/C][/ROW]
[ROW][C]12[/C][C]1436.51312025833[/C][C]606.187968022304[/C][C]1626.7024639[/C][/ROW]
[ROW][C]13[/C][C]1669.334394575[/C][C]666.410564642743[/C][C]1704.3675179[/C][/ROW]
[ROW][C]14[/C][C]1347.45523350833[/C][C]567.987026517009[/C][C]1431.3233417[/C][/ROW]
[ROW][C]15[/C][C]2012.95107651667[/C][C]485.268538061136[/C][C]1452.6584428[/C][/ROW]
[ROW][C]16[/C][C]1243.18234438333[/C][C]399.321770355708[/C][C]1223.7475897[/C][/ROW]
[ROW][C]17[/C][C]2228.11040175[/C][C]347.32728706261[/C][C]1080.039014[/C][/ROW]
[ROW][C]18[/C][C]1430.60173291667[/C][C]393.82759780914[/C][C]1198.882755[/C][/ROW]
[ROW][C]19[/C][C]2431.94791333333[/C][C]404.791537172074[/C][C]1251.294284[/C][/ROW]
[ROW][C]20[/C][C]1615.40362508333[/C][C]580.0211838132[/C][C]1795.409382[/C][/ROW]
[ROW][C]21[/C][C]2474.39596133333[/C][C]704.230875831646[/C][C]1888.237382[/C][/ROW]
[ROW][C]22[/C][C]2068.10929253333[/C][C]795.030760148198[/C][C]2295.8636666[/C][/ROW]
[ROW][C]23[/C][C]2220.97040933333[/C][C]727.222413015652[/C][C]1983.60105[/C][/ROW]
[ROW][C]24[/C][C]2603.89239458333[/C][C]769.357328863915[/C][C]1963.193133[/C][/ROW]
[ROW][C]25[/C][C]2095.78233525[/C][C]745.343719326364[/C][C]2031.766154[/C][/ROW]
[ROW][C]26[/C][C]3082.15134941667[/C][C]665.106301646151[/C][C]1906.322227[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=230744&T=1

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

As an alternative you can also use a QR Code:  

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

Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
1800.653440791667400.8355816842221077.1689468
2938.528898116667411.2024252884751159.7195774
3735.19428515355.597109572321946.2442811
41220.54097975359.5302869929071059.3825976
5721.891798408333268.05408418852892.7225969
61421.678533775280.116692791589943.3613527
7743.946940916667260.21840037403894.5678567
81559.2020925325.5222185109021006.015092
9919.551439133333420.8587251257371275.3758065
101525.07003585833544.5106716537911520.429474
111265.45965556667576.4385644220521532.6493676
121436.51312025833606.1879680223041626.7024639
131669.334394575666.4105646427431704.3675179
141347.45523350833567.9870265170091431.3233417
152012.95107651667485.2685380611361452.6584428
161243.18234438333399.3217703557081223.7475897
172228.11040175347.327287062611080.039014
181430.60173291667393.827597809141198.882755
192431.94791333333404.7915371720741251.294284
201615.40362508333580.02118381321795.409382
212474.39596133333704.2308758316461888.237382
222068.10929253333795.0307601481982295.8636666
232220.97040933333727.2224130156521983.60105
242603.89239458333769.3573288639151963.193133
252095.78233525745.3437193263642031.766154
263082.15134941667665.1063016461511906.322227







Regression: S.E.(k) = alpha + beta * Mean(k)
alpha233.377943473469
beta0.167234434932509
S.D.0.0407746552362585
T-STAT4.10143099833736
p-value0.000407648324467371

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 233.377943473469 \tabularnewline
beta & 0.167234434932509 \tabularnewline
S.D. & 0.0407746552362585 \tabularnewline
T-STAT & 4.10143099833736 \tabularnewline
p-value & 0.000407648324467371 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=230744&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]233.377943473469[/C][/ROW]
[ROW][C]beta[/C][C]0.167234434932509[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0407746552362585[/C][/ROW]
[ROW][C]T-STAT[/C][C]4.10143099833736[/C][/ROW]
[ROW][C]p-value[/C][C]0.000407648324467371[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=230744&T=2

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

As an alternative you can also use a QR Code:  

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

Regression: S.E.(k) = alpha + beta * Mean(k)
alpha233.377943473469
beta0.167234434932509
S.D.0.0407746552362585
T-STAT4.10143099833736
p-value0.000407648324467371







Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha2.26054379979923
beta0.534551041816121
S.D.0.126494685200834
T-STAT4.22587748226277
p-value0.000297321539096121
Lambda0.465448958183879

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & 2.26054379979923 \tabularnewline
beta & 0.534551041816121 \tabularnewline
S.D. & 0.126494685200834 \tabularnewline
T-STAT & 4.22587748226277 \tabularnewline
p-value & 0.000297321539096121 \tabularnewline
Lambda & 0.465448958183879 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=230744&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]2.26054379979923[/C][/ROW]
[ROW][C]beta[/C][C]0.534551041816121[/C][/ROW]
[ROW][C]S.D.[/C][C]0.126494685200834[/C][/ROW]
[ROW][C]T-STAT[/C][C]4.22587748226277[/C][/ROW]
[ROW][C]p-value[/C][C]0.000297321539096121[/C][/ROW]
[ROW][C]Lambda[/C][C]0.465448958183879[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=230744&T=3

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

As an alternative you can also use a QR Code:  

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

Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha2.26054379979923
beta0.534551041816121
S.D.0.126494685200834
T-STAT4.22587748226277
p-value0.000297321539096121
Lambda0.465448958183879



Parameters (Session):
par1 = 12 ;
Parameters (R input):
par1 = 12 ;
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))
j <- 0
k <- 1
for (i in 1:(np*par1))
{
j = j + 1
arr[j,k] <- x[i]
if (j == par1) {
j = 0
k=k+1
}
}
arr
arr.mean <- array(NA,dim=np)
arr.sd <- array(NA,dim=np)
arr.range <- array(NA,dim=np)
for (j in 1:np)
{
arr.mean[j] <- mean(arr[,j],na.rm=TRUE)
arr.sd[j] <- sd(arr[,j],na.rm=TRUE)
arr.range[j] <- max(arr[,j],na.rm=TRUE) - min(arr[,j],na.rm=TRUE)
}
arr.mean
arr.sd
arr.range
(lm1 <- lm(arr.sd~arr.mean))
(lnlm1 <- lm(log(arr.sd)~log(arr.mean)))
(lm2 <- lm(arr.range~arr.mean))
bitmap(file='test1.png')
plot(arr.mean,arr.sd,main='Standard Deviation-Mean Plot',xlab='mean',ylab='standard deviation')
dev.off()
bitmap(file='test2.png')
plot(arr.mean,arr.range,main='Range-Mean Plot',xlab='mean',ylab='range')
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Standard Deviation-Mean Plot',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Section',header=TRUE)
a<-table.element(a,'Mean',header=TRUE)
a<-table.element(a,'Standard Deviation',header=TRUE)
a<-table.element(a,'Range',header=TRUE)
a<-table.row.end(a)
for (j in 1:np) {
a<-table.row.start(a)
a<-table.element(a,j,header=TRUE)
a<-table.element(a,arr.mean[j])
a<-table.element(a,arr.sd[j] )
a<-table.element(a,arr.range[j] )
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Regression: S.E.(k) = alpha + beta * Mean(k)',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'alpha',header=TRUE)
a<-table.element(a,lm1$coefficients[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'beta',header=TRUE)
a<-table.element(a,lm1$coefficients[[2]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'S.D.',header=TRUE)
a<-table.element(a,summary(lm1)$coefficients[2,2])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'T-STAT',header=TRUE)
a<-table.element(a,summary(lm1)$coefficients[2,3])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value',header=TRUE)
a<-table.element(a,summary(lm1)$coefficients[2,4])
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable1.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Regression: ln S.E.(k) = alpha + beta * ln Mean(k)',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'alpha',header=TRUE)
a<-table.element(a,lnlm1$coefficients[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'beta',header=TRUE)
a<-table.element(a,lnlm1$coefficients[[2]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'S.D.',header=TRUE)
a<-table.element(a,summary(lnlm1)$coefficients[2,2])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'T-STAT',header=TRUE)
a<-table.element(a,summary(lnlm1)$coefficients[2,3])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value',header=TRUE)
a<-table.element(a,summary(lnlm1)$coefficients[2,4])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Lambda',header=TRUE)
a<-table.element(a,1-lnlm1$coefficients[[2]])
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable2.tab')