[R] One way ANOVA

2008-08-09 Thread Angelo Scozzarella
Hi, There is a command or a function for One way ANOVA if I know the numbers of subjects for each group, the mean and di standard deviation for each group? Thanks Angelo Scozzarella __ R-help@r-project.org mailing list https://stat.ethz.ch/mai

Re: [R] Reading large datasets and fitting logistic models in R

2008-08-09 Thread Prof Brian Ripley
See also bigglm() in package biglm. On Sat, 9 Aug 2008, Pradheep K E wrote: Hi R-experts, Does anyone have experience using R for handling large scale data (millions of rows, hundreds or thousands of features)? What is the largest size of data that anyone has used with glm? I've used 700,00

Re: [R] instal tar.gz package on windows

2008-08-09 Thread Prof Brian Ripley
In apparently a late response to https://stat.ethz.ch/pipermail/r-help/2008-July/169578.html on Sat, 9 Aug 2008, Andrew Harris wrote: For cross-building a mature package, I stumbled across a convenient utility provided by the R Windows maintainer: _http://win-builder.r-project.org/ It is in

[R] mean-plot (add residuals)

2008-08-09 Thread Jörg Groß
Hi, I want to plot the mean of a variable and add the standard deviation as a line going above and below the mean (like the whiskers in a boxplot) but I don't know how to add these residual-lines. Is there an easy way to do that or do I have to use the line()- function to create them on my

Re: [R] instal tar.gz package on windows

2008-08-09 Thread Andrew Harris
For cross-building a mature package, I stumbled across a convenient utility provided by the R Windows maintainer: _http://win-builder.r-project.org/ _This takes a little while, so for iteration it is much better to install the Rtools, but for a package or two it is great. Best, Andrew Harris,

[R] Reading large datasets and fitting logistic models in R

2008-08-09 Thread Pradheep K E
Hi R-experts, Does anyone have experience using R for handling large scale data (millions of rows, hundreds or thousands of features)? What is the largest size of data that anyone has used with glm? Also, is there a library to read data in sparse data format (like SVMlight format)? Thanks Pradh

[R] Scripting - query

2008-08-09 Thread Gareth Campbell
I have a vector: alleles.present<-c("D3", "D16", ... ) The alleles present changes given the case I'm dealing with - i.e. either all of the alleles I use for my calculations are present, or some of them. Depending on what alleles are present, I need to make matrices and do calculations on those a

Re: [R] Converting nested "for" loops to an "apply" function(s)

2008-08-09 Thread Kurt Newman
Resending. Previous message was truncated. Sorry for possible confusion. > From: [EMAIL PROTECTED] > To: r-help@r-project.org > Date: Sat, 9 Aug 2008 18:25:47 -0400 > Subject: [R] Converting nested "for" loops to an "apply" function(s) > > > Hello, >

Re: [R] box-plot without a box surrounding the plot

2008-08-09 Thread Jörg Groß
par(bty="l") boxplot(count ~ spray, data = InsectSprays, col = "lightgray") #is that what you want? Stephen Thanks! That's what I needed. On Sat, Aug 9, 2008 at 8:12 PM, Jörg Groß <[EMAIL PROTECTED]> wrote: Hi, I want to draw a boxplot and I want to get rid of the surrounding box.

[R] possible problem with rgdal

2008-08-09 Thread Naiara Pinto
Hello all, I recently installed rgdal 0.5.24-1 (kyngchaos framework) and I am having trouble making sense of the row, col information provided. > a = new("GDALDataset", "dummy.tif") > GDALinfo("dummy.tif") rows420 columns 660 bands 1 ll.x-55.5 ll.y-14.5 res.x

Re: [R] box-plot without a box surrounding the plot

2008-08-09 Thread stephen sefick
par(bty="l") boxplot(count ~ spray, data = InsectSprays, col = "lightgray") #is that what you want? Stephen On Sat, Aug 9, 2008 at 8:12 PM, Jörg Groß <[EMAIL PROTECTED]> wrote: > Hi, > > I want to draw a boxplot and I want to get rid of the surrounding box. > So that there are only the axis lines

[R] box-plot without a box surrounding the plot

