[R] multilevel

From: eeecon <121onto_at_gmail.com>
Date: Fri, 20 May 2011 12:54:22 -0700 (PDT)


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_at_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.
Received on Fri 20 May 2011 - 23:04:59 GMT

This quarter's messages: by month, or sorted: [ by date ] [ by thread ] [ by subject ] [ by author ]

All messages

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sat 21 May 2011 - 18:20:08 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive