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(1 + ,1 + ,3 + ,3 + ,1 + ,1 + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,4 + ,2 + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,1 + ,1 + ,3 + ,4 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,1 + ,1 + ,4 + ,5 + ,1 + ,0 + ,1 + ,3 + ,5 + ,2 + ,0 + ,1 + ,2 + ,5 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,5 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,4 + ,1 + ,1 + ,1 + ,5 + ,4 + ,4 + ,0 + ,1 + ,2 + ,3 + ,1 + ,1 + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,1 + ,1 + ,1 + ,3 + ,1 + ,0 + ,1 + ,1 + ,5 + ,1 + ,1 + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,4 + ,1 + ,0 + ,1 + ,2 + ,4 + ,1 + ,1 + ,1 + ,2 + ,3 + ,4 + ,0 + ,1 + ,2 + ,4 + ,1 + ,0 + ,1 + ,2 + ,5 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,1 + ,1 + ,2 + ,3 + ,3 + ,1 + ,1 + ,2 + ,5 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,1 + ,1 + ,2 + ,5 + ,1 + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,4 + ,0 + ,1 + ,0 + ,0 + ,2 + ,5 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,1 + ,1 + ,2 + ,5 + ,4 + ,1 + ,0 + ,4 + ,4 + ,3 + ,0 + ,1 + ,2 + ,4 + ,0 + ,NA + ,1 + ,NA + ,NA + ,NA + ,1 + ,1 + ,2 + ,3 + ,1 + ,0 + ,0 + ,3 + ,3 + ,1 + ,1 + ,1 + ,3 + ,5 + ,1 + ,0 + ,1 + ,1 + ,3 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,1 + ,5 + ,3 + ,NA + ,1 + ,NA + ,NA + ,NA + ,1 + ,1 + ,3 + ,4 + ,3 + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,4 + ,3 + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,4 + ,5 + ,3 + ,0 + ,1 + ,3 + ,4 + ,1 + ,0 + ,0 + ,3 + ,3 + ,2 + ,1 + ,0 + ,2 + ,5 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,4 + ,2 + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,1 + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,0 + ,1 + ,0 + ,1 + ,NA + ,NA + ,NA + ,1 + ,1 + ,3 + ,4 + ,1 + ,0 + ,1 + ,NA + ,NA + ,NA + ,1 + ,1 + ,2 + ,2 + ,2 + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,5 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,1 + ,1 + ,3 + ,5 + ,1 + ,1 + ,1 + ,3 + ,4 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,1 + ,1 + ,4 + ,4 + ,2 + ,0 + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,5 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,3 + ,3 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,4 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,5 + ,1 + ,1 + ,1 + ,3 + ,5 + ,1 + ,0 + ,0 + ,3 + ,4 + ,1 + ,0 + ,1 + ,2 + ,4 + ,2 + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,5 + ,2 + ,NA + ,1 + ,NA + ,NA + ,NA + ,1 + ,1 + ,3 + ,5 + ,1 + ,0 + ,1 + ,3 + ,5 + ,1 + ,0 + ,1 + ,2 + ,5 + ,1 + ,0 + ,1 + ,2 + ,4 + ,1 + ,0 + ,1 + ,3 + ,5 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,5 + ,1 + ,1 + ,1 + ,5 + ,4 + ,3 + ,0 + ,0 + ,3 + ,4 + ,2 + ,1 + ,1 + ,NA + ,NA + ,NA + ,1 + ,1 + ,3 + ,4 + ,2 + ,0 + ,1 + ,4 + ,4 + ,3 + ,1 + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,4 + ,5 + ,2 + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,0 + ,4 + ,4 + ,1 + ,1 + ,1 + ,2 + ,4 + ,2 + ,0 + ,0 + ,2 + ,4 + ,1 + ,0 + ,0 + ,2 + ,5 + ,1 + ,0 + ,0 + ,4 + ,4 + ,2 + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,1 + ,5 + ,3 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,1 + ,0 + ,2 + ,4 + ,1 + ,0 + ,0 + ,3 + ,4 + ,1 + ,0 + ,1 + ,2 + ,4 + ,1 + ,0 + ,1 + ,5 + ,3 + ,2 + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,1 + ,1 + ,4 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,1 + ,0 + ,4 + ,2 + ,1 + ,0 + ,1 + ,4 + ,4 + ,3 + ,1 + ,1 + ,3 + ,3 + ,1 + ,0 + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,NA + ,NA + ,NA + ,1 + ,0 + ,4 + ,4 + ,1 + ,1 + ,1 + ,NA + ,NA + ,NA + ,0 + ,0 + ,3 + ,5 + ,2 + ,0 + ,0 + ,4 + ,4 + ,1 + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,4 + ,2 + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,1 + ,1 + ,1 + ,3 + ,1 + ,1 + ,1 + ,2 + ,5 + ,4 + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,4 + ,1 + ,NA + ,1 + ,2 + ,3 + ,1 + ,0 + ,0 + ,1 + ,5 + ,2 + ,1 + ,1 + ,3 + ,4 + ,3 + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,4 + ,4 + ,2 + ,1 + ,1 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,1 + ,5 + ,2 + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,NA + ,NA + ,NA + ,1 + ,0 + ,3 + ,4 + ,1 + ,1 + ,1 + ,3 + ,5 + ,1 + ,1 + ,1 + ,2 + ,5 + ,1 + ,1 + ,1 + ,5 + ,4 + ,2 + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,1 + ,0 + ,3 + ,4 + ,1 + ,0 + ,0 + ,1 + ,3 + ,1 + ,0 + ,1 + ,NA + ,NA + ,NA + ,1 + ,1 + ,3 + ,5 + ,2 + ,0 + ,1 + ,0 + ,5 + ,1 + ,1 + ,1 + ,2 + ,4 + ,3 + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,0 + ,3 + ,5 + ,1 + ,1 + ,1 + ,3 + ,4 + ,1 + ,1 + ,0 + ,2 + ,5 + ,1 + ,0 + ,0 + ,2 + ,3 + ,1 + ,1 + ,0 + ,4 + ,5 + ,4 + ,NA + ,1 + ,NA + ,NA + ,NA + ,0 + ,1 + ,1 + ,5 + ,1 + ,0 + ,1 + ,3 + ,4 + ,1 + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,1 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,1 + ,2 + ,4 + ,2 + ,NA + ,1 + ,NA + ,NA + ,NA + ,1 + ,0 + ,2 + ,4 + ,1 + ,NA + ,0 + ,NA + ,NA + ,NA + ,1 + ,0 + ,3 + ,5 + ,1 + ,0 + ,0 + ,3 + ,4 + ,1 + ,1 + ,0 + ,2 + ,5 + ,1 + ,0 + ,0 + ,4 + ,5 + ,2 + ,0 + ,0 + ,1 + ,4 + ,1 + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,1 + ,0 + ,2 + ,5 + ,3 + ,1 + ,0 + ,2 + ,4 + ,2 + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,0 + ,3 + ,4 + ,1 + ,0 + ,0 + ,1 + ,5 + ,1 + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,1 + ,0 + ,4 + ,3 + ,3 + ,NA + ,0 + ,NA + ,NA + ,NA + ,1 + ,0 + ,2 + ,5 + ,1 + ,1 + ,0 + ,1 + ,4 + ,1 + ,1 + ,0 + ,3 + ,4 + ,1 + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,0 + ,1 + ,5 + ,1 + ,0 + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,0 + ,3 + ,5 + ,1 + ,0 + ,0 + ,3 + ,4 + ,2 + ,1 + ,0 + ,3 + ,4 + ,1 + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,0 + ,NA + ,NA + ,NA + ,1 + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,0 + ,2 + ,5 + ,1 + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,0 + ,1 + ,5 + ,1 + ,1 + ,0 + ,2 + ,4 + ,1 + ,0 + ,0 + ,1 + ,4 + ,1 + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,0 + ,2 + ,0 + ,1 + ,0 + ,0 + ,3 + ,4 + ,1 + ,0 + ,0 + ,3 + ,4 + ,1 + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,0 + ,2 + ,4 + ,1 + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,0 + ,NA + ,NA + ,NA + ,0 + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,1 + ,0 + ,5 + ,4 + ,2 + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,1 + ,0 + ,2 + ,4 + ,1 + ,0 + ,0 + ,3 + ,5 + ,1 + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,0 + ,3 + ,4 + ,1 + ,0 + ,0 + ,5 + ,4 + ,1 + ,NA + ,0 + ,NA + ,NA + ,NA + ,1 + ,0 + ,3 + ,5 + ,1 + ,1 + ,0 + ,1 + ,5 + ,3 + ,1 + ,0 + ,4 + ,4 + ,1 + ,0 + ,0 + ,2 + ,5 + ,1 + ,0 + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,2 + ,5 + ,1 + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,1 + ,0 + ,4 + ,4 + ,2 + ,1 + ,0 + ,3 + ,4 + ,2 + ,1 + ,0 + ,3 + ,4 + ,1 + ,1 + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,1 + ,0 + ,4 + ,5 + ,2 + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,0 + ,3 + ,4 + ,2 + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,NA + ,0 + ,NA + ,NA + ,NA + ,0 + ,0 + ,NA + ,NA + ,NA) + ,dim=c(5 + ,289) + ,dimnames=list(c('Gender' + ,'Pop' + ,'Stress' + ,'Belonging' + ,'CESD') + ,1:289)) > y <- array(NA,dim=c(5,289),dimnames=list(c('Gender','Pop','Stress','Belonging','CESD'),1:289)) > 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 = '1' > #'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 > 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 Gender Pop Stress Belonging CESD 1 1 1 3 3 1 2 1 1 NA NA NA 3 0 1 2 4 2 4 NA 1 NA NA NA 5 NA 1 NA NA NA 6 NA 1 NA NA NA 7 1 1 3 4 1 8 NA 1 NA NA NA 9 NA 0 NA NA NA 10 1 1 4 5 1 11 0 1 3 5 2 12 0 1 2 5 1 13 NA 1 NA NA NA 14 0 1 2 5 1 15 NA 1 NA NA NA 16 NA 1 NA NA NA 17 0 1 2 4 1 18 1 1 5 4 4 19 0 1 2 3 1 20 1 1 NA NA NA 21 NA 1 NA NA NA 22 1 1 1 3 1 23 0 1 1 5 1 24 1 1 NA NA NA 25 NA 1 NA NA NA 26 0 1 2 4 1 27 0 1 2 4 1 28 1 1 2 3 4 29 0 1 2 4 1 30 0 1 2 5 1 31 NA 1 NA NA NA 32 1 1 2 3 3 33 1 1 2 5 1 34 NA 1 NA NA NA 35 NA 1 NA NA NA 36 1 1 2 5 1 37 NA 0 NA NA NA 38 NA 1 NA NA NA 39 0 1 4 0 1 40 0 0 2 5 1 41 NA 1 NA NA NA 42 1 1 2 5 4 43 1 0 4 4 3 44 0 1 2 4 0 45 NA 1 NA NA NA 46 1 1 2 3 1 47 0 0 3 3 1 48 1 1 3 5 1 49 0 1 1 3 1 50 NA 1 NA NA NA 51 0 1 1 5 3 52 NA 1 NA NA NA 53 1 1 3 4 3 54 NA 1 NA NA NA 55 NA 1 NA NA NA 56 0 1 NA NA NA 57 0 1 NA NA NA 58 0 1 2 4 3 59 NA 1 NA NA NA 60 NA 1 NA NA NA 61 0 1 4 5 3 62 0 1 3 4 1 63 0 0 3 3 2 64 1 0 2 5 1 65 NA 1 NA NA NA 66 0 1 2 4 2 67 NA 1 NA NA NA 68 NA 1 NA NA NA 69 1 1 NA NA NA 70 NA 1 NA NA NA 71 NA 1 NA NA NA 72 0 1 2 0 1 73 0 1 NA NA NA 74 1 1 3 4 1 75 0 1 NA NA NA 76 1 1 2 2 2 77 NA 1 NA NA NA 78 NA 1 NA NA NA 79 0 1 2 5 1 80 NA 1 NA NA NA 81 1 1 3 5 1 82 1 1 3 4 1 83 NA 1 NA NA NA 84 1 1 4 4 2 85 0 1 NA NA NA 86 NA 1 NA NA NA 87 0 1 2 5 1 88 NA 1 NA NA NA 89 0 1 3 3 1 90 NA 1 NA NA NA 91 0 1 2 4 1 92 NA 1 NA NA NA 93 0 1 2 5 1 94 1 1 3 5 1 95 0 0 3 4 1 96 0 1 2 4 2 97 NA 1 NA NA NA 98 0 1 2 5 2 99 NA 1 NA NA NA 100 1 1 3 5 1 101 0 1 3 5 1 102 0 1 2 5 1 103 0 1 2 4 1 104 0 1 3 5 1 105 NA 1 NA NA NA 106 0 1 2 5 1 107 1 1 5 4 3 108 0 0 3 4 2 109 1 1 NA NA NA 110 1 1 3 4 2 111 0 1 4 4 3 112 1 1 NA NA NA 113 NA 1 NA NA NA 114 NA 1 NA NA NA 115 0 1 4 5 2 116 NA 1 NA NA NA 117 NA 1 NA NA NA 118 0 0 4 4 1 119 1 1 2 4 2 120 0 0 2 4 1 121 0 0 2 5 1 122 0 0 4 4 2 123 NA 0 NA NA NA 124 0 1 5 3 1 125 NA 1 NA NA NA 126 NA 0 NA NA NA 127 1 0 2 4 1 128 0 0 3 4 1 129 0 1 2 4 1 130 0 1 5 3 2 131 NA 0 NA NA NA 132 0 1 1 4 1 133 NA 1 NA NA NA 134 1 0 4 2 1 135 0 1 4 4 3 136 1 1 3 3 1 137 0 1 NA NA NA 138 0 1 NA NA NA 139 1 0 4 4 1 140 1 1 NA NA NA 141 0 0 3 5 2 142 0 0 4 4 1 143 NA 0 NA NA NA 144 0 1 2 4 2 145 NA 1 NA NA NA 146 NA 0 NA NA NA 147 1 1 1 3 1 148 1 1 2 5 4 149 NA 1 NA NA NA 150 NA 1 NA NA NA 151 0 1 2 4 1 152 NA 1 2 3 1 153 0 0 1 5 2 154 1 1 3 4 3 155 NA 1 NA NA NA 156 0 1 4 4 2 157 1 1 NA NA NA 158 NA 0 NA NA NA 159 NA 1 NA NA NA 160 0 1 1 5 2 161 NA 1 NA NA NA 162 0 1 NA NA NA 163 1 0 3 4 1 164 1 1 3 5 1 165 1 1 2 5 1 166 1 1 5 4 2 167 NA 1 NA NA NA 168 NA 1 NA NA NA 169 1 0 3 4 1 170 0 0 1 3 1 171 0 1 NA NA NA 172 1 1 3 5 2 173 0 1 0 5 1 174 1 1 2 4 3 175 NA 1 NA NA NA 176 NA 0 NA NA NA 177 NA 1 NA NA NA 178 0 0 3 5 1 179 1 1 3 4 1 180 1 0 2 5 1 181 0 0 2 3 1 182 1 0 4 5 4 183 NA 1 NA NA NA 184 0 1 1 5 1 185 0 1 3 4 1 186 NA 1 NA NA NA 187 NA 1 NA NA NA 188 NA 1 NA NA NA 189 NA 0 NA NA NA 190 0 1 2 4 2 191 NA 1 NA NA NA 192 1 0 2 4 1 193 NA 0 NA NA NA 194 1 0 3 5 1 195 0 0 3 4 1 196 1 0 2 5 1 197 0 0 4 5 2 198 0 0 1 4 1 199 NA 0 NA NA NA 200 0 0 NA NA NA 201 NA 0 NA NA NA 202 NA 0 NA NA NA 203 1 0 2 5 3 204 1 0 2 4 2 205 NA 0 NA NA NA 206 NA 0 NA NA NA 207 NA 0 NA NA NA 208 0 0 3 4 1 209 0 0 1 5 1 210 NA 0 NA NA NA 211 NA 0 NA NA NA 212 1 0 4 3 3 213 NA 0 NA NA NA 214 1 0 2 5 1 215 1 0 1 4 1 216 1 0 3 4 1 217 NA 0 NA NA NA 218 0 0 1 5 1 219 0 0 NA NA NA 220 NA 0 NA NA NA 221 0 0 3 5 1 222 0 0 3 4 2 223 1 0 3 4 1 224 NA 0 NA NA NA 225 0 0 NA NA NA 226 1 0 NA NA NA 227 NA 0 NA NA NA 228 0 0 2 5 1 229 NA 0 NA NA NA 230 0 0 1 5 1 231 1 0 2 4 1 232 0 0 1 4 1 233 NA 0 NA NA NA 234 NA 0 NA NA NA 235 0 0 2 0 1 236 0 0 3 4 1 237 0 0 3 4 1 238 NA 0 NA NA NA 239 0 0 2 4 1 240 NA 0 NA NA NA 241 0 0 NA NA NA 242 0 0 NA NA NA 243 NA 0 NA NA NA 244 1 0 5 4 2 245 NA 0 NA NA NA 246 NA 0 NA NA NA 247 1 0 2 4 1 248 0 0 3 5 1 249 NA 0 NA NA NA 250 NA 0 NA NA NA 251 0 0 3 4 1 252 0 0 5 4 1 253 NA 0 NA NA NA 254 1 0 3 5 1 255 1 0 1 5 3 256 1 0 4 4 1 257 0 0 2 5 1 258 0 0 NA NA NA 259 NA 0 NA NA NA 260 NA 0 NA NA NA 261 NA 0 NA NA NA 262 NA 0 NA NA NA 263 NA 0 NA NA NA 264 NA 0 2 5 1 265 NA 0 NA NA NA 266 NA 0 NA NA NA 267 1 0 4 4 2 268 1 0 3 4 2 269 1 0 3 4 1 270 1 0 NA NA NA 271 NA 0 NA NA NA 272 NA 0 NA NA NA 273 NA 0 NA NA NA 274 NA 0 NA NA NA 275 NA 0 NA NA NA 276 NA 0 NA NA NA 277 NA 0 NA NA NA 278 0 0 NA NA NA 279 NA 0 NA NA NA 280 1 0 4 5 2 281 NA 0 NA NA NA 282 NA 0 NA NA NA 283 NA 0 NA NA NA 284 NA 0 NA NA NA 285 0 0 3 4 2 286 NA 0 NA NA NA 287 NA 0 NA NA NA 288 NA 0 NA NA NA 289 0 0 NA NA NA > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Pop Stress Belonging CESD 0.05153 -0.04993 0.06582 0.01204 0.12200 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -0.6911 -0.4191 -0.2913 0.5391 0.7745 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.05153 0.23321 0.221 0.8255 Pop -0.04993 0.08248 -0.605 0.5459 Stress 0.06582 0.04044 1.627 0.1059 Belonging 0.01204 0.04409 0.273 0.7852 CESD 0.12200 0.05159 2.365 0.0194 * --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 0.4847 on 140 degrees of freedom (144 observations deleted due to missingness) Multiple R-squared: 0.07306, Adjusted R-squared: 0.04658 F-statistic: 2.759 on 4 and 140 DF, p-value: 0.03019 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } Error in if (gqarr[mypoint - kp3 + 1, 2] < 0.01) numsignificant1 <- numsignificant1 + : missing value where TRUE/FALSE needed In addition: Warning messages: 1: In pf(q, df1, df2, lower.tail, log.p) : NaNs produced 2: In pf(q, df1, df2, lower.tail, log.p) : NaNs produced 3: In pf(q, df1, df2, lower.tail, log.p) : NaNs produced 4: In pf(q, df1, df2, lower.tail, log.p) : NaNs produced Execution halted