Re: [R] time function behavior for ts class objects

2009-06-14 Thread Prof Brian Ripley
Try not reducing options(digit=)! I get (with the default value 7) time(xts) Jan Feb Mar Apr May Jun Jul Aug 1973 1973.000 1973.083 1973.167 1973.250 1973.333 1973.417 1973.500 1973.583 1974 1974.000 1974.083 1974.167 1974.250 1974.333 1974.417

Re: [R] 'Errors' with Ubuntu

2009-06-14 Thread Paul Hiemstra
Hi, As Uwe said, this is not really the place to get support for Ubuntu. They have excellent forums. But the problem seems to be that GRUB was not correctly installed. Run an ubuntu live cd and reinstall GRUB from that. The command might be grub-install, but i'm not sure cheers, Paul

[R] a proposal regarding documentation

2009-06-14 Thread Patrick Burns
Proposal That a new mailing list be established that pertains exclusively to R documentation. The purpose of the list would be to discuss weak sections of the documentation and establish fixes for those weak spots. Pro If it works, there would be better documentation. It would be an

Re: [R] a proposal regarding documentation

2009-06-14 Thread Peter Flom
Patrick Burns pbu...@pburns.seanet.com wrote Proposal That a new mailing list be established that pertains exclusively to R documentation. The purpose of the list would be to discuss weak sections of the documentation and establish fixes for those weak spots. Pro If it works, there would be

Re: [R] a proposal regarding documentation

2009-06-14 Thread John Sorkin
Issues (1) Will people use the new listserver? (2) Will developers respond to postings Question (1) Are there any guidelines for creating documentation? If not should they be developed? It seems to me that every help page should include (a) examples, and when approprate (b) a reference to an

Re: [R] a proposal regarding documentation

2009-06-14 Thread stephen sefick
I have a package StreamMetabolism. I believe that documentation is the toughest part. I find it to be straight foward, but then I wrote the package. Lets try one? I don't m ind helping. Stephen Sefick On Sun, Jun 14, 2009 at 7:23 AM, John Sorkinjsor...@grecc.umaryland.edu wrote: Issues

[R] boxplot size

2009-06-14 Thread luscinia
Hey everybody, I am trying to generate multiple boxplots on one page by using the layout () function. This works really well but the plots are always badly proportioned (3 times wider than high) and much too small to be easily readable. I have tried reducing margins to allow for bigger plots,

Re: [R] Gtk objects disappears

2009-06-14 Thread Olivier Nuñez
It woks fine. Thanks very much. Olivier -- Olivier G. Nuñez Email: onu...@iberstat.es Tel : +34 663 03 69 09 Web: http://www.iberstat.es El 14/06/2009, a las 4:30, Michael Lawrence escribió: I would recommend

Re: [R] Problem in installation of Rgraphviz package

2009-06-14 Thread Takeshi N
Hi,there, I could install Rgraphviz in R version 2.6. I could not install this in R version 2.7 or 2.8. Please try in version2.6. ram basnet wrote: Dear R users,   I am not so used to this R software. I have to use the package Rgraphviz but found some problem in the installation

Re: [R] boxplot size

2009-06-14 Thread David Winsemius
Have you actually looked at the help page for layout()? -- David Winsemius On Jun 14, 2009, at 7:35 AM, luscinia wrote: Hey everybody, I am trying to generate multiple boxplots on one page by using the layout () function. This works really well but the plots are always badly

Re: [R] a proposal regarding documentation

2009-06-14 Thread Patrick Burns
I neglected to say that I don't have the time to be one of the drivers for the proposed list. Pat Patrick Burns wrote: Proposal That a new mailing list be established that pertains exclusively to R documentation. The purpose of the list would be to discuss weak sections of the documentation

Re: [R] a proposal regarding documentation

2009-06-14 Thread Patrick Burns
I had only envisioned the scope of the list to be the packages that R-core are responsible for. Expanding to contributed packages would expand both potential usefulness and certain complexity exponentially. Pat stephen sefick wrote: I have a package StreamMetabolism. I believe that

