similar to: Level sets of factors are different (panel.superpose)

I've been able to do most things I'd like to using the map package, but I can't figure out how to fill a region without drawing a black outline. I can draw the outline in a specified colour, but as soon as I set fill to TRUE, I get a black outline in addition. My map has a lot of activity on the coastline so I don't want a black line obscuring what I do there. A pale colour to

I'm asking on this list because the question will be more easily understood than on a Perl list. Lots of talented people using R also know how to use Perl, so it's easier to ask them. Running a Perl script with the -d switch can do some pretty neat things, but as far as I can tell, every line has to be done individually. Is there a way to specify where to stop in the way browser() does

platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 5.1 year 2002 month 06 day 17 language R I've been using bwplot to plot a dataframe that has a structure like this: str(plot.df) `data.frame': 5952

Hi, Now suppose I have just one list called FOO, which has 25 objects, e.g.: [[1]] 1 2 3 4 5 [[2]] 6 7 8 9 10 . . . And I want to do something like: FRED <- data.frame(cbind(unlist(FOO[[1]]), unlist(FOO[[2]]), # ... for all 25 subsets )) Is it possible to do this, without doing unlist(FOO[[i]]) 25

I'm baffled. When I run make check after installing from source, I get a Error 2. From my understanding of how these things work, it would appear to be coming from this (as at the end of base-Ex.Rout.fail: > has.VR <- require(MASS, quietly = TRUE) Attaching package 'MASS': The following object(s) are masked from package:base : confint confint.lm nclass.FD nclass.scott

I have a square matrix wherein a '*' indicates an HSD between the levels indicated by row name and column name. The '.' is simply marking the diagonal. A blank indicates the same group A B C D E F G H I J K L M N A . B . C * * . D * * . E * * * . F * * * .

I've been looking at the kinship package which looks as though it might be appropriate for my purposes. What I can't find is any reference to the data that is used in the example code. A dataframe called d10 with column names, upn, dadid, momid, sex and affect is required. One can get an idea of what sort of values should be in most columns from the description in the pedigree function,

I'm sure there's a simple way to do this, but I can only think of complicated ones. I have a number of character vectors that look something like this: "12 78 23 9 76 43 2 15 41 81 92 5(92 12) (81 78 5 76 9 41) (23 2 15 43)" I wish to get it into a list of numerical vectors like this: $Group [1] 12 78 23 9 76 43 2 15 41 81 92 5 $Subgroup1 [1] 92 12 $Subgroup2 [1] 81 78 5

Of course, that's as trivial as it gets on the command line, but I can't work out how to get a column of numbers that are entered as "10^5" from its 'character' format into a numeric one? I feel a bit embarrassed asking such a simple question. Too much Easter.... Thanks -- Patrick Connolly HortResearch Mt Albert Auckland New Zealand Ph: +64-9 815 4200 x 7188

I seem to remember this was discussed a year or two ago, but I can't find it in the archives. platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 7.0 year 2003 month 04 day 16 language R

I've been asked about "Canonical Variance Analysis" (CVA). I don't see any reference to it searching the R site. Does it go by other names? Genstat describes it thus: Canonical variates analysis operates on a within-group sums of squares and products matrix, calculated from a set of variates and factor that specifies the grouping of units. It finds linear combinations of the

At the end of the 'Details:' section of the help on Anova (car package), it states: The standard R `anova' function calculates sequential (type-I) tests. These rarely test meaningful hypotheses. So I thought I'd try it. However, I was perplexed to get this message: Browse[1]> Anova.glm(leaf.glm1, type ="II") Error in terms.formula(formula, data = data) :

Some years ago when I used S-PLUS, I seem to remember, there was a discussion about a simple method of removing trailing spaces from vector elements. I'd be fairly sure the same would work in R, but I can't find any mention of anything like it in the R archives or with help.search(). There are ways I could do it with substring(), but I seem to remember there was something more elegant.

V&R describes binomial GLMs with mortality out of 20 budworms. Is it appropriate to use the same approach with mortality out of numbers as low as 3? I feel reticent to do so with data that is not very continuous. There are one continuous and one categorical independent variables. Would it be more appropriate to treat the response as an ordered factor with four levels? If so, what family

Say we have a tab delimited file called bug.txt Part Rep Cage Hb pupae 1 1 S 32 1 M 34 L 42 2 S 36 M 28 L 36 read.delim("bug.txt") Part Rep Cage Hb.pupae 1 1 1 S 32 2 NA 1 M 34 3 NA NA L 42 4 NA NA NA 5 NA NA NA 6 NA 2 S 36 7 NA NA M 28 8 NA NA L 36 >

I'm still having trouble getting to grips with time classes. I wish to calculate the difference in days between events. Browse[1]> insp.j$First [1] "2002-02-19 13:00:00 NZDT" Browse[1]> spray.j$Date [1] "2001-11-29 13:00:00 NZDT" Browse[1]> insp.jk - spray.j$Date Time difference of 82 days If I save insp.jk to a vector, I get a nice useful value of 82.

I can use a postscript device to produce a file which I can open with the Gimp, then save quite a good png file. If I try to use the png() function directly, I get something rather like a graphics screen dump. The help file doesn't tell me whether it should be like that or not. _ platform i686-pc-linux-gnu arch i686 os linux-gnu system

I like the fact that the postscript function enables the possbiility of a more useful title than before. However, I'd prefer the default to be the file name. It's very simple for me to make my own postscript function that does just that simply by setting title = file. I always use onefile = TRUE, so it always works (so far). However, I'm a little reluctant to do that in case some

I like the fact that the postscript function enables the possbiility of a more useful title than before. However, I'd prefer the default to be the file name. It's very simple for me to make my own postscript function that does just that simply by setting title = file. I always use onefile = TRUE, so it always works (so far). However, I'm a little reluctant to do that in case some

I've never had this problem before and can't see what could be different from other times I've used keys with lattice. It appears that auto.key is being taken as TRUE when I specify a key list. The list I specify seems to be ignored. Where can I place a browser to figure out what is going on? Having made a list key.symbol from trellis.par.get, and specified a scales list and a