R version 2.8.0 (2008-10-20) Copyright (C) 2008 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. Natural language support but running in an English locale 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 + ,27 + ,5 + ,26 + ,49 + ,35 + ,4 + ,1 + ,36 + ,4 + ,25 + ,45 + ,34 + ,5 + ,1 + ,25 + ,4 + ,17 + ,54 + ,13 + ,2 + ,1 + ,27 + ,3 + ,37 + ,36 + ,35 + ,3 + ,2 + ,25 + ,3 + ,35 + ,36 + ,28 + ,5 + ,2 + ,44 + ,3 + ,15 + ,53 + ,32 + ,4 + ,1 + ,50 + ,4 + ,27 + ,46 + ,35 + ,4 + ,1 + ,41 + ,4 + ,36 + ,42 + ,36 + ,4 + ,1 + ,48 + ,5 + ,25 + ,41 + ,27 + ,4 + ,2 + ,43 + ,4 + ,30 + ,45 + ,29 + ,5 + ,2 + ,47 + ,2 + ,27 + ,47 + ,27 + ,4 + ,2 + ,41 + ,3 + ,33 + ,42 + ,28 + ,3 + ,1 + ,44 + ,2 + ,29 + ,45 + ,29 + ,4 + ,2 + ,47 + ,5 + ,30 + ,40 + ,28 + ,3 + ,2 + ,40 + ,3 + ,25 + ,45 + ,30 + ,3 + ,2 + ,46 + ,3 + ,23 + ,40 + ,25 + ,4 + ,1 + ,28 + ,3 + ,26 + ,42 + ,15 + ,3 + ,1 + ,56 + ,3 + ,24 + ,45 + ,33 + ,4 + ,2 + ,49 + ,4 + ,35 + ,47 + ,31 + ,2 + ,2 + ,25 + ,4 + ,39 + ,31 + ,37 + ,4 + ,2 + ,41 + ,4 + ,23 + ,46 + ,37 + ,3 + ,2 + ,26 + ,3 + ,32 + ,34 + ,34 + ,4 + ,1 + ,50 + ,5 + ,29 + ,43 + ,32 + ,4 + ,1 + ,47 + ,4 + ,26 + ,45 + ,21 + ,3 + ,1 + ,52 + ,2 + ,21 + ,42 + ,25 + ,3 + ,2 + ,37 + ,5 + ,35 + ,51 + ,32 + ,2 + ,2 + ,41 + ,3 + ,23 + ,44 + ,28 + ,4 + ,1 + ,45 + ,4 + ,21 + ,47 + ,22 + ,5 + ,2 + ,26 + ,4 + ,28 + ,47 + ,25 + ,4 + ,1 + ,NA + ,3 + ,30 + ,41 + ,26 + ,2 + ,1 + ,52 + ,4 + ,21 + ,44 + ,34 + ,5 + ,1 + ,46 + ,2 + ,29 + ,51 + ,34 + ,4 + ,1 + ,58 + ,3 + ,28 + ,46 + ,36 + ,3 + ,1 + ,54 + ,5 + ,19 + ,47 + ,36 + ,4 + ,1 + ,29 + ,3 + ,26 + ,46 + ,26 + ,2 + ,2 + ,50 + ,3 + ,33 + ,38 + ,26 + ,3 + ,1 + ,43 + ,2 + ,34 + ,50 + ,34 + ,3 + ,2 + ,30 + ,3 + ,33 + ,48 + ,33 + ,3 + ,2 + ,47 + ,2 + ,40 + ,36 + ,31 + ,5 + ,1 + ,45 + ,3 + ,24 + ,51 + ,33 + ,NA + ,2 + ,48 + ,1 + ,35 + ,35 + ,22 + ,4 + ,2 + ,48 + ,3 + ,35 + ,49 + ,29 + ,4 + ,2 + ,26 + ,4 + ,32 + ,38 + ,24 + ,4 + ,1 + ,46 + ,5 + ,20 + ,47 + ,37 + ,2 + ,2 + ,NA + ,3 + ,35 + ,36 + ,32 + ,4 + ,2 + ,50 + ,3 + ,35 + ,47 + ,23 + ,3 + ,1 + ,25 + ,4 + ,21 + ,46 + ,29 + ,4 + ,1 + ,47 + ,2 + ,33 + ,43 + ,35 + ,1 + ,2 + ,47 + ,2 + ,40 + ,53 + ,20 + ,2 + ,1 + ,41 + ,3 + ,22 + ,55 + ,28 + ,2 + ,2 + ,45 + ,2 + ,35 + ,39 + ,26 + ,4 + ,2 + ,41 + ,4 + ,20 + ,55 + ,36 + ,3 + ,2 + ,45 + ,5 + ,28 + ,41 + ,26 + ,4 + ,2 + ,40 + ,3 + ,46 + ,33 + ,33 + ,3 + ,1 + ,29 + ,4 + ,18 + ,52 + ,25 + ,3 + ,2 + ,34 + ,5 + ,22 + ,42 + ,29 + ,5 + ,1 + ,45 + ,5 + ,20 + ,56 + ,32 + ,3 + ,2 + ,52 + ,3 + ,25 + ,46 + ,35 + ,2 + ,2 + ,41 + ,4 + ,31 + ,33 + ,24 + ,1 + ,2 + ,48 + ,3 + ,21 + ,51 + ,31 + ,2 + ,2 + ,45 + ,3 + ,23 + ,46 + ,29 + ,5 + ,1 + ,54 + ,2 + ,26 + ,46 + ,27 + ,4 + ,2 + ,25 + ,3 + ,34 + ,50 + ,29 + ,4 + ,2 + ,26 + ,4 + ,31 + ,46 + ,29 + ,3 + ,1 + ,28 + ,4 + ,23 + ,51 + ,27 + ,4 + ,2 + ,50 + ,4 + ,31 + ,48 + ,34 + ,4 + ,2 + ,48 + ,4 + ,26 + ,44 + ,32 + ,2 + ,2 + ,51 + ,3 + ,36 + ,38 + ,31 + ,3 + ,2 + ,53 + ,3 + ,28 + ,42 + ,31 + ,4 + ,1 + ,37 + ,3 + ,34 + ,39 + ,31 + ,3 + ,1 + ,56 + ,2 + ,25 + ,45 + ,16 + ,2 + ,1 + ,43 + ,3 + ,33 + ,31 + ,25 + ,4 + ,1 + ,34 + ,3 + ,46 + ,29 + ,27 + ,4 + ,1 + ,42 + ,3 + ,24 + ,48 + ,32 + ,3 + ,2 + ,32 + ,3 + ,32 + ,38 + ,28 + ,5 + ,2 + ,31 + ,5 + ,33 + ,55 + ,25 + ,1 + ,1 + ,46 + ,3 + ,42 + ,32 + ,25 + ,3 + ,2 + ,30 + ,5 + ,17 + ,51 + ,36 + ,3 + ,2 + ,47 + ,4 + ,36 + ,53 + ,36 + ,5 + ,2 + ,33 + ,4 + ,40 + ,47 + ,36 + ,2 + ,1 + ,25 + ,4 + ,30 + ,45 + ,27 + ,3 + ,1 + ,25 + ,5 + ,19 + ,33 + ,29 + ,3 + ,2 + ,21 + ,4 + ,33 + ,49 + ,32 + ,4 + ,2 + ,36 + ,5 + ,35 + ,46 + ,29 + ,2 + ,2 + ,50 + ,3 + ,23 + ,42 + ,31 + ,4 + ,2 + ,48 + ,3 + ,15 + ,56 + ,34 + ,3 + ,2 + ,48 + ,2 + ,38 + ,35 + ,27 + ,3 + ,1 + ,25 + ,3 + ,37 + ,40 + ,28 + ,3 + ,1 + ,48 + ,4 + ,23 + ,44 + ,32 + ,2 + ,2 + ,49 + ,5 + ,41 + ,46 + ,33 + ,3 + ,1 + ,27 + ,5 + ,34 + ,46 + ,29 + ,2 + ,1 + ,28 + ,3 + ,38 + ,39 + ,32 + ,4 + ,2 + ,43 + ,2 + ,45 + ,35 + ,35 + ,4 + ,2 + ,48 + ,3 + ,27 + ,48 + ,33 + ,2 + ,2 + ,48 + ,4 + ,46 + ,42 + ,27 + ,1 + ,1 + ,25 + ,1 + ,26 + ,39 + ,16 + ,5 + ,2 + ,49 + ,4 + ,44 + ,39 + ,32 + ,4 + ,1 + ,26 + ,3 + ,36 + ,41 + ,26 + ,4 + ,1 + ,51 + ,3 + ,20 + ,52 + ,32 + ,4 + ,2 + ,25 + ,4 + ,44 + ,45 + ,38 + ,3 + ,1 + ,29 + ,3 + ,27 + ,42 + ,24 + ,3 + ,1 + ,29 + ,4 + ,27 + ,44 + ,26 + ,1 + ,1 + ,43 + ,2 + ,41 + ,33 + ,19 + ,5 + ,2 + ,46 + ,3 + ,30 + ,42 + ,37 + ,3 + ,1 + ,44 + ,3 + ,33 + ,46 + ,25 + ,3 + ,1 + ,25 + ,3 + ,37 + ,45 + ,24 + ,2 + ,1 + ,51 + ,2 + ,30 + ,40 + ,23 + ,4 + ,1 + ,42 + ,5 + ,20 + ,48 + ,28 + ,4 + ,2 + ,53 + ,5 + ,44 + ,32 + ,38 + ,3 + ,1 + ,25 + ,4 + ,20 + ,53 + ,28 + ,4 + ,2 + ,49 + ,2 + ,33 + ,39 + ,28 + ,4 + ,1 + ,51 + ,3 + ,31 + ,45 + ,26 + ,2 + ,2 + ,20 + ,3 + ,23 + ,36 + ,21 + ,3 + ,2 + ,44 + ,3 + ,33 + ,38 + ,35 + ,3 + ,2 + ,38 + ,4 + ,33 + ,49 + ,31 + ,3 + ,1 + ,46 + ,5 + ,32 + ,46 + ,34 + ,4 + ,2 + ,42 + ,4 + ,25 + ,43 + ,30 + ,5 + ,1 + ,29 + ,NA + ,22 + ,37 + ,30 + ,3 + ,2 + ,46 + ,4 + ,16 + ,48 + ,24 + ,3 + ,2 + ,49 + ,2 + ,36 + ,45 + ,27 + ,2 + ,2 + ,51 + ,3 + ,35 + ,32 + ,26 + ,3 + ,1 + ,38 + ,3 + ,25 + ,46 + ,30 + ,1 + ,1 + ,41 + ,1 + ,27 + ,20 + ,15 + ,4 + ,2 + ,47 + ,3 + ,32 + ,42 + ,28 + ,4 + ,2 + ,44 + ,3 + ,36 + ,45 + ,34 + ,4 + ,2 + ,47 + ,3 + ,51 + ,29 + ,29 + ,3 + ,2 + ,46 + ,3 + ,30 + ,51 + ,26 + ,5 + ,1 + ,44 + ,4 + ,20 + ,55 + ,31 + ,2 + ,2 + ,28 + ,3 + ,29 + ,50 + ,28 + ,2 + ,2 + ,47 + ,4 + ,26 + ,44 + ,33 + ,3 + ,2 + ,28 + ,4 + ,20 + ,41 + ,32 + ,3 + ,1 + ,41 + ,5 + ,40 + ,40 + ,33 + ,2 + ,2 + ,45 + ,4 + ,29 + ,47 + ,31 + ,1 + ,2 + ,46 + ,4 + ,32 + ,42 + ,37 + ,3 + ,1 + ,46 + ,4 + ,33 + ,40 + ,27 + ,5 + ,2 + ,22 + ,3 + ,32 + ,51 + ,19 + ,4 + ,2 + ,33 + ,3 + ,34 + ,43 + ,27 + ,4 + ,1 + ,41 + ,4 + ,24 + ,45 + ,31 + ,4 + ,2 + ,47 + ,5 + ,25 + ,41 + ,38 + ,3 + ,1 + ,25 + ,3 + ,41 + ,41 + ,22 + ,5 + ,2 + ,42 + ,3 + ,39 + ,37 + ,35 + ,3 + ,2 + ,47 + ,3 + ,21 + ,46 + ,35 + ,3 + ,2 + ,50 + ,3 + ,38 + ,38 + ,30 + ,3 + ,1 + ,55 + ,5 + ,28 + ,39 + ,41 + ,3 + ,1 + ,21 + ,3 + ,37 + ,45 + ,25 + ,4 + ,1 + ,NA + ,3 + ,26 + ,46 + ,28 + ,2 + ,1 + ,52 + ,3 + ,30 + ,39 + ,45 + ,2 + ,2 + ,49 + ,4 + ,25 + ,21 + ,21 + ,4 + ,2 + ,46 + ,4 + ,38 + ,31 + ,33 + ,3 + ,1 + ,NA + ,4 + ,31 + ,35 + ,25 + ,3 + ,2 + ,45 + ,3 + ,31 + ,49 + ,29 + ,2 + ,2 + ,52 + ,3 + ,27 + ,40 + ,31 + ,3 + ,1 + ,NA + ,3 + ,21 + ,45 + ,29 + ,3 + ,2 + ,40 + ,4 + ,26 + ,46 + ,31 + ,4 + ,2 + ,49 + ,4 + ,37 + ,45 + ,31 + ,1 + ,1 + ,38 + ,5 + ,28 + ,34 + ,25 + ,1 + ,1 + ,32 + ,5 + ,29 + ,41 + ,27 + ,5 + ,2 + ,46 + ,4 + ,33 + ,43 + ,26 + ,4 + ,2 + ,32 + ,3 + ,41 + ,45 + ,26 + ,3 + ,2 + ,41 + ,3 + ,19 + ,48 + ,23 + ,3 + ,2 + ,43 + ,3 + ,37 + ,43 + ,27 + ,4 + ,1 + ,44 + ,4 + ,36 + ,45 + ,24 + ,3 + ,1 + ,47 + ,5 + ,27 + ,45 + ,35 + ,2 + ,2 + ,28 + ,3 + ,33 + ,34 + ,24 + ,1 + ,1 + ,52 + ,1 + ,29 + ,40 + ,32 + ,1 + ,1 + ,27 + ,2 + ,42 + ,40 + ,24 + ,5 + ,2 + ,45 + ,5 + ,27 + ,55 + ,24 + ,4 + ,1 + ,27 + ,4 + ,47 + ,44 + ,38 + ,3 + ,1 + ,25 + ,4 + ,17 + ,44 + ,36 + ,4 + ,1 + ,28 + ,4 + ,34 + ,48 + ,24 + ,5 + ,1 + ,25 + ,3 + ,32 + ,51 + ,18 + ,4 + ,1 + ,52 + ,4 + ,25 + ,49 + ,34 + ,4 + ,1 + ,44 + ,3 + ,27 + ,33 + ,23 + ,2 + ,2 + ,43 + ,3 + ,37 + ,43 + ,35 + ,3 + ,2 + ,47 + ,4 + ,34 + ,44 + ,22 + ,4 + ,2 + ,52 + ,4 + ,27 + ,44 + ,34 + ,3 + ,2 + ,40 + ,2 + ,37 + ,41 + ,28 + ,4 + ,1 + ,42 + ,3 + ,32 + ,45 + ,34 + ,3 + ,1 + ,45 + ,5 + ,26 + ,44 + ,32 + ,4 + ,1 + ,45 + ,2 + ,29 + ,44 + ,24 + ,1 + ,1 + ,50 + ,5 + ,28 + ,40 + ,34 + ,2 + ,1 + ,49 + ,3 + ,19 + ,48 + ,33 + ,3 + ,1 + ,52 + ,2 + ,46 + ,49 + ,33 + ,3 + ,2 + ,48 + ,3 + ,31 + ,46 + ,29 + ,5 + ,2 + ,51 + ,3 + ,42 + ,49 + ,38 + ,4 + ,2 + ,49 + ,4 + ,33 + ,55 + ,24 + ,3 + ,2 + ,31 + ,4 + ,39 + ,51 + ,25 + ,3 + ,2 + ,43 + ,3 + ,27 + ,46 + ,37 + ,3 + ,2 + ,31 + ,3 + ,35 + ,37 + ,33 + ,3 + ,2 + ,28 + ,4 + ,23 + ,43 + ,30 + ,4 + ,2 + ,43 + ,4 + ,32 + ,41 + ,22 + ,3 + ,2 + ,31 + ,3 + ,22 + ,45 + ,28 + ,2 + ,2 + ,51 + ,3 + ,17 + ,39 + ,24 + ,4 + ,2 + ,58 + ,4 + ,35 + ,38 + ,33 + ,2 + ,2 + ,25 + ,5 + ,34 + ,41 + ,37) + ,dim=c(7 + ,195) + ,dimnames=list(c('Teamwork33' + ,'geslacht' + ,'leeftijd' + ,'opleiding' + ,'Neuroticisme' + ,'Extraversie' + ,'Openheid') + ,1:195)) > y <- array(NA,dim=c(7,195),dimnames=list(c('Teamwork33','geslacht','leeftijd','opleiding','Neuroticisme','Extraversie','Openheid'),1:195)) > 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 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]) > 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 Teamwork33 geslacht leeftijd opleiding Neuroticisme Extraversie Openheid 1 4 1 27 5 26 49 35 2 4 1 36 4 25 45 34 3 5 1 25 4 17 54 13 4 2 1 27 3 37 36 35 5 3 2 25 3 35 36 28 6 5 2 44 3 15 53 32 7 4 1 50 4 27 46 35 8 4 1 41 4 36 42 36 9 4 1 48 5 25 41 27 10 4 2 43 4 30 45 29 11 5 2 47 2 27 47 27 12 4 2 41 3 33 42 28 13 3 1 44 2 29 45 29 14 4 2 47 5 30 40 28 15 3 2 40 3 25 45 30 16 3 2 46 3 23 40 25 17 4 1 28 3 26 42 15 18 3 1 56 3 24 45 33 19 4 2 49 4 35 47 31 20 2 2 25 4 39 31 37 21 4 2 41 4 23 46 37 22 3 2 26 3 32 34 34 23 4 1 50 5 29 43 32 24 4 1 47 4 26 45 21 25 3 1 52 2 21 42 25 26 3 2 37 5 35 51 32 27 2 2 41 3 23 44 28 28 4 1 45 4 21 47 22 29 5 2 26 4 28 47 25 30 4 1 NA 3 30 41 26 31 2 1 52 4 21 44 34 32 5 1 46 2 29 51 34 33 4 1 58 3 28 46 36 34 3 1 54 5 19 47 36 35 4 1 29 3 26 46 26 36 2 2 50 3 33 38 26 37 3 1 43 2 34 50 34 38 3 2 30 3 33 48 33 39 3 2 47 2 40 36 31 40 5 1 45 3 24 51 33 41 NA 2 48 1 35 35 22 42 4 2 48 3 35 49 29 43 4 2 26 4 32 38 24 44 4 1 46 5 20 47 37 45 2 2 NA 3 35 36 32 46 4 2 50 3 35 47 23 47 3 1 25 4 21 46 29 48 4 1 47 2 33 43 35 49 1 2 47 2 40 53 20 50 2 1 41 3 22 55 28 51 2 2 45 2 35 39 26 52 4 2 41 4 20 55 36 53 3 2 45 5 28 41 26 54 4 2 40 3 46 33 33 55 3 1 29 4 18 52 25 56 3 2 34 5 22 42 29 57 5 1 45 5 20 56 32 58 3 2 52 3 25 46 35 59 2 2 41 4 31 33 24 60 1 2 48 3 21 51 31 61 2 2 45 3 23 46 29 62 5 1 54 2 26 46 27 63 4 2 25 3 34 50 29 64 4 2 26 4 31 46 29 65 3 1 28 4 23 51 27 66 4 2 50 4 31 48 34 67 4 2 48 4 26 44 32 68 2 2 51 3 36 38 31 69 3 2 53 3 28 42 31 70 4 1 37 3 34 39 31 71 3 1 56 2 25 45 16 72 2 1 43 3 33 31 25 73 4 1 34 3 46 29 27 74 4 1 42 3 24 48 32 75 3 2 32 3 32 38 28 76 5 2 31 5 33 55 25 77 1 1 46 3 42 32 25 78 3 2 30 5 17 51 36 79 3 2 47 4 36 53 36 80 5 2 33 4 40 47 36 81 2 1 25 4 30 45 27 82 3 1 25 5 19 33 29 83 3 2 21 4 33 49 32 84 4 2 36 5 35 46 29 85 2 2 50 3 23 42 31 86 4 2 48 3 15 56 34 87 3 2 48 2 38 35 27 88 3 1 25 3 37 40 28 89 3 1 48 4 23 44 32 90 2 2 49 5 41 46 33 91 3 1 27 5 34 46 29 92 2 1 28 3 38 39 32 93 4 2 43 2 45 35 35 94 4 2 48 3 27 48 33 95 2 2 48 4 46 42 27 96 1 1 25 1 26 39 16 97 5 2 49 4 44 39 32 98 4 1 26 3 36 41 26 99 4 1 51 3 20 52 32 100 4 2 25 4 44 45 38 101 3 1 29 3 27 42 24 102 3 1 29 4 27 44 26 103 1 1 43 2 41 33 19 104 5 2 46 3 30 42 37 105 3 1 44 3 33 46 25 106 3 1 25 3 37 45 24 107 2 1 51 2 30 40 23 108 4 1 42 5 20 48 28 109 4 2 53 5 44 32 38 110 3 1 25 4 20 53 28 111 4 2 49 2 33 39 28 112 4 1 51 3 31 45 26 113 2 2 20 3 23 36 21 114 3 2 44 3 33 38 35 115 3 2 38 4 33 49 31 116 3 1 46 5 32 46 34 117 4 2 42 4 25 43 30 118 5 1 29 NA 22 37 30 119 3 2 46 4 16 48 24 120 3 2 49 2 36 45 27 121 2 2 51 3 35 32 26 122 3 1 38 3 25 46 30 123 1 1 41 1 27 20 15 124 4 2 47 3 32 42 28 125 4 2 44 3 36 45 34 126 4 2 47 3 51 29 29 127 3 2 46 3 30 51 26 128 5 1 44 4 20 55 31 129 2 2 28 3 29 50 28 130 2 2 47 4 26 44 33 131 3 2 28 4 20 41 32 132 3 1 41 5 40 40 33 133 2 2 45 4 29 47 31 134 1 2 46 4 32 42 37 135 3 1 46 4 33 40 27 136 5 2 22 3 32 51 19 137 4 2 33 3 34 43 27 138 4 1 41 4 24 45 31 139 4 2 47 5 25 41 38 140 3 1 25 3 41 41 22 141 5 2 42 3 39 37 35 142 3 2 47 3 21 46 35 143 3 2 50 3 38 38 30 144 3 1 55 5 28 39 41 145 3 1 21 3 37 45 25 146 4 1 NA 3 26 46 28 147 2 1 52 3 30 39 45 148 2 2 49 4 25 21 21 149 4 2 46 4 38 31 33 150 3 1 NA 4 31 35 25 151 3 2 45 3 31 49 29 152 2 2 52 3 27 40 31 153 3 1 NA 3 21 45 29 154 3 2 40 4 26 46 31 155 4 2 49 4 37 45 31 156 1 1 38 5 28 34 25 157 1 1 32 5 29 41 27 158 5 2 46 4 33 43 26 159 4 2 32 3 41 45 26 160 3 2 41 3 19 48 23 161 3 2 43 3 37 43 27 162 4 1 44 4 36 45 24 163 3 1 47 5 27 45 35 164 2 2 28 3 33 34 24 165 1 1 52 1 29 40 32 166 1 1 27 2 42 40 24 167 5 2 45 5 27 55 24 168 4 1 27 4 47 44 38 169 3 1 25 4 17 44 36 170 4 1 28 4 34 48 24 171 5 1 25 3 32 51 18 172 4 1 52 4 25 49 34 173 4 1 44 3 27 33 23 174 2 2 43 3 37 43 35 175 3 2 47 4 34 44 22 176 4 2 52 4 27 44 34 177 3 2 40 2 37 41 28 178 4 1 42 3 32 45 34 179 3 1 45 5 26 44 32 180 4 1 45 2 29 44 24 181 1 1 50 5 28 40 34 182 2 1 49 3 19 48 33 183 3 1 52 2 46 49 33 184 3 2 48 3 31 46 29 185 5 2 51 3 42 49 38 186 4 2 49 4 33 55 24 187 3 2 31 4 39 51 25 188 3 2 43 3 27 46 37 189 3 2 31 3 35 37 33 190 3 2 28 4 23 43 30 191 4 2 43 4 32 41 22 192 3 2 31 3 22 45 28 193 2 2 51 3 17 39 24 194 4 2 58 4 35 38 33 195 2 2 25 5 34 41 37 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) geslacht leeftijd opleiding Neuroticisme -1.088337 0.069710 0.004534 0.138422 0.021597 Extraversie Openheid 0.064594 0.003175 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -2.89188 -0.59446 0.02314 0.70870 1.99984 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -1.088337 0.858587 -1.268 0.2066 geslacht 0.069710 0.146932 0.474 0.6358 leeftijd 0.004534 0.007785 0.582 0.5610 opleiding 0.138422 0.082040 1.687 0.0933 . Neuroticisme 0.021597 0.010790 2.002 0.0468 * Extraversie 0.064594 0.012356 5.228 4.7e-07 *** Openheid 0.003175 0.014610 0.217 0.8282 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 0.9745 on 181 degrees of freedom (7 observations deleted due to missingness) Multiple R-squared: 0.1605, Adjusted R-squared: 0.1326 F-statistic: 5.765 on 6 and 181 DF, p-value: 1.633e-05 > 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