Re: [R] Monotonic interpolation

RSiteSearch(monotone, restr=func) will give you several packages and functions for monotone smoothing, including the isoreg() function in the standard stats package. You can determine if any of these does what you want. Bert Gunter Genetech Nonclinical Statistics -Original Message

Re: [R] computing distance in miles or km between 2 street addre

There is a well-known (greeedy) algorithm due to Dijkstra for choosing the shortest path = minimum weight path on a weighted digraph between two vertices. I'm sure numerous open source versions of this are available. optim() is not relevant. Bert Gunter Genentech Nonclinical Statistics

Re: [R] Robust linear models and unequal variance

is expected from a specified distribution. Outliers are no longer out when the variance is large. Also look at glm() with the quasi option if you wish to consider fitting a heterogeneous variance structure to initialize a robust method (which could, of course, be distorted by your outliers). Bert

Re: [R] Saving plot into file

?Devices e.g. ?pdf Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of uv Sent: Friday, August 31, 2007 5:41 AM To: r-help@stat.math.ethz.ch Subject: [R] Saving plot into file Hello. I am using R with Mac X11

Re: [R] retrieve p-value from a cox.obj

the summary.coxph object. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of clearsky Sent: Wednesday, August 29, 2007 8:41 AM To: r-help@stat.math.ethz.ch Subject: [R] retrieve p-value from a cox.obj I have a cox.obj

Re: [R] Excel

Erich: This is not a comment either for or against the use of Excel. I only wish to point out that AFAICS, Hadley Wickham's reshape package offers all the pivot table functionality and more. If I am wrong about this, please let me and everyone else know. Bert Gunter Genentech Nonclinical

Re: [R] Experimental Design with R

Please use R's search tools. RSiteSearch(experimental design, restr = funct) finds optBlock() in the AlgDesign package as the 10th hit. Whether this package will have what you want is another issue. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL

Re: [R] subset using noncontiguous variables by name (not index)

, it is not necessary do this. You can construct the call directly as you tried to do. Using the anscombe example, here's how: chooz - c(x1,x3:x4,y2) ## enclose the desired expression in quotes do.call (subset, list( x = anscombe, select = parse(text = chooz))) -- Bert Gunter Genentech Non-Clinical

Re: [R] uneven list to matrix

),dimnames = list(rn,NULL)) ## use dimnames = sort(rn) if you want to sort them for(i in seq(length(alph))){y - alph[[i]]; mx[names(y),i] - y} Bert Gunter Nonclinical Statistics 7-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Christopher Marcum Sent

Re: [R] Regulatory Compliance and Validation Issues

is required. All changes to such code must of course be documented and validated. I believe this is all part of CFR Part 11 requirements. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Frank E Harrell Jr Sent

Re: [R] Convert factor to numeric vector of labels

add that there are arguments for making stringsAsFactors = FALSE; search the archives for discussions why. The memory penalty will have to be paid, of course. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf

Re: [R] Legend on graph

whitespace to the plot region; or to add a grid to the plot without interfering with the legend. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Nguyen Dinh Nguyen Sent: Monday, August 13, 2007 3:42 PM

Re: [R] Mixture of Normals with Large Data

