RE: [R] the huge postscript plot

2003-07-04 Thread Adaikalavan Ramasamy
How about pdf(), jpeg() or png() ? -Original Message- From: Philippe Glaziou [mailto:[EMAIL PROTECTED] Sent: Saturday, July 05, 2003 9:47 AM To: [EMAIL PROTECTED] Subject: Re: [R] the huge postscript plot Yongchao Ge <[EMAIL PROTECTED]> wrote: > I'm just wondering how I can do to make a

Re: [R] the huge postscript plot

2003-07-04 Thread Philippe Glaziou
Yongchao Ge <[EMAIL PROTECTED]> wrote: > I'm just wondering how I can do to make a huge postscript plot be > manageable. My question is that I have to draw around 60,000 points which > makes it painfully slow to print or view in gv or put it into latex > document, though it is very fast to produce

[R] the huge postscript plot

2003-07-04 Thread Yongchao Ge
Hi, I'm just wondering how I can do to make a huge postscript plot be manageable. My question is that I have to draw around 60,000 points which makes it painfully slow to print or view in gv or put it into latex document, though it is very fast to produce the postscript file. A simple example is

Re: [R] Quasi AIC

2003-07-04 Thread Uwe Ligges
Henric Nilsson wrote: > > Dear all, > > Using the quasibinomial and quasipoisson families results in no AIC being > calculated. However, a quasi AIC has actually been defined by Lebreton et al > (1992). In the (in my opinon, at least) very interesting book by Burnham and > Anderson (1998,2002) th

Re: [R] parametros

2003-07-04 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: > > hello, sorry i´m a beginners in R, i use the version 1.6.1 This one is outdated. > for windows and 1.5.5 Such a version was never released, and the below mentioned R-1.5.1 is outdated. > This one for linux (susse), I guess you mean SuSE? > my question is abo

Re: [R] specifying dataframe column names in loops

2003-07-04 Thread Prof Brian Ripley
Sounds like you need to use substitute(). Something like eval(substitute(lf <- locfit(~s, data=age), list(s=s))) Formulae are different: I don't see why density() should be. On Fri, 4 Jul 2003, Simon Goodman wrote: > Dear list members, > > This is probably a naive question (from a new user), b

[R] specifying dataframe column names in loops

2003-07-04 Thread Simon Goodman
Dear list members, This is probably a naive question (from a new user), but I'm having problems getting column names to be recognised in loop that is supposed to apply a set of calculations in turn to each column of a dataframe. The dataframe has 33 columns of output from a simulation. The colum

Re: [R] Problem with fitdistr for beta

2003-07-04 Thread Prof Brian Ripley
In this example shrinking by (1 - 2e-16) leads to a significant change in the distribution: see my probability calculation. And you can't shrink by much less. A beta(0.1, 0.1) is barely a continuous distribution. On Fri, 4 Jul 2003, Spencer Graves wrote: > My standard work-around for the kind

Re: [R] anova.glm() and deviance/df

2003-07-04 Thread Prof Brian Ripley
On Fri, 4 Jul 2003, Henric Nilsson wrote: > Is seems to me that R doesn't allow correct F-tests to be computed in anova.glm > after estimating the dispersion parameter by deviance/df. Or does it? It does when they are correct, e.g. for a Gaussian family. -- Brian D. Ripley, [E

[R] anova.glm() and deviance/df

2003-07-04 Thread Henric Nilsson
Dear R-help, Is seems to me that R doesn't allow correct F-tests to be computed in anova.glm after estimating the dispersion parameter by deviance/df. Or does it? //Henric __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo

[R] Quasi AIC

2003-07-04 Thread Henric Nilsson
Dear all, Using the quasibinomial and quasipoisson families results in no AIC being calculated. However, a quasi AIC has actually been defined by Lebreton et al (1992). In the (in my opinon, at least) very interesting book by Burnham and Anderson (1998,2002) this QAIC (and also QAICc) is covered.

Re: [R] Problem with fitdistr for beta

2003-07-04 Thread Spencer Graves
My standard work-around for the kind of problem you identified is to shrink the numbers just a little towards 0.5. For example: > library(MASS) > a <- rbeta(100,0.1,0.1) > fitdistr(x=a, "beta", start=list(shape1=0.1,shape2=0.1)) Error in optim(start, mylogfn, x = x, hessian = TRUE, ...) :

[R] parametros

2003-07-04 Thread solares
hello, sorry i´m a beginners in R, i use the version 1.6.1 for windows and 1.5.5 for linux (susse), my question is about ¿how i cant to knowledge the parameter of a command tcltk in R?, in Windows i write tt<-tktoplevel() t <- tkmenu(tt) tkconfigure(tt,menu=t) a<- tkmenu(t, tearoff=FALSE) tkcon

[R] parametros

2003-07-04 Thread solares
hello, sorry i´m a beginners in R, i use the version 1.6.1 for windows and 1.5.5 for linux (susse), my question is about ¿how i cant to knowledge the parameter of a command tcltk in R?, in Windows i write tt<-tktoplevel() t <- tkmenu(tt) tkconfigure(tt,menu=t) a<- tkmenu(t, tearoff=FALSE) and

Re: [R] Problem with fitdistr for beta

2003-07-04 Thread Prof Brian Ripley
rbeta(100,0.1,0.1) is generating samples which contain 1, an impossible value for a beta and hence the sample has an infinite log-likelihood. It is clearly documented on the help page that the range is 0 < x < 1. However, that is not so surprising as P(X > 1-1e-16) is about 1% and hence values wi

Re: [R] How to use quasibinomial?

2003-07-04 Thread Henric Nilsson
I think I can now answer my own questions: Henric Nilsson ([EMAIL PROTECTED]) wrote*: >1. I can't seem to get the correct p-values from anova.glm() for the F-tests when >supplying the dispersion argument and having fitted the model using >family=quasibinomial. Using family=quasibinomial does exa

[R] Problem with fitdistr for beta

2003-07-04 Thread Agustin Lobo
I have the following problem: I have a vector x of data (0 a <- rbeta(100,0.1,0.1) > fitdistr(x=a, "beta", start=list(shape1=0.1,shape2=0.1))1) > shape1 shape2 0.09444627 0.12048753 (0.01120670) (0.01550129) but sometimes does not: > a <- rbeta(100,0.1,0.1) > fitdistr(x=a, "beta"

Re: [R] configuration error when installing gtkDevice

2003-07-04 Thread Dimitri Giunchi
Thank you very much to everybody for the help: now it works fine. Dimitri At 08.51 04/07/2003 -0400, you wrote: When compiling the gtkDevice package, we need the gtk+ development libraries, i.e. libgtk.so and header files such as gtk.h. The error message saying 'Cannot find gtk-config' means

Re: [R] configuration error when installing gtkDevice

2003-07-04 Thread Duncan Temple Lang
When compiling the gtkDevice package, we need the gtk+ development libraries, i.e. libgtk.so and header files such as gtk.h. The error message saying 'Cannot find gtk-config' means that we cannot locate the script that is used to tell us where the header files and libraries are located. So this

[R] away from my mail, absent pour le moment

2003-07-04 Thread via the vacation program
I will not be reading my mail from 30_juin_2003 to 09_juillet_2003 Your mail regarding "Re: Application" will be read then. Je ne lirai pas mon courrier du 30_juin_2003 au 09_juillet_2003 Votre message concernant "Re: Application" sera lu a mon retour. __

Re: [R] Looking for graphics/par cheat sheet

2003-07-04 Thread Martin Maechler
> "BDR" == Prof Brian Ripley <[EMAIL PROTECTED]> > on Thu, 3 Jul 2003 16:15:19 +0100 (BST) writes: BDR> There is one in MASS4 p.85. It's a direct descendant BDR> from Bill Venables' `Notes on S'. `An Introduction to BDR> R' is another descendant, and that has two less

Re: [R] configuration error when installing gtkDevice

2003-07-04 Thread Philippe Glaziou
Dimitri Giunchi <[EMAIL PROTECTED]> wrote: > when I try to install gtkDevice on Mandrake Linux 9.1 (R > v. 1.7.1), I get the following comfiguration error (see > below). > > Has anyone else had this problem? > Any hints are greatly appreciated. > > [...] > > trying URL `http://cran.r-project.org/

Re: [R] configuration error when installing gtkDevice

2003-07-04 Thread Prof Brian Ripley
On Fri, 4 Jul 2003, Dimitri Giunchi wrote: > when I try to install gtkDevice on Mandrake Linux 9.1 (R v. 1.7.1), I get > the following comfiguration error (see below). > > Has anyone else had this problem? Yes. Do you have the appropriate gtk development RPMs installed? My current RH8.0 machine

[R] configuration error when installing gtkDevice

2003-07-04 Thread Dimitri Giunchi
Dear all, when I try to install gtkDevice on Mandrake Linux 9.1 (R v. 1.7.1), I get the following comfiguration error (see below). Has anyone else had this problem? Any hints are greatly appreciated. Thank you in advance, Dimitri install.packages("gtkDevice") trying URL `http://cran.r-projec

Re: [R] Help. Import data

2003-07-04 Thread Prof Brian Ripley
On Fri, 4 Jul 2003, Fredrik Karlsson wrote: > Hi Silika, > > This is what I would do: > > >df <- read.spss("datafile", to.data.frame=T) > >for(i in c("age","income")){ #Add more if you like > + df[[i]] <- as.numeric(df[[i]]) > +} Please don't: that does not turn a factor into its numerical equi

Re: [R] Help. Import data

2003-07-04 Thread Fredrik Karlsson
Hi Silika, This is what I would do: >df <- read.spss("datafile", to.data.frame=T) >for(i in c("age","income")){ #Add more if you like + df[[i]] <- as.numeric(df[[i]]) +} /Fredrik On Fri, Jul 04, 2003 at 08:39:07AM +0200, Silika Tereshchenko wrote: > Hoy, > > > I have the following problem.

Re: [R] Help. Import data

2003-07-04 Thread Prof Brian Ripley
How to convert a factor to numeric is one of the questions (with answers) in R's FAQ. Assuming you don't need to do this again, I would just convert what you have imported. Otherwise, explore the max.value.labels argument of read.spss. On Fri, 4 Jul 2003, Silika Tereshchenko wrote: > I have