R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> x <- c(1145.11
+ ,1176.86
+ ,1206.41
+ ,1192.72
+ ,1214.82
+ ,1199.07
+ ,1157.47
+ ,1100.1
+ ,1095.63
+ ,1105.63
+ ,1137.79
+ ,1124.72
+ ,1152.6
+ ,1211.85
+ ,1239.62
+ ,1244.13
+ ,1198.42
+ ,1227.99
+ ,1304.92
+ ,1340.26
+ ,1307.32
+ ,1356.51
+ ,1383.29
+ ,1437.87
+ ,1494.56
+ ,1521.42
+ ,1498.76
+ ,1488.75
+ ,1524.62
+ ,1439.27
+ ,1423.11
+ ,1466.85
+ ,1425.83
+ ,1363.45
+ ,1389.18
+ ,1395.89
+ ,1368.43
+ ,1349.03
+ ,1299.88
+ ,1365.41
+ ,1451.04
+ ,1433.75
+ ,1464.65
+ ,1475.57
+ ,1471.16
+ ,1429.12
+ ,1452.46
+ ,1538.09
+ ,1631.59
+ ,1665.5
+ ,1690.6
+ ,1711.74
+ ,1734.1
+ ,1748.09
+ ,1703.45
+ ,1745.74
+ ,1751.01
+ ,1795.65
+ ,1852.13
+ ,1877.1
+ ,1989.31
+ ,2097.76
+ ,2154.87
+ ,2152.18
+ ,2250.27
+ ,2346.9
+ ,2525.56
+ ,2409.36
+ ,2394.36
+ ,2401.33
+ ,2354.32
+ ,2450.41
+ ,2504.67
+ ,2661.39
+ ,2880.4
+ ,3064.42
+ ,3141.12
+ ,3327.7
+ ,3564.95
+ ,3403.13
+ ,3149.9
+ ,3006.84
+ ,3230.66
+ ,3361.13
+ ,3484.74
+ ,3411.13
+ ,3288.18
+ ,3280.37
+ ,3173.95
+ ,3165.26
+ ,3092.71
+ ,3053.05
+ ,3181.96
+ ,2999.93
+ ,3249.57
+ ,3210.52
+ ,3030.29
+ ,2803.47
+ ,2767.63
+ ,2882.6
+ ,2863.36
+ ,2897.06
+ ,3012.61
+ ,3142.95
+ ,3032.93
+ ,3045.78
+ ,3110.52
+ ,3013.24
+ ,2987.1
+ ,2995.55
+ ,2833.18
+ ,2848.96
+ ,2794.83
+ ,2845.26
+ ,2915.02
+ ,2892.63
+ ,2604.42
+ ,2641.65
+ ,2659.81
+ ,2638.53
+ ,2720.25
+ ,2745.88
+ ,2735.7
+ ,2811.7
+ ,2799.43
+ ,2555.28
+ ,2304.98
+ ,2214.95
+ ,2065.81
+ ,1940.49
+ ,2042
+ ,1995.37
+ ,1946.81
+ ,1765.9
+ ,1635.25
+ ,1833.42
+ ,1910.43
+ ,1959.67
+ ,1969.6
+ ,2061.41
+ ,2093.48
+ ,2120.88
+ ,2174.56
+ ,2196.72
+ ,2350.44
+ ,2440.25
+ ,2408.64
+ ,2472.81
+ ,2407.6
+ ,2454.62
+ ,2448.05
+ ,2497.84
+ ,2645.64
+ ,2756.76
+ ,2849.27
+ ,2921.44
+ ,2981.85
+ ,3080.58
+ ,3106.22
+ ,3119.31
+ ,3061.26
+ ,3097.31
+ ,3161.69
+ ,3257.16
+ ,3277.01
+ ,3295.32
+ ,3363.99
+ ,3494.17
+ ,3667.03
+ ,3813.06
+ ,3917.96
+ ,3895.51
+ ,3801.06
+ ,3570.12
+ ,3701.61
+ ,3862.27
+ ,3970.1
+ ,4138.52
+ ,4199.75
+ ,4290.89
+ ,4443.91
+ ,4502.64
+ ,4356.98
+ ,4591.27
+ ,4696.96
+ ,4621.4
+ ,4562.84
+ ,4202.52
+ ,4296.49
+ ,4435.23
+ ,4105.18
+ ,4116.68
+ ,3844.49
+ ,3720.98
+ ,3674.4
+ ,3857.62
+ ,3801.06
+ ,3504.37
+ ,3032.6
+ ,3047.03
+ ,2962.34
+ ,2197.82
+ ,2014.45
+ ,1862.83
+ ,1905.41
+ ,1810.99
+ ,1670.07
+ ,1864.44
+ ,2052.02
+ ,2029.6
+ ,2070.83
+ ,2293.41
+ ,2443.27
+ ,2513.17
+ ,2466.92
+ ,2502.66)
> par1 = '12'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 216
> (np <- floor(n / par1))
[1] 18
> 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
 [1] 18 18 18 18 18 18 18 18 18 18 18 18