reflexively throwing some silly statistical methodology at them. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Tim Victor Sent: Tuesday, August 07, 2007 3:02 PM To: r-help@stat.math.ethz.ch Subject: Re: [R

Re: [R] Mixture of Normals with Large Data

of the density. You and others are free to disagree. -- Bert On 8/7/07, Bert Gunter [EMAIL PROTECTED] wrote: Why would anyone want to fit a mixture of normals with 110 million observations?? Any questions about the distribution that you would care to ask can be answered directly from the data. Of course

Re: [R] Create vectors form matrices

will always be arcane to those who do not make a serious effort to learn it. It is **not** meant to be intuitive and easy for casual users to just plunge into. It is far too complex and powerful for that. But the rewards are great for serious data analysts who put in the effort. Bert Gunter Genentech

Re: [R] FW: Selecting undefined column of a data frame (was [BioC]read.phenoData vs read.AnnotatedDataFrame)

programming. Bert Gunter Genentech -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Steven McKinney Sent: Friday, August 03, 2007 10:38 AM To: r-help@stat.math.ethz.ch Subject: [R] FW: Selecting undefined column of a data frame (was [BioC]read.phenoData vs

Re: [R] Extracting a website text content using R

Yes, there are. (Please see and follow the posting guide if you wish to obtain something more specific) Bert Gunter Genetech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Am Stat Sent: Wednesday, August 01, 2007 2:19 PM To: r

Re: [R] Slightly OT - use of R

Why? You might receive more useful replies from a relevant subset of users if you specify the purpose you have in mind. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of John Logsdon Sent: Monday, July 30

Re: [R] generating symmetric matrices

it. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Gregory Gentlemen Sent: Friday, July 27, 2007 8:28 PM To: r-help@stat.math.ethz.ch Subject: [R] generating symmetric matrices Greetings, I have a seemingly simple

Re: [R] Constructing correlation matrices

it. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Gregory Gentlemen Sent: Sunday, July 29, 2007 7:32 PM To: r-help@stat.math.ethz.ch Subject: [R] Constructing correlation matrices Greetings, I have a seemingly

Re: [R] doubt about options(graphics.record=T)

Below is an explicit excerpt from the Help file. How, please is this not clear enough? Bert Gunter Genentech Nonclinical Statistics Recorded plot histories are of class SavedPlots. They have a print method, and a subset method. As the individual plots are of class recordedplot they can

Re: [R] Drawing rectangles in multiple panels

) I want to plot for each panel. It is not clear to me how the panel layout information from panel.number(), etc. would be helpful here instead. Am I correct? -- or is there a smarter way to do this that I've missed? Cheers, Bert Bert Gunter Genentech Nonclinical Statistics -Original

Re: [R] why doesn't as.character of this factor create a vector ofcharacters?

storage intensive way (I think) of doing what you wantthat avoids the transpose and as.vector() conversion would be: mx - as.matrix(df) mx[mx[,a]==Abraham,,drop=TRUE] HTH. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED

Re: [R] Salient feature selection

to be replicable/scientifically meaningful. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Andy Weller Sent: Monday, July 02, 2007 8:17 AM To: R-help@stat.math.ethz.ch Subject: [R] Salient feature selection I am

Re: [R] applying max elementwise to two vectors

Please... use and **read** the docs: ?max --- pmax Bert Gunter -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Afshartous, David Sent: Thursday, June 28, 2007 1:20 PM To: r-help@stat.math.ethz.ch Subject: [R] applying max elementwise to two vectors

Re: [R] lme correlation structures

should get if you want to use R as you said) that exposits the ideas at greater length. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Gareth Hughes Sent: Wednesday, June 27, 2007 7:50 AM To: r-help

Re: [R] moving-window (neighborhood) analysis

See the Spatial section under CRAN's Task views Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Carlos Guâno Grohmann Sent: Wednesday, June 27, 2007 8:27 AM To: r-help@stat.math.ethz.ch Subject: [R] moving

Re: [R] exaustive subgrouping or combination

). Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Waverley Sent: Wednesday, June 27, 2007 1:57 PM To: r-help@stat.math.ethz.ch Subject: [R] exaustive subgrouping or combination Dear Colleagues, I am

Re: [R] open .r files with double-click

the RGUI, too, I believe -- but someone may correct me on this. Bert Gunter Nonclinical Statistics 7-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Duncan Murdoch Sent: Saturday, June 09, 2007 4:29 AM To: [EMAIL PROTECTED] Cc: r-help@stat.math.ethz.ch

Re: [R] rlm results on trellis plot

,y,...) panel.abline(lm(y~x),col=blue,lwd=2) panel.abline(rlm(y~x),col = red,lwd=2) }) Note that the coefficients do not need to be explicitly extracted by coef(), as panel.abline will do this automatically. Bert Gunter Genentech Nonclinical Statistics South San

Re: [R] R is not a validated software package..

analysis software infrastructure from SAS to S-Plus, SPSS, or anything else. **DISCLAINMER** My opinions only. They do not in any way represent the view of my company or its employees. Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message

Re: [R] how to find how many modes in 2 dimensions case

Note that the number of modes (local maxima??) is a function of the bandwidth, so I'm not sure your question is even meaningful. Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED

Re: [R] how to input data from the keyboard

