similar to: Automatic splitting/combining nested categorical variable in glm

Dear list, I get some strange results with daply from the plyr package. In the example below, the average age per municipality for employed en unemployed is calculated. If I do this using tapply (see code below) I get the following result: no yes A NA 36.94931 B 51.22505 34.24887 C 48.05759 51.00198 If I do this using daply: municipality no yes

Hi, i am trying to learn something about multilevel analysis using a great "Discovering statistics using R". I constructed some sample data and then tried to fit a model. Generally model fits well, however when trying to fit the same model using z-score (standarizded) variables i got an error: Error in lme.formula(z_wyn ~ z_IQ + Kasa, data = la, random = ~z_IQ | : nlminb problem,

I'm trying to create maps of reptile abundance in different states & counties using data from Herp.net, which provides lists of specimens with the places that they were found. First I would like to parse the list by state using 2-letter abbreviations, since I'm focusing on certain regions. To do this, I've been trying to create a vector (state2) that gives all state names as

cheers , I need help in shapefile manipulations. I have two shapefiles of my country Cameroon. The first contain 10 provinces and each province contains a certain number of administratives units. I dont have ESRI Arc view. I want to add a admistrative unit of the center province to the map with 10 province. Is it possible to to take a peace of shapefile and combining it to another ? Thanks !

I am having trouble getting data labels to display over the provinces in a GADM map of Canada. Specifically, I need the variable "Number" from the data set "by_province", grouped by "region", to appear on the corresponding regions of the map. The data set "by_province" looks like this: long lat order hole piece region

Dear all, I am trying to find the solution for the optimization problem focused on the finding minimum cost. I used the solution proposed by excel solver, but there is a restriction in the number of variables. My data consists of 300 rows represent cities and 6 columns represent the centres. It constitutes a cost matrix, where the cost are distances between each city and each of six centres. ..+

Hi is there a way in R to create a new column vector with the numbers of 2 others columns simply appended? Example : I have a column with provinces codes (1 to 19) I have a second column with districts codes (1 to ##, depending on the province) I want to create a third column with a code that would be unique for each district: Province  District   District_Unique 1             1            11

Hello! I have 2 maps - both created in ssplot and both identical in terms of outline. Is there any way to superimpose Map1 (which has black borders between Canadian provinces) onto Map2 (which is also a map of Canada)? Thanks a lot for your hints! Dimitri ### A. Reading in Canada data at the province and then at the county level: library(raster) getData('ISO3') # Canada's code is

Hi, I've got a question regarding the restructering of a data set. What I have are municipality zip-codes and the names of 5'000 built-up areas within municipalities. The following example shows, what I would like to do: Input (Zip-Codes and Names): # CODE NAME #1 3 aaa #2 3 aab #3 3 aac #4 4 bba #5 4 bbb #6 4

Hi, I''m trying to get rails running on Debian stable with Apache 2, mod_fcgid, Rails 1.0. The application works fine with webrick (*). Let''s start here: @eugene [/var/www/ls2/linuxsi/public] $ ./dispatch.fcgi Status: 500 Internal Server Error Seems to result in the following in the log files: [12/Feb/2006:12:18:40 :: 4105] starting [12/Feb/2006:12:18:40 :: 4105] terminated

Hello helpful R folks, First off, please forgive my English. Second, I'm new with R, I've searched the archives about subsets, and I haven't found quite the help I need. I'm currently analysing a population survey whose data set has about 15000 households (the rows/observations) and 130 variables (the columns). I've managed to import the set into R as a data.frame called

I''m doing a registration form. I saw many registrations forms which asks you for a country and it lists all the countries of the world, and when you select a country the form is actualized with all the states or provinces, and when you select a state/province then the form is actualized with all the cities. I''m looking to do something like that, is there a free service or pay

Dear R-listers, I would like to carry out a very basic descriptive analysis of homicides rates in Italy, taking into account both the spatial dimension (103 provinces) and the temporal dimension (10 years), but no covariates. In practice, what I would like to do is to describe spatio-temporal variation of homicide rates, identifying those combinations of province-year where the homicide rate

Dear R'ers, I am wondering what is the smallest geographicterritorial unit available for formatting in Canada. Provinces? I know that in the US it is the county so that I can color US counties any way I want, for example: ### Example for coloring US counties ### Creating an ARTIFICIAL criterion for coloring US counties: library(maps) allcounties<-data.frame(county=map('county',

I'm using ssplot for drawing a map of Austria and colour the nine provinces regarding their share of employment. Now I wanted to add the figures in each province and failed miserably. Using the locator() and text() function caused the error message "invalid graphics state". I try to show you what I have done below, maybe you can find a general fault in my codes. I know that it's

I have an address model with country_id and province_id fields There is also a full_address method that returns an address that is in a format that the google maps api will be able to return a long-lat coords. Within the full_address method there is a call to obtain the province/state and country name. --------- def full_address ... full_address = [city, self.province.name,

Dear UseRs, I'm using the library maps. I'm drawing maps of Italy. The map available for this library were prepared around *1989*: "This italy database comes from the NUTS III (Tertiary Administrative Units of the European Community) database of the United Nations Environment Programme (UNEP) GRID-Geneva data sets. These were prepared around *1989*" [cited:

Hello: As someone making the move from STATA to R, I'm finding it difficult at times to perform basic tasks in R, so forgive me if I've missed an obvious and easily obtained solution to my problem. I've searched the help guides and the archives and have not been able to find a solution that works. I have a data frame with thousands of observations that looks something like this:

Hi, I have upgraded the openssh (OpenSSH_4.6p1) on Solaris 8 server. After upgrading I am getting the below error messages when we execute the commands in remote systems using ssh. However the ssh works fine without any issue. Please let me know is there is any fix for it. Apr 2 03:05:39 hostname sshd[7294]: [ID 800047 auth.info] Accepted publickey for user from xx.xx.xx.xx port 32808

I have the winbind login working on FC5 but now logins to local accounts cannot authenticate. My config files are here: http://www.pigeonnier.org/nsswitch.conf http://www.pigeonnier.org/pam.d/ http://www.pigeonnier.org/krb.conf Again, if I try to ssh in as a user that exists only as a local account on the remote host, I am rejected. User msh is -not- a AD account and only exists on the FC5