similar to: Use of adj in axis()

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'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'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 want to use question numbers from my survey data (fake data below) as markers in a scatterplot. I'm using "as.character" to convert question numbers to characters. However, plot truncates the 0 (zero) off of question 10. How can I get the ending zero so I can add questions 11 - 20? # Sample code below Question <- c(1,2,3,4,5,6,7,8,9,10) Performance <-

Dear R friends My goal is to eliminate this specific group(1) if the # of NAs in this group greater than 50%(specifically say greater than 3). Would you please show me how to do it. I have a sample data as following: Thanks a lot. Kevin Lin y group f1 f2 f3 30 NA 1 1 1 1 27 NA 1 1 2 2 48 NA 1 2 1 2 40 -0.6066416 1 2 2 1 24

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,

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'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

I investigated why I was getting this error message... Error in Ops.factor(groups[subscripts], vals[1]) : Level sets of factors are different which led me to putting a browser in panel.superpose: Called from: panel.superpose(x, y, subscripts, ...) Browse[1]> vals vals [1] 1 2 3 Levels: 1 2 3 Browse[1]> groups[subscripts] groups[subscripts] [1] 1 2 3 1 2 1 2 3 1 2 3 1 2 3 Levels: 1 2

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 >

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

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

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

Dear All, I would like to perform linear regressions based on Y and all of the combinations of the five predictors, i.e.,(y,x1,x2),(y,x1,x3),....,(y,x1,x2,x4,x5),....,(y,x1,x2,x3,x4,x5). Is there any quick way to do it instead of repeat performing regressions for 31 times? Or, is there any method to manipulate the dataset into the 31 combinations? Thanks for your help!

I've had rpart running on a problem now for a couple of *days*, but I'd expect a decision tree builder to run in minutes if not seconds. Why is rpart slow? Is there anything I can do to make it quicker?

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.

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

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