R version 2.15.2 (2012-10-26) -- "Trick or Treat"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-pc-linux-gnu (32-bit)
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 <- array(list(2001
+ ,100
+ ,95
+ ,102
+ ,103
+ ,91
+ ,99
+ ,101
+ ,91
+ ,114
+ ,101
+ ,103
+ ,85
+ ,2001
+ ,94
+ ,97
+ ,99
+ ,117
+ ,85
+ ,97
+ ,97
+ ,87
+ ,99
+ ,99
+ ,97
+ ,94
+ ,2001
+ ,105
+ ,97
+ ,108
+ ,115
+ ,110
+ ,113
+ ,108
+ ,103
+ ,98
+ ,104
+ ,110
+ ,107
+ ,2001
+ ,95
+ ,97
+ ,92
+ ,74
+ ,90
+ ,100
+ ,95
+ ,97
+ ,91
+ ,99
+ ,97
+ ,98
+ ,2001
+ ,103
+ ,103
+ ,99
+ ,74
+ ,103
+ ,105
+ ,99
+ ,96
+ ,111
+ ,101
+ ,103
+ ,111
+ ,2001
+ ,103
+ ,101
+ ,102
+ ,81
+ ,119
+ ,109
+ ,101
+ ,105
+ ,104
+ ,102
+ ,106
+ ,115
+ ,2001
+ ,100
+ ,96
+ ,87
+ ,86
+ ,76
+ ,91
+ ,92
+ ,74
+ ,100
+ ,93
+ ,89
+ ,76
+ ,2001
+ ,108
+ ,94
+ ,71
+ ,114
+ ,93
+ ,89
+ ,92
+ ,87
+ ,108
+ ,97
+ ,85
+ ,100
+ ,2001
+ ,108
+ ,97
+ ,105
+ ,102
+ ,105
+ ,105
+ ,100
+ ,105
+ ,113
+ ,91
+ ,100
+ ,103
+ ,2001
+ ,120
+ ,101
+ ,115
+ ,85
+ ,92
+ ,120
+ ,106
+ ,118
+ ,113
+ ,97
+ ,106
+ ,117
+ ,2001
+ ,112
+ ,77
+ ,103
+ ,63
+ ,75
+ ,107
+ ,99
+ ,102
+ ,114
+ ,94
+ ,95
+ ,101
+ ,2001
+ ,102
+ ,93
+ ,75
+ ,61
+ ,61
+ ,84
+ ,84
+ ,101
+ ,109
+ ,90
+ ,74
+ ,73
+ ,2002
+ ,105
+ ,45
+ ,97
+ ,87
+ ,80
+ ,101
+ ,106
+ ,86
+ ,116
+ ,105
+ ,94
+ ,84
+ ,2002
+ ,101
+ ,48
+ ,95
+ ,97
+ ,85
+ ,105
+ ,101
+ ,83
+ ,102
+ ,103
+ ,90
+ ,90
+ ,2002
+ ,108
+ ,52
+ ,99
+ ,88
+ ,94
+ ,119
+ ,113
+ ,92
+ ,107
+ ,112
+ ,99
+ ,105
+ ,2002
+ ,107
+ ,49
+ ,100
+ ,67
+ ,78
+ ,114
+ ,110
+ ,87
+ ,111
+ ,114
+ ,100
+ ,111
+ ,2002
+ ,109
+ ,53
+ ,92
+ ,59
+ ,92
+ ,114
+ ,103
+ ,94
+ ,122
+ ,111
+ ,96
+ ,110
+ ,2002
+ ,110
+ ,60
+ ,94
+ ,63
+ ,90
+ ,119
+ ,107
+ ,94
+ ,123
+ ,106
+ ,102
+ ,116
+ ,2002
+ ,111
+ ,51
+ ,89
+ ,86
+ ,72
+ ,99
+ ,98
+ ,75
+ ,108
+ ,112
+ ,88
+ ,85
+ ,2002
+ ,110
+ ,42
+ ,67
+ ,99
+ ,77
+ ,91
+ ,90
+ ,85
+ ,115
+ ,102
+ ,78
+ ,92
+ ,2002
+ ,117
+ ,56
+ ,109
+ ,85
+ ,76
+ ,121
+ ,105
+ ,104
+ ,120
+ ,103
+ ,99
+ ,117
+ ,2002
+ ,130
+ ,51
+ ,113
+ ,74
+ ,89
+ ,128
+ ,116
+ ,109
+ ,117
+ ,105
+ ,107
+ ,119
+ ,2002
+ ,114
+ ,53
+ ,106
+ ,55
+ ,55
+ ,112
+ ,102
+ ,121
+ ,115
+ ,101
+ ,93
+ ,100
+ ,2002
+ ,113
+ ,55
+ ,78
+ ,54
+ ,47
+ ,93
+ ,88
+ ,124
+ ,116
+ ,101
+ ,74
+ ,71
+ ,2003
+ ,110
+ ,44
+ ,102
+ ,81
+ ,91
+ ,108
+ ,114
+ ,88
+ ,118
+ ,117
+ ,96
+ ,82
+ ,2003
+ ,107
+ ,51
+ ,97
+ ,88
+ ,85
+ ,107
+ ,104
+ ,86
+ ,98
+ ,109
+ ,99
+ ,90
+ ,2003
+ ,110
+ ,52
+ ,96
+ ,75
+ ,89
+ ,115
+ ,111
+ ,98
+ ,121
+ ,120
+ ,103
+ ,109
+ ,2003
+ ,113
+ ,54
+ ,99
+ ,55
+ ,90
+ ,121
+ ,111
+ ,94
+ ,118
+ ,115
+ ,102
+ ,112
+ ,2003
+ ,106
+ ,50
+ ,86
+ ,47
+ ,72
+ ,112
+ ,102
+ ,102
+ ,120
+ ,107
+ ,96
+ ,103
+ ,2003
+ ,118
+ ,57
+ ,92
+ ,54
+ ,83
+ ,123
+ ,106
+ ,96
+ ,111
+ ,110
+ ,106
+ ,116
+ ,2003
+ ,118
+ ,49
+ ,86
+ ,71
+ ,72
+ ,101
+ ,104
+ ,79
+ ,117
+ ,110
+ ,95
+ ,89
+ ,2003
+ ,114
+ ,41
+ ,62
+ ,79
+ ,75
+ ,87
+ ,94
+ ,95
+ ,110
+ ,105
+ ,82
+ ,91
+ ,2003
+ ,121
+ ,58
+ ,105
+ ,77
+ ,85
+ ,124
+ ,116
+ ,106
+ ,107
+ ,116
+ ,109
+ ,121
+ ,2003
+ ,130
+ ,63
+ ,108
+ ,57
+ ,81
+ ,125
+ ,118
+ ,116
+ ,115
+ ,116
+ ,114
+ ,123
+ ,2003
+ ,115
+ ,54
+ ,96
+ ,40
+ ,69
+ ,111
+ ,101
+ ,101
+ ,106
+ ,111
+ ,95
+ ,98
+ ,2003
+ ,118
+ ,55
+ ,80
+ ,44
+ ,68
+ ,98
+ ,101
+ ,108
+ ,115
+ ,120
+ ,85
+ ,81
+ ,2004
+ ,111
+ ,56
+ ,95
+ ,67
+ ,94
+ ,102
+ ,109
+ ,92
+ ,112
+ ,111
+ ,98
+ ,84
+ ,2004
+ ,108
+ ,56
+ ,94
+ ,75
+ ,97
+ ,105
+ ,108
+ ,89
+ ,106
+ ,115
+ ,100
+ ,92
+ ,2004
+ ,124
+ ,70
+ ,108
+ ,75
+ ,102
+ ,128
+ ,124
+ ,109
+ ,106
+ ,125
+ ,119
+ ,116
+ ,2004
+ ,115
+ ,69
+ ,97
+ ,49
+ ,94
+ ,125
+ ,117
+ ,97
+ ,114
+ ,116
+ ,109
+ ,112
+ ,2004
+ ,113
+ ,57
+ ,89
+ ,37
+ ,89
+ ,116
+ ,104
+ ,99
+ ,109
+ ,113
+ ,99
+ ,106
+ ,2004
+ ,128
+ ,68
+ ,107
+ ,50
+ ,114
+ ,131
+ ,121
+ ,110
+ ,100
+ ,122
+ ,119
+ ,131
+ ,2004
+ ,117
+ ,53
+ ,87
+ ,63
+ ,82
+ ,98
+ ,101
+ ,76
+ ,105
+ ,123
+ ,94
+ ,83
+ ,2004
+ ,119
+ ,48
+ ,70
+ ,76
+ ,96
+ ,89
+ ,105
+ ,91
+ ,100
+ ,117
+ ,88
+ ,98
+ ,2004
+ ,130
+ ,61
+ ,111
+ ,69
+ ,104
+ ,133
+ ,121
+ ,105
+ ,104
+ ,136
+ ,116
+ ,120
+ ,2004
+ ,126
+ ,62
+ ,105
+ ,49
+ ,88
+ ,114
+ ,116
+ ,103
+ ,112
+ ,121
+ ,109
+ ,121
+ ,2004
+ ,125
+ ,58
+ ,99
+ ,40
+ ,85
+ ,113
+ ,106
+ ,108
+ ,97
+ ,120
+ ,103
+ ,107
+ ,2004
+ ,131
+ ,51
+ ,84
+ ,39
+ ,87
+ ,104
+ ,105
+ ,122
+ ,107
+ ,126
+ ,93
+ ,89
+ ,2005
+ ,116
+ ,51
+ ,87
+ ,54
+ ,86
+ ,108
+ ,107
+ ,92
+ ,104
+ ,116
+ ,100
+ ,81
+ ,2005
+ ,109
+ ,48
+ ,92
+ ,71
+ ,89
+ ,106
+ ,101
+ ,95
+ ,98
+ ,108
+ ,102
+ ,90
+ ,2005
+ ,124
+ ,59
+ ,98
+ ,68
+ ,105
+ ,117
+ ,113
+ ,106
+ ,100
+ ,117
+ ,113
+ ,103
+ ,2005
+ ,119
+ ,54
+ ,95
+ ,43
+ ,83
+ ,123
+ ,109
+ ,98
+ ,97
+ ,113
+ ,112
+ ,117
+ ,2005
+ ,119
+ ,56
+ ,85
+ ,42
+ ,87
+ ,114
+ ,103
+ ,110
+ ,81
+ ,113
+ ,104
+ ,110
+ ,2005
+ ,131
+ ,60
+ ,100
+ ,48
+ ,112
+ ,132
+ ,116
+ ,107
+ ,73
+ ,126
+ ,118
+ ,130
+ ,2005
+ ,111
+ ,51
+ ,79
+ ,58
+ ,97
+ ,92
+ ,98
+ ,69
+ ,89
+ ,114
+ ,94
+ ,79
+ ,2005
+ ,125
+ ,51
+ ,66
+ ,76
+ ,89
+ ,94
+ ,99
+ ,95
+ ,96
+ ,113
+ ,95
+ ,101
+ ,2005
+ ,132
+ ,56
+ ,105
+ ,57
+ ,109
+ ,121
+ ,117
+ ,114
+ ,97
+ ,112
+ ,121
+ ,123
+ ,2005
+ ,127
+ ,53
+ ,96
+ ,44
+ ,88
+ ,114
+ ,107
+ ,104
+ ,98
+ ,113
+ ,114
+ ,111
+ ,2005
+ ,132
+ ,53
+ ,103
+ ,40
+ ,91
+ ,116
+ ,107
+ ,110
+ ,89
+ ,116
+ ,114
+ ,109
+ ,2005
+ ,131
+ ,48
+ ,83
+ ,36
+ ,79
+ ,98
+ ,102
+ ,112
+ ,98
+ ,112
+ ,99
+ ,89
+ ,2006
+ ,122
+ ,50
+ ,91
+ ,60
+ ,115
+ ,112
+ ,103
+ ,92
+ ,91
+ ,119
+ ,112
+ ,87
+ ,2006
+ ,113
+ ,49
+ ,95
+ ,73
+ ,119
+ ,109
+ ,101
+ ,97
+ ,86
+ ,117
+ ,111
+ ,95
+ ,2006
+ ,134
+ ,55
+ ,109
+ ,71
+ ,125
+ ,133
+ ,117
+ ,114
+ ,97
+ ,125
+ ,126
+ ,119
+ ,2006
+ ,119
+ ,50
+ ,92
+ ,45
+ ,96
+ ,118
+ ,103
+ ,93
+ ,102
+ ,113
+ ,112
+ ,110
+ ,2006
+ ,129
+ ,57
+ ,99
+ ,45
+ ,117
+ ,132
+ ,106
+ ,115
+ ,80
+ ,120
+ ,124
+ ,124
+ ,2006
+ ,131
+ ,65
+ ,110
+ ,48
+ ,120
+ ,134
+ ,111
+ ,112
+ ,71
+ ,114
+ ,127
+ ,133
+ ,2006
+ ,117
+ ,53
+ ,88
+ ,60
+ ,104
+ ,97
+ ,94
+ ,76
+ ,91
+ ,114
+ ,101
+ ,84
+ ,2006
+ ,131
+ ,42
+ ,73
+ ,72
+ ,121
+ ,100
+ ,101
+ ,101
+ ,102
+ ,118
+ ,102
+ ,105
+ ,2006
+ ,132
+ ,56
+ ,111
+ ,63
+ ,127
+ ,128
+ ,111
+ ,119
+ ,91
+ ,117
+ ,126
+ ,128
+ ,2006
+ ,141
+ ,58
+ ,112
+ ,32
+ ,118
+ ,135
+ ,114
+ ,118
+ ,94
+ ,121
+ ,129
+ ,127
+ ,2006
+ ,138
+ ,54
+ ,111
+ ,34
+ ,108
+ ,131
+ ,110
+ ,120
+ ,53
+ ,115
+ ,122
+ ,120
+ ,2006
+ ,129
+ ,51
+ ,84
+ ,24
+ ,89
+ ,107
+ ,100
+ ,120
+ ,77
+ ,117
+ ,100
+ ,93
+ ,2007
+ ,127
+ ,59
+ ,102
+ ,65
+ ,137
+ ,122
+ ,104
+ ,99
+ ,70
+ ,119
+ ,122
+ ,98
+ ,2007
+ ,121
+ ,49
+ ,102
+ ,73
+ ,142
+ ,121
+ ,106
+ ,103
+ ,65
+ ,115
+ ,120
+ ,106
+ ,2007
+ ,139
+ ,61
+ ,114
+ ,62
+ ,137
+ ,141
+ ,116
+ ,118
+ ,89
+ ,126
+ ,137
+ ,122
+ ,2007
+ ,129
+ ,52
+ ,99
+ ,32
+ ,123
+ ,125
+ ,104
+ ,103
+ ,70
+ ,118
+ ,124
+ ,116
+ ,2007
+ ,131
+ ,58
+ ,100
+ ,31
+ ,126
+ ,130
+ ,107
+ ,114
+ ,78
+ ,118
+ ,130
+ ,122
+ ,2007
+ ,136
+ ,66
+ ,110
+ ,37
+ ,148
+ ,159
+ ,113
+ ,116
+ ,78
+ ,115
+ ,137
+ ,134
+ ,2007
+ ,129
+ ,62
+ ,93
+ ,48
+ ,116
+ ,111
+ ,104
+ ,84
+ ,73
+ ,122
+ ,114
+ ,88
+ ,2007
+ ,133
+ ,45
+ ,77
+ ,54
+ ,139
+ ,110
+ ,103
+ ,106
+ ,83
+ ,117
+ ,109
+ ,110
+ ,2007
+ ,136
+ ,52
+ ,108
+ ,44
+ ,151
+ ,133
+ ,109
+ ,117
+ ,74
+ ,106
+ ,126
+ ,122
+ ,2007
+ ,151
+ ,59
+ ,120
+ ,41
+ ,124
+ ,135
+ ,123
+ ,125
+ ,102
+ ,111
+ ,141
+ ,135
+ ,2007
+ ,145
+ ,58
+ ,106
+ ,32
+ ,109
+ ,119
+ ,110
+ ,123
+ ,54
+ ,114
+ ,130
+ ,116
+ ,2007
+ ,134
+ ,45
+ ,78
+ ,31
+ ,112
+ ,94
+ ,94
+ ,119
+ ,79
+ ,114
+ ,98
+ ,85
+ ,2008
+ ,136
+ ,65
+ ,100
+ ,49
+ ,136
+ ,118
+ ,114
+ ,100
+ ,86
+ ,125
+ ,130
+ ,106
+ ,2008
+ ,129
+ ,64
+ ,102
+ ,54
+ ,136
+ ,115
+ ,110
+ ,100
+ ,87
+ ,125
+ ,130
+ ,115
+ ,2008
+ ,129
+ ,69
+ ,97
+ ,44
+ ,139
+ ,114
+ ,110
+ ,103
+ ,79
+ ,120
+ ,125
+ ,111
+ ,2008
+ ,139
+ ,71
+ ,101
+ ,31
+ ,138
+ ,131
+ ,113
+ ,104
+ ,64
+ ,121
+ ,136
+ ,133
+ ,2008
+ ,133
+ ,63
+ ,89
+ ,24
+ ,142
+ ,117
+ ,105
+ ,99
+ ,70
+ ,111
+ ,124
+ ,124
+ ,2008
+ ,133
+ ,74
+ ,93
+ ,37
+ ,144
+ ,123
+ ,108
+ ,101
+ ,75
+ ,124
+ ,133
+ ,131
+ ,2008
+ ,137
+ ,63
+ ,89
+ ,38
+ ,147
+ ,106
+ ,101
+ ,73
+ ,72
+ ,120
+ ,121
+ ,97
+ ,2008
+ ,127
+ ,52
+ ,62
+ ,42
+ ,201
+ ,89
+ ,95
+ ,86
+ ,83
+ ,126
+ ,102
+ ,97
+ ,2008
+ ,144
+ ,73
+ ,96
+ ,36
+ ,196
+ ,116
+ ,112
+ ,110
+ ,74
+ ,116
+ ,131
+ ,131
+ ,2008
+ ,150
+ ,67
+ ,95
+ ,31
+ ,170
+ ,116
+ ,113
+ ,115
+ ,82
+ ,117
+ ,130
+ ,127
+ ,2008
+ ,132
+ ,63
+ ,80
+ ,24
+ ,177
+ ,97
+ ,96
+ ,101
+ ,78
+ ,106
+ ,106
+ ,101
+ ,2008
+ ,139
+ ,70
+ ,67
+ ,29
+ ,190
+ ,82
+ ,93
+ ,112
+ ,77
+ ,102
+ ,93
+ ,88
+ ,2009
+ ,123
+ ,66
+ ,71
+ ,38
+ ,138
+ ,92
+ ,91
+ ,89
+ ,77
+ ,106
+ ,100
+ ,76
+ ,2009
+ ,122
+ ,60
+ ,73
+ ,44
+ ,133
+ ,90
+ ,91
+ ,93
+ ,72
+ ,97
+ ,99
+ ,87
+ ,2009
+ ,136
+ ,66
+ ,81
+ ,33
+ ,131
+ ,99
+ ,101
+ ,103
+ ,76
+ ,108
+ ,112
+ ,110
+ ,2009
+ ,133
+ ,68
+ ,77
+ ,23
+ ,110
+ ,99
+ ,98
+ ,91
+ ,75
+ ,99
+ ,109
+ ,102
+ ,2009
+ ,127
+ ,68
+ ,68
+ ,19
+ ,124
+ ,89
+ ,94
+ ,88
+ ,69
+ ,101
+ ,102
+ ,99
+ ,2009
+ ,139
+ ,81
+ ,77
+ ,27
+ ,150
+ ,106
+ ,102
+ ,93
+ ,67
+ ,106
+ ,116
+ ,117
+ ,2009
+ ,131
+ ,75
+ ,73
+ ,29
+ ,163
+ ,84
+ ,96
+ ,65
+ ,68
+ ,105
+ ,103
+ ,83
+ ,2009
+ ,132
+ ,55
+ ,54
+ ,34
+ ,138
+ ,78
+ ,92
+ ,82
+ ,73
+ ,103
+ ,91
+ ,90
+ ,2009
+ ,136
+ ,79
+ ,85
+ ,26
+ ,133
+ ,101
+ ,106
+ ,102
+ ,69
+ ,102
+ ,119
+ ,116
+ ,2009
+ ,142
+ ,52
+ ,86
+ ,28
+ ,123
+ ,100
+ ,105
+ ,102
+ ,76
+ ,107
+ ,117
+ ,117
+ ,2009
+ ,133
+ ,56
+ ,79
+ ,18
+ ,107
+ ,96
+ ,97
+ ,122
+ ,67
+ ,100
+ ,106
+ ,96
+ ,2009
+ ,132
+ ,66
+ ,67
+ ,24
+ ,122
+ ,80
+ ,94
+ ,105
+ ,69
+ ,101
+ ,92
+ ,73
+ ,2010
+ ,121
+ ,66
+ ,72
+ ,29
+ ,141
+ ,87
+ ,95
+ ,83
+ ,68
+ ,105
+ ,102
+ ,66
+ ,2010
+ ,124
+ ,59
+ ,76
+ ,38
+ ,136
+ ,90
+ ,95
+ ,85
+ ,64
+ ,118
+ ,104
+ ,73
+ ,2010
+ ,145
+ ,78
+ ,90
+ ,33
+ ,140
+ ,113
+ ,114
+ ,102
+ ,69
+ ,129
+ ,124
+ ,114
+ ,2010
+ ,135
+ ,70
+ ,84
+ ,22
+ ,109
+ ,105
+ ,107
+ ,86
+ ,67
+ ,124
+ ,118
+ ,107
+ ,2010
+ ,128
+ ,65
+ ,75
+ ,20
+ ,109
+ ,100
+ ,100
+ ,84
+ ,71
+ ,128
+ ,109
+ ,102
+ ,2010
+ ,142
+ ,88
+ ,90
+ ,31
+ ,128
+ ,116
+ ,112
+ ,93
+ ,58
+ ,129
+ ,129
+ ,125
+ ,2010
+ ,130
+ ,75
+ ,77
+ ,27
+ ,162
+ ,89
+ ,101
+ ,64
+ ,57
+ ,128
+ ,105
+ ,80
+ ,2010
+ ,131
+ ,62
+ ,60
+ ,28
+ ,147
+ ,87
+ ,100
+ ,81
+ ,69
+ ,125
+ ,100
+ ,95
+ ,2010
+ ,141
+ ,85
+ ,92
+ ,28
+ ,148
+ ,111
+ ,111
+ ,100
+ ,76
+ ,125
+ ,125
+ ,120
+ ,2010
+ ,140
+ ,82
+ ,88
+ ,25
+ ,103
+ ,110
+ ,107
+ ,96
+ ,74
+ ,130
+ ,116
+ ,117
+ ,2010
+ ,142
+ ,83
+ ,83
+ ,21
+ ,102
+ ,104
+ ,105
+ ,93
+ ,77
+ ,125
+ ,112
+ ,99
+ ,2010
+ ,140
+ ,78
+ ,69
+ ,24
+ ,100
+ ,85
+ ,104
+ ,102
+ ,81
+ ,122
+ ,97
+ ,64
+ ,2011
+ ,132
+ ,81
+ ,73
+ ,28
+ ,117
+ ,96
+ ,106
+ ,78
+ ,77
+ ,129
+ ,107
+ ,82
+ ,2011
+ ,132
+ ,75
+ ,78
+ ,33
+ ,139
+ ,99
+ ,105
+ ,92
+ ,64
+ ,124
+ ,114
+ ,97
+ ,2011
+ ,151
+ ,91
+ ,92
+ ,31
+ ,122
+ ,117
+ ,114
+ ,99
+ ,67
+ ,144
+ ,130
+ ,121)
+ ,dim=c(13
+ ,123)
+ ,dimnames=list(c('Jaar'
+ ,'Voedingsmiddelen'
+ ,'Tabaksproducten'
+ ,'Textiel'
+ ,'Kleding'
+ ,'Leer'
+ ,'Hout'
+ ,'Papier'
+ ,'Uitgeverijen'
+ ,'Cokes'
+ ,'Chemische'
+ ,'Rubber'
+ ,'Nietmetaalhoudende')
+ ,1:123))
> y <- array(NA,dim=c(13,123),dimnames=list(c('Jaar','Voedingsmiddelen','Tabaksproducten','Textiel','Kleding','Leer','Hout','Papier','Uitgeverijen','Cokes','Chemische','Rubber','Nietmetaalhoudende'),1:123))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '4'
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from 'package:base':
as.Date, as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Textiel Jaar Voedingsmiddelen Tabaksproducten Kleding Leer Hout Papier
1 102 2001 100 95 103 91 99 101
2 99 2001 94 97 117 85 97 97
3 108 2001 105 97 115 110 113 108
4 92 2001 95 97 74 90 100 95
5 99 2001 103 103 74 103 105 99
6 102 2001 103 101 81 119 109 101
7 87 2001 100 96 86 76 91 92
8 71 2001 108 94 114 93 89 92
9 105 2001 108 97 102 105 105 100
10 115 2001 120 101 85 92 120 106
11 103 2001 112 77 63 75 107 99
12 75 2001 102 93 61 61 84 84
13 97 2002 105 45 87 80 101 106
14 95 2002 101 48 97 85 105 101
15 99 2002 108 52 88 94 119 113
16 100 2002 107 49 67 78 114 110
17 92 2002 109 53 59 92 114 103
18 94 2002 110 60 63 90 119 107
19 89 2002 111 51 86 72 99 98
20 67 2002 110 42 99 77 91 90
21 109 2002 117 56 85 76 121 105
22 113 2002 130 51 74 89 128 116
23 106 2002 114 53 55 55 112 102
24 78 2002 113 55 54 47 93 88
25 102 2003 110 44 81 91 108 114
26 97 2003 107 51 88 85 107 104
27 96 2003 110 52 75 89 115 111
28 99 2003 113 54 55 90 121 111
29 86 2003 106 50 47 72 112 102
30 92 2003 118 57 54 83 123 106
31 86 2003 118 49 71 72 101 104
32 62 2003 114 41 79 75 87 94
33 105 2003 121 58 77 85 124 116
34 108 2003 130 63 57 81 125 118
35 96 2003 115 54 40 69 111 101
36 80 2003 118 55 44 68 98 101
37 95 2004 111 56 67 94 102 109
38 94 2004 108 56 75 97 105 108
39 108 2004 124 70 75 102 128 124
40 97 2004 115 69 49 94 125 117
41 89 2004 113 57 37 89 116 104
42 107 2004 128 68 50 114 131 121
43 87 2004 117 53 63 82 98 101
44 70 2004 119 48 76 96 89 105
45 111 2004 130 61 69 104 133 121
46 105 2004 126 62 49 88 114 116
47 99 2004 125 58 40 85 113 106
48 84 2004 131 51 39 87 104 105
49 87 2005 116 51 54 86 108 107
50 92 2005 109 48 71 89 106 101
51 98 2005 124 59 68 105 117 113
52 95 2005 119 54 43 83 123 109
53 85 2005 119 56 42 87 114 103
54 100 2005 131 60 48 112 132 116
55 79 2005 111 51 58 97 92 98
56 66 2005 125 51 76 89 94 99
57 105 2005 132 56 57 109 121 117
58 96 2005 127 53 44 88 114 107
59 103 2005 132 53 40 91 116 107
60 83 2005 131 48 36 79 98 102
61 91 2006 122 50 60 115 112 103
62 95 2006 113 49 73 119 109 101
63 109 2006 134 55 71 125 133 117
64 92 2006 119 50 45 96 118 103
65 99 2006 129 57 45 117 132 106
66 110 2006 131 65 48 120 134 111
67 88 2006 117 53 60 104 97 94
68 73 2006 131 42 72 121 100 101
69 111 2006 132 56 63 127 128 111
70 112 2006 141 58 32 118 135 114
71 111 2006 138 54 34 108 131 110
72 84 2006 129 51 24 89 107 100
73 102 2007 127 59 65 137 122 104
74 102 2007 121 49 73 142 121 106
75 114 2007 139 61 62 137 141 116
76 99 2007 129 52 32 123 125 104
77 100 2007 131 58 31 126 130 107
78 110 2007 136 66 37 148 159 113
79 93 2007 129 62 48 116 111 104
80 77 2007 133 45 54 139 110 103
81 108 2007 136 52 44 151 133 109
82 120 2007 151 59 41 124 135 123
83 106 2007 145 58 32 109 119 110
84 78 2007 134 45 31 112 94 94
85 100 2008 136 65 49 136 118 114
86 102 2008 129 64 54 136 115 110
87 97 2008 129 69 44 139 114 110
88 101 2008 139 71 31 138 131 113
89 89 2008 133 63 24 142 117 105
90 93 2008 133 74 37 144 123 108
91 89 2008 137 63 38 147 106 101
92 62 2008 127 52 42 201 89 95
93 96 2008 144 73 36 196 116 112
94 95 2008 150 67 31 170 116 113
95 80 2008 132 63 24 177 97 96
96 67 2008 139 70 29 190 82 93
97 71 2009 123 66 38 138 92 91
98 73 2009 122 60 44 133 90 91
99 81 2009 136 66 33 131 99 101
100 77 2009 133 68 23 110 99 98
101 68 2009 127 68 19 124 89 94
102 77 2009 139 81 27 150 106 102
103 73 2009 131 75 29 163 84 96
104 54 2009 132 55 34 138 78 92
105 85 2009 136 79 26 133 101 106
106 86 2009 142 52 28 123 100 105
107 79 2009 133 56 18 107 96 97
108 67 2009 132 66 24 122 80 94
109 72 2010 121 66 29 141 87 95
110 76 2010 124 59 38 136 90 95
111 90 2010 145 78 33 140 113 114
112 84 2010 135 70 22 109 105 107
113 75 2010 128 65 20 109 100 100
114 90 2010 142 88 31 128 116 112
115 77 2010 130 75 27 162 89 101
116 60 2010 131 62 28 147 87 100
117 92 2010 141 85 28 148 111 111
118 88 2010 140 82 25 103 110 107
119 83 2010 142 83 21 102 104 105
120 69 2010 140 78 24 100 85 104
121 73 2011 132 81 28 117 96 106
122 78 2011 132 75 33 139 99 105
123 92 2011 151 91 31 122 117 114
Uitgeverijen Cokes Chemische Rubber Nietmetaalhoudende t
1 91 114 101 103 85 1
2 87 99 99 97 94 2
3 103 98 104 110 107 3
4 97 91 99 97 98 4
5 96 111 101 103 111 5
6 105 104 102 106 115 6
7 74 100 93 89 76 7
8 87 108 97 85 100 8
9 105 113 91 100 103 9
10 118 113 97 106 117 10
11 102 114 94 95 101 11
12 101 109 90 74 73 12
13 86 116 105 94 84 13
14 83 102 103 90 90 14
15 92 107 112 99 105 15
16 87 111 114 100 111 16
17 94 122 111 96 110 17
18 94 123 106 102 116 18
19 75 108 112 88 85 19
20 85 115 102 78 92 20
21 104 120 103 99 117 21
22 109 117 105 107 119 22
23 121 115 101 93 100 23
24 124 116 101 74 71 24
25 88 118 117 96 82 25
26 86 98 109 99 90 26
27 98 121 120 103 109 27
28 94 118 115 102 112 28
29 102 120 107 96 103 29
30 96 111 110 106 116 30
31 79 117 110 95 89 31
32 95 110 105 82 91 32
33 106 107 116 109 121 33
34 116 115 116 114 123 34
35 101 106 111 95 98 35
36 108 115 120 85 81 36
37 92 112 111 98 84 37
38 89 106 115 100 92 38
39 109 106 125 119 116 39
40 97 114 116 109 112 40
41 99 109 113 99 106 41
42 110 100 122 119 131 42
43 76 105 123 94 83 43
44 91 100 117 88 98 44
45 105 104 136 116 120 45
46 103 112 121 109 121 46
47 108 97 120 103 107 47
48 122 107 126 93 89 48
49 92 104 116 100 81 49
50 95 98 108 102 90 50
51 106 100 117 113 103 51
52 98 97 113 112 117 52
53 110 81 113 104 110 53
54 107 73 126 118 130 54
55 69 89 114 94 79 55
56 95 96 113 95 101 56
57 114 97 112 121 123 57
58 104 98 113 114 111 58
59 110 89 116 114 109 59
60 112 98 112 99 89 60
61 92 91 119 112 87 61
62 97 86 117 111 95 62
63 114 97 125 126 119 63
64 93 102 113 112 110 64
65 115 80 120 124 124 65
66 112 71 114 127 133 66
67 76 91 114 101 84 67
68 101 102 118 102 105 68
69 119 91 117 126 128 69
70 118 94 121 129 127 70
71 120 53 115 122 120 71
72 120 77 117 100 93 72
73 99 70 119 122 98 73
74 103 65 115 120 106 74
75 118 89 126 137 122 75
76 103 70 118 124 116 76
77 114 78 118 130 122 77
78 116 78 115 137 134 78
79 84 73 122 114 88 79
80 106 83 117 109 110 80
81 117 74 106 126 122 81
82 125 102 111 141 135 82
83 123 54 114 130 116 83
84 119 79 114 98 85 84
85 100 86 125 130 106 85
86 100 87 125 130 115 86
87 103 79 120 125 111 87
88 104 64 121 136 133 88
89 99 70 111 124 124 89
90 101 75 124 133 131 90
91 73 72 120 121 97 91
92 86 83 126 102 97 92
93 110 74 116 131 131 93
94 115 82 117 130 127 94
95 101 78 106 106 101 95
96 112 77 102 93 88 96
97 89 77 106 100 76 97
98 93 72 97 99 87 98
99 103 76 108 112 110 99
100 91 75 99 109 102 100
101 88 69 101 102 99 101
102 93 67 106 116 117 102
103 65 68 105 103 83 103
104 82 73 103 91 90 104
105 102 69 102 119 116 105
106 102 76 107 117 117 106
107 122 67 100 106 96 107
108 105 69 101 92 73 108
109 83 68 105 102 66 109
110 85 64 118 104 73 110
111 102 69 129 124 114 111
112 86 67 124 118 107 112
113 84 71 128 109 102 113
114 93 58 129 129 125 114
115 64 57 128 105 80 115
116 81 69 125 100 95 116
117 100 76 125 125 120 117
118 96 74 130 116 117 118
119 93 77 125 112 99 119
120 102 81 122 97 64 120
121 78 77 129 107 82 121
122 92 64 124 114 97 122
123 99 67 144 130 121 123
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Jaar Voedingsmiddelen Tabaksproducten
2.293e+04 -1.146e+01 -3.154e-01 4.050e-02
Kleding Leer Hout Papier
4.960e-02 -3.973e-02 3.919e-01 6.072e-01
Uitgeverijen Cokes Chemische Rubber
3.097e-02 -1.757e-02 -2.258e-01 7.955e-01
Nietmetaalhoudende t
-3.341e-01 8.411e-01
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-10.630 -2.980 0.003 2.824 10.867
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.293e+04 4.876e+03 4.704 7.54e-06 ***
Jaar -1.146e+01 2.436e+00 -4.705 7.48e-06 ***
Voedingsmiddelen -3.154e-01 1.154e-01 -2.734 0.007312 **
Tabaksproducten 4.050e-02 3.106e-02 1.304 0.194925
Kleding 4.960e-02 3.576e-02 1.387 0.168305
Leer -3.973e-02 2.255e-02 -1.762 0.080863 .
Hout 3.919e-01 8.431e-02 4.649 9.41e-06 ***
Papier 6.072e-01 1.180e-01 5.146 1.19e-06 ***
Uitgeverijen 3.097e-02 4.532e-02 0.683 0.495875
Cokes -1.757e-02 4.846e-02 -0.363 0.717626
Chemische -2.258e-01 6.264e-02 -3.604 0.000473 ***
Rubber 7.955e-01 1.044e-01 7.619 1.00e-11 ***
Nietmetaalhoudende -3.342e-01 4.898e-02 -6.822 5.24e-10 ***
t 8.411e-01 2.299e-01 3.659 0.000392 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 4.23 on 109 degrees of freedom
Multiple R-squared: 0.9181, Adjusted R-squared: 0.9083
F-statistic: 94.01 on 13 and 109 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 0.2921235 0.58424694 0.70787653
[2,] 0.8793700 0.24125998 0.12062999
[3,] 0.8360713 0.32785737 0.16392869
[4,] 0.8743623 0.25127543 0.12563771
[5,] 0.8640909 0.27181821 0.13590910
[6,] 0.8376799 0.32464010 0.16232005
[7,] 0.8285252 0.34294960 0.17147480
[8,] 0.7920656 0.41586874 0.20793437
[9,] 0.8473979 0.30520412 0.15260206
[10,] 0.8272841 0.34543180 0.17271590
[11,] 0.8126458 0.37470842 0.18735421
[12,] 0.7717260 0.45654796 0.22827398
[13,] 0.8580268 0.28394646 0.14197323
[14,] 0.8857830 0.22843391 0.11421696
[15,] 0.8497302 0.30053956 0.15026978
[16,] 0.8499787 0.30004266 0.15002133
[17,] 0.8242968 0.35140645 0.17570323
[18,] 0.7766872 0.44662569 0.22331285
[19,] 0.7868241 0.42635170 0.21317585
[20,] 0.7456152 0.50876954 0.25438477
[21,] 0.8588299 0.28234010 0.14117005
[22,] 0.8496136 0.30077279 0.15038640
[23,] 0.8380553 0.32388950 0.16194475
[24,] 0.8644866 0.27102683 0.13551342
[25,] 0.8456833 0.30863344 0.15431672
[26,] 0.8218258 0.35634839 0.17817419
[27,] 0.8740583 0.25188333 0.12594167
[28,] 0.8690840 0.26183205 0.13091603
[29,] 0.8417504 0.31649915 0.15824957
[30,] 0.9470742 0.10585160 0.05292580
[31,] 0.9699943 0.06001130 0.03000565
[32,] 0.9626391 0.07472172 0.03736086
[33,] 0.9496839 0.10063221 0.05031610
[34,] 0.9362732 0.12745357 0.06372678
[35,] 0.9185072 0.16298561 0.08149280
[36,] 0.8945715 0.21085707 0.10542853
[37,] 0.8679636 0.26407274 0.13203637
[38,] 0.8357622 0.32847561 0.16423780
[39,] 0.8037862 0.39242753 0.19621376
[40,] 0.8970654 0.20586927 0.10293464
[41,] 0.8853192 0.22936158 0.11468079
[42,] 0.8751212 0.24975760 0.12487880
[43,] 0.8684386 0.26312275 0.13156138
[44,] 0.8452783 0.30944346 0.15472173
[45,] 0.8192507 0.36149854 0.18074927
[46,] 0.7915831 0.41683387 0.20841694
[47,] 0.7534343 0.49313148 0.24656574
[48,] 0.7125322 0.57493557 0.28746778
[49,] 0.6863556 0.62728879 0.31364439
[50,] 0.6899039 0.62019222 0.31009611
[51,] 0.7651732 0.46965358 0.23482679
[52,] 0.7846572 0.43068566 0.21534283
[53,] 0.7854585 0.42908305 0.21454152
[54,] 0.8044764 0.39104718 0.19552359
[55,] 0.8421009 0.31579828 0.15789914
[56,] 0.8271252 0.34574965 0.17287483
[57,] 0.7949419 0.41011615 0.20505808
[58,] 0.7845411 0.43091773 0.21545886
[59,] 0.7448598 0.51028035 0.25514018
[60,] 0.7409995 0.51800107 0.25900053
[61,] 0.7218330 0.55633410 0.27816705
[62,] 0.8724952 0.25500961 0.12750480
[63,] 0.8406538 0.31869246 0.15934623
[64,] 0.8752375 0.24952506 0.12476253
[65,] 0.8425329 0.31493428 0.15746714
[66,] 0.8267480 0.34650406 0.17325203
[67,] 0.7844952 0.43100967 0.21550484
[68,] 0.7844251 0.43114984 0.21557492
[69,] 0.7451740 0.50965198 0.25482599
[70,] 0.7916561 0.41668790 0.20834395
[71,] 0.8728212 0.25435752 0.12717876
[72,] 0.8333240 0.33335209 0.16667604
[73,] 0.7913437 0.41731260 0.20865630
[74,] 0.7901423 0.41971542 0.20985771
[75,] 0.7515992 0.49680167 0.24840083
[76,] 0.8327566 0.33448677 0.16724339
[77,] 0.7937589 0.41248214 0.20624107
[78,] 0.8510742 0.29785169 0.14892585
[79,] 0.8163911 0.36721780 0.18360890
[80,] 0.8202565 0.35948698 0.17974349
[81,] 0.8526200 0.29476004 0.14738002
[82,] 0.8347256 0.33054884 0.16527442
[83,] 0.8608063 0.27838746 0.13919373
[84,] 0.7947833 0.41043336 0.20521668
[85,] 0.7507975 0.49840505 0.24920253
[86,] 0.7717112 0.45657767 0.22828883
[87,] 0.6682246 0.66355072 0.33177536
[88,] 0.5735990 0.85280193 0.42640097
[89,] 0.4390621 0.87812423 0.56093788
[90,] 0.3570045 0.71400901 0.64299549
> postscript(file="/var/wessaorg/rcomp/tmp/1ki461353238880.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/wessaorg/rcomp/tmp/2mb9q1353238880.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/wessaorg/rcomp/tmp/3fnbq1353238880.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/wessaorg/rcomp/tmp/4szlt1353238880.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/wessaorg/rcomp/tmp/5o1ga1353238880.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> qqnorm(mysum$resid, main='Residual Normal Q-Q Plot')
> qqline(mysum$resid)
> grid()
> dev.off()
null device
1
> (myerror <- as.ts(mysum$resid))
Time Series:
Start = 1
End = 123
Frequency = 1
1 2 3 4 5
-0.990321917 2.663998492 -2.946816455 -2.446177309 2.526076307
6 7 8 9 10
1.046513058 -5.625858534 -7.949504117 2.277068241 6.720882737
11 12 13 14 15
5.329583336 -3.259028071 3.436546848 4.965514239 -2.105260996
16 17 18 19 20
3.933212096 2.911550925 -4.443419639 4.300925932 -3.094725751
21 22 23 24 25
10.867365047 4.516820091 9.389689998 1.180371370 4.980704231
26 27 28 29 30
1.979536117 0.145010381 2.589523045 -4.876754819 -5.772341610
31 32 33 34 35
-3.372591583 -8.607574158 -1.656348838 -1.112587757 3.793783887
36 37 38 39 40
-3.036833302 2.874610334 1.213992589 -5.134740401 -8.277657918
41 42 43 44 45
-0.322821488 -0.754066300 3.919384723 -4.200621324 1.610580424
46 47 48 49 50
7.024588715 6.272236547 -0.092118699 -4.391007689 0.796409135
51 52 53 54 55
-3.252081842 -3.258965161 -3.429714253 -1.466918544 -2.938188168
56 57 58 59 60
-9.317248872 -3.055967373 -3.614636952 3.320299941 -3.378919691
61 62 63 64 65
-2.464193839 2.575143187 0.891052662 -0.803107321 -2.602784868
66 67 68 69 70
3.214238072 5.812190759 -3.969830358 4.193382072 1.991298084
71 72 73 74 75
3.958340227 -2.028463602 2.352872954 2.152774486 -0.796450097
76 77 78 79 80
1.722942259 -4.310911233 -10.629662935 -2.647065901 -6.228538865
81 82 83 84 85
0.412296504 -0.404977842 -0.785295630 2.773067934 -1.517947615
86 87 88 89 90
3.855147909 0.097048796 -3.017874592 -2.035742284 -5.919974305
91 92 93 94 95
0.025153151 -2.012027443 -0.014637251 -1.441917598 3.879492823
96 97 98 99 100
4.667866054 1.224839981 4.826193337 6.613040707 0.263286367
101 102 103 104 105
0.637470064 -3.012567782 2.188163491 -1.993613616 -3.363468675
106 107 108 109 110
3.461299035 -1.719577418 -2.651517005 -2.070439685 4.048629571
111 112 113 114 115
2.752677406 1.546707617 2.533624412 -1.676948593 4.743763725
116 117 118 119 120
-3.513206668 1.823669460 5.345127820 0.003030483 -8.089634858
121 122 123
-0.848281066 1.264075998 6.315792280
> postscript(file="/var/wessaorg/rcomp/tmp/6xpsk1353238880.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> dum <- cbind(lag(myerror,k=1),myerror)
> dum
Time Series:
Start = 0
End = 123
Frequency = 1
lag(myerror, k = 1) myerror
0 -0.990321917 NA
1 2.663998492 -0.990321917
2 -2.946816455 2.663998492
3 -2.446177309 -2.946816455
4 2.526076307 -2.446177309
5 1.046513058 2.526076307
6 -5.625858534 1.046513058
7 -7.949504117 -5.625858534
8 2.277068241 -7.949504117
9 6.720882737 2.277068241
10 5.329583336 6.720882737
11 -3.259028071 5.329583336
12 3.436546848 -3.259028071
13 4.965514239 3.436546848
14 -2.105260996 4.965514239
15 3.933212096 -2.105260996
16 2.911550925 3.933212096
17 -4.443419639 2.911550925
18 4.300925932 -4.443419639
19 -3.094725751 4.300925932
20 10.867365047 -3.094725751
21 4.516820091 10.867365047
22 9.389689998 4.516820091
23 1.180371370 9.389689998
24 4.980704231 1.180371370
25 1.979536117 4.980704231
26 0.145010381 1.979536117
27 2.589523045 0.145010381
28 -4.876754819 2.589523045
29 -5.772341610 -4.876754819
30 -3.372591583 -5.772341610
31 -8.607574158 -3.372591583
32 -1.656348838 -8.607574158
33 -1.112587757 -1.656348838
34 3.793783887 -1.112587757
35 -3.036833302 3.793783887
36 2.874610334 -3.036833302
37 1.213992589 2.874610334
38 -5.134740401 1.213992589
39 -8.277657918 -5.134740401
40 -0.322821488 -8.277657918
41 -0.754066300 -0.322821488
42 3.919384723 -0.754066300
43 -4.200621324 3.919384723
44 1.610580424 -4.200621324
45 7.024588715 1.610580424
46 6.272236547 7.024588715
47 -0.092118699 6.272236547
48 -4.391007689 -0.092118699
49 0.796409135 -4.391007689
50 -3.252081842 0.796409135
51 -3.258965161 -3.252081842
52 -3.429714253 -3.258965161
53 -1.466918544 -3.429714253
54 -2.938188168 -1.466918544
55 -9.317248872 -2.938188168
56 -3.055967373 -9.317248872
57 -3.614636952 -3.055967373
58 3.320299941 -3.614636952
59 -3.378919691 3.320299941
60 -2.464193839 -3.378919691
61 2.575143187 -2.464193839
62 0.891052662 2.575143187
63 -0.803107321 0.891052662
64 -2.602784868 -0.803107321
65 3.214238072 -2.602784868
66 5.812190759 3.214238072
67 -3.969830358 5.812190759
68 4.193382072 -3.969830358
69 1.991298084 4.193382072
70 3.958340227 1.991298084
71 -2.028463602 3.958340227
72 2.352872954 -2.028463602
73 2.152774486 2.352872954
74 -0.796450097 2.152774486
75 1.722942259 -0.796450097
76 -4.310911233 1.722942259
77 -10.629662935 -4.310911233
78 -2.647065901 -10.629662935
79 -6.228538865 -2.647065901
80 0.412296504 -6.228538865
81 -0.404977842 0.412296504
82 -0.785295630 -0.404977842
83 2.773067934 -0.785295630
84 -1.517947615 2.773067934
85 3.855147909 -1.517947615
86 0.097048796 3.855147909
87 -3.017874592 0.097048796
88 -2.035742284 -3.017874592
89 -5.919974305 -2.035742284
90 0.025153151 -5.919974305
91 -2.012027443 0.025153151
92 -0.014637251 -2.012027443
93 -1.441917598 -0.014637251
94 3.879492823 -1.441917598
95 4.667866054 3.879492823
96 1.224839981 4.667866054
97 4.826193337 1.224839981
98 6.613040707 4.826193337
99 0.263286367 6.613040707
100 0.637470064 0.263286367
101 -3.012567782 0.637470064
102 2.188163491 -3.012567782
103 -1.993613616 2.188163491
104 -3.363468675 -1.993613616
105 3.461299035 -3.363468675
106 -1.719577418 3.461299035
107 -2.651517005 -1.719577418
108 -2.070439685 -2.651517005
109 4.048629571 -2.070439685
110 2.752677406 4.048629571
111 1.546707617 2.752677406
112 2.533624412 1.546707617
113 -1.676948593 2.533624412
114 4.743763725 -1.676948593
115 -3.513206668 4.743763725
116 1.823669460 -3.513206668
117 5.345127820 1.823669460
118 0.003030483 5.345127820
119 -8.089634858 0.003030483
120 -0.848281066 -8.089634858
121 1.264075998 -0.848281066
122 6.315792280 1.264075998
123 NA 6.315792280
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2.663998492 -0.990321917
[2,] -2.946816455 2.663998492
[3,] -2.446177309 -2.946816455
[4,] 2.526076307 -2.446177309
[5,] 1.046513058 2.526076307
[6,] -5.625858534 1.046513058
[7,] -7.949504117 -5.625858534
[8,] 2.277068241 -7.949504117
[9,] 6.720882737 2.277068241
[10,] 5.329583336 6.720882737
[11,] -3.259028071 5.329583336
[12,] 3.436546848 -3.259028071
[13,] 4.965514239 3.436546848
[14,] -2.105260996 4.965514239
[15,] 3.933212096 -2.105260996
[16,] 2.911550925 3.933212096
[17,] -4.443419639 2.911550925
[18,] 4.300925932 -4.443419639
[19,] -3.094725751 4.300925932
[20,] 10.867365047 -3.094725751
[21,] 4.516820091 10.867365047
[22,] 9.389689998 4.516820091
[23,] 1.180371370 9.389689998
[24,] 4.980704231 1.180371370
[25,] 1.979536117 4.980704231
[26,] 0.145010381 1.979536117
[27,] 2.589523045 0.145010381
[28,] -4.876754819 2.589523045
[29,] -5.772341610 -4.876754819
[30,] -3.372591583 -5.772341610
[31,] -8.607574158 -3.372591583
[32,] -1.656348838 -8.607574158
[33,] -1.112587757 -1.656348838
[34,] 3.793783887 -1.112587757
[35,] -3.036833302 3.793783887
[36,] 2.874610334 -3.036833302
[37,] 1.213992589 2.874610334
[38,] -5.134740401 1.213992589
[39,] -8.277657918 -5.134740401
[40,] -0.322821488 -8.277657918
[41,] -0.754066300 -0.322821488
[42,] 3.919384723 -0.754066300
[43,] -4.200621324 3.919384723
[44,] 1.610580424 -4.200621324
[45,] 7.024588715 1.610580424
[46,] 6.272236547 7.024588715
[47,] -0.092118699 6.272236547
[48,] -4.391007689 -0.092118699
[49,] 0.796409135 -4.391007689
[50,] -3.252081842 0.796409135
[51,] -3.258965161 -3.252081842
[52,] -3.429714253 -3.258965161
[53,] -1.466918544 -3.429714253
[54,] -2.938188168 -1.466918544
[55,] -9.317248872 -2.938188168
[56,] -3.055967373 -9.317248872
[57,] -3.614636952 -3.055967373
[58,] 3.320299941 -3.614636952
[59,] -3.378919691 3.320299941
[60,] -2.464193839 -3.378919691
[61,] 2.575143187 -2.464193839
[62,] 0.891052662 2.575143187
[63,] -0.803107321 0.891052662
[64,] -2.602784868 -0.803107321
[65,] 3.214238072 -2.602784868
[66,] 5.812190759 3.214238072
[67,] -3.969830358 5.812190759
[68,] 4.193382072 -3.969830358
[69,] 1.991298084 4.193382072
[70,] 3.958340227 1.991298084
[71,] -2.028463602 3.958340227
[72,] 2.352872954 -2.028463602
[73,] 2.152774486 2.352872954
[74,] -0.796450097 2.152774486
[75,] 1.722942259 -0.796450097
[76,] -4.310911233 1.722942259
[77,] -10.629662935 -4.310911233
[78,] -2.647065901 -10.629662935
[79,] -6.228538865 -2.647065901
[80,] 0.412296504 -6.228538865
[81,] -0.404977842 0.412296504
[82,] -0.785295630 -0.404977842
[83,] 2.773067934 -0.785295630
[84,] -1.517947615 2.773067934
[85,] 3.855147909 -1.517947615
[86,] 0.097048796 3.855147909
[87,] -3.017874592 0.097048796
[88,] -2.035742284 -3.017874592
[89,] -5.919974305 -2.035742284
[90,] 0.025153151 -5.919974305
[91,] -2.012027443 0.025153151
[92,] -0.014637251 -2.012027443
[93,] -1.441917598 -0.014637251
[94,] 3.879492823 -1.441917598
[95,] 4.667866054 3.879492823
[96,] 1.224839981 4.667866054
[97,] 4.826193337 1.224839981
[98,] 6.613040707 4.826193337
[99,] 0.263286367 6.613040707
[100,] 0.637470064 0.263286367
[101,] -3.012567782 0.637470064
[102,] 2.188163491 -3.012567782
[103,] -1.993613616 2.188163491
[104,] -3.363468675 -1.993613616
[105,] 3.461299035 -3.363468675
[106,] -1.719577418 3.461299035
[107,] -2.651517005 -1.719577418
[108,] -2.070439685 -2.651517005
[109,] 4.048629571 -2.070439685
[110,] 2.752677406 4.048629571
[111,] 1.546707617 2.752677406
[112,] 2.533624412 1.546707617
[113,] -1.676948593 2.533624412
[114,] 4.743763725 -1.676948593
[115,] -3.513206668 4.743763725
[116,] 1.823669460 -3.513206668
[117,] 5.345127820 1.823669460
[118,] 0.003030483 5.345127820
[119,] -8.089634858 0.003030483
[120,] -0.848281066 -8.089634858
[121,] 1.264075998 -0.848281066
[122,] 6.315792280 1.264075998
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2.663998492 -0.990321917
2 -2.946816455 2.663998492
3 -2.446177309 -2.946816455
4 2.526076307 -2.446177309
5 1.046513058 2.526076307
6 -5.625858534 1.046513058
7 -7.949504117 -5.625858534
8 2.277068241 -7.949504117
9 6.720882737 2.277068241
10 5.329583336 6.720882737
11 -3.259028071 5.329583336
12 3.436546848 -3.259028071
13 4.965514239 3.436546848
14 -2.105260996 4.965514239
15 3.933212096 -2.105260996
16 2.911550925 3.933212096
17 -4.443419639 2.911550925
18 4.300925932 -4.443419639
19 -3.094725751 4.300925932
20 10.867365047 -3.094725751
21 4.516820091 10.867365047
22 9.389689998 4.516820091
23 1.180371370 9.389689998
24 4.980704231 1.180371370
25 1.979536117 4.980704231
26 0.145010381 1.979536117
27 2.589523045 0.145010381
28 -4.876754819 2.589523045
29 -5.772341610 -4.876754819
30 -3.372591583 -5.772341610
31 -8.607574158 -3.372591583
32 -1.656348838 -8.607574158
33 -1.112587757 -1.656348838
34 3.793783887 -1.112587757
35 -3.036833302 3.793783887
36 2.874610334 -3.036833302
37 1.213992589 2.874610334
38 -5.134740401 1.213992589
39 -8.277657918 -5.134740401
40 -0.322821488 -8.277657918
41 -0.754066300 -0.322821488
42 3.919384723 -0.754066300
43 -4.200621324 3.919384723
44 1.610580424 -4.200621324
45 7.024588715 1.610580424
46 6.272236547 7.024588715
47 -0.092118699 6.272236547
48 -4.391007689 -0.092118699
49 0.796409135 -4.391007689
50 -3.252081842 0.796409135
51 -3.258965161 -3.252081842
52 -3.429714253 -3.258965161
53 -1.466918544 -3.429714253
54 -2.938188168 -1.466918544
55 -9.317248872 -2.938188168
56 -3.055967373 -9.317248872
57 -3.614636952 -3.055967373
58 3.320299941 -3.614636952
59 -3.378919691 3.320299941
60 -2.464193839 -3.378919691
61 2.575143187 -2.464193839
62 0.891052662 2.575143187
63 -0.803107321 0.891052662
64 -2.602784868 -0.803107321
65 3.214238072 -2.602784868
66 5.812190759 3.214238072
67 -3.969830358 5.812190759
68 4.193382072 -3.969830358
69 1.991298084 4.193382072
70 3.958340227 1.991298084
71 -2.028463602 3.958340227
72 2.352872954 -2.028463602
73 2.152774486 2.352872954
74 -0.796450097 2.152774486
75 1.722942259 -0.796450097
76 -4.310911233 1.722942259
77 -10.629662935 -4.310911233
78 -2.647065901 -10.629662935
79 -6.228538865 -2.647065901
80 0.412296504 -6.228538865
81 -0.404977842 0.412296504
82 -0.785295630 -0.404977842
83 2.773067934 -0.785295630
84 -1.517947615 2.773067934
85 3.855147909 -1.517947615
86 0.097048796 3.855147909
87 -3.017874592 0.097048796
88 -2.035742284 -3.017874592
89 -5.919974305 -2.035742284
90 0.025153151 -5.919974305
91 -2.012027443 0.025153151
92 -0.014637251 -2.012027443
93 -1.441917598 -0.014637251
94 3.879492823 -1.441917598
95 4.667866054 3.879492823
96 1.224839981 4.667866054
97 4.826193337 1.224839981
98 6.613040707 4.826193337
99 0.263286367 6.613040707
100 0.637470064 0.263286367
101 -3.012567782 0.637470064
102 2.188163491 -3.012567782
103 -1.993613616 2.188163491
104 -3.363468675 -1.993613616
105 3.461299035 -3.363468675
106 -1.719577418 3.461299035
107 -2.651517005 -1.719577418
108 -2.070439685 -2.651517005
109 4.048629571 -2.070439685
110 2.752677406 4.048629571
111 1.546707617 2.752677406
112 2.533624412 1.546707617
113 -1.676948593 2.533624412
114 4.743763725 -1.676948593
115 -3.513206668 4.743763725
116 1.823669460 -3.513206668
117 5.345127820 1.823669460
118 0.003030483 5.345127820
119 -8.089634858 0.003030483
120 -0.848281066 -8.089634858
121 1.264075998 -0.848281066
122 6.315792280 1.264075998
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/wessaorg/rcomp/tmp/7u1651353238880.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/wessaorg/rcomp/tmp/8ehiq1353238880.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/wessaorg/rcomp/tmp/94zen1353238880.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/wessaorg/rcomp/tmp/10zbm51353238880.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/11x32f1353238880.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/123ext1353238880.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/13o87x1353238880.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/14fx911353238880.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/wessaorg/rcomp/tmp/158dbu1353238880.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/wessaorg/rcomp/tmp/16ulgd1353238880.tab")
+ }
>
> try(system("convert tmp/1ki461353238880.ps tmp/1ki461353238880.png",intern=TRUE))
character(0)
> try(system("convert tmp/2mb9q1353238880.ps tmp/2mb9q1353238880.png",intern=TRUE))
character(0)
> try(system("convert tmp/3fnbq1353238880.ps tmp/3fnbq1353238880.png",intern=TRUE))
character(0)
> try(system("convert tmp/4szlt1353238880.ps tmp/4szlt1353238880.png",intern=TRUE))
character(0)
> try(system("convert tmp/5o1ga1353238880.ps tmp/5o1ga1353238880.png",intern=TRUE))
character(0)
> try(system("convert tmp/6xpsk1353238880.ps tmp/6xpsk1353238880.png",intern=TRUE))
character(0)
> try(system("convert tmp/7u1651353238880.ps tmp/7u1651353238880.png",intern=TRUE))
character(0)
> try(system("convert tmp/8ehiq1353238880.ps tmp/8ehiq1353238880.png",intern=TRUE))
character(0)
> try(system("convert tmp/94zen1353238880.ps tmp/94zen1353238880.png",intern=TRUE))
character(0)
> try(system("convert tmp/10zbm51353238880.ps tmp/10zbm51353238880.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
7.698 1.206 9.109