2008-08-09 Thread Jörg Groß
Hi, I want to draw a boxplot and I want to get rid of the surrounding box. So that there are only the axis lines left on the bottom and on the left. I tried a litte bit with box() but I dont get it the way I want. Can somebody help me out? __ R-he

Re: [R] How to repress the annoying complains from X window system

2008-08-09 Thread dusa.adrian
Dear Marc, Marc Schwartz wrote: > > On Fri, 2007-05-18 at 11:25 -0400, Hao Liu wrote: >> Dear All: >> >> I am running some GUI functions in linux environment, they runs fine, >> however I constantly get this kind of message in R console: >> >> Warning: X11 protocol error: BadWindow (invalid

Re: [R] effective matrix subset

2008-08-09 Thread jgarcia
It seems that this solution provided by Dan (and also available in SPoetry; I'm sorry I didn't notice it) is the fastest and simplest. I was using a more standard approach: V <- t(A)[(0:(nrow(A)-1))*ncol(A)+X], That wasn't bad, but I was confident that you, R gurus, could outperform this. This is

Re: [R] how to interpret t.test output

2008-08-09 Thread Peter Dalgaard
Felipe Carrillo wrote: Hi Ted: Thanks for your prompt reply and explanation. That's what I was wondering, why would one need to test mu=0 ,which is the t.test default. But reading Peter Dalgaard's book and looking at some examples online, I saw t.test being used like that; t.test(datasetname) wi

[R] Converting nested "for" loops to an "apply" function(s)

2008-08-09 Thread Kurt Newman
Hello, I would like to know more about how to use the "apply" family and have attempted to convert nested "for" loops in example code from Contributed Documentation ("The Friendly Beginners' R Course” by Toby Marthews (ZIP, 2007-03-01)") to an "apply" function(s). The relevant code is: dista

[R] how to interpret t.test output

2008-08-09 Thread Felipe Carrillo
I take it back, Peter Dalgaard's book uses t.test with mu=7725 and no mu=0. I got the script online. Hi Ted: Thanks for your prompt reply and explanation. That's what I was wondering, why would one need to test mu=0 ,which is the t.test default. But reading Peter Dalgaard's book and looking at

Re: [R] how to interpret t.test output

2008-08-09 Thread Felipe Carrillo
Hi Ted: Thanks for your prompt reply and explanation. That's what I was wondering, why would one need to test mu=0 ,which is the t.test default. But reading Peter Dalgaard's book and looking at some examples online, I saw t.test being used like that; t.test(datasetname) with no other arguments.

Re: [R] Multivariate Regression with Weights

2008-08-09 Thread Arne Henningsen
On Saturday 09 August 2008 23:25, Arne Henningsen wrote: > > -Original Message- > > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] > > Sent: Monday, August 04, 2008 5:15 PM > > To: Zhang Yanwei - Princeton-MRAm > > Subject: RE: [R] Multivariate Regression with Weights > > > > the syste

Re: [R] how to interpret t.test output

2008-08-09 Thread Ted Harding
On 09-Aug-08 20:31:33, Felipe Carrillo wrote: ># Hi all: > #I got a vector with fish lengths(mm) > # Can someone help me interpret the output of > # a t.test in plain english? > # Based on the t.test below I can say that > # I reject the null hypothesis because > # th

Re: [R] Multivariate Regression with Weights

2008-08-09 Thread Arne Henningsen
> -Original Message- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] > Sent: Monday, August 04, 2008 5:15 PM > To: Zhang Yanwei - Princeton-MRAm > Subject: RE: [R] Multivariate Regression with Weights > > the systemfit package can do that and the documentation is quite nice > also. T

[R] how to interpret t.test output

2008-08-09 Thread Felipe Carrillo
# Hi all: #I got a vector with fish lengths(mm) # Can someone help me interpret the output of # a t.test in plain english? # Based on the t.test below I can say that # I reject the null hypothesis because # the p-value is smaller than the the significance # leve

Re: [R] Box.test degrees of freedom

2008-08-09 Thread David Stoffer
I believe tsdiag() uses the correct degrees of freedom in applying Box.test, but the graphic shows "lag" on the horizontal axis when it should display "degrees of freedom". raf.rossignol wrote: > > Hello, > > Prof Brian Ripley wrote: >> >> I think you are referring to its application to t