Please do your homework: help.search(input) Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Miguel Caro Sent: Thursday, June 07, 2007 11:01 AM To: r-help@stat.math.ethz.ch Subject: [R] how to input data

Re: [R] Comparing multiple distributions

what the package is for. But plotting the data should always be at least the first thing you do anyway. I often find it to be sufficient, too. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of jiho Sent

Re: [R] runif with weights

replacement Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Ken Knoblauch Sent: Wednesday, May 30, 2007 5:59 PM To: [EMAIL PROTECTED] Cc: r-help@stat.math.ethz.ch Subject: [R] runif with weights Not sure why you have

Re: [R] normality tests [Broadcast]

on hierarchical models (e.g. repeated measures/mixed effect models) has also dealt with lack of independence. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of wssecn Sent: Friday, May 25, 2007 2:59 PM To: r-help

Re: [R] trouble understanding why ...==NaN isn't true

1. NaN is a character string, **not** NaN; hence is.nan(NaN) yields FALSE. 2. Please read the docs! ?NaN explicitly says: Do not test equality to NaN, or even use identical, since systems typically have many different NaN values. Bert Gunter Genentech Nonclinical Statistics -Original

Re: [R] Is it possible to print a data.frame without the row names?

?write.table Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Bos, Roger Sent: Thursday, May 24, 2007 7:17 AM To: r-help@stat.math.ethz.ch Subject: [R] Is it possible to print a data.frame without the row

Re: [R] R-help with apply and ccf

(a[i,],b[i,])) ## Note: value is a vector, not an array Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Michael Andric Sent: Tuesday, May 22, 2007 8:35 AM To: r-help@stat.math.ethz.ch Subject: [R] R-help

Re: [R] Simple programming question

. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Lauri Nikkinen Sent: Friday, May 18, 2007 8:02 AM To: Gabor Grothendieck Cc: r-help@stat.math.ethz.ch Subject: Re: [R] Simple programming question Thank you all

Re: [R] using lm() with variable formula

... and note that if a matrix of responses is on the left of ~ , separate regressions will be simultaneously fit to each of the columns of the matrix. Note that this **is** in TFM -- ?lm. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto

Re: [R] repeated measures regression

also be found on CRAN. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of John Christie Sent: Thursday, May 17, 2007 10:06 AM To: R-help@stat.math.ethz.ch Subject: [R] repeated measures regression How does one go

Re: [R] bug or feature?

... but it **is** explicitly documented in ?subset: For data frames, the subset argument works on the rows. Note that subset will be evaluated in the data frame, so columns can be referred to (by name) as variables in the expression (see the examples). Bert Gunter Genentech Nonclinical

Re: [R] Testing for existence inside a function

parent.frame() is the frame of the caller of foo1, which is what is wanted if foo1 is to work when called within a function. Note that parent.frame() would also work when foo1 is called at the command line. Further corrections/clarifications welcome, of course. Bert Gunter Genentech Nonclinical Statistics

Re: [R] confidence intervals on multiple comparisons

from this, but you must pay attention to the details here. Bert Gunter Genentech Nonclinical statistics Salvatore Enrico Indiogine

Re: [R] Looking for a cleaner way to implement a setting certainindices of a matrix to 1 function

Suggestion: You might make it easier for folks to help if you explained in clear and simple terms what you are trying to do. Code is hard to deconstruct. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf

Re: [R] Neural Nets (nnet) - evaluating success rate of predictions

on it is now also a biased estimate. You need yet new validation data for that. Of course, there are all sort of cross validation schemes one can use to avoid -- or maybe mitigate -- these issues: most books on statistical classification/machine learning discuss this in detail. Bert Gunter Genentech

Re: [R] Multiple scatterplots

maybe something like a trellis plot of the different distributions is what is wanted?? -- but I may well misunderstand. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of John Kane Sent: Wednesday, May 02, 2007

Re: [R] thousand separator (was RE: weight)

regular expressions alone, so I, too, await a clever reply. -- Bert Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Marc Schwartz Sent: Monday, April 30, 2007 10:02 AM To: Liaw, Andy Cc: r-help

Re: [R] thousand separator (was RE: weight)

