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('Treatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'Treatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'Treatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'No' + ,'Treatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'No' + ,'Treatment' + ,NA + ,'UsedStats' + ,'No' + ,'Treatment' + ,NA + ,'UsedStats' + ,'Yes' + ,'Treatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'Treatment' + ,NA + ,'UsedStats' + ,'Yes' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'Treatment' + ,NA + ,'UsedStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'Treatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'Treatment' + ,NA + ,'UsedStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'Treatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'Yes' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'Treatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'Treatment' + ,NA + ,'UsedStats' + ,'No' + ,'Treatment' + ,NA + ,'UsedStats' + ,'Yes' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'Yes' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'Treatment' + ,NA + ,'UsedStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'Treatment' + ,NA + ,'UsedStats' + ,'Yes' + ,'Treatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'Treatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'Treatment' + ,NA + ,'UsedStats' + ,'Yes' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'Treatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'No' + ,'Treatment' + ,NA + ,'UsedStats' + ,'Yes' + ,'Treatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'UsedStats' + ,'Yes' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,'NoTreatment' + ,NA + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'Treatment' + ,'UsedStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'Treatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'Treatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'Treatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'Treatment' + ,'UsedStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'Treatment' + ,'UsedStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'Treatment' + ,'UsedStats' + ,'No' + ,NA + ,'Treatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'UsedStats' + ,'No' + ,NA + ,'Treatment' + ,'UsedStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'Treatment' + ,'UsedStats' + ,'No' + ,NA + ,'NoTreatment' + ,'UsedStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'Treatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'UsedStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'UsedStats' + ,'No' + ,NA + ,'Treatment' + ,'UsedStats' + ,'No' + ,NA + ,'Treatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'UsedStats' + ,'Yes' + ,NA + ,'Treatment' + ,'UsedStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'Treatment' + ,'NoStats' + ,'No' + ,NA + ,'Treatment' + ,'UsedStats' + ,'No' + ,NA + ,'Treatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'NoStats' + ,'No' + ,NA + ,'NoTreatment' + ,'UsedStats' + ,'Yes' + ,NA + ,'NoTreatment' + ,'UsedStats' + ,'Yes' + ,NA + ,'NoTreatment' + ,'UsedStats' + ,'No') + ,dim=c(4 + ,154) + ,dimnames=list(c('T40' + ,'T20' + ,'Used' + ,'CorrectAnalysis') + ,1:154)) > y <- array(NA,dim=c(4,154),dimnames=list(c('T40','T20','Used','CorrectAnalysis'),1:154)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } There were 50 or more warnings (use warnings() to see the first 50) > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '4' > par3 <- 'No Linear Trend' > par2 <- 'Do not include Seasonal Dummies' > par1 <- '4' > #'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, 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 CorrectAnalysis T40 T20 Used 1 NA NA NA NA 2 NA NA NA NA 3 NA NA NA NA 4 NA NA NA NA 5 NA NA NA NA 6 NA NA NA NA 7 NA NA NA NA 8 NA NA NA NA 9 NA NA NA NA 10 NA NA NA NA 11 NA NA NA NA 12 NA NA NA NA 13 NA NA NA NA 14 NA NA NA NA 15 NA NA NA NA 16 NA NA NA NA 17 NA NA NA NA 18 NA NA NA NA 19 NA NA NA NA 20 NA NA NA NA 21 NA NA NA NA 22 NA NA NA NA 23 NA NA NA NA 24 NA NA NA NA 25 NA NA NA NA 26 NA NA NA NA 27 NA NA NA NA 28 NA NA NA NA 29 NA NA NA NA 30 NA NA NA NA 31 NA NA NA NA 32 NA NA NA NA 33 NA NA NA NA 34 NA NA NA NA 35 NA NA NA NA 36 NA NA NA NA 37 NA NA NA NA 38 NA NA NA NA 39 NA NA NA NA 40 NA NA NA NA 41 NA NA NA NA 42 NA NA NA NA 43 NA NA NA NA 44 NA NA NA NA 45 NA NA NA NA 46 NA NA NA NA 47 NA NA NA NA 48 NA NA NA NA 49 NA NA NA NA 50 NA NA NA NA 51 NA NA NA NA 52 NA NA NA NA 53 NA NA NA NA 54 NA NA NA NA 55 NA NA NA NA 56 NA NA NA NA 57 NA NA NA NA 58 NA NA NA NA 59 NA NA NA NA 60 NA NA NA NA 61 NA NA NA NA 62 NA NA NA NA 63 NA NA NA NA 64 NA NA NA NA 65 NA NA NA NA 66 NA NA NA NA 67 NA NA NA NA 68 NA NA NA NA 69 NA NA NA NA 70 NA NA NA NA 71 NA NA NA NA 72 NA NA NA NA 73 NA NA NA NA 74 NA NA NA NA 75 NA NA NA NA 76 NA NA NA NA 77 NA NA NA NA 78 NA NA NA NA 79 NA NA NA NA 80 NA NA NA NA 81 NA NA NA NA 82 NA NA NA NA 83 NA NA NA NA 84 NA NA NA NA 85 NA NA NA NA 86 NA NA NA NA 87 NA NA NA NA 88 NA NA NA NA 89 NA NA NA NA 90 NA NA NA NA 91 NA NA NA NA 92 NA NA NA NA 93 NA NA NA NA 94 NA NA NA NA 95 NA NA NA NA 96 NA NA NA NA 97 NA NA NA NA 98 NA NA NA NA 99 NA NA NA NA 100 NA NA NA NA 101 NA NA NA NA 102 NA NA NA NA 103 NA NA NA NA 104 NA NA NA NA 105 NA NA NA NA 106 NA NA NA NA 107 NA NA NA NA 108 NA NA NA NA 109 NA NA NA NA 110 NA NA NA NA 111 NA NA NA NA 112 NA NA NA NA 113 NA NA NA NA 114 NA NA NA NA 115 NA NA NA NA 116 NA NA NA NA 117 NA NA NA NA 118 NA NA NA NA 119 NA NA NA NA 120 NA NA NA NA 121 NA NA NA NA 122 NA NA NA NA 123 NA NA NA NA 124 NA NA NA NA 125 NA NA NA NA 126 NA NA NA NA 127 NA NA NA NA 128 NA NA NA NA 129 NA NA NA NA 130 NA NA NA NA 131 NA NA NA NA 132 NA NA NA NA 133 NA NA NA NA 134 NA NA NA NA 135 NA NA NA NA 136 NA NA NA NA 137 NA NA NA NA 138 NA NA NA NA 139 NA NA NA NA 140 NA NA NA NA 141 NA NA NA NA 142 NA NA NA NA 143 NA NA NA NA 144 NA NA NA NA 145 NA NA NA NA 146 NA NA NA NA 147 NA NA NA NA 148 NA NA NA NA 149 NA NA NA NA 150 NA NA NA NA 151 NA NA NA NA 152 NA NA NA NA 153 NA NA NA NA 154 NA NA NA NA > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : 0 (non-NA) cases Calls: lm -> lm.fit Execution halted