Re: [R] a proposal regarding documentation

2009-06-14 Thread stephen sefick
That make sense to just do the R-core documentation. I was suggesting using my package because I understand it and could have a dialog about the design choices for the documentation, which would also improve the package itself. I am starting my PhD program tomorrow and will probably not have

Re: [R] a proposal regarding documentation

2009-06-14 Thread spencerg
The documentation for the nlme package was improved a few years ago by an informal process of this nature. When I first started following r-help, I answered many questions suggesting in part the the person read some portion of Pinheiro and Bates (2000). At that time, none of the help

Re: [R] a proposal regarding documentation

2009-06-14 Thread John Sorkin
Perhaps help pages should have links to relevant portions of the WIKI. John John David Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics University of Maryland School of Medicine Division of Gerontology Baltimore VA Medical Center 10 North Greene Street GRECC (BT/18/GR) Baltimore, MD

[R] Centering graphics by default in Sweave

2009-06-14 Thread Frank E Harrell Jr
Is there an elegant way to get Sweave to center graphics by default? I'd like to use \centerline{\includegraphics{}} etc. to save some vertical space that \begin{center} ... \end{center} uses, and I'd like to avoid centering with each fig=T= Thanks Frank -- Frank E Harrell Jr Professor and

Re: [R] a proposal regarding documentation

2009-06-14 Thread Peter Flom
Maybe once the list is going, people who are willing to submit their packages to such scrutiny and critique could join in? Peter -Original Message- From: Patrick Burns pbu...@pburns.seanet.com Sent: Jun 14, 2009 12:30 PM To: stephen sefick ssef...@gmail.com Cc: r-h...@stat.math.ethz.ch,

[R] learning about panel functions in lattice

2009-06-14 Thread Juliet Hannah
Hi All, I am trying to understand panel functions. Let's use this example. library(lattice) time-c(rep(1:10,5)) y -time+rnorm(50,5,2) group-c(rep('A',30),rep('B',20)) subject-c(rep('a',10),rep('b',10),rep('c',10),rep('d',10),rep('e',10)) myData -data.frame(subject,group,time,y) head(myData)

Re: [R] learning about panel functions in lattice

2009-06-14 Thread Gabor Grothendieck
2009/6/14 Juliet Hannah juliet.han...@gmail.com: Hi All, I am trying to understand panel functions. Let's use this example. library(lattice) time-c(rep(1:10,5)) y -time+rnorm(50,5,2) group-c(rep('A',30),rep('B',20)) subject-c(rep('a',10),rep('b',10),rep('c',10),rep('d',10),rep('e',10))

[R] if + is.na

2009-06-14 Thread Grześ
Hello! I wont to use a function is.na() I have two vectors: a=c(1,NA,3,3,3) b=c(0,0,0,0,0) and when I use is.na function it's ok: is.na(a) [1] FALSE TRUE FALSE FALSE FALSE but I would create sth like this: for i in 1:length(a){ if (wsp[i] == is.na(a)) {b=43} } or like this

Re: [R] if + is.na

2009-06-14 Thread Ben Bolker
Grześ wrote: Hello! I wont to use a function is.na() I have two vectors: a=c(1,NA,3,3,3) b=c(0,0,0,0,0) and when I use is.na function it's ok: is.na(a) [1] FALSE TRUE FALSE FALSE FALSE but I would create sth like this: for i in 1:length(a){ if (wsp[i] == is.na(a))

Re: [R] if + is.na

2009-06-14 Thread Rolf Turner
On 15/06/2009, at 7:48 AM, Grześ wrote: Hello! I wont to use a function is.na() I have two vectors: a=c(1,NA,3,3,3) b=c(0,0,0,0,0) and when I use is.na function it's ok: is.na(a) [1] FALSE TRUE FALSE FALSE FALSE but I would create sth like this: for i in 1:length(a){ if (wsp[i]

Re: [R] External signal in ODE written in C (using deSolve and approx1?)