Re: [R] RPro

2008-08-09 Thread Marc Schwartz
on 08/09/2008 12:45 PM Marc Schwartz wrote: on 08/09/2008 12:13 PM Charles C. Berry wrote: On Fri, 8 Aug 2008, Kenn Konstabel wrote: There's more to this trend: SPSS and Statistica now advertise "R language support" : http://www.statsoft.com/industries/Rlanguage.htm http://www.spss.com/spssd

Re: [R] RPro

2008-08-09 Thread Marc Schwartz
on 08/09/2008 12:13 PM Charles C. Berry wrote: On Fri, 8 Aug 2008, Kenn Konstabel wrote: There's more to this trend: SPSS and Statistica now advertise "R language support" : http://www.statsoft.com/industries/Rlanguage.htm http://www.spss.com/spssdirections/na/sessions.cfm?sessionType=2 If

Re: [R] RPro

2008-08-09 Thread Charles C. Berry
On Fri, 8 Aug 2008, Kenn Konstabel wrote: There's more to this trend: SPSS and Statistica now advertise "R language support" : http://www.statsoft.com/industries/Rlanguage.htm http://www.spss.com/spssdirections/na/sessions.cfm?sessionType=2 If you can't beat R, join R. --- Should someone s

[R] Reshape set operations?

