[R] finding paired values from common subjects -- vector operation?

2007-10-20 Thread gallon li
I have two vectors for values collected from a group of subjects, say a=c(100,200,150,120,140,180) b=c(200,300,420,130) I also have two vectors which indicate the corresponding subjects for a and b, say for a, the subjects are suba=c(1,2,3,4,5,6) for b, the subjects are subb=c(1,3,5,6) Then

[R] building a copy of R with certain packages in Linux

2007-10-20 Thread Edna Bell
Dear R Gurus: I would like to have a version of R with certain packages only to be put together in Linux. Does something like that exist, please? I'm thinking it might be a sort of R CMD build but I don't know how to proceed. Thanks for any help. Sincerely, Edna Bell

Re: [R] Secondary Y axis title

2007-10-20 Thread Gabor Grothendieck
mtext does not support srt -- use text instead. There is a comment about that in the example although perhaps the wording could be clearer. On 10/20/07, Dave Hewitt <[EMAIL PROTECTED]> wrote: > The controls for the margin sizes, which set the amount of space allocated > for the labels, are the 'm

Re: [R] formatting a list

2007-10-20 Thread jim holtman
another choice is: x <- scan('temp.txt', what=c(rep(list(0), 19))) On 10/20/07, Tomas Vaisar <[EMAIL PROTECTED]> wrote: > Hi, > > I am new to R and need to read in a file with 19 columns and 7000 rows > and make it into a list of 7000 lists with 19 items each. For a > simpler case of 10 by 10 ta

[R] mixed model with factor