2009-06-14 Thread Thomas Petzoldt
Glenn Woodart wrote: Dear Thomas Thank you for your excellent help. I agree method 3 would be the best one, and this is the one I have been trying to get working (like you described). So far my skills with pointers, what is returned etc in pure C code is limited, so my attempts have failed so

Re: [R] [R-help] how to install own R withour root?

2009-06-14 Thread Patrick Connolly
Or you could do what I suggested below. I had assumed you had a ~/bin/ directory and knew how to make a symbolic link. If that assumption is incorrect, let us know and we can tell you how to do it. As Bill suggests, removing the present old version would be preferable, but I can understand you

[R] estimate the reliability of a scale with dichotomous items

2009-06-14 Thread Martin Batholdy
hi, How can I compute a reliability score of a scale consisting only of dichotomous items? thanks for any help! __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] a proposal regarding documentation

2009-06-14 Thread Philippe Grosjean
John Sorkin wrote: Perhaps help pages should have links to relevant portions of the WIKI. John John David Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics University of Maryland School of Medicine Division of Gerontology Baltimore VA Medical Center 10 North Greene Street GRECC

Re: [R] a proposal regarding documentation

2009-06-14 Thread Peter Flom
I certainly don't have anything against the WIKI, but I think that the documentation is where the action is, especially for newbies. It's the natural first step when you want to learn about a function or when you get an error message you don't understand. Peter Peter L. Flom, PhD

Re: [R] estimate the reliability of a scale with dichotomous items

