similar to: lmList and missing values

Dear all, I am trying to plot some data with differing range in y-values with type="b", adding error bars and break the y-axis into two parts, one lower part from 12 to 20, and one upper part from 34 to 40. I have tried to follow the basic ideas from the script provided here by Jim Lemon: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/56487.html My attempt looks like this:

------- Forwarded message follows ------- From: Petr Pikal <petr.pikal at precheza.cz> To: Henrik Parn <henrik.parn at bio.ntnu.no> Subject: Re: [R] paired t-test. Need to rearrange data? Date sent: Tue, 08 Aug 2006 16:13:47 +0200 Hi Uff, it takes me a bit headache but this shall do it ?unstack ?table ?t new.list<-unstack(test.data,y~id)

Dear R users, I successfully installed SciViews the other day. However, when I try to run it now, the command/script window does not appear. Strange. Well, actually I see a tendency to a script window (in the lower part of the sciview window where it is suppose to be) during start up, but when the program is entirely open, the script window is gone. I have tried to uninstall and reinstall

Dear all, I have some csv-files (originating from Excel-files) containing empty cells. In my example file I have four variables of different classes, each with some empty cells in the original csv-file: > test <- read.csv2("test.csv", dec=".") > test id id2 x y 1 a 1 NA 2 b e NA 2.2 3 f 3 3.3 4 c g 4 4.4 > class(test$id) [1]

Dear all, I run a linear model with three significant explanatory variabels x1: a factor with 4 levels x2 and x3: factors with two levels each x4: continuous model <- lm(y ~ x1 + x2 * x3 + x4) <> The data is not perfectly balanced between the different factor-combinations and I use treatment contrasts. <> With an aov.object, I assume I could have used model.tables(aov.object,

Dear all, During my pre-R era I tried (yes, tried) to understand mixed models by working through the 'rat example' in Sokal and Rohlfs Biometry (2000) 3ed p 288-292. The same example was later used by Crawley (2002) in his Statistical Computing p 363-373 and I have seen the same data being used elsewhere in the litterature. Because this example is so thoroughly described, I thought

Dear all, I wish to thank Christoph Buser and John Wilkinson for their input, and especially John his examples and for for pointing me to the thread 'Doubt about nested aov output' where the rat-example was hiding...:

Dear R developers, I find the which.min function very handy. However, maybe you could consider a tiny addition to the example data in the the help text. By adding another zero to the vector x... x <- c(1:4, 0, 0:5,11) ...the example would more explicitly show that which.min actually 'giv[es] the index of the /first/ minimum or maximum respectively of x'. This will also more

Dear all, I wish to create a POSIXct variable from date and time variables using the ymd_hm function in package lubridate. In some cases data for time is missing, which causes a problem for ymd_hm. I wish to find a smooth way to handle this. # Some example data: x <- data.frame(date = c("2011-09-22", "2011-07-28"), time = c("15:00", NA)) x # paste date and

Dear all, I have received some data on birds that looks sth like this: # a unique id for each individual id <- c(1,1,1,2,2,2,3,3,3,4,4,5,6) # the year the bird was measured year <- c(1995, 1996, 1997, 1995, 1996, 1997, 1996, 1997, 1998, 1996, 1997, 1997, 1998) # the year the bird was hatched year.hatch <- c(1995, 1995, 1995, 1995, 1995, 1995, 1996, 1996, 1996, 1996, 1996, 1997, 1998)

Dear R-developers, A student asked me today of how to specify the limits of the x-axis. I knew that he should use xlim, but I tried to encourage him to have a try himself with the various help functions. I do not judge if he used the correct search strategy or right key words, but anyway here is what he tried: he looked at ?plot. There xlim is not mentioned. He checked ?par. There you find

Dear all, # 1. I use tapply to calculate various things based on a grouping variable. # simple example using the warpbreaks dataset: tmp <- with(warpbreaks, tapply(breaks, tension, range)) tmp dimnames(tmp) # 2. I wish to plot the result of each element of tmp, using the names of the elements in tmp as 'main' in each plot # My attempt: # par(mfrow = c(3, 1)) # my preferred

Dear R developers, I find the format of descriptions of the arguments in the read.table help text slightly inconsistent. For example, the logical arguments comes in seven different formats, more or less explicit about the consequences of a TRUE (or FALSE): 1. check.names logical. If TRUE then the names... 2. blank.lines.skip logical: if TRUE blank lines... 3. flush logical: if TRUE,

Dear R developers, In the help text for ave() you can read in Arguments section: "... Grouping variables, typically factors, all of the same length as x" ...and in Value: "If ... is g1,g2,". I think a description in the Arguments section of how two or more grouping variables is specified together with an Example would be very helpful. In absence of a help text I just

Dear R developers, ?help.search ...opens a tab in my browser (Firefox 3.5.6, Win XP): http://127.0.0.1:31642/library/utils/html/help.search.html When I click on the link for "?", (http://127.0.0.1:31642/library/utils/help/?) in the last line of "Details", I get an error message: "Server error: invalid response from R". Similarly, The link to "??"

Hi R-helpers. #I start with the reproducible example: firm<-c(rep(1,10),rep(2,10),rep(3,10),rep(4,10),rep(5,10)) year<-c(rep(1998:2007,5)) industry<-c(rep(1,20),rep(5,10),rep(7,10),rep(9,10)) X1<-rnorm(50) X2<-rnorm(50,mean=0.5,sd=0.1) Y<-rnorm(50,mean=0,sd=0.5) data<-data.frame(firm, industry,year,X1,X2,Y) data #I need to calculate for all the industries the following

Attached is a suggested patch for src/main/character.c. It does not fix a bug, but rather provides an improvement on the recent extension of do_grep. In essence, instead of four occurences of 'invert ^ LOGICAL(ind)[i]' that accommodate for the option 'invert' added to grep, there is one occurence of 'LOGICAL(ind)[i] ^= invert' that modifies in-place the logical vector of

given what ?identical says, i find the following odd: x = 1:10 y = 1:10 all.equal(x,y) [1] TRUE identical(x,y) [1] TRUE y[11] = 11 y = y[1:10] all.equal(x,y) [1] TRUE identical(x,y) [1] FALSE y [1] 1 2 3 4 5 6 7 8 9 10 length(y) [1] 10 looks like a bug. platform i686-pc-linux-gnu arch i686 os linux-gnu system

The attached patch provides a modification to the recent fix/improvement to do_grep already included in the most recent development version. The original fix added new functionality to the grep function by adding a new parameter, 'invert'. In the source code for the underlying do_grep, the value of the parameter is used to invert the logical match-no match flag vector ind. The

is there any way i can use R as a tool for scripts in unix, as wotan[hrue]$ cat example.R #!/usr/bin/R x=1 havard -- Havard Rue Department of Mathematical Sciences Norwegian University of Science and Technology N-7491 Trondheim, Norway Voice: +47-7359-3533 URL : http://www.math.ntnu.no/~hrue Fax : +47-7359-3524 Email: havard.rue at math.ntnu.no