search for: stahel

...ue where logical needed For any other equal from and to , it works: > seq(1,1,1) [1] 1 The error occurs in the statement dd <- abs(del)/max(abs(to), abs(from)) of seq.default -- for obvious reasons. This is certainly easy to debug... Thank you for your efforts. Werner Stahel ----------------- This message was sent by --------------------------- Werner Stahel http://stat.ethz.ch/~stahel Seminar fuer Statistik phone : +41 1 632 34 30 ETH-Zentrum, LEO D8 fax : +41 1 632 12 28 CH-8092 Zurich...

>>>>> "WSt" == Werner Stahel <stahel@stat.math.ethz.ch> writes: [to me privately] WSt> I do not want to call this a bug, but it bugged me. Please try the WSt> following: WSt> tform <- sqrt(RADAI) ~ sqrt(RADAI.e) + TAGE.ej + SPITAL + ARZT + DAS28 + WSt> SJC + TJC + DGA + HAQ + PGA + PAIN +...

...lm model containing a bs smooth term. I asked Martin how I could see the estimated function. He said that such a possibility -- which is available in S-Plus' plot.gam -- was on your list of desirable functions. This is a vote to give it high priority. (I appreciate all your efforts!) Werner Stahel ------------------ This message was sent by --------------------------------- Werner Stahel http://stat.ethz.ch/~stahel Seminar fuer Statistik phone : +41 1 632 34 30 ETH-Zentrum, LEO fax : +4...

...dit the information below-- Version: platform = i686-pc-linux-gnu arch = i686 os = linux-gnu system = i686, linux-gnu status = major = 1 minor = 1.1 year = 2000 month = August day = 15 language = R Search Path: .GlobalEnv, d.complexity, package:SfS, Autoloads, package:base -- Werner Stahel ------------------ This message was sent by --------------------------------- Werner Stahel http://stat.ethz.ch/~stahel Seminar fuer Statistik phone : +41 1 632 34 30 ETH-Zentrum, LEO fax : +4...

Hallo, I'm looking for a "sunflower" plot. I noticed the Splus function "p.sunflowers" (Andreas Ruckstuhl, Werner Stahel, Martin Maechler, Tim Hesterberg) don't work... If its not available I would port this function... Thanks for 1 answer P *************************************************************************** P.Malewski Tel.: 0531 500965 Maschplatz 8 Email: P.Malewski at tu-bs.de **************...

As some of you R-devel readers may know, the plot() method for "lm" objects is based in large parts on contributions by John Maindonald, subsequently "massaged" by me and other R-core members. In the statistics litterature on applied regression, people have had diverse oppinions on what (and how many!) plots should be used for goodness-of-fit / residual diagnostics, and to my

>>>>> On Tue, 31 Aug 1999 14:57, Werner Stahel <stahel@stat.math.ethz.ch> said: WSt> Here is a suggestion. It seems that qqplots, comparing a sample WSt> to a distribution other than the normal, are not explicitly WSt> available in S or R. I found (in S-plus / Trellis it is, see below) WSt> qqplot(y, rt(...

Martin Maechler <maechler at stat.math.ethz.ch> wrote: > Werner Stahel (in our stat group) has been using hacked versions of lm > and some hacked lm methods which exactly address this, > i.e. they follow the "b)" approach; however I think that it's still a > hack that only works in some (most used) cases. > One would have to change quite a...

Dear All, When I run rlm to obtain robust standard errors, my output does not include p-values. Is there any reason p-values should not be used in this case? Is there an argument I could use in rlm so that the output does include p-values? Thanks in advance, Celso [[alternative HTML version deleted]]

I'd like to get interaction terms in a model to be in another form. Namely, suppose I had variables age and group, the latter a factor with levels A, B, C, with age * group in the model. What I would like are the variables "age:group=A", "age:group=B" and "age:group=C" (and group itself of course). The coefficients of the model will then be the age effect

...ow return a P-value of NaN in cases with df=0 for test="Chisq" (they gave NaN until pchisq was altered in R 2.6.2). o The default methods of drop1() and add1() now re-evaluate in the object's formula environment, rather than simply in the parent (informal bug report by Werner Stahel). o The undocumented limit of 8191 bytes on comments (from # to the end of the line) in R code has been removed. o qgamma() warns if used with shape < 1e-10, when it may be unreliable. (PR#12324) o Installing a source package was not using the package's declared encoding...

...ow return a P-value of NaN in cases with df=0 for test="Chisq" (they gave NaN until pchisq was altered in R 2.6.2). o The default methods of drop1() and add1() now re-evaluate in the object's formula environment, rather than simply in the parent (informal bug report by Werner Stahel). o The undocumented limit of 8191 bytes on comments (from # to the end of the line) in R code has been removed. o qgamma() warns if used with shape < 1e-10, when it may be unreliable. (PR#12324) o Installing a source package was not using the package's declared encoding...

I'm sorry if this is a basic question, but I'm stumped. I'm just trying to plot the residuals from a linear model against another variable in the data frame. Here are the lines I'm trying to execute: size <- read.table(file="/u67/abasl70/surveys/annenberg/mega/smschl.dat", header=T) sizef <- data.frame(size, row.names=size$unit) attach(sizef) mschmod <- lm