Re: [R] Silhouette question

2003-06-20 Thread Prof Brian Ripley
?try is your friend here. On Sat, 21 Jun 2003, Jonck van der Kogel wrote: > Hi all, > I am momentarily experimenting with Silhouette from the cluster library > but I am getting some errors. Since Silhouette can be seen as a quality > measure for a clustering what I want to do is run a series of

Re: [R] combining mathematical notation and value substitution

2003-06-20 Thread Faheem Mitha
> > "Monotonic Multigamma run (" * n == len * ", " * theta == t1 * > > > > Is this a valid expression? My understanding of an expression is that it > > contains one more more statements. > > That's only part of the expression. This is the full expression > "Monotonic Multigamma run (" * n == l

Re: [R] combining mathematical notation and value substitution

2003-06-20 Thread Thomas Lumley
On Fri, 20 Jun 2003, Faheem Mitha wrote: > > > On Fri, 20 Jun 2003, Uwe Ligges wrote: > > > > > t1 <- theta # you cannot use theta as variable and math. symbol > > plot(1:10, main = > >substitute("Monotonic Multigamma run (" * n == len * ", " * > > theta == t1 * ").", list(len = len,

[R] [OFF] stepwise using REML???

2003-06-20 Thread Ronaldo Reis Jr.
Hi, I know that is not possible make a stepwise procedure using REML in R, I can use ML for this. For nested design it may be very dangerous due the difference in variance structure, mainly in a splitplot design. ML make significative variables that REML dont make. I read an article that is m

[R] negative binomial regression

2003-06-20 Thread Ross Nelson
I have been using the negative binomial regression function in the MASS library to generate incidence rate ratios and confidence intervals. Can this function also generate a gradient value that would estimate the slope between the dependent and independent variables? or can it be generated by

Re: [R] combining mathematical notation and value substitution

2003-06-20 Thread Faheem Mitha
On Fri, 20 Jun 2003, Uwe Ligges wrote: > Faheem Mitha wrote: [snip] > > main=paste("Monotonic Multigamma run (n=", > > deparse(substitute(len)),", ", > > expression(theta),"=", deparse(substitute(theta)),").") > > No, it won't work that way, because you have to specify an S expression > in order

[R] Silhouette question

2003-06-20 Thread Jonck van der Kogel
Hi all, I am momentarily experimenting with Silhouette from the cluster library but I am getting some errors. Since Silhouette can be seen as a quality measure for a clustering what I want to do is run a series of different clusterings and store the one with the highest Silhouette value. In tha

Re: [R] stepwise regression

2003-06-20 Thread Spencer Graves
Consider the following: > library(MASS) > Dat <- data.frame(x1=rnorm(9), x2=rnorm(9), y1=rnorm(9), y2=rnorm(9)) > mdl <- "~x1+x2+I(x1^2)+I(x2^2)+x1:x2" > y <- names(Dat)[3:4] > fits <- list() > for(i in 1:2){ + fit0 <- lm(formula(paste(y[i], "~1")), Dat) + fits[[i]] <- stepAIC(fit0, mdl) + } Sta

[R] glmm and overall goodness of fit

2003-06-20 Thread Ronaldo Reis Jr.
Hi, exist in R any glmm function that have any tools for test for overall goodness of fit? Thanks Ronaldo -- O papel da impressora é sempre mais forte na parte picotada. -- | // | \\ [***] |> ( õ õ ) [Ronaldo Reis Júnior] | V [UF

Re: [R] stepwise regression

2003-06-20 Thread Roger D. Peng
Try, help.search("stepwise") It brings up the functions step() and stepAIC() from MASS. Andrew Taylor wrote: Hi, S-PLUS includes the function "stepwise" which can use a variety of methods to conduct stepwise multiple linear regression on a set of predictors. Does a similar function exist in R?

Re: [R] stepwise regression

2003-06-20 Thread Spencer Graves
Have you considered stepAIC in library(MASS) or step? hth. spencer graves Andrew Taylor wrote: Hi, S-PLUS includes the function "stepwise" which can use a variety of methods to conduct stepwise multiple linear regression on a set of predictors. Does a similar function exist in R? I'm having d

[R] stepwise regression

2003-06-20 Thread Andrew Taylor
Hi, S-PLUS includes the function "stepwise" which can use a variety of methods to conduct stepwise multiple linear regression on a set of predictors. Does a similar function exist in R? I'm having difficulty finding one. If there is one it must be under a different name because I get an error m

[R] RE: Application

2003-06-20 Thread kemasse
There are no attached files. >= Original Message From [EMAIL PROTECTED] = >See the attached file for details. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

RE: [R] Spedd: R vs S-plus

2003-06-20 Thread Ted Harding
On 20-Jun-03 Ted Harding wrote: > Sorry to raise what has probably been discussed before, > but I an repeatedly struck by the comparative slowness > of S-plus for Windows compared with R for Linux when doing > much the same thing. Thanks to all who so promptly responded with comments and informati

Re: [R] namespaces not available

