R version 2.6.0 (2007-10-03) Copyright (C) 2007 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(1.0137 + ,89.97 + ,0.9834 + ,99.8 + ,0.9643 + ,112.99 + ,0.947 + ,93.69 + ,0.906 + ,108.02 + ,0.9492 + ,99.11 + ,0.9397 + ,94.33 + ,0.9041 + ,83.75 + ,0.8721 + ,106.37 + ,0.8552 + ,109.63 + ,0.8564 + ,105.5 + ,0.8973 + ,96.13 + ,0.9383 + ,102.48 + ,0.9217 + ,101.37 + ,0.9095 + ,112.76 + ,0.892 + ,95.57 + ,0.8742 + ,102.81 + ,0.8532 + ,104.13 + ,0.8607 + ,97.52 + ,0.9005 + ,85.29 + ,0.9111 + ,101.01 + ,0.9059 + ,108.48 + ,0.8883 + ,101.33 + ,0.8924 + ,87.57 + ,0.8833 + ,97.44 + ,0.87 + ,96.06 + ,0.8758 + ,106.67 + ,0.8858 + ,102.67 + ,0.917 + ,104.54 + ,0.9554 + ,102.46 + ,0.9922 + ,103.35 + ,0.9778 + ,83.27 + ,0.9808 + ,108.22 + ,0.9811 + ,115.23 + ,1.0014 + ,103.7 + ,1.0183 + ,93.61 + ,1.0622 + ,100.25 + ,1.0773 + ,100.56 + ,1.0807 + ,108.86 + ,1.0848 + ,105.43 + ,1.1582 + ,104.77 + ,1.1663 + ,109.13 + ,1.1372 + ,106.13 + ,1.1139 + ,82.27 + ,1.1222 + ,113.6 + ,1.1692 + ,117.73 + ,1.1702 + ,104.83 + ,1.2286 + ,104.61 + ,1.2613 + ,102.93 + ,1.2646 + ,106.95 + ,1.2262 + ,123.45 + ,1.1985 + ,111.99 + ,1.2007 + ,103.95 + ,1.2138 + ,122.05 + ,1.2266 + ,108.04 + ,1.2176 + ,93.72 + ,1.2218 + ,119.61 + ,1.249 + ,118.29 + ,1.2991 + ,117.14 + ,1.3408 + ,112.76 + ,1.3119 + ,105.97 + ,1.3014 + ,107.96 + ,1.3201 + ,122.27 + ,1.2938 + ,114.54 + ,1.2694 + ,110.15 + ,1.2165 + ,120.02 + ,1.2037 + ,103.94 + ,1.2292 + ,96.18 + ,1.2256 + ,121.01 + ,1.2015 + ,110.55 + ,1.1786 + ,120.04 + ,1.1856 + ,114.19) + ,dim=c(2 + ,72) + ,dimnames=list(c('wk' + ,'uit') + ,1:72)) > y <- array(NA,dim=c(2,72),dimnames=list(c('wk','uit'),1:72)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'First Differences' > par2 = 'Include Monthly Dummies' > par1 = '2' > #'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) > 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) + } Error in cbind(x, x2) : number of rows of matrices must match (see arg 2) Execution halted