> arr
         [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]    [,8]    [,9]
 [1,] 1145.11 1152.60 1494.56 1368.43 1631.59 1989.31 2504.67 3484.74 3030.29
 [2,] 1176.86 1211.85 1521.42 1349.03 1665.50 2097.76 2661.39 3411.13 2803.47
 [3,] 1206.41 1239.62 1498.76 1299.88 1690.60 2154.87 2880.40 3288.18 2767.63
 [4,] 1192.72 1244.13 1488.75 1365.41 1711.74 2152.18 3064.42 3280.37 2882.60
 [5,] 1214.82 1198.42 1524.62 1451.04 1734.10 2250.27 3141.12 3173.95 2863.36
 [6,] 1199.07 1227.99 1439.27 1433.75 1748.09 2346.90 3327.70 3165.26 2897.06
 [7,] 1157.47 1304.92 1423.11 1464.65 1703.45 2525.56 3564.95 3092.71 3012.61
 [8,] 1100.10 1340.26 1466.85 1475.57 1745.74 2409.36 3403.13 3053.05 3142.95
 [9,] 1095.63 1307.32 1425.83 1471.16 1751.01 2394.36 3149.90 3181.96 3032.93
[10,] 1105.63 1356.51 1363.45 1429.12 1795.65 2401.33 3006.84 2999.93 3045.78
[11,] 1137.79 1383.29 1389.18 1452.46 1852.13 2354.32 3230.66 3249.57 3110.52
[12,] 1124.72 1437.87 1395.89 1538.09 1877.10 2450.41 3361.13 3210.52 3013.24
        [,10]   [,11]   [,12]   [,13]   [,14]   [,15]   [,16]   [,17]   [,18]
 [1,] 2987.10 2720.25 1946.81 2350.44 2981.85 3667.03 4443.91 3844.49 1905.41
 [2,] 2995.55 2745.88 1765.90 2440.25 3080.58 3813.06 4502.64 3720.98 1810.99
 [3,] 2833.18 2735.70 1635.25 2408.64 3106.22 3917.96 4356.98 3674.40 1670.07
 [4,] 2848.96 2811.70 1833.42 2472.81 3119.31 3895.51 4591.27 3857.62 1864.44
 [5,] 2794.83 2799.43 1910.43 2407.60 3061.26 3801.06 4696.96 3801.06 2052.02
 [6,] 2845.26 2555.28 1959.67 2454.62 3097.31 3570.12 4621.40 3504.37 2029.60
 [7,] 2915.02 2304.98 1969.60 2448.05 3161.69 3701.61 4562.84 3032.60 2070.83
 [8,] 2892.63 2214.95 2061.41 2497.84 3257.16 3862.27 4202.52 3047.03 2293.41
 [9,] 2604.42 2065.81 2093.48 2645.64 3277.01 3970.10 4296.49 2962.34 2443.27
