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