2007-10-20 Thread Irene Mantzouni
Dear all I am trying to fit a mixed model with a factor and a random effect on a slope: y~F*x+...,random=~x where F is a factor with 2 levels and x the covariate. the random effects for the 2 levels of F should be equal so I am fitting the model like: ex.lme=lme(y~x+F+z*x+F*x-1, random=list(g

Re: [R] Secondary Y axis title

2007-10-20 Thread Dave Hewitt
The controls for the margin sizes, which set the amount of space allocated for the labels, are the 'mai' numbers in the first 'par' call. The fourth number changed to 0.9 or greater should give you more room like you're looking for. You can use the alternative 'mtext' way of adding the label (i

Re: [R] Secondary Y axis title

2007-10-20 Thread Gabor Grothendieck
In the zoo example you can control how far it lies from the axis by adjusting the text() command and you can alter the margins using the values for the mai parameter in the par() command. On 10/20/07, John Theal <[EMAIL PROTECTED]> wrote: > I still haven't resolved this problem. However, when I r

Re: [R] formatting a list

2007-10-20 Thread Charles C. Berry
On Sat, 20 Oct 2007, Tomas Vaisar wrote: > Hi, > > I am new to R and need to read in a file with 19 columns and 7000 rows > and make it into a list of 7000 lists with 19 items each. For a > simpler case of 10 by 10 table I used x <-scan("file", > list(0,0,0,0,0,0,0,0,0,0)), perhaps clumsy, but it

Re: [R] Secondary Y axis title

2007-10-20 Thread John Theal
I still haven't resolved this problem. However, when I run the plot.zoo example the label appears, it also appears if I enter the commands manually. The problem that persists is that the label has to lie so close to the secondary y axis, that (in my case) it becomes almost indistinguishable f

[R] formatting a list

2007-10-20 Thread Tomas Vaisar
Hi, I am new to R and need to read in a file with 19 columns and 7000 rows and make it into a list of 7000 lists with 19 items each. For a simpler case of 10 by 10 table I used x <-scan("file", list(0,0,0,0,0,0,0,0,0,0)), perhaps clumsy, but it did the job. However with the large 19x7000 (which n

Re: [R] mussy data

2007-10-20 Thread Gavin Simpson
On Sat, 2007-10-20 at 18:15 +0700, Sandra Stattman wrote: > Hi, I am writing my B.SC thesis and I have to do something in R, my data are > > iindexidx > 1110.2 > 2210.4 > 3310.3 > 4120.5 > 5220.3 > 6320.1 > 71

Re: [R] Secondary Y axis title

2007-10-20 Thread Gabor Grothendieck
That's strange. When I run example(plot.zoo) its also cut off but if I run the very same code by copying and pasting it into my session then its not cut off: > packageDescription("zoo")$Version [1] "1.4-0" > R.version.string [1] "R version 2.6.0 Patched (2007-10-08 r43124)" library(zoo) examp

Re: [R] beginner's tutorial, books, etc re: time-series analysis, ARMA/ARIMA models...

2007-10-20 Thread Gabor Grothendieck
There are two vignettes in the zoo package that discuss basic manipulations of time series. library(zoo) vignette("zoo") vignette("zoo-quickref") On 10/19/07, Thomas Pujol <[EMAIL PROTECTED]> wrote: > Can anyone point me towards a beginner's tutorial, books, etc for someone > that on time-serie

[R] saturation binding in nlme

2007-10-20 Thread David Airey
To estimate saturation binding parameters Bmax and Kd in a receptor saturation binding experiment, I use the following nonlinear equation and the nls() function: bmax*X*dummy + ns*X + background = total binding kd+X where X is concentration, and dummy is an indicator to allow sh

[R] beginner's tutorial, books, etc re: time-series analysis, ARMA/ARIMA models...

2007-10-20 Thread Thomas Pujol
Can anyone point me towards a beginner's tutorial, books, etc for someone that on time-series analysis, ARMA/ARIMA models, etc? keywords: time-series, autoregression, autoregressive, tutorial, beginner, ARMA, ARIMA __ [[alternative HT

Re: [R] Secondary Y axis title

2007-10-20 Thread Dave Hewitt
I was curious about the exact same question that John Theal posed - how to get a second y-axis label for a plot of two data series against a common x (in my case, time). I tried two different methods - one in lattice and one in plot. Both times the y-axis on the right side appears but the label

Re: [R] plot.Design

2007-10-20 Thread Frank E Harrell Jr
John Smith wrote: > Dear R-users: > > I am trying to use the following code to reproduce the figures on page 140 > of Prof. Frank Harrell's book 'Regression Modeling Strategies': > > rm(list=ls()) > options(width=128) > library(Hmisc) > library(Design) > > getHdata(counties) > counties$older <-

[R] mussy data

2007-10-20 Thread Sandra Stattman
Hi, I am writing my B.SC thesis and I have to do something in R, my data are iindexidx 1110.2 2210.4 3310.3 4120.5 5220.3 6320.1 7130.2 8230.1 9330.4 and for the analysis they h

Re: [R] modeling large data

2007-10-20 Thread Benilton Carvalho
And it's likely you want to combine biglm with a database approach, either by using your own dbs (RSQLite or equivalent) or using tools like the 'sqlitedf' pkg. B On Oct 20, 2007, at 12:26 PM, Benilton Carvalho <[EMAIL PROTECTED]> wrote: > And I suggested to use the 'biglm' package. Didn't

Re: [R] modeling large data

2007-10-20 Thread Benilton Carvalho
And I suggested to use the 'biglm' package. Didn't it work? It has examples that I found useful. B On Oct 20, 2007, at 11:41 AM, "Wensui Liu" <[EMAIL PROTECTED]> wrote: > Hi, Dear Listers, > Several days ago, I posted a question regarding modeling large dataset > in R. Could anyone with such e

Re: [R] unable to interactively label curves on a plot

2007-10-20 Thread Frank E Harrell Jr
Quin Wills wrote: > Hello all > > Simple question for the gurus. I'm trying to interactively label curves > on a single plot. The labcurve() function from Hmisc seems like the way > to do this (?). I just can't seem to get it to work. Toy example: > > > x <- 1:10 > > y1 <- x^2 > > y2 <- 2*x

Re: [R] r achives

2007-10-20 Thread Peter Dalgaard
John Sorkin wrote: > Might I suggest that the Email deamon include the URL of the archives to the > footer it attaches to Email messages. > Johh > > Er, it is one click from what is in there already. Is it really not possible for people to guess that the archives for a mailing list might b

Re: [R] path to libgfortran 'hardcoded' in R?

2007-10-20 Thread Charles C. Berry
On Sat, 20 Oct 2007, Rainer Hurling wrote: > I am using R-2.6.0 on FreeBSD 8.0-CURRENT (i386). In the last days I had > problems when building packages SparseM, lme4 and randomForest. > > The below message shows for randomForest, that 'libgfortran' was not > found. The same error appeared with Spa

[R] pairs, par("plt")

2007-10-20 Thread Oliver Soong
I'm having some confusion over the coordinate system after using pairs. I'm not interested in the content of the actual pairs plot, although the number of pairs seems to matter a bit. I'm purely interested in knowing where my points will be plotted on the device. However, after using pairs, the p

Re: [R] modeling large data

2007-10-20 Thread Wensui Liu
Hi, Dear Listers, Several days ago, I posted a question regarding modeling large dataset in R. Could anyone with such experience shed some light on it? I truly appreciate it. wensui On 10/17/07, Wensui Liu <[EMAIL PROTECTED]> wrote: > Hi, Dear Listers, > I am just curious if R is able to model a

Re: [R] r achives [follow-up]

2007-10-20 Thread Ted Harding
On 20-Oct-07 15:23:25, Ted Harding wrote: > On 20-Oct-07 13:53:11, Alan Zaslavsky wrote: >> I think that should be >> >> https://stat.ethz.ch/pipermail/r-help/ >> >> (note https prefix) > > Well, I thought so too: My "bookmark" has https. > But, in replying to Raymond Chiruka, I wondered wh

Re: [R] r achives

2007-10-20 Thread Ted Harding
On 20-Oct-07 13:53:11, Alan Zaslavsky wrote: > I think that should be > > https://stat.ethz.ch/pipermail/r-help/ > > (note https prefix) Well, I thought so too: My "bookmark" has https. But, in replying to Raymond Chiruka, I wondered whether the https was necessary. So I tried http, and it

Re: [R] r achives

2007-10-20 Thread John Sorkin
Might I suggest that the Email deamon include the URL of the archives to the footer it attaches to Email messages. Johh John Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics Baltimore VA Medical Center GRECC, University of Maryland School of Medicine Claude D. Pepper OAIC, University of Ma

Re: [R] Getting at what a named object represents in a function...

2007-10-20 Thread Prof Brian Ripley
Because of lazy evaluation, your function can find the expression for 'param' by deparse(substitute(param)), provided that is called before 'param' is actually used. However, although that would allow you to change 'list' in the calling environment, that is regarded as really bad programming p

Re: [R] extending zoo (or ts) objects

2007-10-20 Thread Gabor Grothendieck
On 10/20/07, Felix Andrews <[EMAIL PROTECTED]> wrote: > list(...) > > I would like to extend zoo objects for use in a hydrological > application. By that I mean I want to attach various attributes, such > as id numbers, coordinates of the observation location, a units > specification, etc. I might

Re: [R] r achives

2007-10-20 Thread Alan Zaslavsky
I think that should be https://stat.ethz.ch/pipermail/r-help/ (note https prefix) > From: (Ted Harding) <[EMAIL PROTECTED]> > Subject: Re: [R] r achives > > On 19-Oct-07 09:36:52, raymond chiruka wrote: > > sorry but how do i accsess r archives > > If you mean the archive of postings to t

[R] Getting at what a named object represents in a function...

2007-10-20 Thread Martin Waller
Hi, I'm pretty new to R. I have an object (say a list) and I I have a function that I call on various columns in that list (excuse terminology if it's wrong/ambiguous). Imagine its like this (actual values are unimportant) and called mylist: >mylist AB 15 25 36 48 50

Re: [R] please help me

2007-10-20 Thread Ben Bolker
azadeh sadeghian wrote: > > dear list > I am student M.S. statistics in department statistics . I am working in > the function "nls" in the [R 2.3.1] with 246 data and want to fit a model > to vectors( v and u ) but I have > a problem to use it > u > > > PS -- I just looked back at

Re: [R] please help me

2007-10-20 Thread Ben Bolker
azadeh sadeghian wrote: > > dear list > I am student M.S. statistics in department statistics . I am working in > the function "nls" in the [R 2.3.1] with 246 data and want to fit a model > to vectors( v and u ) but I have > a problem to use it > > Thanks for providing a reproducible e

[R] path to libgfortran 'hardcoded' in R?

2007-10-20 Thread Rainer Hurling
I am using R-2.6.0 on FreeBSD 8.0-CURRENT (i386). In the last days I had problems when building packages SparseM, lme4 and randomForest. The below message shows for randomForest, that 'libgfortran' was not found. The same error appeared with SparseM and lme4. -

[R] extending zoo (or ts) objects

2007-10-20 Thread Felix Andrews
list(...) I would like to extend zoo objects for use in a hydrological application. By that I mean I want to attach various attributes, such as id numbers, coordinates of the observation location, a units specification, etc. I might also want to add to the "class" attribute (inheriting from "zoo")

[R] How to plot with txt or symbol instead of number as x-axis scale ?

2007-10-20 Thread Yong Wang
Dear R-list My question is exactly as the subject, I did search the archieve but failed to find relevant message, in more details: The simplest example, 4 factors and 4 value like below: 1020 30 40 A BC D if I plot the value, the x-axis will be scaled as 1,2,3,4 w

[R] please help me

2007-10-20 Thread azadeh sadeghian
dear list I am student M.S. statistics in department statistics . I am working in the function "nls" in the [R 2.3.1] with 246 data and want to fit a model to vectors( v and u ) but I have a problem to use it u 5.00e-13 2.179057e+03 6.537171e+03 1.089529e+04 1.525340e+04 1.961151e+04

[R] please help me

2007-10-20 Thread azadeh sadeghian
dear list I am student M.S. statistics in department statistics . I am working in the function "nls" in the [R 2.3.1] with 246 data and want to fit a model to vectors( v and u ) but I have a problem to use it u 5.00e-13 2.179057e+03 6.537171e+03 1.089529e+04 1.525340e+04 1.961151e+04