Hi, My code indicates there may be a bug in multilevel. I doubt this is actually the case, can anyone tell me what is wrong with my code?
The data file for this code can be downloaded here: http://cameron.econ.ucdavis.edu/mmabook/mma15p4gev.asc Here is the code that generates the bug: rm(list = ls()) library("multilevel") #v1 fishing.data <- scan("mn.dat") fishing.data.array <- array(fishing.data, c(11, 4, 1182)) fishing.data.array <- aperm(fishing.data.array, c(3,2,1)) means <- aggregate(fishing.data.array[,,3:5], as.data.frame(cbind(fishing.data.array[,1:4,2])), mean) ## Table 1 ## > means ## V1 V2 V3 V4 V1 V2 V3 V4 V5 V6 ## 1 1 0 0 0 35.69949 35.69949 97.80913 125.00316 0.2791948 0.2190015 ## 2 0 1 0 0 30.57133 30.57133 82.42908 109.76335 0.2614444 0.2025348 ## 3 0 0 1 0 137.52710 137.52710 41.60681 70.58409 0.2082868 0.1297646 ## 4 0 0 0 1 120.64835 120.64835 44.56376 75.09694 0.2519077 0.1595341 ## V7 V8 V9 V10 V11 V12 ## 1 0.1593985 0.5176090 4.051617 4.051617 4.051617 4.051617 ## 2 0.1501489 0.4980798 3.387172 3.387172 3.387172 3.387172 ## 3 0.1775411 0.6539167 4.654107 4.654107 4.654107 4.654107 ## 4 0.1771628 0.6914998 3.880900 3.880900 3.880900 3.880900 # ## Note V1 and V2 above are the same. This is suspicious. ## In fact (1,V1) from Table 1 should agree with (1,V1) from Table 2 below: # #v2 fishing.data <- scan("mn.dat") dim(fishing.data) <- c(11,1182*4) fishing.data <- t(fishing.data) means <- aggregate(fishing.data[,3:5], as.data.frame(fishing.data[,2]*fishing.data[,6:8]), mean) ## Table 2 ## > means ## V1 V2 V3 V1 V2 V3 ## 1 0 0 0 95.85924 0.2717639 4.097600 ## 2 1 0 0 30.57133 0.2025348 3.387172 ## 3 0 1 0 41.60681 0.1775411 4.654107 ## 4 0 0 1 75.09694 0.6914998 3.880900 -- View this message in context: http://r.789695.n4.nabble.com/multilevel-tp3539421p3539421.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