2008-08-09 Thread rkevinburton
I have mange to use the library reshape to give me data structures that I want. Specifically: m2008 <- melt(t2008, id.var=c("DayOfYear","Category","SubCategory","Sku"), measure.var=c("Quantity")) m2007 <- melt(t2007, id.var=c("DayOfYear","Category","SubCategory","Sku"), measure.var=c("Quantity

Re: [R] levels values of cut()

2008-08-09 Thread Prof Brian Ripley
On Sat, 9 Aug 2008, baptiste auguie wrote: Thank you all for the precious tips. For memory I've made the following wrapper function for this. I wonder whether a short note on these regular expressions could be useful on the help page of cut(). Already there in R-devel cutIntervals <

Re: [R] Color of box frame in Legend (Was: Matrix barplot)

2008-08-09 Thread Andreas Tille
On Mon, 28 Jul 2008, Nutter, Benjamin wrote: Try sourcing in the 'new.legend' function below. It's the legend function with a new argument called 'box.col'. The argument will change the color of the box surrounding the legend. If I understand what it is you are looking for, this should work.

Re: [R] effective matrix subset

2008-08-09 Thread Marc Schwartz
on 08/09/2008 06:52 AM Dan Davison wrote: On Sat, Aug 09, 2008 at 06:29:59AM -0500, Marc Schwartz wrote: on 08/09/2008 06:01 AM [EMAIL PROTECTED] wrote: Hi; If we have a matrix A, and a vector X, where length(X)=nrow(A), and X contains a wanted column for each row in A, in row ascending order.

Re: [R] levels values of cut()

2008-08-09 Thread baptiste auguie
Thank you all for the precious tips. For memory I've made the following wrapper function for this. I wonder whether a short note on these regular expressions could be useful on the help page of cut(). cutIntervals <- function(x, ...){ dotArgs <- unlist(c(...)) if( any(names(dotArgs

Re: [R] Tick marks that correspond with bars on barplot

2008-08-09 Thread Jim Lemon
On Fri, 2008-08-08 at 16:04 -0300, Megan J Bellamy wrote: > Hello all, > > I have created a barplot that shows change in hardwood/softwood density from > 1965 to 2005 in 5 year periods (1965,1970, etc). I would like to have an > X-axis where the labels for the years line up after every two bars

Re: [R] [lme4]Coef output with binomial lmer

2008-08-09 Thread Mark Difford
Hi Tom, >> 1|ass%in%pop%in%fam This is "non-standard," but as you have found, it works. The correct translation is in fact 1|fam/pop/ass and not 1|ass/pop/fam as suggested by Harold Doran. Dropping %, ass%in%pop%in%fam reads [means] as: nest ass in pop [= pop/ass], and then nest this in fam ==

Re: [R] effective matrix subset

2008-08-09 Thread Dan Davison
On Sat, Aug 09, 2008 at 06:29:59AM -0500, Marc Schwartz wrote: > on 08/09/2008 06:01 AM [EMAIL PROTECTED] wrote: >> Hi; >> If we have a matrix A, and a vector X, where length(X)=nrow(A), and X >> contains a wanted column for each row in A, in row ascending order. How >> would be the most effective

Re: [R] effective matrix subset

2008-08-09 Thread Marc Schwartz
on 08/09/2008 06:01 AM [EMAIL PROTECTED] wrote: Hi; If we have a matrix A, and a vector X, where length(X)=nrow(A), and X contains a wanted column for each row in A, in row ascending order. How would be the most effective way to extract the desired vector V (with length(V)=nrow(A))? A <- matri

[R] effective matrix subset

2008-08-09 Thread jgarcia
Hi; If we have a matrix A, and a vector X, where length(X)=nrow(A), and X contains a wanted column for each row in A, in row ascending order. How would be the most effective way to extract the desired vector V (with length(V)=nrow(A))? Wishes, Javier __

Re: [R] levels values of cut()

2008-08-09 Thread Prof Brian Ripley
On Sat, 9 Aug 2008, baptiste auguie wrote: Dear list, I have the following example, from which I am hoping to retrieve numeric values of the factor levels (that is, without the brackets): x <- seq(1, 15, length=100) y <- sin(x) my.cuts <- cut(which(abs(y) < 1e-1), 3) levels(my.cuts) his

Re: [R] RPro

2008-08-09 Thread Ajay ohri
I think a bit of standardization and commercial exposure would help R, given the fact it is much much better than any other expensive stats package yet has no gotten the market share it deserves because everyone is either building it or using it, hence no one has actually taken the effort to promot

Re: [R] re cursive root finding

2008-08-09 Thread baptiste auguie
Hi all, with a lot more fiddling around, I've come up with this solution: x <- seq(1, 15, length=100) y <- jitter(cos(x), a=0.2) plot(x, y) findZero <- function(x, y, guess = mean(range(x)), exclusion= 5*diff(range(x))/length(x)) { interval <- c(guess - exclusion, guess+exc

Re: [R] levels values of cut()

2008-08-09 Thread Stephen Tucker
Not sure what you're looking for, but does this help? Extending your code, > library(gsubfn) > t(strapply(levels(my.cuts),"([0-9.]+),([0-9.]+)", + function(...) as.numeric(c(...)),backref=-2,simplify=TRUE)) [,1] [,2] [1,] 15.9 38.3 [2,] 38.3 60.7 [3,] 60.7 83.1 - Original Mes

[R] levels values of cut()

2008-08-09 Thread baptiste auguie
Dear list, I have the following example, from which I am hoping to retrieve numeric values of the factor levels (that is, without the brackets): x <- seq(1, 15, length=100) y <- sin(x) my.cuts <- cut(which(abs(y) < 1e-1), 3) levels(my.cuts) hist() does not suit me for this, as it does n

Re: [R] read.table question

2008-08-09 Thread Dan Davison
On Fri, Aug 08, 2008 at 07:27:13PM -0700, Alessandro wrote: > Hi All. > > > > I have a file txt with 3 columns (X, Y and Z). every rows has 4 decimal > place (i.e. x.). I use read.table to import the data in R, but with > summary(), I don't see the decimal place after the dot. Is there any

Re: [R] import/export txt file

2008-08-09 Thread Dan Davison
On Fri, Aug 08, 2008 at 04:44:13PM -0700, Alessandro wrote: > Hi All, > > > > I have 2 questions: > > 1. Import: when I import my txt file (X,Y and Z) in R with "testground > <- read.table(file="c:/work_LIDAR_USA/R_kriging/ground26841492694149.txt", > header=T)", I lost the 4 number afte

Re: [R] Lattice: regression lines within grouped xyplot panels

2008-08-09 Thread Rainer Hurling
On 09.08.2008 01:05 (UTC+1), Deepayan Sarkar wrote: On Fri, Aug 8, 2008 at 2:38 PM, Rainer Hurling <[EMAIL PROTECTED]> wrote: Dear community, I am looking for a possibility to draw 'regression lines' instead of 'smooth' lines in grouped xyplots. The following code should give you a small exampl