Re: [R] help

2006-04-26 Thread Dieter Menne
Anamika Chaudhuri yahoo.com> writes: > I was trying to calculate an age adjusted rate for which we have the formula as > sum(weight*rate)/sum(weight). I am using the framingham dataset called fram here.I cant figure out how I > could change the logical statement I have for the standard dat

[R] plot cdf

2006-04-26 Thread Nongluck Klibbua
hi, I would like to know what is the function to plot cdf. Thanks, Luck __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] Memory usage and limit

2006-04-26 Thread Paul Roebuck
On Wed, 26 Apr 2006, Min Shao wrote: > I recently made a 64-bit build of R-2.2.1 under Solaris 9 using gcc v.3.4.2. > The server has 12GB memory, 6 Sparc CPUs and plenty of swap space. I was the > only user at the time of the following experiment. > > I wanted to benchmark R's capability to read l

[R] What are the differences between ACF and PACF in time series analysis?

2006-04-26 Thread Michael
Hi all, I am desperately looking for answer to my previous question: what are the differences between ACF and PACF in time series and their applications? I got confused a lot by these two functions in R... Already having ACF, why do people decide to create PACF? Thanks a lot [[alternativ

Re: [R] Incomplete Trio in TDT analysis

2006-04-26 Thread Steve Buyske
Using only informative incomplete trios and dropping others turns out to bias the TDT (Curtis & Sham 95). There are a variety of modifications to handle case-parent trios, starting (I think) with Weinberg's 1999 paper that uses the EM algorithm and a likelihood-ratio test. I don't know of any o

Re: [R] problem with get command [Broadcast]

2006-04-26 Thread Liaw, Andy
Sorry, couldn't resist... From: Gabor Grothendieck > > ov$vn1 is not a variable. It is the result of applying the $ > function to the ov and vn1 arguments. ... which would suggest: get("$")(get("ov"), "vn1"). 8-) Andy > For example, using BOD which is a data frame that comes with > R, ra

[R] help

2006-04-26 Thread Anamika Chaudhuri
Hi, I was trying to calculate an age adjusted rate for which we have the formula as sum(weight*rate)/sum(weight). I am using the framingham dataset called fram here.I cant figure out how I could change the logical statement I have for the standard dataset (in red) to calculate weight an

[R] Incomplete Trio in TDT analysis

2006-04-26 Thread Farrel Buchinsky
I am involved in a study where, as in most of life, men demonstrate themselves to be recalcitrant. So while we have many probands and most of their mothers we only have about 50% of the trios being complete. I have been running tdt and trio.types. It appears as if it is ignoring the duos. Sometime

Re: [R] stl function

2006-04-26 Thread Gabor Grothendieck
stl does this internally: x <- na.action(as.ts(x)) so stl(x, s.window, na.action = f) is the same as stl(f(as.ts(x)), s.window) e.g. nottem[25] <- NA # nottem is a built in data set in R stl(nottem, "per") # error stl(nottem, "per", na.action = na.contiguous) library(zoo) stl(nottem, "per",

Re: [R] problem with get command

2006-04-26 Thread Gabor Grothendieck
ov$vn1 is not a variable. It is the result of applying the $ function to the ov and vn1 arguments. For example, using BOD which is a data frame that comes with R, rather than get("BOD$Time") use get("BOD")[["Time"]] On 4/26/06, Thomas Davidoff <[EMAIL PROTECTED]> wrote: > I don't understand what

Re: [R] help in R

2006-04-26 Thread Anupam Tyagi
Anamika Chaudhuri yahoo.com> writes: > > Hi, > > I cant understand where I am going wrong.Below is my code.I would really > appreciate your help. Dear Anamika, to improve your chances to get a response from this list, and to help in people responde quickly and accurately to you, please run

[R] problem with get command

2006-04-26 Thread Thomas Davidoff
I don't understand what my error is in the following: I need to use the get command on a series of variables, but can't for some reason that I don't understand. Why am I told no such variable as ov$vn1 after getting a summary report on that very variable? > summary(ov$vn1) Min. 1st Qu. Med

[R] Sign of loadings and scores from PCA in cross validation

2006-04-26 Thread Anne S Jacobson
Hi, The help file on princomp says that 'The signs of the columns of the loadings and scores are arbitrary, and so may differ between differnt programs for PCA, and even between different builds of R'. Does anyone know 1. Whether this applys to all functions like svd, eigen ... etc? 2. Duri

Re: [R] Polygon-like interactive selection of plotted points

2006-04-26 Thread Marc Schwartz (via MN)
boundary. I tried the followingm using identify() instead: library(splancs) set.seed(20060426) xy <- cbind(x=runif(100), y=runif(100)) plot(xy) # Select the SAME 7 vertices here pts <- identify(xy, n = 7) # These are the row indices into 'xy' # for the c

[R] Memory usage and limit

2006-04-26 Thread Min Shao
Hello everyone, I recently made a 64-bit build of R-2.2.1 under Solaris 9 using gcc v.3.4.2. The server has 12GB memory, 6 Sparc CPUs and plenty of swap space. I was the only user at the time of the following experiment. I wanted to benchmark R's capability to read large data files and used a dat

[R] left-truncation in survreg

2006-04-26 Thread singyee ling
Dear R-users, I know that a few people have asked whether survreg handles left-truncation data and the reply that i have seen so far is that it does. However, when I try to use survreg on left-truncated data, I got the following error message. > survcs3<-survreg(Surv(start,end,status)~AG, data=

Re: [R] Polygon-like interactive selection of plotted points

2006-04-26 Thread Roger Bivand
this (though > > > it is always possible). > > > > > > However, using: > > > > > > RSiteSearch("points inside polygon") > > > > > > brings up several function hits that, if put together with the above > > > interactive f

Re: [R] Were to find appropriate functions for a given task in R

2006-04-26 Thread Erik Iverson
> Then, people tend to > define its own functions (I'm doing this too), and a lack of > standardization makes it difficult to keep everything into control. If you think of R as more of a language rather than a pre-packaged statistical program, I feel that helps. In the C++ world, there are p

Re: [R] Polygon-like interactive selection of plotted points

2006-04-26 Thread Clint Bowman
if put together with the above > > interactive functions, could be used to do what you wish. That is, input > > the matrix of x,y coords of the interactively selected polygon and the > > x,y coords of the underlying points set to return the points inside or > > outside the polyg

Re: [R] copying previously installed libraries to R 2.3.0

2006-04-26 Thread Gabor Grothendieck
Note that the FAQ is really no different from batchfiles since the FAQ does not address how to move or copy the packages and the batchfiles scripts just do that (with the safeguard that they will not overwrite any packages that are already there so you could, for example, manually install a few pac

[R] comparing mantel tests for three different species

2006-04-26 Thread Patrick Kuss
Hi, I have calculated mantel statistics for three plant species to check for isolation by distance (thus, comparing genetic distance to geographic distance). Now, I would like to test whether the species have significantly different isolation by distance behaviour. Is there a R-function that all

[R] Fwd: Re: Help needed

2006-04-26 Thread Anamika Chaudhuri
Note: forwarded message attached. - --- Begin Message --- Subject: Re: [R] Help needed To: Petr Pikal <[EMAIL PROTECTED]> In-Reply-To: <[EMAIL PROTECTED]> MIME-Version: 1.0 Content-Type: multipart/alternative; boundary="0-1025507389-1146080390=:929

[R] Course - May 2006, R/Splus Fundamentals and Programming Techniques***San Francisco, Washington DC

2006-04-26 Thread elvis
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce 2-day "R/S-plus Fundamentals and Programming Techniques" in San Francisco: www.xlsolutions-corp.com/Rfund.htm San Francisco *** May 25-26 Washington, DC *** May TBD Reserve your seat now at the early

Re: [R] www.r-project.org

2006-04-26 Thread Adaikalavan Ramasamy
I am coming late into the discussion, so apologies if the following points are redundant. 1) IMHO, the most important feature that would make life a lot easier for everyone is having search engines on the main webpage. I know you can click on the "Search" on the left hand side pane but putting it

Re: [R] cdf of weibull distribution

2006-04-26 Thread J. Hosking
Peter Dalgaard wrote: > Sachin J <[EMAIL PROTECTED]> writes: > > >>Hi, >> >> I have a data set which is assumed to follow weibull distr'. How can I find of cdf for this data. For example, for normal data I used (package - lmomco) >> >> >cdfnor(15,parnor(lmom.ub(c(df$V1 If X is a We

[R] help in R

2006-04-26 Thread Anamika Chaudhuri
Hi, I cant understand where I am going wrong.Below is my code.I would really appreciate your help. Thanks. > genfile<-read.table("c:/tina/phd/bs871/hw/genfile.txt",skip=1) > > #read in SNP data > snp.dat <- as.matrix(genfile) > snp.name <- scan("c:/tina/phd/bs871/hw/genfile.

Re: [R] copying previously installed libraries to R 2.3.0

2006-04-26 Thread Thomas Harte
the windoze faq that you refer to doesn't quite address the question that i asked, but thanks all the same. 2.8 What's the best way to upgrade? That's a matter of taste. For most people the best thing to do is to uninstall R (see the previous Q), install the new version, copy any installed pa

Re: [R] Code for "Screenshots" graphics (following on from ease-of-use issues on www.r-project.org)

2006-04-26 Thread Thomas Lumley
On Wed, 26 Apr 2006, Marc Schwartz (via MN) wrote: > On Wed, 2006-04-26 at 11:05 -0700, John McHenry wrote: >> Does anyone know where the code for the graphics on: >> >> http://www.r-project.org/screenshots/screenshots.html >> >> lives? > > > demo(graphics) > > demo(image) > > demo(persp) > >

Re: [R] accessing a SOAP based web service from R

2006-04-26 Thread Seth Falcon
Rajarshi Guha <[EMAIL PROTECTED]> writes: > Hi, > I have a number of web services that run on a server and can be > accessed via SOAP as well as by creating a specific URL. > > Now, for certain services, the argument that I need to send is about 10 > to 20 characters and the URL approach works f

Re: [R] Polygon-like interactive selection of plotted points

2006-04-26 Thread Roger Bivand
polygon and the > x,y coords of the underlying points set to return the points inside or > outside the polygon boundaries. This sequence from the splancs package should do it: library(splancs) set.seed(20060426) xy <- cbind(x=runif(100), y=runif(100)) plot(xy) poly <- getpoly() # this

Re: [R] Code for "Screenshots" graphics (following on from ease-of-use issues on www.r-project.org)

2006-04-26 Thread Marc Schwartz (via MN)
On Wed, 2006-04-26 at 11:05 -0700, John McHenry wrote: > Does anyone know where the code for the graphics on: > > http://www.r-project.org/screenshots/screenshots.html > > lives? demo(graphics) demo(image) demo(persp) These should cover each of the screen shots and then some. I

Re: [R] copying previously installed libraries to R 2.3.0

2006-04-26 Thread Christos Hatzis
See Windows FAQ 2.8 - works well. -Christos -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Thomas Harte Sent: Wednesday, April 26, 2006 2:54 PM To: r-help@stat.math.ethz.ch Subject: [R] copying previously installed libraries to R 2.3.0 hi all, is the

[R] copying previously installed libraries to R 2.3.0

2006-04-26 Thread Thomas Harte
hi all, is there a new mechanism in R 2.3.0 for copying libraries from, say, R 2.2.1 to R 2.3.0? i ask because gabor grothendieck comments in his copydir.bat (from gabor's batchfiles at: http://cran.r-project.org/contrib/extra/batchfiles/batchfiles_0.2-5.zip ): ``:: I personally upgraded my

[R] accessing a SOAP based web service from R

2006-04-26 Thread Rajarshi Guha
Hi, I have a number of web services that run on a server and can be accessed via SOAP as well as by creating a specific URL. Now, for certain services, the argument that I need to send is about 10 to 20 characters and the URL approach works fine. However if the argument is a few kilobytes I woul

Re: [R] Polygon-like interactive selection of plotted points

2006-04-26 Thread roger koenker
?in.convex.hull in the package tripack. url:www.econ.uiuc.edu/~rogerRoger Koenker email[EMAIL PROTECTED]Department of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678Champaign, IL 61820 On Apr 26, 2006

Re: [R] help using tapply

2006-04-26 Thread Gabor Grothendieck
On 4/26/06, Dimitri Szerman <[EMAIL PROTECTED]> wrote: > Dear R-mates, > > # Here's what I am trying to do. I have a dataset like this: > > id = c(rep(1,8), rep(2,8)) > dur1 <- c( 17,18,19,18,24,19,24,24 ) > est1 <- c( rep(1,5), rep(2,3) ) > dur2 <- c(1,1,3,4,8,12,13,14) > est2 <- rep(1,8) > > myda

Re: [R] Polygon-like interactive selection of plotted points

2006-04-26 Thread Marc Schwartz (via MN)
On Wed, 2006-04-26 at 18:13 +0100, Florian Nigsch wrote: > [Please CC me for all replies, since I am not currently subscribed to > the list.] > > Hi all, > > I have the following problem/question: Imagine you have a two- > dimensional plot, and you want to select a number of points, around >

[R] Code for "Screenshots" graphics (following on from ease-of-use issues on www.r-project.org)

2006-04-26 Thread John McHenry
Does anyone know where the code for the graphics on: http://www.r-project.org/screenshots/screenshots.html lives? - [[alternative HTML version deleted]] __ R-help@stat.math.ethz

[R] Generalized SVD

2006-04-26 Thread Philippe Hupé
Hi, I need to computed the GSVD of two matrices. For doing so I used in my C code the lapack routine dggsvd. But when I source my file gsvd.R I have the following error: Error in eval.with.vis(expr, envir, enclos) : BLAS/LAPACK routine 'DGGSVD' gave error code -1 Is there a problem w

Re: [R] Linker problem in installing 64-bit R

2006-04-26 Thread Min Shao
It turned out that I do not have the 64-bit readline library installed. After disabling it using the "-without-readline" option for ./configure, I successfully made a 64-bit build under Solaris 9. Thanks for your help. On 4/20/06, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > > Note: > > > conf

Re: [R] The beautiful math plot

2006-04-26 Thread Anupam Tyagi
Gesmann, Markus lloyds.com> writes: > > See: > ?image I would like to use this image in a publication. How do I get permission from the creater of this image? Anupam Tyagi. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinf

Re: [R] www.r-project.org

2006-04-26 Thread Anupam Tyagi
Gabor Grothendieck gmail.com> writes: > > Maybe a separate web site that shows R off or maybe just > a pointer to the R Graph Gallery. I think the idea of a pointer to the R Graph Gallery is a better one than changing the current R website. I think that the main R web site should be accessible

[R] Polygon-like interactive selection of plotted points

2006-04-26 Thread Florian Nigsch
[Please CC me for all replies, since I am not currently subscribed to the list.] Hi all, I have the following problem/question: Imagine you have a two- dimensional plot, and you want to select a number of points, around which you could draw a polygon. The points of the polygon are defined

Re: [R] A question about is.interger function

2006-04-26 Thread Michael H. Prager
Not to belabor the point, but if you really need to distinguish integers from other numerical constants, you could try something like this: > x <- as.integer(5) > typeof(x) [1] "integer" > is.integer(x) [1] TRUE > on 4/26/2006 1:20 AM Leon said the following: > Hi, All > I am using is.integ

Re: [R] Were to find appropriate functions for a given task in R

2006-04-26 Thread Xiaohua Dai
Hi Albert, Maybe you can summarize what you find during your teaching and put them on one website. Everyone gives a piece of firewood, then we can build up a strong fire. Thanks Xiaohua On 4/26/06, Albert Sorribas <[EMAIL PROTECTED]> wrote: > > Thanks for the indications. All these suggestio

[R] help using tapply

2006-04-26 Thread Dimitri Szerman
Dear R-mates, # Here's what I am trying to do. I have a dataset like this: id = c(rep(1,8), rep(2,8)) dur1 <- c( 17,18,19,18,24,19,24,24 ) est1 <- c( rep(1,5), rep(2,3) ) dur2 <- c(1,1,3,4,8,12,13,14) est2 <- rep(1,8) mydata = data.frame(id, estat=c(est1, est2),

Re: [R] Were to find appropriate functions for a given task in R

2006-04-26 Thread Albert Sorribas
Thanks for the indications. All these suggestion will help a lot. The TINN-R application is impressive. However, I think something more PROBLEM directed is needed. I agree that Google and the ?RSiteSearch helps a lot, but yet the explanations (like this discussion) are not clear enough in many cas

[R] stl function

2006-04-26 Thread Andrea Toreti
Hi, I have a monthly time series with missing values and I would use stl function to identify seasonality. I tried all settings of na.action but the result is the same: stl(tm245,s.window=11, na.action=na.pass) Error in stl(tm245, s.window = 11, na.action = na.pass) : NA/NaN/Inf in for

[R] garch in tseries

2006-04-26 Thread Erin Hodgess
Hello again! Is there a way to include a mean in the garch function in the library(tseries), please? I tried include.mean=T in the function statement but it didn't work thanks in advance! R Version 2.2.1 Windows Sincerely, Erin mailto: [EMAIL PROTECTED] ___

Re: [R] www.r-project.org

2006-04-26 Thread Gabor Grothendieck
Maybe a separate web site that shows R off or maybe just a pointer to the R Graph Gallery. On 4/26/06, Duncan Murdoch <[EMAIL PROTECTED]> wrote: > Romain Francois wrote: > > Dear R users and developpers, > > > > My question is adressed to both of you, so I choose R-help to post it. > > > > Are the

[R] garchFit from fSeries

2006-04-26 Thread Erin Hodgess
Dear R People: I'm trying to use the garchFit function from the library(fSeries) However, R freezes every time that I use it. Is anyone else having this problem, please? Thanks in advance! R Version 2.2.1 Windows. Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathemat

Re: [R] Were to find appropriate functions for a given task in R

2006-04-26 Thread Gabor Grothendieck
There is a reference sheet here: http://www.rpad.org/Rpad/R-refcard.pdf a function finder here: http://biostat.mc.vanderbilt.edu/s/finder/find.html and task views here: http://cran.r-project.org/src/contrib/Views/ Also use of RSiteSearch and help.search from within R can be helpful. On 4

Re: [R] Were to find appropriate functions for a given task in R

2006-04-26 Thread Berton Gunter
1. Check out the R reference card at http://www.rpad.org/Rpad/Rpad-refcard.pdf . There are also several others available from the CRAN website. 2. Check out TINN-R, a Windows text/R code editor that integrates the above and provides function "tips" inline to give the syntax of many R functions. 3

Re: [R] Were to find appropriate functions for a given task in R

2006-04-26 Thread Dimitris Rizopoulos
you may find RSiteSearch() and help.search() functions useful for this purpose, considering usually different key words for the same topic. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 3

[R] Were to find appropriate functions for a given task in R

2006-04-26 Thread Albert Sorribas
This is a generic request concerning were to look for finding appropriate information on a precise procedure in R. I’m using R for teaching introductory statistics and my students are learning how to deal with it. However, I find it difficult to locate some of the procedures. For instance, for basi

Re: [R] [O/T] undergrads and R

2006-04-26 Thread John Fox
Dear Ales, > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of Ales Ziberna > Sent: Wednesday, April 26, 2006 5:21 AM > To: John Fox > Cc: r-help@stat.math.ethz.ch > Subject: Re: [R] [O/T] undergrads and R > > Dear John Fox, > > thank you for your rep

Re: [R] cdf of weibull distribution

2006-04-26 Thread Peter Dalgaard
Sachin J <[EMAIL PROTECTED]> writes: > Hi, > > I have a data set which is assumed to follow weibull distr'. How can I find > of cdf for this data. For example, for normal data I used (package - lmomco) > > >cdfnor(15,parnor(lmom.ub(c(df$V1 > > Also, lmomco package does not

[R] cdf of weibull distribution

2006-04-26 Thread Sachin J
Hi, I have a data set which is assumed to follow weibull distr'. How can I find of cdf for this data. For example, for normal data I used (package - lmomco) >cdfnor(15,parnor(lmom.ub(c(df$V1 Also, lmomco package does not have functions for finding cdf for some of the distri

Re: [R] environment

2006-04-26 Thread Thomas Lumley
On Wed, 26 Apr 2006, Romain Francois wrote: > Hi, > > Consider the code : > > g <- function(){ > print(.x) > .x <- 3 > } > > f <- function(){ > environment(g) <- environment() > .x <- 2 > g() > .x > } > > > f() > [1] 2 > [1] 2 > > > I would like f() to return 3. How can I do that ? Am I compl

Re: [R] A question about is.interger function

2006-04-26 Thread Thomas Lumley
On Wed, 26 Apr 2006, Leon wrote: > Hi, All > I am using is.integer function to examine whether an object is an integer or > not, but I get such results, As the help page says, is.integer() Creates or tests for objects of type '"integer"'. That is, is.integer() tests to see whether the var

Re: [R] lwd - Windows

2006-04-26 Thread Francisco J. Zagmutt
Dear Prof Brian Ripley Thanks for replying to my question. I am sorry but it is still not clear to me if there is a way to handle an "intermediate" line width. In my windows system when I create a figure with 4 plots par(mfrow=c(2,2),lwd=1) the lines look a bit too thin and using lwd=2 makes

[R] installation problem with 2.3.0

2006-04-26 Thread Eric Evans
Hello, I'm in the process of trying to upgrade our R software from version 2.2.3 to 2.3.0 on our Sun workstations. In the installation process I'm encountering some puzzling errors whe I do the 'make install.' After installing the doc files successfully the installation process crashes while

Re: [R] [Rd] R 2.3.0: Use of NULL as an environment is deprecated

2006-04-26 Thread Antonio, Fabio Di Narzo
2006/4/26, Uwe Ligges <[EMAIL PROTECTED]>: > > Pfaff, Bernhard Dr. wrote: > >>"BernPf" == Pfaff, Bernhard Dr <[EMAIL PROTECTED]> > >>on Wed, 26 Apr 2006 08:28:34 +0100 writes: > > > > > > BernPf> Dear Andy, Uwe and remaining list subscribers, > > BernPf> now, coming closer to th

Re: [R] www.r-project.org

2006-04-26 Thread Jacques Veslot
what about improving the "Links" page, with for instance a list of links - such as thoses in task views - grouped by theme, instead of modifying r-project site ? -- --- [EMAIL PROTECTED] CNRS UMR 8090 - http://www-good.ibl.fr Gén

Re: [R] www.r-project.org

2006-04-26 Thread Duncan Murdoch
Romain Francois wrote: > Dear R users and developpers, > > My question is adressed to both of you, so I choose R-help to post it. > > Are there any plans to jazz up the main R website : http://www.r-project.org > The look it have now is the same for a long time and kind of sad > compared to othe

Re: [R] Generalized linear mixed models

2006-04-26 Thread mike waters
Peter, The types of model you are looking to fit are *nonlinear* models. In fact, to be more precise, they are examples of nonlinear ordinal regression models. Check out the nordr package in R. a link to documentation is: http://rss.acs.unt.edu/Rdoc/library/gnlm/html/nordr.html Regards, Mike

Re: [R] Running R on Windows 2000 Terminal Services

2006-04-26 Thread Duncan Murdoch
Barry Rowlingson wrote: > Gavin Simpson wrote: > >>Dear list, >> >>My employer uses a Windows 2000 Terminal Server-based system for its >>college-wide managed computer service - computers connect directly to >>the WTS servers for their sessions, using a Citrix ICA client. When I >>asked them to in

Re: [R] nls and factor

2006-04-26 Thread Christian Ritz
Hi Manuel, an alternative to the approach pointed out by Prof. Ripley is to use the package 'drc' which allows one or more parameters in a non-linear regression model to depend on a factor. You will need the latest version available at www.bioassay.dk (an older version is available on CRAN).

Re: [R] About regression and plot

2006-04-26 Thread Gabor Grothendieck
On 4/26/06, Daniel Yang <[EMAIL PROTECTED]> wrote: > --- I changed to format to text instead of html --- > > Dear R-help, > > This is my first R day. I want to ask some more beginner's questions. Read the posting guide at the bottom of each post to r-help. > > Q1. How can I obtain the covariance

Re: [R] About regression and plot

2006-04-26 Thread Romain Francois
Le 26.04.2006 15:16, Daniel Yang a écrit : > Dear R-help, > > This is my first R day. I want to ask some more beginner's questions. > > Q1. How can I obtain the covariance matrix for parameter estimates of a > multiple regression? I checked ?lm but didn't get the information. > > Q2. How can I s

[R] About regression and plot

2006-04-26 Thread Daniel Yang
--- I changed to format to text instead of html --- Dear R-help, This is my first R day. I want to ask some more beginner's questions. Q1. How can I obtain the covariance matrix for parameter estimates of a multiple regression? I checked ?lm but didn't get the information. Q2. How can I see t

[R] About regression and plot

2006-04-26 Thread Daniel Yang
Dear R-help, This is my first R day. I want to ask some more beginner's questions. Q1. How can I obtain the covariance matrix for parameter estimates of a multiple regression? I checked ?lm but didn't get the information. Q2. How can I see the old graphs in the graph window? Q3. Can R plot a

Re: [R] The beautiful math plot

2006-04-26 Thread Gesmann, Markus
See: ?image Markus Gesmann FPMA Lloyd's Market Analysis Lloyd's * One Lime Street * London * EC3M 7HA Telephone +44 (0)20 7327 6472 Facsimile +44 (0)20 7327 5718 http://www.lloyds.com -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Daniel Yang Sent: 26

[R] The beautiful math plot

2006-04-26 Thread Daniel Yang
Dear R-help, How can I replicate the beautiful math plot found in the right-hand side of http://www.r-project.org/screenshots/desktop.jpg? I tried the following code but didn't obtain something as beautiful. r <- seq(-10, 10, len=100) y <- cos(r^2)*exp(-r/6) par(pty="s") plot(r,y,type="l") T

[R] The beautiful math plot

2006-04-26 Thread Daniel Yang
Dear R-help, How can I replicate the beautiful math plot found in the right-hand side of http://www.r-project.org/screenshots/desktop.jpg? I tried the following code but didn't obtain something as beautiful. r <- seq(-10, 10, len=100) y <- cos(r^2)*exp(-r/6) par(pty="s") plot(r,y,type="l") Tha

Re: [R] Is there a way....

2006-04-26 Thread François Pinard
[Levent TERLEMEZ] >I would like to get rid of counting lines in fix() when i made >a mistake in coding? Is there an easy way to an line numbers to editor? You may configure R to use the editor you want, the way you want. For example, if you want fix() to start Vim in graphical mode, highlighti

Re: [R] Is there a way....

2006-04-26 Thread Pfaff, Bernhard Dr.
Hello Levent, there has recently been a thread on the ESS list by using setnu in EMACS: https://stat.ethz.ch/pipermail/ess-help/2006-April/003450.html Bernhard -Ursprüngliche Nachricht- Von: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Im Auftrag von Uwe Ligges Gesendet: Mittwoch, 26.

Re: [R] environment

2006-04-26 Thread Gabor Grothendieck
Also, if you don;t need to create child objects which override .x, and I don't think you do here, p could be further simplified to this (only the print statement has been changed): p <- proto(.x = 2, g = function(.) { print(.x); .$.x <- 3 }) On 4/26/06, Gabor Grothendieck <[EMAIL PROTECTED]> wro

Re: [R] Is there a way....

2006-04-26 Thread Uwe Ligges
Levent TERLEMEZ wrote: > Hello, > > I would like to get rid of counting lines in fix() when i made a > mistake in coding? Is there an easy way to an line numbers to editor? I think the feature you are requesting is not available. Beside, it is a good idea to use some appropriate editor and

Re: [R] Partially crossed and nested random factors in lme/lmer

2006-04-26 Thread Corsin Müller, Zoologisches Inst.
I do have Pinheiro and Bates 2000 and went through it several times – unfortunately without success for my particular issue. I was not aware of the article in R news or the mlmRev package you referred to, though. I will try with that. The example I gave was not self-contained. Sorry, my fault.

Re: [R] environment

2006-04-26 Thread Gabor Grothendieck
A third possibility is to using the proto package and define a proto object (an environment with special meaning for $) containing the two components .x and g like this: library(proto) p <- proto(.x = 2, g = function(.) { print(.$.x); .$.x <- 3 }) p$.x # 2 print(p$g()) # 2, 3 p$.x # 3 or you ca

Re: [R] new.frame()

2006-04-26 Thread Gabor Grothendieck
Also, the proto and R.oo packages provide object oriented ways of working with environments. On 4/26/06, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > ?new.env > ?local > > should help you. R works with environments, basically a frame plus an > enclosure. > > On Wed, 26 Apr 2006, Anna Whitfield

Re: [R] www.r-project.org

2006-04-26 Thread Romain Francois
Le 26.04.2006 12:14, Prof Brian Ripley a écrit : > The webmasters have not AFAICS replied in this thread, and perhaps you > should wait for them. > > In particular, making cran.r-project.org an autodirector to the > nearest mirror and replacing frames have I believe already been agreed > and jus

Re: [R] lwd - Windows

2006-04-26 Thread Prof Brian Ripley
On Wed, 26 Apr 2006, Francisco J. Zagmutt wrote: > Dear all > > Is there a way or trick in windows to plot a line width that is not an > integer i.e 1.5? > I am aware that the documentation for window devices states "Line widths as > controlled by par(lwd=) are in multiples of the pixel size, and

Re: [R] www.r-project.org

2006-04-26 Thread Barry Rowlingson
Romain Francois wrote: > > What about > - place somewhere a div called 'Focus on a package' where we could have > a short presentation of a package, etc ... or for a CRAN task view (to > do that, php would be great, but we can do I don't know perl scripts to > generate static html pages) > - a

Re: [R] regression modeling

2006-04-26 Thread John Maindonald
An interesting example of this is the forest cover data set that is available from http://www.ics.uci.edu/~mlearn The proportions of the different cover types change systematically as one moves through the file. It seems that distance through the file is a proxy for the geographical co-ordinates.

Re: [R] new.frame()

2006-04-26 Thread Prof Brian Ripley
?new.env ?local should help you. R works with environments, basically a frame plus an enclosure. On Wed, 26 Apr 2006, Anna Whitfield wrote: > Hello, > > I would like to know whether R has a homogeneous function > of S-plus's new.frame(), which create explicit frames in > the evaluator and prov

Re: [R] [O/T] undergrads and R

2006-04-26 Thread Ales Ziberna
Dear John Fox, thank you for your reply. I appologize for not finding the "sort.names" option myself. As for the option of choosing several variables at a time, I would like to see this possibility everywhere where possible. Let me give a few examples: 1. Calculating the mean of several variab

[R] new.frame()

2006-04-26 Thread Anna Whitfield
Hello, I would like to know whether R has a homogeneous function of S-plus's new.frame(), which create explicit frames in the evaluator and provide a locale for computations that can be shared among various functions. new.frame() in S-plus: http://www.uni-muenster.de/ZIV/Mitarbeiter/BennoSueselbe

Re: [R] www.r-project.org

2006-04-26 Thread Prof Brian Ripley
The webmasters have not AFAICS replied in this thread, and perhaps you should wait for them. In particular, making cran.r-project.org an autodirector to the nearest mirror and replacing frames have I believe already been agreed and just need resources to implement (and first www. and cran.r-pro

Re: [R] Partially crossed and nested random factors in lme/lmer

2006-04-26 Thread Spencer Graves
1. Are you familiar with Pinheiro and Bates (2000) Mixed-Effects Models in S and S-Plus (Springer)? I've learned a lot from this book both about mixed models generally and about "lme" in particular -- and even indirectly about "lmer". Examples in this book would, I believe, help yo

Re: [R] Minimum Volume Ellipsoid Estimator ("cov.mve"??)

2006-04-26 Thread Prof Brian Ripley
On Wed, 26 Apr 2006, Öhagen Patrik wrote: Dear All, The Minimum Volume Ellipsoid estimator ("cov.mve") was part of S+ but I cannot find it in R. Where can I find it? help.search("cov.mve") gives cov.rob(MASS) Resistant Estimation of Multivariate Location a

Re: [R] Generalized linear mixed models

2006-04-26 Thread Prof Brian Ripley
On Tue, 25 Apr 2006, Peter Tait wrote: > Hi, > > I would like to fit a generalized linear mixed model (glmm) with a 3 > level response. > > My data is from a longitudinal study with multiple observations/patient > and multiple patients / country. > > Is there an R package that will fit a proportio

Re: [R] environment

2006-04-26 Thread Romain Francois
Le 26.04.2006 11:02, Peter Dalgaard a écrit : > Romain Francois <[EMAIL PROTECTED]> writes: > > >> Hi, >> >> Consider the code : >> >> g <- function(){ >> print(.x) >> .x <- 3 >> } >> >> f <- function(){ >> environment(g) <- environment() >> .x <- 2 >> g() >> .x >> } >> >> > f() >> [

[R] Minimum Volume Ellipsoid Estimator ("cov.mve"??)

2006-04-26 Thread Öhagen Patrik
Dear All, The Minimum Volume Ellipsoid estimator ("cov.mve") was part of S+ but I cannot find it in R. Where can I find it? Why was it excluded? Cheers, Patrik Sweden __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-hel

Re: [R] environment

2006-04-26 Thread Peter Dalgaard
Romain Francois <[EMAIL PROTECTED]> writes: > Hi, > > Consider the code : > > g <- function(){ > print(.x) > .x <- 3 > } > > f <- function(){ > environment(g) <- environment() > .x <- 2 > g() > .x > } > > > f() > [1] 2 > [1] 2 > > > I would like f() to return 3. How can I do that

Re: [R] program and comment

2006-04-26 Thread Hans-Peter
2006/4/26, Guojun Zhu <[EMAIL PROTECTED]>: > Hi. I am a newbie to R. I need to do a courses *do* have a look at e.g.: - http://cran.r-project.org/manuals.html (An Introduction to R) - http://www.ms.unimelb.edu.au/~andrewpr/r-users/icebreakeR.pdf (contributed manual) - http://cran.r-project.org/d

Re: [R] program and comment

2006-04-26 Thread Dirk De Becker
Guojun Zhu wrote: >Hi. I am a newbie to R. I need to do a courses >projects with some manipulation lots of data record >and some back-and-forth linear regression. I do not >want to spend $100 for a SAS lincence, which professor >used. So I decide to start to learn R to finish it. >I am using R

[R] program and comment

2006-04-26 Thread Guojun Zhu
Hi. I am a newbie to R. I need to do a courses projects with some manipulation lots of data record and some back-and-forth linear regression. I do not want to spend $100 for a SAS lincence, which professor used. So I decide to start to learn R to finish it. I am using R in Windows. I feel it i

[R] environment

2006-04-26 Thread Romain Francois
Hi, Consider the code : g <- function(){ print(.x) .x <- 3 } f <- function(){ environment(g) <- environment() .x <- 2 g() .x } > f() [1] 2 [1] 2 I would like f() to return 3. How can I do that ? Am I completely out of place ? Doing that, I want to avoid to pass .x as a parameter

  1   2   >