2003-06-20 Thread Göran Broström
On Fri, 20 Jun 2003, Prof Brian Ripley wrote: > Might it help to note that the version of MASS which shipped with 1.7.0 > didn't get the import of biplot from mva right? I think this was what > happened with that version of VR. I upgraded to R-1.7.1 and now the old .RData file loads when I star

Re: [R] Spedd: R vs S-plus

2003-06-20 Thread Prof Brian Ripley
On Fri, 20 Jun 2003, Uwe Ligges wrote: > (Ted Harding) wrote: > > Hi Folks, > > > > Sorry to raise what has probably been discussed before, > > but I an repeatedly struck by the comparative slowness > > of S-plus for Windows compared with R for Linux when doing > > much the same thing. > > > > I

Re: [R] Spedd: R vs S-plus

2003-06-20 Thread A.J. Rossini
(Ted Harding) <[EMAIL PROTECTED]> writes: > My instincts say that "WIndows" is probably a significant > factor in the comnparison, but still ... > > Ideally, to compare R with S-plus, one should look at them > both on the same OS (Unix or Windows) on the same machine. > > Can anyone give me clean

[R] Fitting particular repeated measures model with lme()

2003-06-20 Thread J.R. Lockwood
Hello, I have a simulated data structure in which students are nested within teachers, and with each student are associated two test scores. There are 20 classrooms and 25 students per classroom, for a total of 500 students and two scores per student. Here are the first 10 lines of my dataframe

RE: [R] Spedd: R vs S-plus

2003-06-20 Thread Philippe Grosjean
You will find a speed comparison of S-PLUS and R, but also Matlab, Octave, Scilab, Ox and O-Matrix under Windows at http://www.sciviews.org/other/benchmark.htm. It seems that, at least under Windows your impressions are confirmed. Also, I measured that R is a little bit faster under Linux Mandrake

Re: [R] Power Law Exponents

2003-06-20 Thread Spencer Graves
Have you looked at coef(lm (log10(vouty) ~ log10(voutx)))? This the slope / coefficient of log10(voutx) should be an estimate of the power. Why do yo think the exponent is wrong? If the plot looks like a rough ellipse, then this is a standard way to fit this. This assumes that log10(

Re: [R] Spedd: R vs S-plus

2003-06-20 Thread Uwe Ligges
(Ted Harding) wrote: Hi Folks, Sorry to raise what has probably been discussed before, but I an repeatedly struck by the comparative slowness of S-plus for Windows compared with R for Linux when doing much the same thing. I don't have a direct comparison, because they're not running on the same ma

[R] Spedd: R vs S-plus

2003-06-20 Thread Ted Harding
Hi Folks, Sorry to raise what has probably been discussed before, but I an repeatedly struck by the comparative slowness of S-plus for Windows compared with R for Linux when doing much the same thing. I don't have a direct comparison, because they're not running on the same machine; but machine W

Re: [R] namespaces not available

2003-06-20 Thread Luke Tierney
On Fri, 20 Jun 2003, Göran Broström wrote: > This happened to me with R-1.7.0, Linux RH9: > > -- > > q() > Save workspace image? [y/n/c]: y > Warning messages: > 1: namespaces may not be available when loading > 2: names in persist

[R] Power Law Exponents

2003-06-20 Thread Martin H. Robinson
I am having difficulty with the calculation of the power law exponent for set of nodes within a graph. Specifically, I am interested in the distribution of in-degree and out-degree among communities of web pages where the web pages are the nodes of the graph and the hyperlinks the edges. Accordi

Re: [R] namespaces not available

2003-06-20 Thread Prof Brian Ripley
I think you have saved in your workspace an object whose environment is a namespace, probably mva, that is not available when .RData is loaded but is once the session is running. On Fri, 20 Jun 2003, Göran Broström wrote: > This happened to me with R-1.7.0, Linux RH9: > > --

[R] namespaces not available

2003-06-20 Thread Göran Broström
This happened to me with R-1.7.0, Linux RH9: -- > q() Save workspace image? [y/n/c]: y Warning messages: 1: namespaces may not be available when loading 2: names in persistent strings are currently ignored [EMAIL PROTECTED]:~/R/t

Re: [R] combining mathematical notation and value substitution

2003-06-20 Thread Uwe Ligges
Faheem Mitha wrote: Dear People, I need to make a label which both contains math notation as well as substitutes a value for an object. In the following label, len and theta are one dim variables, and I am substituting their values appropriately. This label looks fine except that I want the greek

[R] Re: combining mathematical notation and value substitution

2003-06-20 Thread Faheem Mitha
On Fri, 20 Jun 2003, Faheem Mitha wrote: > Hmm. Tried this, didn't work either. Inspired by pg 32 of "R for > Beginners". > > main=paste("Monotonic Multigamma run", >as.expression(substitute(n==length,list(length=len))), >as.expression(substitute(theta==th,list(th=theta)

[R] Re: combining mathematical notation and value substitution

2003-06-20 Thread Faheem Mitha
On Fri, 20 Jun 2003, Faheem Mitha wrote: > > Dear People, > > I need to make a label which both contains math notation as well as > substitutes a value for an object. In the following label, len and theta > are one dim variables, and I am substituting their values appropriately. > This label loo