Nothing! My mistake! gsub -- not sub -- is what you want to get 'em all. -- Bert Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Marc Schwartz Sent: Monday, April 30, 2007 10:18 AM To: Bert Gunter Cc: r

Re: [R] exclude the unfit data from the iteration

?try Wrap each iteration in a try() call Also ?tryCatch if you want to get fancy -- and can understand the rather arcane docs. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Mohammad Ehsanul Karim Sent

Re: [R] fitting mixed models to censored data?

, or even impossible, to check from the observed data. This was ca 1994, I believe, so I don't know whether this view is still held among experts (which I am not). But if it is, you may do well to be careful of whatever SAS does even if you do have to go running off to it. Cheers, Bert Gunter Genentech

Re: [R] summary and min max

I believe it is fair to say that this is where (S3 to keep it simple) classes come in handy: Class the sorts of objects you're working with, say MyClass, and then write your own summary.MyClass() method. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL

Re: [R] inconsistent output using 'round'

of significant digits to be printed in values. The default, NULL, uses getOption(digits) And, lo and behold, your output shows a minimum of 3 **significant** digits with more being used in tables to line up values that are both greater and less than 1. -- Bert Gunter Genentech Non-Clinical

[R] bquote in plot.default vs plot.formula ?

, Bert Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 Version info: R.Version() $platform [1] i386-pc-mingw32 $arch [1] i386 $os [1] mingw32 $system [1] i386, mingw32 $status [1] $major [1] 2 $minor [1] 4.1 $year [1] 2006 $month [1] 12 $day [1

Re: [R] nlm() and optim()

havereproduced results from complex modelling that rely on numerical optimization for model fitting can be difficult. Careful and parsimonious modelling is vital. Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED

Re: [R] creating a data frame from a list

makeFrame(lst) a b c A 1 2 NA B 8 3 2 C NA 0 NA D NA NA 0 Cheers, Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Dimitri Szerman Sent: Thursday, April

Re: [R] Wikibooks

culture. Disclaimer 2: Others have repeatedly made similar comments (asking us to refer people to the docs rather than providing explicit answers, I mean). Cheers, Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL

[R] Completely off topic, but amusing?

does validated software really mean? (Rhetorical question -- no reply sought). Cheers to all, Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] Colored boxes with values in the box

Sounds like ?image what you are looking for, perhaps? Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Pappu, Kartik Sent: Thursday, March 22, 2007 3:42 PM To: r

Re: [R] Bad points in regression

(Robust) will bring up numerous possibilities.rrcov and robustbase are at least two packages devoted to this, but the functionality is available in many others (e.g. rlm() in MASS). Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message

Re: [R] Connecting R-help and Google Groups?

and able to host and maintain our lists as part of the CRAN infrastructure, CRAN maintains control. I think this is wise and prudent. I am happy to be educated to the contrary if I misunderstand how this would work. Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650

Re: [R] How to modify a column of a matrix

?cut ## if you have several bins, where ifelse becomes messy Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Marc Schwartz Sent: Monday, March 12, 2007 11:25 AM

Re: [R] hwo can i get a vector that...

apply(yourMatrix,1,which.max) Bert Gunter Nonclinical Statistics 7-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of bunny , lautloscrew.com Sent: Wednesday, March 07, 2007 2:12 PM To: r-help@stat.math.ethz.ch Subject: [R] hwo can i get a vector

Re: [R] how to edit my R codes into a efficient way

that one uses is something like: ContentsOfFiles - lapply(filenameVector, functionThatReadsFile,additionalParametersto Function) More specifically, ContentsOfFiles - lapply(filenameVector, read.csv, header=TRUE, quote=,fill=TRUE) see ?lapply Bert Gunter Genentech Nonclinical Statistics South San

[R] Off topic:Spam on R-help increase?

Folks: In the past 2 days I have seen a large increase of spam getting into R-help. Are others experiencing this problem? If so, has there been some change to the spam filters on the R-servers? If not, is the problem on my end? Feel free to reply privately. Thanks. Bert Gunter Genentech

Re: [R] Recalling and printing multiple graphs. Is there somethingin the HISTORY menu that will help?

know if this approach works if you don't get a definitive answer elsewhere. Cheers, Bert Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of John Sorkin Sent: Tuesday

Re: [R] from function to its name?

, Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Seth Falcon Sent: Friday, March 02, 2007 9:18 AM To: r-help@stat.math.ethz.ch Subject: Re: [R] from function to its name

Re: [R] R code for Statistical Models in S ?

of stuff not even mentioned there.See also the Green book. Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Charilaos Skiadas Sent: Thursday, March 01, 2007 12:56 PM

Re: [R] Default par() options

par() or X11() using do.call() or something similar; but that seems like more than you need here. Anyway, HTH. -- Bert Gunter Genentech Non-Clinical Statistics South San Francisco, CA The business of the statistician is to catalyze the scientific learning process. - George E. P. Box

Re: [R] Packages in R for least median squares regression and computingoutliers (thompson tau technique etc.)

Packages MASS and robustbase both have this functionality. There may also be others. Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of lalitha viswanath Sent

Re: [R] What is a expression good for?

. There are also special tools available for when you really might want to do this sort of thing -- eg. ?formula, ?terms for altering model specifications. But it is tricky to do right and in full generality -- e.g. ?eval and the above references for some of the issues. Bert Gunter Genentech

Re: [R] fitting of all possible models

... Below -- Bert Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Frank E Harrell Jr Sent: Tuesday, February 27, 2007 5:14 AM To: Indermaur Lukas Cc: r-help

Re: [R] looping

. Bert Gunter Nonclinical Statistics 7-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Neil Hepburn Sent: Monday, February 26, 2007 4:11 PM To: r-help@stat.math.ethz.ch Subject: [R] looping Greetings: I am looking for some help (probably really

Re: [R] Repeated measures in Classification and Regresssion Trees

, Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Andrew Park Sent: Friday, February 23, 2007 7:51 AM To: r-help@stat.math.ethz.ch Subject: [R] Repeated measures

Re: [R] investigating interactions with mixed models

this. As usual, statistical significance and scientific relevance are not equivalent, and the degree of overlap between the two is often difficult to judge. Cheers, Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED

Re: [R] convert to binary to decimal

why not simply: sum(x * 2^(rev(seq_along(x)) - 1)) ? Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 Bert Gunter Nonclinical Statistics 7-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Roland Rau Sent

Re: [R] font size in plots

In general, most methods for R's generic plot command (try: getAnywhere(plot.hclust)) in R's base graphics system accept further arguments in the (...) portion that provide these sorts of capabilities. ?par will tell you about these further graphical parameters. Bert Gunter Genentech

Re: [R] Putting splom in a function

the cast as.expression() gives. Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Roberto Perdisci Sent: Wednesday, February 14, 2007 1:05 PM To: r-help@stat.math.ethz.ch

Re: [R] Problem with subsets and xyplot

of the respective variables in x. so meansbymsa does not have the same number of rows as your original data frame, which it must for subsetting to work properly (meansbymsa[,2] was recycled to be of the right length by default, which produces the nonsense you got. See ?xyplot) Bert Gunter Genentech

Re: [R] R in Industry

... two main drawbacks of R at our firm (as viewed by our IT dept) are lack of guaranteed support as well as the difficulty in finding candidates. -- Just an aside: lack of guaranteed support -- absolutely true in theory, absolutely false in practice. I doubt that the voluntary support found on

Re: [R] lme in R and Splus-7

for the function, I got quite different estimation results from these two software. Anyone has this experience before? Why don't you try searching the archives yourself to see? Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374

Re: [R] strange error in robust package

). Bert Gunter Nonclinical Statistics 7-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Monica Pisica Sent: Monday, February 05, 2007 2:34 PM To: r-help@stat.math.ethz.ch Subject: [R] strange error in robust package Importance: High Hi everybody, I

Re: [R] Loop with string variable AND customizable summary output

