R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-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(41
+ ,38
+ ,39
+ ,32
+ ,30
+ ,35
+ ,31
+ ,33
+ ,34
+ ,37
+ ,35
+ ,29
+ ,39
+ ,31
+ ,34
+ ,36
+ ,36
+ ,35
+ ,37
+ ,38
+ ,38
+ ,31
+ ,36
+ ,34
+ ,38
+ ,35
+ ,39
+ ,38
+ ,33
+ ,37
+ ,32
+ ,33
+ ,36
+ ,32
+ ,38
+ ,38
+ ,39
+ ,38
+ ,32
+ ,32
+ ,32
+ ,33
+ ,31
+ ,31
+ ,39
+ ,38
+ ,37
+ ,39
+ ,39
+ ,32
+ ,41
+ ,32
+ ,36
+ ,35
+ ,33
+ ,37
+ ,33
+ ,33
+ ,34
+ ,33
+ ,31
+ ,28
+ ,27
+ ,32
+ ,37
+ ,31
+ ,34
+ ,37
+ ,34
+ ,30
+ ,32
+ ,33
+ ,29
+ ,31
+ ,36
+ ,33
+ ,29
+ ,31
+ ,35
+ ,33
+ ,37
+ ,32
+ ,34
+ ,33
+ ,38
+ ,32
+ ,35
+ ,33
+ ,38
+ ,28
+ ,37
+ ,35
+ ,38
+ ,39
+ ,33
+ ,34
+ ,36
+ ,38
+ ,38
+ ,32
+ ,32
+ ,38
+ ,32
+ ,30
+ ,32
+ ,33
+ ,34
+ ,38
+ ,32
+ ,32
+ ,37
+ ,32
+ ,39
+ ,34
+ ,29
+ ,34
+ ,37
+ ,36
+ ,35
+ ,34
+ ,30
+ ,28
+ ,38
+ ,34
+ ,34
+ ,35
+ ,31
+ ,35
+ ,34
+ ,31
+ ,35
+ ,37
+ ,36
+ ,35
+ ,30
+ ,27
+ ,39
+ ,40
+ ,35
+ ,37
+ ,38
+ ,36
+ ,31
+ ,38
+ ,34
+ ,39
+ ,38
+ ,41
+ ,34
+ ,27
+ ,39
+ ,30
+ ,37
+ ,37
+ ,34
+ ,31
+ ,28
+ ,31
+ ,37
+ ,27
+ ,33
+ ,36
+ ,37
+ ,38
+ ,35
+ ,37
+ ,37
+ ,33
+ ,32
+ ,34
+ ,33
+ ,31
+ ,38
+ ,39
+ ,33
+ ,34
+ ,29
+ ,32
+ ,33
+ ,33
+ ,31
+ ,36
+ ,36
+ ,32
+ ,35
+ ,41
+ ,32
+ ,28
+ ,29
+ ,30
+ ,39
+ ,36
+ ,37
+ ,35
+ ,35
+ ,31
+ ,37
+ ,34
+ ,32
+ ,36
+ ,38
+ ,36
+ ,37
+ ,35
+ ,36
+ ,37
+ ,32
+ ,28
+ ,33
+ ,39
+ ,40
+ ,32
+ ,38
+ ,35
+ ,41
+ ,39
+ ,36
+ ,35
+ ,43
+ ,42
+ ,30
+ ,34
+ ,31
+ ,33
+ ,32
+ ,41
+ ,32
+ ,33
+ ,37
+ ,34
+ ,37
+ ,32
+ ,33
+ ,40
+ ,34
+ ,40
+ ,33
+ ,35
+ ,38
+ ,36
+ ,33
+ ,37
+ ,31
+ ,27
+ ,38
+ ,39
+ ,37
+ ,38
+ ,33
+ ,31
+ ,31
+ ,33
+ ,39
+ ,32
+ ,44
+ ,39
+ ,33
+ ,36
+ ,35
+ ,33
+ ,32
+ ,33
+ ,28
+ ,32
+ ,40
+ ,37
+ ,27
+ ,30
+ ,37
+ ,38
+ ,32
+ ,29
+ ,28
+ ,22
+ ,34
+ ,35
+ ,30
+ ,35
+ ,35
+ ,34
+ ,31
+ ,35
+ ,32
+ ,34
+ ,30
+ ,34
+ ,30
+ ,35
+ ,31
+ ,23
+ ,40
+ ,31
+ ,32
+ ,27
+ ,36
+ ,36
+ ,32
+ ,31
+ ,35
+ ,32
+ ,38
+ ,39
+ ,42
+ ,37
+ ,34
+ ,38
+ ,35
+ ,39
+ ,35
+ ,34
+ ,33
+ ,31
+ ,36
+ ,32
+ ,32
+ ,37
+ ,33
+ ,36
+ ,34
+ ,32
+ ,32
+ ,35
+ ,34
+ ,36)
+ ,dim=c(2
+ ,162)
+ ,dimnames=list(c('Connected'
+ ,'Separate')
+ ,1:162))
>  y <- array(NA,dim=c(2,162),dimnames=list(c('Connected','Separate'),1:162))
>  for (i in 1:dim(x)[1])
+  {
+  	for (j in 1:dim(x)[2])
+  	{
+  		y[i,j] <- as.numeric(x[i,j])
+  	}
+  }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '3'
> library(lattice)
> library(lmtest)
Loading required package: zoo
> 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])
Error in cbind(x[, par1], x[, 1:k != par1]) : subscript out of bounds
Execution halted