[10,] 2641.65 1940.49 2120.88 2756.76 3295.32 4138.52 4435.23 2197.82 2513.17
[11,] 2659.81 2042.00 2174.56 2849.27 3363.99 4199.75 4105.18 2014.45 2466.92
[12,] 2638.53 1995.37 2196.72 2921.44 3494.17 4290.89 4116.68 1862.83 2502.66
      [,19]
 [1,]    NA
 [2,]    NA
 [3,]    NA
 [4,]    NA
 [5,]    NA
 [6,]    NA
 [7,]    NA
 [8,]    NA
 [9,]    NA
[10,]    NA
[11,]    NA
[12,]    NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]    [,6]    [,7]    [,8]    [,9]
 [1,]  31.75  59.25  26.86 -19.40  33.91  108.45  156.72  -73.61 -226.82
 [2,]  29.55  27.77 -22.66 -49.15  25.10   57.11  219.01 -122.95  -35.84
 [3,] -13.69   4.51 -10.01  65.53  21.14   -2.69  184.02   -7.81  114.97
 [4,]  22.10 -45.71  35.87  85.63  22.36   98.09   76.70 -106.42  -19.24
 [5,] -15.75  29.57 -85.35 -17.29  13.99   96.63  186.58   -8.69   33.70
 [6,] -41.60  76.93 -16.16  30.90 -44.64  178.66  237.25  -72.55  115.55
 [7,] -57.37  35.34  43.74  10.92  42.29 -116.20 -161.82  -39.66  130.34
 [8,]  -4.47 -32.94 -41.02  -4.41   5.27  -15.00 -253.23  128.91 -110.02
 [9,]  10.00  49.19 -62.38 -42.04  44.64    6.97 -143.06 -182.03   12.85
[10,]  32.16  26.78  25.73  23.34  56.48  -47.01  223.82  249.64   64.74
[11,] -13.07  54.58   6.71  85.63  24.97   96.09  130.47  -39.05  -97.28
[12,]  27.88  56.69 -27.46  93.50 112.21   54.26  123.61 -180.23  -26.14
        [,10]   [,11]   [,12]  [,13]  [,14]   [,15]   [,16]   [,17]   [,18]
 [1,]    8.45   25.63 -180.91  89.81  98.73  146.03   58.73 -123.51  -94.42
 [2,] -162.37  -10.18 -130.65 -31.61  25.64  104.90 -145.66  -46.58 -140.92
 [3,]   15.78   76.00  198.17  64.17  13.09  -22.45  234.29  183.22  194.37
 [4,]  -54.13  -12.27   77.01 -65.21 -58.05  -94.45  105.69  -56.56  187.58
 [5,]   50.43 -244.15   49.24  47.02  36.05 -230.94  -75.56 -296.69  -22.42
 [6,]   69.76 -250.30    9.93  -6.57  64.38  131.49  -58.56 -471.77   41.23
 [7,]  -22.39  -90.03   91.81  49.79  95.47  160.66 -360.32   14.43  222.58
 [8,] -288.21 -149.14   32.07 147.80  19.85  107.83   93.97  -84.69  149.86
 [9,]   37.23 -125.32   27.40 111.12  18.31  168.42  138.74 -764.52   69.90
[10,]   18.16  101.51   53.68  92.51  68.67   61.23 -330.05 -183.37  -46.25
[11,]  -21.28  -46.63   22.16  72.17 130.18   91.14   11.50 -151.62   35.74
[12,]   81.72  -48.56  153.72  60.41 172.86  153.02 -272.19   42.58      NA
      [,19]
 [1,]    NA
 [2,]    NA
 [3,]    NA
 [4,]    NA
 [5,]    NA
 [6,]    NA
 [7,]    NA
 [8,]    NA
 [9,]    NA