,function(z)summary(lm(y~x,data=get(z should work. Note: both of these are untested for the general case where they might be used within a function and may not find the right z unless you pay attention to scope, especially in the get() construction. Bert Gunter Genentech Nonclinical Statistics

Re: [R] strange behaviour with equality after simple subtraction

FAQ on R 7.31. ?all.equal ?identical Have you read these? Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Mike Prager Sent: Friday, January 26, 2007 8:41 AM To: r-help

Re: [R] Robust PCA?

You seem not to have received a reply. You can use cov.rob in MASS or cov.Mcd in robustbase or undoubtedly others to obtain a robust covariance matrix and then use that for PCA. -- Bert Bert Gunter Nonclinical Statistics 7-7374 -Original Message- From: [EMAIL PROTECTED] [mailto

Re: [R] Effect size in GLIM models

by well-designed studies and experiments, hindered by data mining of observational data. I don't think any of this is profound, just sometimes forgotten; however, I would welcome public or private reaction to this comment, and especially refinement/corrections. Bert Gunter Genentech Nonclinical

Re: [R] R editor vs. Tinn-R

, analyses, etc. that I have used on the data. I do not know what is more orderly and useful than that! So would you care to elaborate? Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 __ R-help@stat.math.ethz.ch mailing list https

Re: [R] eval(parse(text vs. get when accessing a function

of Thomas's axiom . Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Ramon Diaz-Uriarte Sent: Friday, January 05, 2007 10:02 AM To: r-help; [EMAIL PROTECTED] Subject: [R] eval(parse(text

[R] Some Windows code for GUI-izing workspace loading

use the recently used files options on all my software all the time and wanted to emulate this for R. Suggestions for improvement (or better yet, code!) or information about bugs or other stupidities gratefully appreciated. Cheers, Bert Gunter Genentech Nonclinical Statistics South San Francisco

Re: [R] na.action and simultaneous regressions

Gunter -Original Message- From: Ravi Varadhan [mailto:[EMAIL PROTECTED] Sent: Wednesday, January 03, 2007 2:15 PM To: 'Bert Gunter'; 'Talbot Katz'; r-help@stat.math.ethz.ch Subject: RE: [R] na.action and simultaneous regressions No, Bert, lm doesn't produce a list each of whose

Re: [R] sorting by name

This is trivial. help([) and An Introduction to R will tell you how. P.S. As earlier posts today have mentioned, stepwise variable selection is generally a bad idea. Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 -Original Message- From: [EMAIL PROTECTED

Re: [R] 2 questions

the questions unworthy of my time to respond to. Whether right or wrong, queries posted in a way that conveys that impression are less likely to elicit good replies. I guess the moral is that on this list anyway, good behavior is rewarded, and bad behavior is ignored. Cheers, Bert Gunter

Re: [R] combinations of m objects into r groups

This issue has come up before: RSiteSearch(nkpartitions) will find references for you on CRAN. You might also try http://ranau.cs.ui.ac.id/book/AlgDesignManual/BOOK/BOOK4/NODE153.HTM for some background, or google on set partitions. Bottom line: it ain't trivial. Cheers, Bert Gunter

Re: [R] Remove from a string

or something similar would be the better approach. Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Marc Schwartz Sent: Friday, December 08, 2006 6:14 AM To: Katharina Vedovelli Cc: r

Re: [R] Summary shows wrong maximum

Folks: Is So this is at best a matter of opinion, and credentials do matter for opinions. -- Brian Ripley an R fortunes candidate? -- Bert Gunter On Tue, 5 Dec 2006, Oliver Czoske wrote: On Mon, 4 Dec 2006, Uwe Ligges wrote: Sebastian Spaeth wrote: Hi all, I have a list

Re: [R] Usage of apply

prefer them. (They are also written to do their looping as efficiently as possible, which explicit looping in user code may not.) Of course, vectorized calculations (colMeans() in the example below) **are** much faster and usually clearer than explicit loops. Bert Gunter Genentech Nonclinical

Re: [R] Summary shows wrong maximum

... -- Bert Bert Gunter Nonclinical Statistics 7-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Mike Prager Sent: Wednesday, December 06, 2006 11:46 AM To: r-help@stat.math.ethz.ch Subject: Re: [R] Summary shows wrong maximum I don't know about candidacy

Re: [R] stat question - not R question so ignore if not interested

Monday night beer and statistics sessions: after you're done and perhaps have written up and presented your (intricate!) analysis, you're always worried that someone might come along and say, Well, did you consider...? Cheers, Bert Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA

Re: [R] tests for NULL objects

standard conventions, but even if you disagree, they are certainly not worth making a fuss over and certainly cannot be changed without breaking a lot of code, I'm sure. Bert Gunter Nonclinical Statistics 7-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED

  1   2   >