similar to: lmer function in R

Hi, I recently found the MNP package. Out of curiosity, I tried to reproduce results from Greene (Econometric Analysis, fourth edition) on page 874. The signs of the estimates are all opposite those of Greene's table. Might anyone be able to tell me what I am doing wrong? I have attached the function call, the coefficients, and a few rows of the data. The dataset has 210 observations.

Thanks John and Chris for your commentary. I''m constantly amazed at what is out there -- projects I didn''t know existed. I think Chris and I are looking for something similiar -- fault tolerant/high-availability over slow/fast links with ease of administration and (ideally) zero downtime. I''d like to virtualize my storage along with my Xen machines. After looking at

Hello all. I noticed that the default setting for breaks in the construction of histograms in hist() is ?right = TRUE?. I think ?right=FALSE? would be more consistent with usual definitions of lower and upper limits for bins in applied statistics, and I suggest that you consider making it the default for hist(). For example, I generated the following frequency distribution for duration of

Dear R users, I have a city map in shape file (polygon). I also have some points that I hope to plot them to the city map. The only information I have about those points are their relative longitude and latitude to the city center by miles. Is there a way that R can help me to do this? Thanks. Gary [[alternative HTML version deleted]]

Hi all, Could you please help me? I am trying to understand why this line works: lm1x = lm(y~X-1, tmp) Here it seems that I was combining the design matrix and the data frame... And X below is not a single column, in fact, it's a bunch of columns in matrix form... I don't understand why this line works... Is it just luck, i.e. if we change the data-set and/or formulas to something

Hi, I am having difficuly deleting and renaming files with users from a foreign domain using acls. My setup is as follows : I have two Samba (3.0.14a)/LDAP domains connected via a VPN (OpenVPN) with a bi-directional trust relationship established. The trust relationship appears to be working correctly. I can log on onto PC's at either end on either Domain :) and the browse lists of both

Dear R community, I'm running "a mlogit" function with a very large data set. Since the data size is so large, I got the Error: Error: cannot allocate vector of size 5.4 Mb In addition: Warning messages: 1: In mapply("*", X, P, SIMPLIFY = FALSE) : Reached total allocation of 16340Mb: see help(memory.size) 2: In mapply("*", X, P, SIMPLIFY = FALSE) : Reached

Dear R users, I have created a Loess surface in R, in which x is relative longitude by miles, y is relative latitude by miles, and z is population density at the neighborhood level. The purpose is to identify some population centers in the region. I'm wondering if there is a way to determine the coordinates (x,y) of each center, so I can know exactly where they are. Let me use the

Dear all, I was wondering whether anyone knows about a program which converts (part of) a data definition file (such as a .dct file for stata, or .sas file for sas) into an R-script. Here is an example with a .sas file: <-- snipp VALUE HTYPE (default=32) 0 = "Not in household" 1 = "Married couple family household" 2 =

Dear all, sorry, this is for sure really basic, but I searched a lot in the internet, and just couldn't find a solution. The problem is to create new variables from a data frame which contains both individual and group variables, such as mean age for an household. My data frame: df hhid h.age 1 10010020 23 2 10010020 23 3 10010126 42 4 10010126 60 5 10010142

Dear R users, I am trying to convert a dataset from one format to several rectangular datasets. A consultant helped design the data entry program for our survey using Delphi/Pascal and for each household the information is stored in a file called "EA-HM-HH.TXT" where EA is the enumeration area number, HM is the homestead number and HH is the household number. Within this file the

Hi R-users!! I have some trouble with the survey pakage and i would be very glad if you can give me an advice. I have a sample from a survey where household were interviewed. The sample has 4 criteria on which the stratification was based: REGION, SIZE OF HOUSEHOLD, SIZE OF LOCALITY, AGE OF HEAD OF HOUSEHOLD. Since i don't have the whole information in each cell of the cross

Dear All, I'm analysing a negative binomial dataset from a population-based study. Many covariates were determined on household level, so all members of a household have the same value for those covariates. In STATA, there seems to be an option for 'clustered analysis' for neg-bin regression. Does an equivalent exist for R(MASS)'s glm.nb or a comparable function? Many thanks for

Hello, I need some help in creating a new variable. I need to create a 'couple identifier', which gives a unique code for every couple/triple/... in a household. So, I can identify couples. To do this, I should use 4 variables: * SERIAL = a unique numeric code for each household * PERNUM = a unique numeric code for each person * SPLOC = the numeric code of the spouse in the

I have a people table with four types of people: clients, spouses, children, and others all setup using single-table inheritance with a foreign key back to a household record. A Household has_one client and spouse, and has_many children and others. I want to use a single "Household.find(@session [:household_id], :include => [:client, :spouse, :children, :others])"

Hi all, I have a dataset of individuals where the variable ID corresponds to the identification of the household where the individual lives. rel.head stands for the relationship with the household head. so rel.head=1 is the household head, rel.head=2 is the spouse, rel.head=3 is the children. Here is an example to see how it looks like: df<-data.frame(ID=c("17100",

Dear R-users, Recently I had to analyze a dataset from household survey. The sample design ensured, that each household in the population has the same probability of being sampled. However the data were gathered from only one adult individual in each household, who was randomly choosen by an interviewer (via "Kish grid"). To equalize the probabilities for each INDIVIDUAL a casewise

I've written a function that takes some input data output from a simulation model and creates some graphs. It's not very complicated code, and it works perfectly fine if I just run the code as is. But I have converted it into a function so we call it externally, and when I try to source the code to test the function, I get the error message "variable names are limited to 256

I would appreciate help in knowing how to repeat categorical variable code given in column=A, by the number in a matching column=B. For example, I have a categorical variable code attributed to a household=A and want to replicate the code for all member of the household, as given in column=B. I would like to have one sequence of categorical variable codes for individuals in column C. I have ~9000

Hi, I am trying to run a logistic regression to look at the risk of malaria infection in individuals. I want to account for intra household correlation and so want to include a household level random effect. I have been using the lmer command in lme4 package but am getting some strange results that are completely different to those I get using STATA. Can I just check that this is the correct