similar to: Overimposing one map in ssplot onto another

I'm working on some functions for generalized canonical discriminant analysis in conjunction with the heplots package. I've written a candisc.mlm function that takes an mlm object and computes a candisc object containing canonical scores, coeficients, etc. But I'm stumped on how to construct a mlm for the canonical scores, in a function using the *same* right-hand-side of the model

Hello, I have a file with two dependent variables (three and five) and one independent variable. I do i.mod <- lm(cbind(three, five) ~ species, data=i.txt) and get the following output: Coefficients: three five (Intercept) 9.949 9.586 species -1.166 -1.156 I do a" i.can<-candisc(i.mod,data=i): and get the following output: Canonical Discriminant Analysis

Hello, I have a question regarding the candisc package. My data are: species three five 1 2.95 6.63 1 2.53 7.79 1 3.57 5.65 1 3.16 5.47 2 2.58 4.46 2 2.16 6.22 2 3.27 3.52 I put these in a table and then a linear model >newdata <- lm(cbind(three, five) ~ species, data=rawdata) and then do a candisc on them >candata<-candisc(newdata)

I am trying to create a frequency distribution and I am a bit confused. Here are the commands I have entered: > data <- read.csv(file="40609_sortedfinal.csv",head=TRUE,sep=",") > NumberOfActionsByStatus = data$STATUS > NumberOfActionsByUser = data$ETS_LOGIN > NumberOfBidOffer = data$BID_OFFER > NumberOfActionsByUser.freq = table(NumberOfActionsByUser) >

Hi all, I have made a plot with ssplot, using a SpatialPointsDataFrame. The content is quite simple, as I have 9 grid points with lon/lat coordinates and 9 values attached to these coordinates. They are in a square area of 3 by 3 gridboxes. I would like to lay a map from maps() over these values, but when I try this, the grids of the maps (smaller) do not overlap with the grids I have

Dear Rers, is there a way to color counties on a full US map based on a criterion one establishes (i.e., all counties I assign the same number should be the same color)? I explored a bit and looks like the package "maps" might be of help. library(maps) One could get a map of the US: map('usa') One could get countries within a US state: map('county', 'iowa', fill

Hi all, I''ve got the latest win32-mmap code checked into CVS. Unfortunately, it seems that I''m not able to open an existing mapping and retrieve set data. Below is a simple example that seems like it ought to work but doesn''t. Any ideas? # map1.rb require ''win32/mmap'' include Win32 mmap = MMap.new(:name => ''alpha'', :size

I'm having difficulty overlaying maps when writing to a file graphics device. My command sequence has the structure plot(map1) par(new = T) plot(map2) On the screen device, it works fine. When I attempt something like png(file = "map.png") plot(map1) par(new = T) plot(map2) dev.off() only the last map appears, the previous ones having been cleared. Can someone clarify? 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 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

Dear All, Please consider the snippet at the end of the email. I often download some maps (in the R format) from http://www.gadm.org/ However, when I run (typically more than once) a variation of the script below (based on http://bit.ly/1b3W0Aa ), I often get Error in load(url(paste("http://gadm.org/data/rda/", fileName, "_adm", : cannot open the connection In

Dear All, Please see the short script at the end of the email, which I assembled looking for bits and pieces on the web. It essentially does what I need: it plots several countries as a color-coded map. I just would like to fine-tune a bit the final image, in particular (1) Select my own colors for "painting" the countries (i.e. associate manually a color to every level) (2) Be

Hello all. I am trying to run PAN, multilevel multiple imputation program, in R to impute missing data in a longitudinal dataset. I could successfully run the multiple imputation when I only imputed one variable. However, when I tried to impute a time-varying covariate as well as a response variable, I received an error message, “Error: subscript out of bounds.” Can anyone tell if my commands

Dear All, I need to do something relatively simple: generate a map of Europe and paint the various states with different colors (only 3-4 are needed) according to a rule. I would like to keep it as simple as possible and the script resorting to the sp package that I found at http://bit.ly/Oc71ub is exactly what I am looking for (I need to repeat the exercise with Europe instead of

Dear all, in GADM map there are three levels (nation, province and precinct) for each country of the world but for all of them you are never able to plot only one part of a chosen country. To be sure, I am trying to plot only one region of “Italy” and colour the different precincts in it. So far I am able to colour only the part of my interest but the programme still plot the whole country. Is

Dear alltogether, I tried pan() to impute NAs for longitudinal data. The terminology in the following output follows the pan manpage. No data are attached to this script as this may be too huge. y = 15 responses pred = at first just intercept was tried (later on covariates should follow) subj = 168 different subjects with 4 to 6 observations for each subject at time points t1, t2, ..., t6 #

Hi all, I realized the current implmentation has a problem - you can only get the last value set? I realized, after looking at the old C code, that it actually stores values in a hash and marshals the hash, not the values themselves. That seemed clunky to me, though. I thought it would be more interesting if we just marshalled the entire mmap object and passed that back and forth.

Hi All, I am unable to get the plot geomap in googleVis package. data is as follows > head(index.ret) country ytd 1 Argentina -10.18 2 Australia -3.42 3 Austria -2.70 4 Belgium 1.94 5 Brazil -7.16 6 Canada 0.56 > map1 = gvisGeoMap(index.ret,locationvar = 'country', numvar = 'ytd') > plot(map1) But it just displays a blank page, showing an

I have csv data that extend beyond the area I want for an existing map. I want using the boundaries of the polygon shape file as a cookie cutter so that I can overlay the csv data on map without including anything outside the map boundaries and create a dbf file or shapefile of the clipped data .  The reproducible example: ############################################### library(RColorBrewer)