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(3 + ,13 + ,14 + ,13 + ,3 + ,25 + ,55 + ,147 + ,3 + ,12 + ,8 + ,13 + ,5 + ,158 + ,7 + ,71 + ,4 + ,15 + ,12 + ,16 + ,6 + ,NA + ,NA + ,NA + ,3 + ,12 + ,7 + ,12 + ,6 + ,143 + ,10 + ,NA + ,2 + ,10 + ,10 + ,11 + ,5 + ,67 + ,74 + ,43 + ,3 + ,12 + ,7 + ,12 + ,3 + ,NA + ,NA + ,NA + ,4 + ,15 + ,16 + ,18 + ,8 + ,148 + ,138 + ,8 + ,2 + ,9 + ,11 + ,11 + ,4 + ,28 + ,NA + ,NA + ,3 + ,12 + ,14 + ,14 + ,4 + ,114 + ,113 + ,34 + ,4 + ,11 + ,6 + ,9 + ,4 + ,NA + ,NA + ,NA + ,3 + ,11 + ,16 + ,14 + ,6 + ,123 + ,115 + ,103 + ,3 + ,11 + ,11 + ,12 + ,6 + ,145 + ,9 + ,NA + ,4 + ,15 + ,16 + ,11 + ,5 + ,113 + ,114 + ,73 + ,2 + ,7 + ,12 + ,12 + ,4 + ,152 + ,59 + ,159 + ,2 + ,11 + ,7 + ,13 + ,6 + ,NA + ,NA + ,NA + ,3 + ,11 + ,13 + ,11 + ,4 + ,36 + ,114 + ,113 + ,2 + ,10 + ,11 + ,12 + ,6 + ,NA + ,NA + ,NA + ,4 + ,14 + ,15 + ,16 + ,6 + ,8 + ,102 + ,44 + ,2 + ,10 + ,7 + ,9 + ,4 + ,108 + ,NA + ,NA + ,1 + ,6 + ,9 + ,11 + ,4 + ,112 + ,86 + ,NA + ,3 + ,11 + ,7 + ,13 + ,2 + ,51 + ,17 + ,41 + ,4 + ,15 + ,14 + ,15 + ,7 + ,43 + ,45 + ,74 + ,2 + ,11 + ,15 + ,10 + ,5 + ,120 + ,123 + ,NA + ,3 + ,12 + ,7 + ,11 + ,4 + ,13 + ,24 + ,NA + ,3 + ,14 + ,15 + ,13 + ,6 + ,55 + ,5 + ,NA + ,4 + ,15 + ,17 + ,16 + ,6 + ,103 + ,123 + ,32 + ,3 + ,9 + ,15 + ,15 + ,7 + ,127 + ,136 + ,126 + ,3 + ,13 + ,14 + ,14 + ,5 + ,14 + ,4 + ,154 + ,4 + ,13 + ,14 + ,14 + ,6 + ,135 + ,76 + ,129 + ,4 + ,16 + ,8 + ,14 + ,4 + ,38 + ,99 + ,98 + ,4 + ,13 + ,8 + ,8 + ,4 + ,11 + ,98 + ,82 + ,3 + ,12 + ,14 + ,13 + ,7 + ,43 + ,67 + ,45 + ,4 + ,14 + ,14 + ,15 + ,7 + ,141 + ,92 + ,8 + ,3 + ,11 + ,8 + ,13 + ,4 + ,62 + ,13 + ,NA + ,2 + ,9 + ,11 + ,11 + ,4 + ,62 + ,24 + ,129 + ,4 + ,16 + ,16 + ,15 + ,6 + ,135 + ,129 + ,31 + ,3 + ,12 + ,10 + ,15 + ,6 + ,117 + ,117 + ,117 + ,2 + ,10 + ,8 + ,9 + ,5 + ,82 + ,11 + ,99 + ,4 + ,13 + ,14 + ,13 + ,6 + ,145 + ,20 + ,55 + ,4 + ,16 + ,16 + ,16 + ,7 + ,87 + ,91 + ,132 + ,4 + ,14 + ,13 + ,13 + ,6 + ,76 + ,111 + ,58 + ,4 + ,15 + ,5 + ,11 + ,3 + ,124 + ,NA + ,NA + ,2 + ,5 + ,8 + ,12 + ,3 + ,151 + ,58 + ,NA + ,2 + ,8 + ,10 + ,12 + ,4 + ,131 + ,NA + ,NA + ,3 + ,11 + ,8 + ,12 + ,6 + ,127 + ,146 + ,101 + ,4 + ,16 + ,13 + ,14 + ,7 + ,76 + ,129 + ,31 + ,5 + ,17 + ,15 + ,14 + ,5 + ,25 + ,48 + ,147 + ,2 + ,9 + ,6 + ,8 + ,4 + ,NA + ,NA + ,NA + ,3 + ,9 + ,12 + ,13 + ,5 + ,58 + ,111 + ,132 + ,3 + ,13 + ,16 + ,16 + ,6 + ,115 + ,32 + ,123 + ,2 + ,10 + ,5 + ,13 + ,6 + ,130 + ,112 + ,39 + ,2 + ,6 + ,15 + ,11 + ,6 + ,17 + ,51 + ,136 + ,3 + ,12 + ,12 + ,14 + ,5 + ,102 + ,53 + ,141 + ,2 + ,8 + ,8 + ,13 + ,4 + ,21 + ,131 + ,NA + ,4 + ,14 + ,13 + ,13 + ,5 + ,NA + ,NA + ,NA + ,3 + ,12 + ,14 + ,13 + ,5 + ,14 + ,76 + ,135 + ,3 + ,11 + ,12 + ,12 + ,4 + ,110 + ,106 + ,118 + ,4 + ,16 + ,16 + ,16 + ,6 + ,133 + ,26 + ,154 + ,1 + ,8 + ,10 + ,15 + ,2 + ,83 + ,44 + ,NA + ,4 + ,15 + ,15 + ,15 + ,8 + ,56 + ,63 + ,116 + ,2 + ,7 + ,8 + ,12 + ,3 + ,NA + ,NA + ,NA + ,4 + ,16 + ,16 + ,14 + ,6 + ,44 + ,116 + ,88 + ,4 + ,14 + ,19 + ,12 + ,6 + ,70 + ,119 + ,25 + ,4 + ,16 + ,14 + ,15 + ,6 + ,36 + ,18 + ,113 + ,2 + ,9 + ,6 + ,12 + ,5 + ,5 + ,134 + ,157 + ,4 + ,14 + ,13 + ,13 + ,5 + ,118 + ,138 + ,26 + ,2 + ,11 + ,15 + ,12 + ,6 + ,17 + ,41 + ,38 + ,4 + ,13 + ,7 + ,12 + ,5 + ,79 + ,NA + ,NA + ,4 + ,15 + ,13 + ,13 + ,6 + ,122 + ,57 + ,53 + ,1 + ,5 + ,4 + ,5 + ,2 + ,119 + ,101 + ,NA + ,4 + ,15 + ,14 + ,13 + ,5 + ,36 + ,114 + ,106 + ,4 + ,13 + ,13 + ,13 + ,5 + ,36 + ,113 + ,106 + ,3 + ,11 + ,11 + ,14 + ,5 + ,141 + ,122 + ,102 + ,3 + ,11 + ,14 + ,17 + ,6 + ,NA + ,14 + ,138 + ,3 + ,12 + ,12 + ,13 + ,6 + ,37 + ,10 + ,142 + ,3 + ,12 + ,15 + ,13 + ,6 + ,110 + ,27 + ,73 + ,3 + ,12 + ,14 + ,12 + ,5 + ,10 + ,39 + ,130 + ,3 + ,12 + ,13 + ,13 + ,5 + ,14 + ,133 + ,86 + ,4 + ,14 + ,8 + ,14 + ,4 + ,157 + ,42 + ,78 + ,1 + ,6 + ,6 + ,11 + ,2 + ,59 + ,NA + ,NA + ,3 + ,7 + ,7 + ,12 + ,4 + ,77 + ,58 + ,NA + ,4 + ,14 + ,13 + ,12 + ,6 + ,129 + ,133 + ,4 + ,3 + ,14 + ,13 + ,16 + ,6 + ,125 + ,151 + ,91 + ,3 + ,10 + ,11 + ,12 + ,5 + ,87 + ,111 + ,132 + ,4 + ,13 + ,5 + ,12 + ,3 + ,61 + ,139 + ,NA + ,3 + ,12 + ,12 + ,12 + ,6 + ,146 + ,126 + ,NA + ,2 + ,9 + ,8 + ,10 + ,4 + ,96 + ,139 + ,NA + ,3 + ,12 + ,11 + ,15 + ,5 + ,133 + ,138 + ,14 + ,4 + ,16 + ,14 + ,15 + ,8 + ,47 + ,52 + ,97 + ,3 + ,10 + ,9 + ,12 + ,4 + ,74 + ,67 + ,45 + ,4 + ,14 + ,10 + ,16 + ,6 + ,109 + ,97 + ,NA + ,3 + ,10 + ,13 + ,15 + ,6 + ,30 + ,137 + ,149 + ,4 + ,16 + ,16 + ,16 + ,7 + ,116 + ,56 + ,57 + ,4 + ,15 + ,16 + ,13 + ,6 + ,149 + ,3 + ,105 + ,3 + ,12 + ,11 + ,12 + ,5 + ,19 + ,78 + ,NA + ,2 + ,10 + ,8 + ,11 + ,4 + ,96 + ,NA + ,NA + ,1 + ,8 + ,4 + ,13 + ,6 + ,NA + ,NA + ,NA + ,2 + ,8 + ,7 + ,10 + ,3 + ,21 + ,NA + ,NA + ,2 + ,11 + ,14 + ,15 + ,5 + ,26 + ,118 + ,128 + ,3 + ,13 + ,11 + ,13 + ,6 + ,156 + ,39 + ,29 + ,4 + ,16 + ,17 + ,16 + ,7 + ,53 + ,63 + ,148 + ,4 + ,16 + ,15 + ,15 + ,7 + ,72 + ,78 + ,93 + ,3 + ,14 + ,17 + ,18 + ,6 + ,27 + ,26 + ,4 + ,3 + ,11 + ,5 + ,13 + ,3 + ,66 + ,50 + ,NA + ,1 + ,4 + ,4 + ,10 + ,2 + ,71 + ,104 + ,158 + ,3 + ,14 + ,10 + ,16 + ,8 + ,66 + ,54 + ,144 + ,3 + ,9 + ,11 + ,13 + ,3 + ,40 + ,104 + ,NA + ,4 + ,14 + ,15 + ,15 + ,8 + ,57 + ,148 + ,122 + ,2 + ,8 + ,10 + ,14 + ,3 + ,3 + ,30 + ,149 + ,2 + ,8 + ,9 + ,15 + ,4 + ,12 + ,38 + ,17 + ,3 + ,11 + ,12 + ,14 + ,5 + ,107 + ,132 + ,91 + ,3 + ,12 + ,15 + ,13 + ,7 + ,80 + ,132 + ,111 + ,3 + ,11 + ,7 + ,13 + ,6 + ,98 + ,84 + ,99 + ,4 + ,14 + ,13 + ,15 + ,6 + ,155 + ,71 + ,40 + ,4 + ,15 + ,12 + ,16 + ,7 + ,111 + ,125 + ,132 + ,4 + ,16 + ,14 + ,14 + ,6 + ,81 + ,25 + ,123 + ,4 + ,16 + ,14 + ,14 + ,6 + ,50 + ,66 + ,54 + ,2 + ,11 + ,8 + ,16 + ,6 + ,49 + ,86 + ,90 + ,3 + ,14 + ,15 + ,14 + ,6 + ,96 + ,61 + ,86 + ,4 + ,14 + ,12 + ,12 + ,4 + ,2 + ,60 + ,152 + ,3 + ,12 + ,12 + ,13 + ,4 + ,1 + ,144 + ,152 + ,4 + ,14 + ,16 + ,12 + ,5 + ,22 + ,120 + ,123 + ,2 + ,8 + ,9 + ,12 + ,4 + ,64 + ,139 + ,100 + ,4 + ,13 + ,15 + ,14 + ,6 + ,56 + ,131 + ,116 + ,4 + ,16 + ,15 + ,14 + ,6 + ,144 + ,159 + ,59 + ,3 + ,12 + ,6 + ,14 + ,5 + ,NA + ,NA + ,NA + ,4 + ,16 + ,14 + ,16 + ,8 + ,94 + ,18 + ,5 + ,3 + ,12 + ,15 + ,13 + ,6 + ,25 + ,123 + ,147 + ,3 + ,11 + ,10 + ,14 + ,5 + ,93 + ,18 + ,139 + ,1 + ,4 + ,6 + ,4 + ,4 + ,NA + ,NA + ,NA + ,4 + ,16 + ,14 + ,16 + ,8 + ,48 + ,123 + ,81 + ,4 + ,15 + ,12 + ,13 + ,6 + ,30 + ,105 + ,3 + ,2 + ,10 + ,8 + ,16 + ,4 + ,19 + ,NA + ,NA + ,3 + ,13 + ,11 + ,15 + ,6 + ,NA + ,NA + ,NA + ,4 + ,15 + ,13 + ,14 + ,6 + ,10 + ,68 + ,37 + ,3 + ,12 + ,9 + ,13 + ,4 + ,78 + ,157 + ,5 + ,4 + ,14 + ,15 + ,14 + ,6 + ,93 + ,94 + ,69 + ,2 + ,7 + ,13 + ,12 + ,3 + ,NA + ,NA + ,NA + ,5 + ,19 + ,15 + ,15 + ,6 + ,95 + ,87 + ,NA + ,3 + ,12 + ,14 + ,14 + ,5 + ,50 + ,156 + ,142 + ,4 + ,12 + ,16 + ,13 + ,4 + ,86 + ,139 + ,17 + ,3 + ,13 + ,14 + ,14 + ,6 + ,33 + ,145 + ,100 + ,4 + ,15 + ,14 + ,16 + ,4 + ,152 + ,55 + ,70 + ,2 + ,8 + ,10 + ,6 + ,4 + ,51 + ,41 + ,NA + ,3 + ,12 + ,10 + ,13 + ,4 + ,48 + ,25 + ,123 + ,3 + ,10 + ,4 + ,13 + ,6 + ,97 + ,47 + ,109 + ,2 + ,8 + ,8 + ,14 + ,5 + ,77 + ,NA + ,NA + ,10 + ,15 + ,15 + ,6 + ,130 + ,143 + ,37 + ,4 + ,15 + ,16 + ,14 + ,6 + ,8 + ,102 + ,44 + ,4 + ,16 + ,12 + ,15 + ,8 + ,84 + ,148 + ,98 + ,3 + ,13 + ,12 + ,13 + ,7 + ,51 + ,153 + ,11 + ,4 + ,16 + ,15 + ,16 + ,7 + ,33 + ,32 + ,9 + ,2 + ,9 + ,9 + ,12 + ,4 + ,6 + ,106 + ,NA + ,4 + ,14 + ,12 + ,15 + ,6 + ,116 + ,63 + ,57 + ,4 + ,14 + ,14 + ,12 + ,6 + ,88 + ,56 + ,63 + ,3 + ,12 + ,11 + ,14 + ,2 + ,142 + ,39 + ,66) + ,dim=c(8 + ,156) + ,dimnames=list(c('ManyFriends' + ,'Popularity' + ,'KnowingPeople' + ,'Liked' + ,'Celebrity' + ,'FirstFriend' + ,'SecondFriend' + ,'ThirdFriend') + ,1:156)) > y <- array(NA,dim=c(8,156),dimnames=list(c('ManyFriends','Popularity','KnowingPeople','Liked','Celebrity','FirstFriend','SecondFriend','ThirdFriend'),1:156)) > 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 ManyFriends Popularity KnowingPeople Liked Celebrity FirstFriend 1 3 13 14 13 3 25 2 3 12 8 13 5 158 3 4 15 12 16 6 NA 4 3 12 7 12 6 143 5 2 10 10 11 5 67 6 3 12 7 12 3 NA 7 4 15 16 18 8 148 8 2 9 11 11 4 28 9 3 12 14 14 4 114 10 4 11 6 9 4 NA 11 3 11 16 14 6 123 12 3 11 11 12 6 145 13 4 15 16 11 5 113 14 2 7 12 12 4 152 15 2 11 7 13 6 NA 16 3 11 13 11 4 36 17 2 10 11 12 6 NA 18 4 14 15 16 6 8 19 2 10 7 9 4 108 20 1 6 9 11 4 112 21 3 11 7 13 2 51 22 4 15 14 15 7 43 23 2 11 15 10 5 120 24 3 12 7 11 4 13 25 3 14 15 13 6 55 26 4 15 17 16 6 103 27 3 9 15 15 7 127 28 3 13 14 14 5 14 29 4 13 14 14 6 135 30 4 16 8 14 4 38 31 4 13 8 8 4 11 32 3 12 14 13 7 43 33 4 14 14 15 7 141 34 3 11 8 13 4 62 35 2 9 11 11 4 62 36 4 16 16 15 6 135 37 3 12 10 15 6 117 38 2 10 8 9 5 82 39 4 13 14 13 6 145 40 4 16 16 16 7 87 41 4 14 13 13 6 76 42 4 15 5 11 3 124 43 2 5 8 12 3 151 44 2 8 10 12 4 131 45 3 11 8 12 6 127 46 4 16 13 14 7 76 47 5 17 15 14 5 25 48 2 9 6 8 4 NA 49 3 9 12 13 5 58 50 3 13 16 16 6 115 51 2 10 5 13 6 130 52 2 6 15 11 6 17 53 3 12 12 14 5 102 54 2 8 8 13 4 21 55 4 14 13 13 5 NA 56 3 12 14 13 5 14 57 3 11 12 12 4 110 58 4 16 16 16 6 133 59 1 8 10 15 2 83 60 4 15 15 15 8 56 61 2 7 8 12 3 NA 62 4 16 16 14 6 44 63 4 14 19 12 6 70 64 4 16 14 15 6 36 65 2 9 6 12 5 5 66 4 14 13 13 5 118 67 2 11 15 12 6 17 68 4 13 7 12 5 79 69 4 15 13 13 6 122 70 1 5 4 5 2 119 71 4 15 14 13 5 36 72 4 13 13 13 5 36 73 3 11 11 14 5 141 74 3 11 14 17 6 NA 75 3 12 12 13 6 37 76 3 12 15 13 6 110 77 3 12 14 12 5 10 78 3 12 13 13 5 14 79 4 14 8 14 4 157 80 1 6 6 11 2 59 81 3 7 7 12 4 77 82 4 14 13 12 6 129 83 3 14 13 16 6 125 84 3 10 11 12 5 87 85 4 13 5 12 3 61 86 3 12 12 12 6 146 87 2 9 8 10 4 96 88 3 12 11 15 5 133 89 4 16 14 15 8 47 90 3 10 9 12 4 74 91 4 14 10 16 6 109 92 3 10 13 15 6 30 93 4 16 16 16 7 116 94 4 15 16 13 6 149 95 3 12 11 12 5 19 96 2 10 8 11 4 96 97 1 8 4 13 6 NA 98 2 8 7 10 3 21 99 2 11 14 15 5 26 100 3 13 11 13 6 156 101 4 16 17 16 7 53 102 4 16 15 15 7 72 103 3 14 17 18 6 27 104 3 11 5 13 3 66 105 1 4 4 10 2 71 106 3 14 10 16 8 66 107 3 9 11 13 3 40 108 4 14 15 15 8 57 109 2 8 10 14 3 3 110 2 8 9 15 4 12 111 3 11 12 14 5 107 112 3 12 15 13 7 80 113 3 11 7 13 6 98 114 4 14 13 15 6 155 115 4 15 12 16 7 111 116 4 16 14 14 6 81 117 4 16 14 14 6 50 118 2 11 8 16 6 49 119 3 14 15 14 6 96 120 4 14 12 12 4 2 121 3 12 12 13 4 1 122 4 14 16 12 5 22 123 2 8 9 12 4 64 124 4 13 15 14 6 56 125 4 16 15 14 6 144 126 3 12 6 14 5 NA 127 4 16 14 16 8 94 128 3 12 15 13 6 25 129 3 11 10 14 5 93 130 1 4 6 4 4 NA 131 4 16 14 16 8 48 132 4 15 12 13 6 30 133 2 10 8 16 4 19 134 3 13 11 15 6 NA 135 4 15 13 14 6 10 136 3 12 9 13 4 78 137 4 14 15 14 6 93 138 2 7 13 12 3 NA 139 5 19 15 15 6 95 140 3 12 14 14 5 50 141 4 12 16 13 4 86 142 3 13 14 14 6 33 143 4 15 14 16 4 152 144 2 8 10 6 4 51 145 3 12 10 13 4 48 146 3 10 4 13 6 97 147 2 8 8 14 5 77 148 10 15 15 6 130 143 149 15 16 14 6 8 102 150 16 12 15 8 84 148 151 13 12 13 7 51 153 152 16 15 16 7 33 32 153 9 9 12 4 6 106 154 14 12 15 6 116 63 155 14 14 12 6 88 56 156 12 11 14 2 142 39 SecondFriend ThirdFriend 1 55 147 2 7 71 3 NA NA 4 10 NA 5 74 43 6 NA NA 7 138 8 8 NA NA 9 113 34 10 NA NA 11 115 103 12 9 NA 13 114 73 14 59 159 15 NA NA 16 114 113 17 NA NA 18 102 44 19 NA NA 20 86 NA 21 17 41 22 45 74 23 123 NA 24 24 NA 25 5 NA 26 123 32 27 136 126 28 4 154 29 76 129 30 99 98 31 98 82 32 67 45 33 92 8 34 13 NA 35 24 129 36 129 31 37 117 117 38 11 99 39 20 55 40 91 132 41 111 58 42 NA NA 43 58 NA 44 NA NA 45 146 101 46 129 31 47 48 147 48 NA NA 49 111 132 50 32 123 51 112 39 52 51 136 53 53 141 54 131 NA 55 NA NA 56 76 135 57 106 118 58 26 154 59 44 NA 60 63 116 61 NA NA 62 116 88 63 119 25 64 18 113 65 134 157 66 138 26 67 41 38 68 NA NA 69 57 53 70 101 NA 71 114 106 72 113 106 73 122 102 74 14 138 75 10 142 76 27 73 77 39 130 78 133 86 79 42 78 80 NA NA 81 58 NA 82 133 4 83 151 91 84 111 132 85 139 NA 86 126 NA 87 139 NA 88 138 14 89 52 97 90 67 45 91 97 NA 92 137 149 93 56 57 94 3 105 95 78 NA 96 NA NA 97 NA NA 98 NA NA 99 118 128 100 39 29 101 63 148 102 78 93 103 26 4 104 50 NA 105 104 158 106 54 144 107 104 NA 108 148 122 109 30 149 110 38 17 111 132 91 112 132 111 113 84 99 114 71 40 115 125 132 116 25 123 117 66 54 118 86 90 119 61 86 120 60 152 121 144 152 122 120 123 123 139 100 124 131 116 125 159 59 126 NA NA 127 18 5 128 123 147 129 18 139 130 NA NA 131 123 81 132 105 3 133 NA NA 134 NA NA 135 68 37 136 157 5 137 94 69 138 NA NA 139 87 NA 140 156 142 141 139 17 142 145 100 143 55 70 144 41 NA 145 25 123 146 47 109 147 NA NA 148 37 4 149 44 4 150 98 3 151 11 4 152 9 2 153 NA 4 154 57 4 155 63 3 156 66 3 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Popularity KnowingPeople Liked Celebrity 5.128224 0.345179 0.100532 -0.513301 0.043978 FirstFriend SecondFriend ThirdFriend 0.003250 -0.002802 -0.006297 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -4.78709 -0.59833 -0.09234 0.44689 6.16123 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 5.128224 1.202098 4.266 4.46e-05 *** Popularity 0.345179 0.072418 4.766 6.24e-06 *** KnowingPeople 0.100532 0.057342 1.753 0.0826 . Liked -0.513301 0.081513 -6.297 7.74e-09 *** Celebrity 0.043978 0.009037 4.867 4.15e-06 *** FirstFriend 0.003250 0.003172 1.024 0.3080 SecondFriend -0.002802 0.003212 -0.872 0.3850 ThirdFriend -0.006297 0.003335 -1.888 0.0619 . --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 1.458 on 102 degrees of freedom (46 observations deleted due to missingness) Multiple R-squared: 0.757, Adjusted R-squared: 0.7403 F-statistic: 45.39 on 7 and 102 DF, p-value: < 2.2e-16 > 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 Execution halted