[10,]    NA
[11,]    NA
[12,]    NA
> 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
> postscript(file="/var/www/html/rcomp/tmp/16doi1292691010.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(arr.mean,type='b',ylab='mean',main='Mean Plot',xlab='Periodic Index')
> mtext(paste('#blocks = ',np))
> abline(overall.mean,0)
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/26doi1292691010.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(arr.median,type='b',ylab='median',main='Median Plot',xlab='Periodic Index')
> mtext(paste('#blocks = ',np))
> abline(overall.median,0)
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/3z4631292691010.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(arr.midrange,type='b',ylab='midrange',main='Midrange Plot',xlab='Periodic Index')
> mtext(paste('#blocks = ',np))
> abline(overall.midrange,0)
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/4z4631292691010.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> 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'))
$stats
         [,1]    [,2]    [,3]     [,4]     [,5]    [,6]     [,7]    [,8]
[1,] 1145.110 1176.86 1206.41 1192.720 1198.420 1199.07 1157.470 1100.10
[2,] 1631.590 1665.50 1635.25 1711.740 1734.100 1748.09 1703.450 1745.74
[3,] 2427.555 2550.82 2572.17 2642.255 2601.215 2504.95 2486.805 2453.60
[4,] 3030.290 3080.58 3106.22 3119.310 3141.120 3165.26 3092.710 3142.95
[5,] 4443.910 4502.64 4356.98 4591.270 4696.960 4621.40 4562.840 4202.52
         [,9]   [,10]   [,11]    [,12]
[1,] 1095.630 1105.63 1137.79 1124.720
[2,] 1751.010 1795.65 1852.13 1862.830
[3,] 2523.845 2457.25 2410.62 2476.535
[4,] 3149.900 3006.84 3230.66 3210.520
[5,] 4296.490 4435.23 4199.75 4290.890

$n
 [1] 18 18 18 18 18 18 18 18 18 18 18 18

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 1906.666 2023.831 2024.367 2118.062 2077.227 1977.182 1969.431 1933.266
[2,] 2948.444 3077.809 3119.973 3166.448 3125.203 3032.718 3004.179 2973.934
         [,9]    [,10]    [,11]    [,12]
[1,] 2002.885 2006.191 1897.242 1974.642
[2,] 3044.805 2908.309 2923.998 2978.428

$out
numeric(0)

$group
numeric(0)

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12"

Warning message:
In bxp(list(stats = c(1145.11, 1631.59, 2427.555, 3030.29, 4443.91,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5z4631292691010.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> 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'))
$stats
         [,1]     [,2]    [,3]     [,4]    [,5]    [,6]     [,7]    [,8]
[1,] -226.820 -162.370 -22.450 -106.420 -244.15 -72.550 -161.820 -288.21
[2,]  -73.610 -122.950  -2.690  -56.560  -75.56 -44.640  -57.370  -84.69
[3,]   29.305  -27.135  42.655    4.915    2.65  20.415   24.885   -4.44
[4,]   89.810   27.770 183.220   77.010   47.02  76.930   91.810   93.97
[5,]  156.720  219.010 234.290  187.580  186.58 237.250  222.580  149.86
        [,9]  [,10]    [,11]   [,12]
[1,] -182.03 -47.01 -151.620 -180.23
[2,]  -62.38  18.16  -21.280  -26.14
[3,]   15.58  42.92   23.565   56.69
[4,]   49.19  68.67   85.630  112.21
[5,]  168.42 101.51  130.470  172.86

$n
 [1] 18 18 18 18 18 18 18 18 18 18 18 17

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,] -31.55417 -83.26457 -26.57966 -44.82775 -42.99997 -24.85883 -30.67106
[2,]  90.16417  28.99457 111.88966  54.65775  48.29997  65.68883  80.44106
         [,8]      [,9]    [,10]     [,11]     [,12]
[1,] -70.9747 -25.96974 24.10959 -16.24931   3.67341
[2,]  62.0947  57.12974 61.73041  63.37931 109.70659

$out
 [1] -296.69 -250.30 -471.77 -360.32 -764.52  223.82  249.64 -330.05 -183.37
[10] -272.19

$group
 [1]  5  6  6  7  9 10 10 10 10 12

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12"

Warning message:
In bxp(list(stats = c(-226.82, -73.6099999999997, 29.3050000000001,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6sdn61292691010.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> 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'))
$stats
         [,1]     [,2]    [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 1095.630 1152.600 1363.45 1299.880 1631.590 1989.310 2504.670 2999.930
[2,] 1115.175 1219.920 1409.50 1366.920 1697.025 2153.525 2943.620 3128.985
[3,] 1151.290 1274.525 1453.06 1442.395 1739.920 2350.610 3145.510 3196.240
[4,] 1195.895 1348.385 1496.66 1467.905 1773.330 2405.345 3344.415 3284.275
[5,] 1214.820 1437.870 1524.62 1538.090 1877.100 2525.560 3564.950 3484.740
         [,9]    [,10]    [,11]    [,12]    [,13]    [,14]    [,15]    [,16]
[1,] 2767.630 2604.420 1940.490 1635.250 2350.440 2981.850 3570.120 4105.180
[2,] 2872.980 2650.730 2053.905 1871.925 2424.445 3088.945 3751.335 4249.505
[3,] 3012.925 2839.220 2430.130 1964.635 2463.715 3140.500 3878.890 4439.570
[4,] 3039.355 2903.825 2740.790 2107.180 2701.200 3286.165 4054.310 4577.055
[5,] 3142.950 2995.550 2811.700 2196.720 2921.440 3494.170 4290.890 4696.960
       [,17]    [,18] [,19]
[1,] 1862.83 1670.070    NA
[2,] 2580.08 1884.925    NA
[3,] 3275.70 2061.425    NA
[4,] 3761.02 2455.095    NA
[5,] 3857.62 2513.170    NA

$n
 [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 1114.473 1215.931 1413.306 1396.335 1705.117 2235.753 2962.705 3125.411
[2,] 1188.107 1333.119 1492.814 1488.455 1774.723 2465.467 3328.315 3267.069
        [,9]    [,10]    [,11]    [,12]    [,13]    [,14]    [,15]    [,16]
[1,] 2937.04 2723.782 2116.837 1857.334 2337.485 3050.547 3740.701 4290.172
[2,] 3088.81 2954.658 2743.423 2071.936 2589.945 3230.453 4017.079 4588.968
        [,17]    [,18] [,19]
[1,] 2737.065 1801.367    NA
[2,] 3814.335 2321.483    NA

$out
numeric(0)

$group
numeric(0)

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12" "13" "14" "15"
[16] "16" "17" "18" NA  

Warning message:
In bxp(list(stats = c(1095.63, 1115.175, 1151.29, 1195.895, 1214.82,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7sdn61292691010.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> 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'))
$stats
         [,1]     [,2]     [,3]
[1,] 2464.708 2410.620 2346.860
[2,] 2483.967 2455.425 2396.724
[3,] 2516.259 2495.878 2436.771
[4,] 2524.847 2561.495 2474.906
[5,] 2548.686 2642.255 2546.668

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 2497.613 2447.498 2401.112
[2,] 2534.904 2544.257 2472.431

$out
numeric(0)

$group
numeric(0)

$names
[1] "mean"     "median"   "midrange"

Warning message:
In bxp(list(stats = c(2464.70833333333, 2483.9675, 2516.25861111111,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/16doi1292691010.ps tmp/16doi1292691010.png",intern=TRUE))
character(0)
> try(system("convert tmp/26doi1292691010.ps tmp/26doi1292691010.png",intern=TRUE))
character(0)
> try(system("convert tmp/3z4631292691010.ps tmp/3z4631292691010.png",intern=TRUE))
character(0)
> try(system("convert tmp/4z4631292691010.ps tmp/4z4631292691010.png",intern=TRUE))
character(0)
> try(system("convert tmp/5z4631292691010.ps tmp/5z4631292691010.png",intern=TRUE))
character(0)
> try(system("convert tmp/6sdn61292691010.ps tmp/6sdn61292691010.png",intern=TRUE))
character(0)
> try(system("convert tmp/7sdn61292691010.ps tmp/7sdn61292691010.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.253   0.998   4.012