RE: [R] ANOVA with both discreet and continuous variable

2005-04-22 Thread Bill.Venables
An anonymous enquirer asks: : -Original Message- : From: [EMAIL PROTECTED] : [mailto:[EMAIL PROTECTED] On Behalf Of array chip : Sent: Saturday, 23 April 2005 2:32 PM : To: r-help@stat.math.ethz.ch : Subject: [R] ANOVA with both discreet and continuous variable : : : Hi all, : : I have

[R] ANOVA with both discreet and continuous variable

2005-04-22 Thread array chip
Hi all, I have dataset with 2 independent variable, one (x1) is continuous, the other (x2) is a categorical variable with 2 levels. The dependent variable (y) is continuous. When I run linear regression y~x1*x2, I found that the p value for the continuous independent variable x1 changes when diffe

Re: [R] Problem with R-2.1.0: install.packages() doesn't work

2005-04-22 Thread Marc Schwartz
Hi all, In follow up to my message earlier, I spent some time applying the modification to the menu() function that I referenced earlier. I also recompiled a local copy of R on my FC3 system using: ./configure --without-tcltk This results in: > capabilities() jpeg pngtcltk X11

RE: [R] Anova - interpretation of the interaction term

2005-04-22 Thread Bill.Venables
: -Original Message- : From: [EMAIL PROTECTED] : [mailto:[EMAIL PROTECTED] On Behalf Of : michael watson (IAH-C) : Sent: Friday, 22 April 2005 7:47 PM : To: r-help@stat.math.ethz.ch : Subject: [R] Anova - interpretation of the interaction term : : : Hi : : So carrying on my use of an

[R] Cyrillic R: how can I help? (or Re: R 2.1.0, Windows, non-Latin locales)

2005-04-22 Thread iit512
Dear R-experts, I am using R in Russia since 2000. One can find some Russian documentation for R (including the reference card) on: http://herba.msu.ru/shipunov/software/r/r-en.htm All your criticism would be appreciated. > We really do need help from the users in the beta-test period, > espe

Re: [R] Expression in panel.text

2005-04-22 Thread Deepayan Sarkar
On Friday 22 April 2005 10:12, David Orme wrote: > Hi, > > I've got a lattice xyplot and I want to superimpose correlation > coefficients and p values on each panel. I've been trying to get this > to work using something of the form: > > x <- rnorm(400) > y <- rnorm(400) > a <- gl(4, 100) > > xyplo

Re: [R] Problem with R-2.1.0: install.packages() doesn't work

2005-04-22 Thread Roger D. Peng
What happens if you don't set the 'CRAN' option via 'options()' first and just run 'install.packages()'? -roger Waichler, Scott R wrote: I installed R-2.1.0 from source on a Linux box running Red Hat Enterprise Linux WS release 4 but install.packages() wouldn't work (see below). When I install R

[R] Expression in panel.text

2005-04-22 Thread David Orme
Hi, I've got a lattice xyplot and I want to superimpose correlation coefficients and p values on each panel. I've been trying to get this to work using something of the form: x <- rnorm(400) y <- rnorm(400) a <- gl(4, 100) xyplot(y~x | a, panel=function(x,y, ...){

Re: [R] Need help arranging the plot in different fashion than the default format

2005-04-22 Thread Marc Schwartz
On Fri, 2005-04-22 at 16:19 -0500, Ghosh, Sandeep wrote: > I want to change the way the plot is appearing in the barchart .. > > testdata <- as.data.frame(t(structure(c > (1,2004,"LV1",3.8,2,87,2,2004,"LV1",3.2,3,28,3,2004,"LV1",3.4,3,88,4,2004, >"LV1",3,2,26,5,2004,"LV1",3.8,2,87,6,2004,"LV1",3.2

Re: [R] hi: I am newbie for R too...

2005-04-22 Thread ronggui
i think the log-linear model can deal with this. and the glm with poisson family,loglin,loglm from MASS packages can estimate log-linear model. hope help:) On Fri, 22 Apr 2005 22:15:23 +0800 hau chung chan <[EMAIL PROTECTED]> wrote: > I have so much questions about using R...but i will try my b

RE: [R] Tool for update

2005-04-22 Thread Andy Bunn
Hi Ronaldo: First, this script was discussed a few days ago. Note that under 2.1 update.packages() has changed for the better. Look at the NEWS file. So, if you want to keep updating on Tuesdays, despite sensible advice to the contrary offered in this thread http://finzi.psych.upenn.edu/R/Rhelp02

RE: [R] if(foo == TRUE) .. etc

2005-04-22 Thread Thomas Lumley
On Fri, 22 Apr 2005, bogdan romocea wrote: Great suggestion; it made me change all my Ts/Fs to TRUE/FALSE. Given F <- TRUE T <- FALSE is it possible to forbid T to stand for TRUE, and F for FALSE in function(...,something=T)? Or, alternatively, never allow F <- whatever and T <- whatever? All

Re: [R] Problem with R-2.1.0: install.packages() doesn't work

2005-04-22 Thread Marc Schwartz
On Fri, 2005-04-22 at 18:00 -0400, Manuel Morales wrote: > On Fri, 2005-04-22 at 15:42 -0500, Marc Schwartz wrote: > > On Fri, 2005-04-22 at 15:44 -0400, Manuel Morales wrote: > > > On Fri, 2005-04-22 at 17:48 +0200, Uwe Ligges wrote: > > > > Waichler, Scott R wrote: > > > > > > > > > I installed

Re: [R] help with POSIX

2005-04-22 Thread Peter Dalgaard
"Ghosh, Sandeep" <[EMAIL PROTECTED]> writes: > For the r script below > > >datestr <- "01/01/2004" > >as.POSIXct(as.Date(datestr, "%d/%m/%Y")) > I get the following output > "2003-12-31 18:00:00 Central Standard Time" > > Why is the date a day before. I guess its something to do with the > tim

[R] Hoaglin Outlier Method

2005-04-22 Thread khobson
I am a new user of R so please bear with me. I have reviewed some R books, FAQs and such but the volume of material is great. I am in the process of porting my current SAS and SVS Script code to Lotus Approach, R and WordPerfect. My question is, can you help me determine the best R method to

[R] How to override coerion error in 'scan'

2005-04-22 Thread Kittler, Richard
I am using 'read.csv' in V2.0.1 to read in a CSV file with the colClasses option and am getting an error from 'scan' when it encounters a non-numeric value for a 'numeric' column, i.e. > ds <- read.csv(in_file, nrows=irow, row.names=NULL, colClasses=zclass, comment.char="")

Re: [R] when can we expect Prof Tierney's compiled R?

2005-04-22 Thread Luke Tierney
On Wed, 20 Apr 2005, Peter Dalgaard wrote: Luke Tierney <[EMAIL PROTECTED]> writes: Vectorized operations in R are also as fast as compiled C (because that is what they are :-)). A compiler such as the one I'm working on will be able to make most difference for non-vectorizable or not very vectori

[R] density

2005-04-22 Thread Hui Han
Hi, I used the density function in the R package, and got the following results. I just wonder how to explain them. What is Min, 1st Qu, Median, and so on? I could not find an explanation from help(density). The plot doesn't seem to match the x and y value either. Thanks in advance for any help

Re: [R] Tool for update

2005-04-22 Thread Ronaldo Reis-Jr.
Em Sex 25 Mar 2005 11:16, abunn escreveu: > I edited Rprofile to update everything on Tuesdays. I've been doing this > since 2.0 and I think I've had R running almost every Tuesday, which begs > the question of what I would be doing if R hadn't come into existence. > > In any case, It works pretty

[R] algorithm used in k-mean clustering

2005-04-22 Thread Asha Jayanthi
Hi, I have used the kmean fucntion in R to produce some results for my analysis. I like to know the specific underlying algorithm used for the implementation of the function kmean in R. I tried looking for some documents but could not find any. I obtained the kmean result for k ranging from 2 to

Re: [R] if(foo == TRUE) .. etc

2005-04-22 Thread Vincent Goulet
Le 22 Avril 2005 13:41, Liaw, Andy a écrit : > > From: bogdan romocea > > > > Great suggestion; it made me change all my Ts/Fs to TRUE/FALSE. > > Given > >    F <- TRUE > >    T <- FALSE > > is it possible to forbid T to stand for TRUE, and F for FALSE in > >    function(...,something=T)? > > Or, a

[R] defining custom survreg distributions

2005-04-22 Thread Michael Townsley
Dear R Helpers, I am trying to carry out survival analysis in the presence of long term survivors (immunes/cureds). This involves using a split population model (some call it a mixture model) where the assumption of eventual failure is relaxed. I am following closely the formulation by Malle

[R] Alternative to aggregate() for averaging

2005-04-22 Thread Sebastien Durand
Hello, I am working with very large matrices 1 by 2000 and for each row I have a grouping vector and I would like to average each row by groups... Presently I am using this command and it behaves correctly even with NA values is the groupvector: data<-aggregate(mymatrix,list(groupvector),me

[R] Problem with R-2.1.0: install.packages() doesn't work

2005-04-22 Thread Waichler, Scott R
I installed R-2.1.0 from source on a Linux box running Red Hat Enterprise Linux WS release 4 but install.packages() wouldn't work (see below). When I install R-2.0.1 from RPM on the same system, everything is fine. Version 2.1.0 (2005-04-18), ISBN 3-900051-07-0 . . . > options(CRAN = "http://c

Re: [R] algorithm used in k-mean clustering

2005-04-22 Thread Gavin Simpson
Asha Jayanthi wrote: Hi, I have used the kmean fucntion in R to produce some results for my analysis. I like to know the specific underlying algorithm used for the implementation of the function kmean in R. I tried looking for some documents but could not find any. I obtained the kmean result

[R] help with POSIX

2005-04-22 Thread Ghosh, Sandeep
For the r script below >datestr <- "01/01/2004" >as.POSIXct(as.Date(datestr, "%d/%m/%Y")) I get the following output "2003-12-31 18:00:00 Central Standard Time" Why is the date a day before. I guess its something to do with the time, but is there a way to get it to return 2004-01-01 instead?

[R] Problems with "help.start()"

2005-04-22 Thread Kenneth Roy Cabrera Torres
Hi R users and developers: I just install the new R version 2.1.0 in a linux platform. I get this error when I call the function help.start() Making links in per-session dir ... Error in gsub(pattern, replacement, x, ignore.case, extended, fixed) : input string 28 is invalid in this locale

[R] Need help arranging the plot in different fashion than the default format

2005-04-22 Thread Ghosh, Sandeep
I want to change the way the plot is appearing in the barchart .. testdata <- as.data.frame(t(structure(c(1,2004,"LV1",3.8,2,87,2,2004,"LV1",3.2,3,28,3,2004,"LV1",3.4,3,88,4,2004,"LV1",3,2,26,5,2004,"LV1",3.8,2,87,6,2004,"LV1",3.2,3,28,7,2004,"LV1",3.4,3,88,8,2004,"LV1",3,2,26,9,2004,"LV1",3.8,2,

RE: [R] if(foo == TRUE) .. etc

2005-04-22 Thread Liaw, Andy
> From: bogdan romocea > > Great suggestion; it made me change all my Ts/Fs to TRUE/FALSE. > Given >F <- TRUE >T <- FALSE > is it possible to forbid T to stand for TRUE, and F for FALSE in >function(...,something=T)? > Or, alternatively, never allow F <- whatever and T <- whatever?

[R] dr ()

2005-04-22 Thread Jessica Higgs
Hi all-- A quick question about the dr () function. I am using this function to reduce the dimensions of a data set I have that involves 14 predictor variables and one predictant or response. The goal is to discover which variables play the most important role in determining the response and, t

[R] shared library configuration; Gnome GUI

2005-04-22 Thread Paul Johnson
Hello, everybody: On a Fedora Core 3 Linux system, I built R-2.1 using an updated version of the spec file that was used to make the RPMs for version 2.0.1 on the CRAN system. The build was fine, and packages updates perfectly. Thanks! Then I got curious about the package gnomeGUI. While tryin

Re: [R] Help needed with lattice graph!

2005-04-22 Thread Paul Murrell
Hi Deepayan Sarkar wrote: On Friday 22 April 2005 03:29, Sander Oom wrote: Dear R users, If I manage to sort out this graph, it is certainly a candidate for the new R graph gallery (http://addictedtor.free.fr/graphiques/displayGallery.php)! I created the following lattice graph: library(lattice) tm

RE: [R] Hoaglin Outlier Method

2005-04-22 Thread Liaw, Andy
That looks just like how `outliers' are determined in boxplots. You can use the output of boxplot.stats() to compute the limits. [EDA purists would tell you that those shound be letter values (or `F' for fourths), not quartiles.] HTH, Andy > -Original Message- > From: [EMAIL PROTECTED]

Re: [R] help with POSIX

2005-04-22 Thread Thomas Lumley
On Fri, 22 Apr 2005, Ghosh, Sandeep wrote: For the r script below datestr <- "01/01/2004" as.POSIXct(as.Date(datestr, "%d/%m/%Y")) I get the following output "2003-12-31 18:00:00 Central Standard Time" Why is the date a day before. I guess its something to do with the time, but is there a way to g

Re: [R] Installing packages from source code

2005-04-22 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: Hi everybody, I have trouble in installing packages from source code by following Section 5.1 in manual R-admin.pdf . I am using R 2.1.0 and Win NT. Following the Windows toolset section in the manual, I download the tool set package from: http://www.murdoch

Re: [R] R 2.1.0 for Windows installation error? atanh not in R.dll?

2005-04-22 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: Could someone please tell me what I did wrong to create this message or what I should do to correct this problem? I downloaded 2.1.0 Windows binary and installed into C:/R/rw2010, using the installer. I ran md5check.exe in C:/R/rw2010/bin/ and got "No errors." The problem

Re: [R] Help needed with lattice graph!

2005-04-22 Thread Deepayan Sarkar
On Friday 22 April 2005 03:29, Sander Oom wrote: > Dear R users, > > If I manage to sort out this graph, it is certainly a candidate for the > new R graph gallery > (http://addictedtor.free.fr/graphiques/displayGallery.php)! > > I created the following lattice graph: > > library(lattice) > tmp <- e

[R] Re: an interesting qqnorm question

2005-04-22 Thread WeiWei Shi
hope it is not b/c some central limit therory, otherwise my initial plan will fail :) On 4/22/05, WeiWei Shi <[EMAIL PROTECTED]> wrote: > Hi, r-gurus: > > I happened to have a question in my work: > > I have a dataset, which has only one dimention, like > 0.99037297527605 > 0.991179836732708 > 0

RE: [R] Problem with R-2.1.0: install.packages() doesn't work

2005-04-22 Thread Liaw, Andy
> From: Uwe Ligges > > Waichler, Scott R wrote: > > > I installed R-2.1.0 from source on a Linux box running Red Hat > > Enterprise Linux WS release 4 but install.packages() > wouldn't work (see > > below). When I install R-2.0.1 from RPM on the same > system, everything > > is fine. > > > >

Re: [R] Problem with R-2.1.0: install.packages() doesn't work

2005-04-22 Thread Manuel Morales
On Fri, 2005-04-22 at 17:48 +0200, Uwe Ligges wrote: > Waichler, Scott R wrote: > > > I installed R-2.1.0 from source on a Linux box running Red Hat > > Enterprise Linux WS release 4 but install.packages() wouldn't work (see > > below). > >>install.packages("rgenoud") > > > > --- Please select a

Re: [R] Problem with R-2.1.0: install.packages() doesn't work

2005-04-22 Thread Manuel Morales
On Fri, 2005-04-22 at 15:42 -0500, Marc Schwartz wrote: > On Fri, 2005-04-22 at 15:44 -0400, Manuel Morales wrote: > > On Fri, 2005-04-22 at 17:48 +0200, Uwe Ligges wrote: > > > Waichler, Scott R wrote: > > > > > > > I installed R-2.1.0 from source on a Linux box running Red Hat > > > > Enterprise

Re: [R] Error when downloading and installing ALL R packages

2005-04-22 Thread Uwe Ligges
Bernd Weiss wrote: Hi, after updating to 2.1 (see below) I am no longer able to install all R packages as mentioned at . After finishing the download, I received the following error: [...] trying URL 'http://www.stats.ox.ac.uk/pub/RWin/bin/w

Re: [R] Help needed with lattice graph!

2005-04-22 Thread Romain Francois
Well done ! The result of your work has been uploaded to the gallery (graph 48). Romain Le 22.04.2005 22:39, Paul Murrell a écrit : Hi Deepayan Sarkar wrote: On Friday 22 April 2005 03:29, Sander Oom wrote: Dear R users, If I manage to sort out this graph, it is certainly a candidate for the new R

RE: [R] if(foo == TRUE) .. etc

2005-04-22 Thread bogdan romocea
Great suggestion; it made me change all my Ts/Fs to TRUE/FALSE. Given F <- TRUE T <- FALSE is it possible to forbid T to stand for TRUE, and F for FALSE in function(...,something=T)? Or, alternatively, never allow F <- whatever and T <- whatever? I don't know what the technical side is

[R] an interesting qqnorm question

2005-04-22 Thread WeiWei Shi
Hi, r-gurus: I happened to have a question in my work: I have a dataset, which has only one dimention, like 0.99037297527605 0.991179836732708 0.995635340631367 0.997186769599305 0.991632565640424 0.984047197106486 0.99225943762649 1.00555642128421 0.993725402926564 the data is saved in a f

Re: [R] Problem with R-2.1.0: install.packages() doesn't work

2005-04-22 Thread Uwe Ligges
Roger D. Peng wrote: What happens if you don't set the 'CRAN' option via 'options()' first and just run 'install.packages()'? Should be the same as long as no CRAN mirror has been choosen and X11 is not accessible (see my other message). The code needs to be a bit more defensive here - or presen

Re: [R] Problem with R-2.1.0: install.packages() doesn't work

2005-04-22 Thread Marc Schwartz
On Fri, 2005-04-22 at 15:44 -0400, Manuel Morales wrote: > On Fri, 2005-04-22 at 17:48 +0200, Uwe Ligges wrote: > > Waichler, Scott R wrote: > > > > > I installed R-2.1.0 from source on a Linux box running Red Hat > > > Enterprise Linux WS release 4 but install.packages() wouldn't work (see > > >

Re: [R] Problem with R-2.1.0: install.packages() doesn't work

2005-04-22 Thread Peter Dalgaard
Manuel Morales <[EMAIL PROTECTED]> writes: > On Fri, 2005-04-22 at 17:48 +0200, Uwe Ligges wrote: > > Waichler, Scott R wrote: > > > > > I installed R-2.1.0 from source on a Linux box running Red Hat > > > Enterprise Linux WS release 4 but install.packages() wouldn't work (see > > > below). > >

Re: [R] Problem with R-2.1.0: install.packages() doesn't work

2005-04-22 Thread Uwe Ligges
Waichler, Scott R wrote: I installed R-2.1.0 from source on a Linux box running Red Hat Enterprise Linux WS release 4 but install.packages() wouldn't work (see below). When I install R-2.0.1 from RPM on the same system, everything is fine. Version 2.1.0 (2005-04-18), ISBN 3-900051-07-0 . . . op

Re: [R] Factor Analysis functions...

2005-04-22 Thread Pierre BADY
hi all, In the library ‘ade4’, there are two eigenanalysis which enable the ordination of the categorical variables. 1- Multiple Correspondence Analysis (MCA, Tenenhaus & Young 1985) performs the multiple correspondence analysis of a factor table (see the function ‘dudi.acm’). this function is eq

Re: [R] Installing packages from source on WindowsXP

2005-04-22 Thread Marc Schwartz
On Fri, 2005-04-22 at 09:46 -0400, Tyler Smith wrote: > Ok, I uninstalled my copy of HHW (and deleted the extra file from the > Rtools/bin folder), and reinstalled from the murdoch-sutherland site. > From what I read there I think it's the same version, and I actually > get a warning that I alr

[R] hi: I am newbie for R too...

2005-04-22 Thread hau chung chan
I have so much questions about using R...but i will try my best to solve myself before bothering you guys firstly...yes no maleold 13070 young34

Re: [R] Installing packages from source on WindowsXP

2005-04-22 Thread Uwe Ligges
Tyler Smith wrote: Ok, I uninstalled my copy of HHW (and deleted the extra file from the Rtools/bin folder), and reinstalled from the murdoch-sutherland site. From what I read there I think it's the same version, and I actually get a warning that I already have a newer version installed on my

Re: [R] ugly loop

2005-04-22 Thread Bill Simpson
On Fri, 22 Apr 2005, Marc Schwartz wrote: > Thus, I just need to use t(ver) instead of ver: > > > x <- pt$x[t(ver)] > > y <- pt$y[t(ver)] > > z <- pt$z[t(ver)] > That should do it? Yep!! Thanks very much Marc and others who suggested this. Cheers Bill __

[R] density estimation

2005-04-22 Thread Bernard Palagos
hello sorry for my english I would like estimate density for multivariate variable,( f(x,y) , f(x,y ,z) for example) ; for calculate mutual information how is posible with R? thanks Bernard Bernard Palagos Unité Mixte de Recherche Cemagref - Agro.M - CIRAD Information et Technologie pour les Ag

[R] Installing packages from source on WindowsXP

2005-04-22 Thread Tyler Smith
Ok, I uninstalled my copy of HHW (and deleted the extra file from the Rtools/bin folder), and reinstalled from the murdoch-sutherland site. From what I read there I think it's the same version, and I actually get a warning that I already have a newer version installed on my machine (1.3.1) anyw

Re: [R] ugly loop

2005-04-22 Thread Marc Schwartz
On Fri, 2005-04-22 at 09:31 -0400, Bill Simpson wrote: > Thanks Marc for your help. > > > > The following code is slow and ugly: > > > > > > count<-0 > > > for(i in 1:nrow(ver)) > > > for(j in 1:ncol(ver)) > > > { > > > count<-count+1 > > > x[count]<-pt$x[ver[i,j]] > > > y[count

Re: [R] ugly loop

2005-04-22 Thread Woodrow Setzer
Almost there; you need the transpose of v, since Bill originally had columns changing faster: e.g. x <- pt$x[t(ver)] -Original Message- From: Marc Schwartz <[EMAIL PROTECTED]> Sent: Apr 22, 2005 9:17 AM To: Bill Simpson <[EMAIL PROTECTED]> Cc: R-Help Subject: Re: [R] ugly loop On Fri, 2

Re: [R] ugly loop

2005-04-22 Thread Bill Simpson
Thanks Marc for your help. > > The following code is slow and ugly: > > > > count<-0 > > for(i in 1:nrow(ver)) > > for(j in 1:ncol(ver)) > > { > > count<-count+1 > > x[count]<-pt$x[ver[i,j]] > > y[count]<-pt$y[ver[i,j]] > > z[count]<-pt$z[ver[i,j]] > > } > > > > Please

Re: [R] ugly loop

2005-04-22 Thread Dimitris Rizopoulos
then use x <- c(t(ver)) Best, --D Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/16/336899 Fax: +32/16/337015 Web: http://www.med.kuleuven.ac.be/biostat/ http://www.student.k

RE: [R] lsfit result - how to compute t-values for coefficients

2005-04-22 Thread Liaw, Andy
You can try something like: fit <- lsfit(X, array(Y, dim = c(t.length, prod(d.dim[2:4] rsq <- sapply(ls.print(fit, print.it=FALSE)$coef.table, function(x) x[-1,c(1, 3)]) dim(rsq) <- c(ncol(X), 2, dim(rsq)[2]) dimnames(rsq) <- list(NULL, c("Estimate", "t-value"), prod(d.dim[2:4])

Re: [R] ugly loop

2005-04-22 Thread Bill Simpson
To clarify: I want to get rid of the loop over i,j Here is a simpler example. ver is a 2D matrix count<-0 for(i in 1:nrow(ver)) for(j in 1:ncol(ver)) { count<-count+1 x[count]<-ver[i,j] } Bill __ R-help@stat.math.ethz.ch mailing list

Re: [R] ugly loop

2005-04-22 Thread Marc Schwartz
On Fri, 2005-04-22 at 08:58 -0400, Bill Simpson wrote: > The following code is slow and ugly: > > count<-0 > for(i in 1:nrow(ver)) > for(j in 1:ncol(ver)) > { > count<-count+1 > x[count]<-pt$x[ver[i,j]] > y[count]<-pt$y[ver[i,j]] > z[count]<-pt$z[ver[i,j]] > } > > Please

Re: [R] ugly loop

2005-04-22 Thread Dimitris Rizopoulos
maybe you want something like this: x <- c(t(pt$x)) Best, --D Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/16/336899 Fax: +32/16/337015 Web: http://www.med.kuleuven.ac.be/biost

[R] ugly loop

2005-04-22 Thread Bill Simpson
The following code is slow and ugly: count<-0 for(i in 1:nrow(ver)) for(j in 1:ncol(ver)) { count<-count+1 x[count]<-pt$x[ver[i,j]] y[count]<-pt$y[ver[i,j]] z[count]<-pt$z[ver[i,j]] } Please help me make it better. Thanks! Bill

Re: [R] Error when downloading and installing ALL R packages

2005-04-22 Thread roger bos
I have no idea how to solve that error, but I use the following bit of code (which I picked up from someone else on this list) to update ALL packages on CRAN and it works on 2.1.0: x <- packageStatus(repositories="http://cran.r-project.org/src/contrib";) st <- x$avai["Status"] install.packages(row

[R] lsfit result - how to compute t-values for coefficients

2005-04-22 Thread Christoph Lehmann
Hi I used lsfit instead of lm since I have a huge Y data-set (X being constant for all Y). Since I require the t-values for all coefficients: which would be the fastest way to compute them, eg for the example: ## using lsfit with a matrix response: t.length <- 5 d.dim <- c(t.length,7,8,9) # dim

[R] fBasics question: Get dates corresponding to maximum values

2005-04-22 Thread Neuro LeSuperHéros
fBasics question: Get maximum dates Hello, This two series data set has been created with the timeSeries() function from fBasics str(total) Formal class 'timeSeries' [package "fBasics"] with 7 slots ..@ Data : num [1:262, 1:2] 8703 8603 8573 8680 8668 ... .. ..- attr(*, "dimnames")=List

RE: [R] Assign factor and levels inside function

2005-04-22 Thread Tim Howard
Aha! You've just opened the door to another level for this blundering R user. I even went back to my well-used copy of "An Introduction to R" to see where I missed this standard approach for processing new data. Nothing clear but certainly alluded to in many of the function examples. I don't

RE: [R] Using R to illustrate the Central Limit Theorem

2005-04-22 Thread Ted Harding
On 21-Apr-05 [EMAIL PROTECTED] wrote: > Here's a bit of a refinement on Ted's first suggestion. > [ corrected from runif(M*k), N, k) to runif(N*k), N, k) ] > > N <- 1 > graphics.off() > par(mfrow = c(1,2), pty = "s") > for(k in 1:20) { > m <- (rowMeans(matrix(runif(N*k), N, k)) - 0.5)

Re: [R] How to know if a classification tree is predicitve or not?

2005-04-22 Thread Uwe Ligges
Laure Maton wrote: Hello, I would like to know how to know if a classification tree is predictive or not ? Is it sufficient to analyse results of cross validation? Depends on your interpretation of "predictive", maybe you want to look at stuff like sensitivity and specificity as well. Uwe Ligges

RE: [R] Anova - interpretation of the interaction term

2005-04-22 Thread Ted Harding
On 22-Apr-05 michael watson \(IAH-C\) wrote: > [...] > So from a statistical point of view, I can see exactly why the > interaction term is significant, but what my colleage wants to know is > that WITHIN the vaccinated group, does infection decrease expression > significantly? And within the unva

Re: [R] Installing packages from source on WindowsXP

2005-04-22 Thread Uwe Ligges
Tyler Smith wrote: Hi, I'm having some problems installing packages from the source files on Windows, using R CMD INSTALL pkg. I'm running WindowsXP, and I've followed the instructions as per the README.packages file from the R installation. I ran into a hitch, with the install failing following

Re: [R] pointer to comments re Paul Murrell's new book, R, & SAS on Andrew Gelman's blog

2005-04-22 Thread Peter Dalgaard
Douglas Bates <[EMAIL PROTECTED]> writes: > Tony Plate wrote: > > There are some interesting comments re Paul Murrell's new book, R, & > > SAS on Andrew Gelman's blog: > > http://www.stat.columbia.edu/~cook/movabletype/archives/2005/04/a_new_book_on_r.html > > I found the comments quite interesti

RE: [R] Required Packages etiquette

2005-04-22 Thread Liaw, Andy
[This is probably more appropriate for R-devel rather than R-help.] IMHO the criterion for requiring a package (say `X') is that, if some function(s) in your package can not perform its tasks correctly without something in `X', then you should list `X' as required. You can not assume the user wil

[R] Anova - interpretation of the interaction term

2005-04-22 Thread michael watson \(IAH-C\)
Hi So carrying on my use of analysis of variance to check for the effects of two factors. It's made simpler by the fact that both my factors have only two levels each, creating four unique groups. I have a highly significant interaction term. In the context of the experiment, this makes sense.

Re: [R] How to tell if R is running in batch mode

2005-04-22 Thread Peter Dalgaard
"Khabie-Zeitoune, David" <[EMAIL PROTECTED]> writes: > Is there a way to programmatically tell whether R is running in batch or > GUI mode? interactive() should do (it's not distinguishing between GUI and terminal mode though). -- O__ Peter Dalgaard Blegdamsvej 3 c/ /'_

Re: [R] How to tell if R is running in batch mode

2005-04-22 Thread Duncan Murdoch
Khabie-Zeitoune, David wrote: Hi Is there a way to programmatically tell whether R is running in batch or GUI mode? ?interactive __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! htt

Re: [R] Beginner in R

2005-04-22 Thread Peter Dalgaard
Norbert Billet <[EMAIL PROTECTED]> writes: > hello ( and sorry for my poor english ... ) > > I'm a newbie on R software and I need to obtain this kind of system : > > a structure, like a liste : > > my_struct <- list() > my_struct$a <- a_value > my_struct$b <- another_value > my_struct$c <- one

[R] How to tell if R is running in batch mode

2005-04-22 Thread Khabie-Zeitoune, David
Hi Is there a way to programmatically tell whether R is running in batch or GUI mode? Thanks Dave ___ David Khabie-Zeitoune Quantitative Arbitrage Brevan Howard Asset Management direct: +44 (0)20 7022 6167 mobile:+44 (0)7

[R] Intership: R programmer, London

2005-04-22 Thread Vincent_Zoonekynd
State Street Global Advisors (one of the world's leading assets management companies)'s alternative strategies group is looking for interns able to bring a significant contribution to our R&D projects. Here are some of the directions we are currently investigating: - Digital Signal Processing

[R] Help needed with lattice graph!

2005-04-22 Thread Sander Oom
Dear R users, If I manage to sort out this graph, it is certainly a candidate for the new R graph gallery (http://addictedtor.free.fr/graphiques/displayGallery.php)! I created the following lattice graph: library(lattice) tmp <- expand.grid(geology = c("Sand","Clay","Silt","Rock"), species =

Re: [R] pointer to comments re Paul Murrell's new book, R, & SAS on Andrew Gelman's blog

2005-04-22 Thread Douglas Bates
Tony Plate wrote: There are some interesting comments re Paul Murrell's new book, R, & SAS on Andrew Gelman's blog: http://www.stat.columbia.edu/~cook/movabletype/archives/2005/04/a_new_book_on_r.html I found the comments quite interesting. I also find it somewhat amusing when people speculate

Re: [R] lme4: apparently different results between 0.8-2 and 0.95-6

2005-04-22 Thread Douglas Bates
Steve Cumming wrote: I've been using lme4 to fit Poisson GLMMs with crossed random effects. The data are counts(y) sampled at 55 sites over 4 (n=12) or 5 (n=43) years. Most models use three fixed effects: x1 is a two level factor; x2 and x3 are continuous. We are including random intercepts for YEA

Re: [R] pointer to comments re Paul Murrell's new book, R, & SAS on Andrew Gelman's blog

2005-04-22 Thread Campbell
Brief comment on blog. S Plus is mentioned only in passing. I've always thought comparisons between SAS and R are invalid as they do different things. If I wasn't using R I'd be coding in either C or Fortran. SAS appears to be a database with statistical functionality attached. I've always th

[R] Beginner in R

2005-04-22 Thread Norbert Billet
hello ( and sorry for my poor english ... ) I'm a newbie on R software and I need to obtain this kind of system : a structure, like a liste : my_struct <- list() my_struct$a <- a_value my_struct$b <- another_value my_struct$c <- one_more_value and a function with two args : the first is a instance

[R] Error when downloading and installing ALL R packages

2005-04-22 Thread Bernd Weiss
Hi, after updating to 2.1 (see below) I am no longer able to install all R packages as mentioned at . After finishing the download, I received the following error: [...] trying URL 'http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2