Home » date » 2010 » Dec » 14 »

*The author of this computation has been verified* R Software Module: /rwasp_spectrum.wasp (opens new window with default values) Title produced by software: Spectral Analysis Date of computation: Tue, 14 Dec 2010 15:25:02 +0000   Cite this page as follows: Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2010/Dec/14/t1292340293g4wua3y0eidsp5j.htm/, Retrieved Thu, 23 May 2013 06:42:59 +0000   Original text written by user:   IsPrivate? No (this computation is public)   User-defined keywords:   System-generated keywords (parent): t122876443062pvsyb6eh8ar26 (pk = 30802) Estimated Impact 44   Dataseries X: » Textfile « » CSV « » Stem and Leaf « » Histogram « » Kernel Density « » Harrell-Davis Quantiles « » Central Tendency « » Variability « 655362 873127 1107897 1555964 1671159 1493308 2957796 2638691 1305669 1280496 921900 867888 652586 913831 1108544 1555827 1699283 1509458 3268975 2425016 1312703 1365498 934453 775019 651142 843192 1146766 1652601 1465906 1652734 2922334 2702805 1458956 1410363 1019279 936574 708917 885295 1099663 1576220 1487870 1488635 2882530 2677026 1404398 1344370 936865 872705 628151 953712 1160384 1400618 1661511 1495347 2918786 2775677 1407026 1370199 964526 850851 683118 847224 1073256 1514326 1503734 1507712 2865698 2788128 1391596 1366378 946295 859626   Output produced by software: Summary of computational transaction Raw Input view raw input (R code) Raw Output view raw output of R engine Computing time 3 seconds R Server 'Gwilym Jenkins' @ 72.249.127.135 Raw Periodogram Parameter Value Box-Cox transformation parameter (lambda) 1 Degree of non-seasonal differencing (d) 0 Degree of seasonal differencing (D) 0 Seasonal Period (s) 12 Frequency (Period) Spectrum 0.0139 (72) 62478617825.3127 0.0278 (36) 89128095696.3262 0.0417 (24) 84731853984.6067 0.0556 (18) 175777510725.099 0.0694 (14.4) 167248478466.740 0.0833 (12) 11234666087755.3 0.0972 (10.2857) 41322945855.9152 0.1111 (9) 85785134330.5371 0.125 (8) 49184534770.8798 0.1389 (7.2) 25578088059.7966 0.1528 (6.5455) 15034543597.4116 0.1667 (6) 1413150054626.46 0.1806 (5.5385) 17977231321.1271 0.1944 (5.1429) 1968168767.08022 0.2083 (4.8) 1498439370.98786 0.2222 (4.5) 3036206160.63415 0.2361 (4.2353) 13601757309.5936 0.25 (4) 1953482768642.24 0.2639 (3.7895) 2382725250.52715 0.2778 (3.6) 5929811889.50213 0.2917 (3.4286) 5632615327.15129 0.3056 (3.2727) 959114923.157202 0.3194 (3.1304) 8310138999.01496 0.3333 (3) 1181621223631.19 0.3472 (2.88) 17727381454.3681 0.3611 (2.7692) 16459815226.2094 0.375 (2.6667) 9237112245.75524 0.3889 (2.5714) 16784887855.0472 0.4028 (2.4828) 2896581629.83148 0.4167 (2.4) 457764000973.218 0.4306 (2.3226) 7058804467.61872 0.4444 (2.25) 27571483230.1647 0.4583 (2.1818) 10444018421.0633 0.4722 (2.1176) 4143590730.50703 0.4861 (2.0571) 11938751493.2789 0.5 (2) 3798311494.85855   Charts produced by software: http://www.freestatistics.org/blog/date/2010/Dec/14/t1292340293g4wua3y0eidsp5j/1x61s1292340298.png (opens in new window) http://www.freestatistics.org/blog/date/2010/Dec/14/t1292340293g4wua3y0eidsp5j/1x61s1292340298.ps (opens in new window) Click here to open pdf file. http://www.freestatistics.org/blog/date/2010/Dec/14/t1292340293g4wua3y0eidsp5j/2x61s1292340298.png (opens in new window) http://www.freestatistics.org/blog/date/2010/Dec/14/t1292340293g4wua3y0eidsp5j/2x61s1292340298.ps (opens in new window) Click here to open pdf file.   Parameters (Session): par1 = 1 ; par2 = 0 ; par3 = 0 ; par4 = 12 ;   Parameters (R input): par1 = 1 ; par2 = 0 ; par3 = 0 ; par4 = 12 ;   R code (references can be found in the software module): par1 <- as.numeric(par1) par2 <- as.numeric(par2) par3 <- as.numeric(par3) par4 <- as.numeric(par4) if (par1 == 0) { x <- log(x) } else { x <- (x ^ par1 - 1) / par1 } if (par2 > 0) x <- diff(x,lag=1,difference=par2) if (par3 > 0) x <- diff(x,lag=par4,difference=par3) bitmap(file='test1.png') r <- spectrum(x,main='Raw Periodogram') dev.off() bitmap(file='test2.png') cpgram(x,main='Cumulative Periodogram') dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Raw Periodogram',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Parameter',header=TRUE) a<-table.element(a,'Value',header=TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Box-Cox transformation parameter (lambda)',header=TRUE) a<-table.element(a,par1) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Degree of non-seasonal differencing (d)',header=TRUE) a<-table.element(a,par2) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Degree of seasonal differencing (D)',header=TRUE) a<-table.element(a,par3) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Seasonal Period (s)',header=TRUE) a<-table.element(a,par4) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Frequency (Period)',header=TRUE) a<-table.element(a,'Spectrum',header=TRUE) a<-table.row.end(a) for (i in 1:length(r$freq)) { a<-table.row.start(a) mylab <- round(r$freq[i],4) mylab <- paste(mylab,' (',sep='') mylab <- paste(mylab,round(1/r$freq[i],4),sep='') mylab <- paste(mylab,')',sep='') a<-table.element(a,mylab,header=TRUE) a<-table.element(a,round(r$spec[i],6)) a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable.tab')