search for: rauberta

...lt;- factor(rep(1:4, each=25)) plot(c(0,4), c(-3,3), type="n", xaxt="n", yaxt="n") boxplot(x ~ f, add=TRUE) My expectation is that a high-level plotting function will not mess with the axes when called with add=TRUE. Is this a bug or is my expectation wrong? Richard Raubertas Biometrics Research, RY84-16 Merck & Co. ------------------------------------------------------------------------------

...nited States.1252 [2] LC_CTYPE=English_United States.1252 [3] LC_MONETARY=English_United States.1252 [4] LC_NUMERIC=C [5] LC_TIME=English_United States.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base Richard Raubertas Merck & Co. Notice: This e-mail message, together with any attachme...{{dropped:11}}

...sd' (based on 'var') has similar documentation. However, I get: var(numeric(0)) [1] NA rather than an error. Personally I prefer that basic summary functions like 'var' not throw errors even in corner cases. But either way, the result and the docs are inconsistent. Richard Raubertas > sessionInfo() R version 3.5.1 (2018-07-02) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows 7 x64 (build 7601) Service Pack 1 Matrix products: default locale: [1] LC_COLLATE=English_United States.1252 [2] LC_CTYPE=English_United States.1252 [3] LC_MONETARY=English_Unite...

Hi, Is it possible to add labels to variables in R (so as to have a better description of what the variables represent)? Thanks, Neil

...39;s y[40] <- 50 # add an outlier plot(lm(y ~ x, na.action=na.omit)) # outlier correctly labeled in all # four plots plot(lm(y ~ x, na.action=na.exclude)) # labels attached to wrong points # in the QQ plot (only) Rich Raubertas Biometrics Research Merck & Co. ------------------------------------------------------------------------------ Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (Whitehouse Station, New Jersey, USA), and/or its affiliates (which may be...

...r 1 minor 4.1 year 2002 month 01 day 30 language R > Following is a diagnosis of the problem by one of my colleagues: > From: Huntsinger, Reid > Sent: Monday, April 22, 2002 5:15 PM > To: Raubertas, Richard > Subject: RE: R problem > > [ ... ] This happens > because of the attempt to convert. The function type.convert > is called on columns with unspecified types; this tries to > convert to "logical, integer, numeric or complex" in that > order (first one...

...ite=FALSE' argument appears to be ignored in file.copy(): > file.exists(c("file1.txt", "file2.txt")) [1] TRUE TRUE > file.copy(from="file2.txt", to="file1.txt", overwrite=FALSE) [1] TRUE > and sure enough, file1.txt has been overwritten. Rich Raubertas Biometrics Research, RY33-300 Merck & Co. ------------------------------------------------------------------------------ Notice: This e-mail message, together with any attachments, ...{{dropped}}

Dear R-helpers, I have found a slightly annoying problem when trying to plot lines on graphs. I first created my data using tapply, thus:- y1=as.vector(fit1$coef$random$id) x1=tapply(o1,id,median,na.rm=T) x2=tapply(o2,id,median,na.rm=T) #then I plot the data, thus:- plot(x1[x2==0],y[x2==0]) #if I now try to fit the linear regression, R 'hangs up'

...cause each panel to adjust its axis to match the set of 'a' values observed in that panel, but unfortunately it seems 'relation' is ignored for factors. Would it make sense to have 'relation' work for factors? Is there some other way to accomplish the same effect? Rich Raubertas Merck & Co.

..."b", temp) [1] "b" "Nb" It seems to me that paste(temp[1], temp[2]) should return a single character NA: it asks to concatenate "a" with some unknown string, so the result should be unknown as well. This is certainly how numeric NA's are handled. Rich Raubertas Biometrics Research, RY33-300 Merck & Co. ------------------------------------------------------------------------------ Notice: This e-mail message, together with any attachments, ...{{dropped}}

...if (!na.rm && any(is.na(x))) NA else sqrt(var(x, na.rm = na.rm)) } else if (is.data.frame(x)) { sapply(x, sd, na.rm = na.rm) } else { x <- as.vector(x) if (!na.rm && any(is.na(x))) NA else sqrt(var(x, na.rm = na.rm)) } } Rich Raubertas Merck & Co. ------------------------------------------------------------------------------Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates (which may...

Dear R-users, Does anybody know how to construct a block-diagonal matrix (with the blocks being different matrixs, concerning the dimension and the values) ? I would like to know also if there is any package that generates values from a multivariate skew normal distribution. Thanks all, Caio --------------------------------- [[alternative HTML version deleted]]

...it is no longer missing) but does not set 'noZeroColL'. In the similar function 'anova.lme' the problem is avoided by the=20 first line Lmiss <- missing(L) and then testing whether 'Lmiss' is TRUE in the rest of the=20 function, rather than 'missing(L)'. Rich Raubertas Merck & Co. =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D= =3D=3D=3D=3D=3D=3D=3D=3D=3D > sessionInfo() R version 2.8.1 (2008-12-22)=20 i386-pc-mingw32=20 locale: LC_COLLATE=3DEnglish_United States.1252;LC_CTYPE=3DEnglish_United States.1252;LC_MONETARY=3DEngl...

I know nothing about your package, but "list" is a terrible name for it, as "list" is also the name of a data type in R. -- Jeff

I''m trying to deploy a rails app for development on slicehost and I can''t seem to get mongrel to work correctly. When I run mongrel_rails start I get: ** Starting Mongrel listening at 0.0.0.0:3000 ** Starting Rails with development environment... ** Rails loaded. ** Loading any Rails specific GemPlugins ** Signals ready. TERM => stop. USR2 => restart. INT => stop

Hi, I have noticed the following: > a <- array(1:4, c(2, 2)) > A <- array(1:4, c(2,2,2)) > A + a Error in A + a : non-conformable arrays It works with a matrix + a vector, why doesn't it work with arrays? Am I missing something? How would you do the above operation efficiently (ie I need to add a matrix to each "plane" of 3-dim array)? At the moment I am using

...quot; is in the environment of f1, which was # made available to f2 above. x[i] <<- i^2 NULL } f1() # prints NA 4 NA NA Of course this requires that f1 anticipate that f2 will want to modify its objects; i.e., f1 and f2 must be written as a pair, to know about each other. Richard Raubertas Biometrics Research Merck & Co. > -----Original Message----- > From: Thomas Lumley [mailto:tlumley at u.washington.edu] > Sent: Friday, May 17, 2002 11:44 AM > To: oehl_list at gmx.de > Cc: r-help at stat.math.ethz.ch; Martin.Schlather at uni-bayreuth.de > Subject: Re: [R] W...

In general, how can I increase a vector of length m (< n) to length n by padding it with m - n missing values, without losing attributes? The two approaches I've tried, using length<- and adding missings with c, do not work in general: > a <- as.Date("2008-01-01") > c(a, NA) [1] "2008-01-01" NA > length(a) <- 2 > a [1] 13879 NA > b <-

Dear all, I would like to inquire, if it is a desired feature that the combination with rbind() of two data frames with factors columns does not sort the factors levels of the combined data frame. > str(rbind(data.frame(a = factor(c(4, 3))), data.frame(a = factor(c(2, 1))))) 'data.frame': 4 obs. of 1 variable: $ a: Factor w/ 4 levels

Dear List: My question is more statistical than R oriented (although it originates from my work with nlme). I know statistical questions are occasionally posted, so I hope my question is relevant to the list as I cannot turn up a solution anywhere else. I will frame it in the context of an R related issue. To illustrate the problem, consider student achievement test score data with multiple