2009-06-14 Thread William Revelle
At 11:02 PM +0200 6/14/09, Martin Batholdy wrote: hi, How can I compute a reliability score of a scale consisting only of dichotomous items? RSiteSearch(cronbach's alpha) suggests the ltm, psych, psy, psycho, epicalc, packages (among others). There are a number of ways of estimating

Re: [R] a proposal regarding documentation

2009-06-14 Thread spencerg
I agree that the documentation is a primary source, made now more accessible with the availability of the RSiteSearch package that allows more structures searches of the help pages in contributed packages than previously available. On the other hand, it can sometimes be difficult to get

[R] How to fix my nested conditional IF ELSE code?

2009-06-14 Thread Mark Na
Hi, I've been struggling most of the morning with an IF ELSE problem, and I wonder if someone might be able to sort me out. Here's what I need to do (dummy example, my data are more complicated): If type = A or B or C and status = a then count = 1 and status = b then count = 2 and

Re: [R] a proposal regarding documentation

2009-06-14 Thread Gabor Grothendieck
In PHP and also in MySQL the manual has a wiki capability so that users can add notes at the end of each page, e.g. http://www.php.net/manual/en/functions.variable-functions.php http://dev.mysql.com/doc/refman/4.1/en/update.html That would combine documentation and wiki into one. Here it would

Re: [R] How to fix my nested conditional IF ELSE code?

2009-06-14 Thread Chuck Cleland
On 6/14/2009 6:18 PM, Mark Na wrote: Hi, I've been struggling most of the morning with an IF ELSE problem, and I wonder if someone might be able to sort me out. Here's what I need to do (dummy example, my data are more complicated): If type = A or B or C and status = a then count =

Re: [R] a proposal regarding documentation

2009-06-14 Thread Ted Harding
On 14-Jun-09 22:23:24, Gabor Grothendieck wrote: In PHP and also in MySQL the manual has a wiki capability so that users can add notes at the end of each page, e.g. http://www.php.net/manual/en/functions.variable-functions.php http://dev.mysql.com/doc/refman/4.1/en/update.html That

[R] Leave One Out Cross Validation

2009-06-14 Thread muddz
Hi All, I have been trying to get this LOO-Cross Validation method to work on R for the past 3 weeks but have had no luck. I am hoping someone would kindly help me. Essentially I want to code the LOO-Cross Validation for the 'Local Constant' and 'Local Linear' constant estimators. I want to

Re: [R] a proposal regarding documentation

2009-06-14 Thread Gabor Grothendieck
On Sun, Jun 14, 2009 at 6:41 PM, Ted Hardingted.hard...@manchester.ac.uk wrote: That is an interesting and potentially very fruitful idea. I'd like to propose a quasi-readonly mode in which wiki users could interpolate the comments at chosen points within the original text of a help-page (but

Re: [R] a proposal regarding documentation

2009-06-14 Thread Peter Flom
Well, suppose I wanted to suggest changes to some documentation, or write an alternate help file for some function. Where would I put it? Let's say I type, in R, ?median. Now suppose I have suggestions. If I look at http://wiki.r-project.org/rwiki/doku.php?id=idx=rdoc:base I don't see the

[R] Comparative evaluation studies of R packages and functions. Was: Re: a proposal regarding documentation

2009-06-14 Thread Ravi Varadhan
Hi Spencer, You have brought up an important topic. I agree with you that a comparative evaluation of R packages providing similar capabilities is a worthwhile task. Your example of gam is a very good one. Let me give you a couple of other examples of comparative evaluation that I am

[R] Identifying clusters of size n

2009-06-14 Thread Nathan S. Watson-Haigh
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Is there a library which is capable of identifying distinct clusters of size n from a series of XY coordinates? Failing this, I'd like to be able to to something like: Using a sliding window of size n along the x-axis I'd like to determine the

Re: [R] How to fix my nested conditional IF ELSE code?

2009-06-14 Thread David Winsemius
On Jun 14, 2009, at 6:33 PM, Chuck Cleland wrote: On 6/14/2009 6:18 PM, Mark Na wrote: Here's what I need to do (dummy example, my data are more complicated): If type = A or B or C and status = a then count = 1 and status = b then count = 2 and status = c then count = 3 Else

Re: [R] Centering graphics by default in Sweave

2009-06-14 Thread Marc Schwartz
On Jun 14, 2009, at 2:28 PM, Frank E Harrell Jr wrote: Is there an elegant way to get Sweave to center graphics by default? I'd like to use \centerline{\includegraphics{}} etc. to save some vertical space that \begin{center} ... \end{center} uses, and I'd like to avoid centering with each

Re: [R] How to fix my nested conditional IF ELSE code?

2009-06-14 Thread Gabor Grothendieck
Note that TRUE and FALSE become 1 and 0 when used in arithmetic formulae so: result - with(DF, (type %in% c(A, B, C)) * (1 * (status == a) + 2 * (status == b) + 3 * (status == c)) + (type %in% c(D, E, F)) * (9 * (status == a) + 8 * (status == b) +

[R] lack of memory for logistic regression in R?

2009-06-14 Thread Michael
Hi all, I am getting the following error message: mymodel = glm(response ~ . , family=binomial, data=C); Error: cannot allocate vector of size 734.2 Mb In addition: Warning messages: 1: In array(0, c(n, n), list(levs, levs)) : Reached total allocation of 1535Mb: see help(memory.size) 2: In

Re: [R] FBI Homicide data?

2009-06-14 Thread spencerg
Hi, Mario: The following suggest to me that these data may not be available in any contributed R package: library(RSiteSearch) fbi - RSiteSearch.function('FBI Homicide') hits(fbi) # 0 hom - RSiteSearch.function('Homicide') hits(hom) # 4 HTML(hom) The 4 items exposed via the

[R] books on Time series

2009-06-14 Thread Antonio Olinto
Dear list fellows, I want to study time series and use R to analyse time series of fishing data from several species (landings and cpue) investigating the correlation between them and with environmental factors (water temperature, wind, etc.). Searching at Amazon I found three books with

Re: [R] FBI Homicide data?

2009-06-14 Thread David Winsemius
That graphic was created as part of a Ph D thesis in the mid 1990's. I sent the OP off-list a link to that work and the FBI webpage that describes access to such data. -- David On Jun 14, 2009, at 9:11 PM, spencerg wrote: Hi, Mario: The following suggest to me that these data may

Re: [R] lack of memory for logistic regression in R?

2009-06-14 Thread David Winsemius
On Jun 14, 2009, at 9:06 PM, Michael wrote: Hi all, I am getting the following error message: mymodel = glm(response ~ . , family=binomial, data=C); Error: cannot allocate vector of size 734.2 Mb In addition: Warning messages: 1: In array(0, c(n, n), list(levs, levs)) : Reached total

Re: [R] Identifying clusters of size n

2009-06-14 Thread Dylan Beaudette
On Sun, Jun 14, 2009 at 4:39 PM, Nathan S. Watson-Haighnathan.watson-ha...@csiro.au wrote: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Is there a library which is capable of identifying distinct clusters of size n from a series of XY coordinates? Failing this, I'd like to be able to to

Re: [R] if + is.na

2009-06-14 Thread Jeremiah Rounds
Your error message is because if wants a single value and you are giving it a vector. Typically you want to use functions all or any to correct this error message (look them up ?all ?any) and eg if(any(is.na(...))) But in this case to accomplish the task you're after I don't even think

Re: [R] Identifying clusters of size n

2009-06-14 Thread Dylan Beaudette
On Sun, Jun 14, 2009 at 7:26 PM, Nathan S. Watson-Haighnathan.watson-ha...@csiro.au wrote: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Dylan Beaudette wrote: On Sun, Jun 14, 2009 at 4:39 PM, Nathan S. Watson-Haighnathan.watson-ha...@csiro.au wrote: -BEGIN PGP SIGNED MESSAGE-

Re: [R] Centering graphics by default in Sweave

2009-06-14 Thread Frank E Harrell Jr
Thanks very much Marc for the very useful ideas. I do hope there will be an even more concise way to do this. At present I'm not using figure environments in Sweave as in stat reports I don't usually need floats. I tried re-defining \includegraphics to include \centerline but got into

Re: [R] lack of memory for logistic regression in R?

2009-06-14 Thread Michael
I have to use logistic regression... On Sun, Jun 14, 2009 at 8:04 PM, Frank E Harrell Jrf.harr...@vanderbilt.edu wrote: Also it would be useful to compare glm with the lrm function in the Design package, for speed and memory use. Frank David Winsemius wrote: On Jun 14, 2009, at 9:06 PM,

Re: [R] How to fix my nested conditional IF ELSE code?

2009-06-14 Thread Mark Na
Thanks Gabor, this is quite clever and it's nice to see another way of doing it (without ifelse). Mark On Sun, Jun 14, 2009 at 6:51 PM, Gabor Grothendieck ggrothendi...@gmail.com wrote: Note that TRUE and FALSE become 1 and 0 when used in arithmetic formulae so: result - with(DF,

Re: [R] lack of memory for logistic regression in R?

2009-06-14 Thread Frank E Harrell Jr
Michael wrote: I have to use logistic regression... Good. That's what lrm is dedicated to. Frank On Sun, Jun 14, 2009 at 8:04 PM, Frank E Harrell Jrf.harr...@vanderbilt.edu wrote: Also it would be useful to compare glm with the lrm function in the Design package, for speed and memory use.

[R] Sweave template

2009-06-14 Thread Frank E Harrell Jr
Dear Group, I have made significant improvements to our Sweave template, have made the template self-contained (i.e., you can run it yourself and it will find the datasets it needs), and have included the output pdf file. This is at http://biostat.mc.vanderbilt.edu/SweaveTemplate . You will

Re: [R] Identifying clusters of size n

2009-06-14 Thread Nathan S. Watson-Haigh
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Dylan Beaudette wrote: On Sun, Jun 14, 2009 at 4:39 PM, Nathan S. Watson-Haighnathan.watson-ha...@csiro.au wrote: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Is there a library which is capable of identifying distinct clusters of size n

Re: [R] Identifying clusters of size n

2009-06-14 Thread Mose
Hey Nathan, You might like the DBSCAN algorithm. http://en.wikipedia.org/wiki/DBSCAN There's an implementation in the 'fpc' package. http://cran.r-project.org/web/packages/fpc/index.html -Mose On Sun, Jun 14, 2009 at 7:36 PM, Dylan Beaudettedylan.beaude...@gmail.com wrote: On Sun, Jun 14,