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 <- array(list(14 + ,12 + ,18 + ,11 + ,11 + ,14 + ,12 + ,12 + ,16 + ,21 + ,18 + ,12 + ,14 + ,22 + ,14 + ,11 + ,15 + ,10 + ,15 + ,13 + ,17 + ,10 + ,19 + ,8 + ,10 + ,15 + ,16 + ,14 + ,18 + ,10 + ,14 + ,14 + ,14 + ,14 + ,17 + ,11 + ,14 + ,10 + ,16 + ,13 + ,18 + ,7 + ,11 + ,14 + ,14 + ,12 + ,12 + ,14 + ,17 + ,11 + ,9 + ,9 + ,16 + ,11 + ,14 + ,15 + ,15 + ,14 + ,11 + ,13 + ,16 + ,9 + ,13 + ,15 + ,17 + ,10 + ,15 + ,11 + ,14 + ,13 + ,16 + ,8 + ,9 + ,20 + ,15 + ,12 + ,17 + ,10 + ,13 + ,10 + ,15 + ,9 + ,16 + ,14 + ,16 + ,8 + ,12 + ,14 + ,12 + ,11 + ,11 + ,13 + ,15 + ,9 + ,15 + ,11 + ,17 + ,15 + ,13 + ,11 + ,16 + ,10 + ,14 + ,14 + ,11 + ,18 + ,12 + ,14 + ,12 + ,11 + ,15 + ,12 + ,16 + ,13 + ,15 + ,9 + ,12 + ,10 + ,12 + ,15 + ,8 + ,20 + ,13 + ,12 + ,11 + ,12 + ,14 + ,14 + ,15 + ,13 + ,10 + ,11 + ,11 + ,17 + ,12 + ,12 + ,15 + ,13 + ,15 + ,14 + ,14 + ,13 + ,16 + ,15 + ,15 + ,13 + ,15 + ,10 + ,13 + ,11 + ,12 + ,19 + ,17 + ,13 + ,13 + ,17 + ,15 + ,13 + ,13 + ,9 + ,15 + ,11 + ,16 + ,10 + ,15 + ,9 + ,16 + ,12 + ,15 + ,12 + ,14 + ,13 + ,15 + ,13 + ,14 + ,12 + ,13 + ,15 + ,7 + ,22 + ,17 + ,13 + ,13 + ,15 + ,15 + ,13 + ,14 + ,15 + ,13 + ,10 + ,16 + ,11 + ,12 + ,16 + ,14 + ,11 + ,17 + ,11 + ,15 + ,10 + ,17 + ,10 + ,12 + ,16 + ,16 + ,12 + ,11 + ,11 + ,15 + ,16 + ,9 + ,19 + ,16 + ,11 + ,15 + ,16 + ,10 + ,15 + ,10 + ,24 + ,15 + ,14 + ,11 + ,15 + ,13 + ,11 + ,14 + ,15 + ,18 + ,12 + ,16 + ,10 + ,14 + ,14 + ,14 + ,13 + ,14 + ,9 + ,14 + ,15 + ,12 + ,15 + ,14 + ,14 + ,15 + ,11 + ,15 + ,8 + ,15 + ,11 + ,13 + ,11 + ,17 + ,8 + ,17 + ,10 + ,19 + ,11 + ,15 + ,13 + ,13 + ,11 + ,9 + ,20 + ,15 + ,10 + ,15 + ,15 + ,15 + ,12 + ,16 + ,14 + ,11 + ,23 + ,14 + ,14 + ,11 + ,16 + ,15 + ,11 + ,13 + ,12 + ,15 + ,10 + ,16 + ,14 + ,14 + ,12 + ,15 + ,12 + ,16 + ,11 + ,16 + ,12 + ,11 + ,13 + ,12 + ,11 + ,9 + ,19 + ,16 + ,12 + ,13 + ,17 + ,16 + ,9 + ,12 + ,12 + ,9 + ,19 + ,13 + ,18 + ,13 + ,15 + ,14 + ,14 + ,19 + ,11 + ,13 + ,9 + ,12 + ,18 + ,13 + ,16) + ,dim=c(2 + ,162) + ,dimnames=list(c('Happ' + ,'Depr') + ,1:162)) > y <- array(NA,dim=c(2,162),dimnames=list(c('Happ','Depr'),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' > #'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!) > library(lattice) > library(lmtest) Loading required package: zoo Attaching package: 'zoo' The following object(s) are masked from package:base : 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]) Error in cbind(x[, par1], x[, 1:k != par1]) : subscript out of bounds Execution halted