Re: [R] ESS: Viewing Restored Workspace At Startup

2011-09-13 Thread Patrick Connolly
On Tue, 13-Sep-2011 at 06:53AM -0700, Rich Shepard wrote: |> I use ESS and have no problems invoking it when I fire up emacs. When it |> loads it tells me, [Previously saved workspace restored]. However, I do not |> know how to view that history of commands so I can continue what I was |> doing.

Re: [R] ESS: Viewing Restored Workspace At Startup

2011-09-13 Thread Rich Shepard
On Wed, 14 Sep 2011, Patrick Connolly wrote: The objects will be available already. See what you have with the ls() command. They are stored in .RData so the fact that you can't read .RData directly matters not. Patrick, I'm aware of .Rdata and use ls() regularly outside of R. What I did n

Re: [R] barplot in hexagram layout

2011-09-13 Thread Schatzi
I decided to go with circles instead of rectangles. Thank you for your help. Here is the new code: dev.new(width=2.5, height=3,mar=c(0,0,0,0)) par(mfrow=c(1,1),mar=c(0,0,0,0)) x=c(-1,1,1,-1,-3,-3,5) y=c(1.2,0.6,-.7,-1.3,-.7,0.6,5) plot(0,xlim=c(-4,2),ylim=c(-2,2),type="n",axes=FALSE,xlab="",ylab="

Re: [R] writing a loop using several list()-objects

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 10:26 AM, travellina wrote: Hi there, I am looking for some help with a loop which uses three list()- objects: "sequence" has got numbers, "start" and "end" have got the begin/end time points which should be applied to the numbers in "sequence". length() of all three o

Re: [R] space in directory name

2011-09-13 Thread Mikkel Grum
You are missing \\ between Documents and settings and Administrator. - Original Message - From: dbonneau To: r-help@r-project.org Cc: Sent: Tuesday, September 13, 2011 10:59 AM Subject: [R] space in directory name Hi, I am trying to read a text file located in following paths. I am get

[R] S4 vs Reference Classes

2011-09-13 Thread Joseph Park
Hi, I'm looking for some guidance on whether to use S4 or Reference Classes for an analysis application I'm developing. I'm a C++/Python developer, and like to 'think' in OOD. I started my app with S4, thinking that was the best set of OO features in R. However, it appears that o

[R] count regression zero count comparison

2011-09-13 Thread jekang
Hello, It might be more of a statistical question than an R question. I was reading http://cran.r-project.org/web/packages/pscl/vignettes/countreg.pdf, and I was wondering why the following functions were used to compare zero counts (observed and predicted), instead of just using hist(fitted(fm_

Re: [R] How to calculate the power of Wilcoxon signed rank test

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 3:34 PM, tn85 wrote: Hello All, I perform a Wilcoxon signed rank test for two sets of data to test whether they two have significantly different means. I would also like to know the power of this test. Given that none of the various "Wilcoxon tests" are for differen

Re: [R] Force regression line to a 1:1 relationship

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 11:56 AM, David Winsemius wrote: On Sep 13, 2011, at 11:44 AM, David Winsemius wrote: On Sep 13, 2011, at 9:43 AM, RCulloch wrote: Dear John, Thank you for that, and for explaining why the abline() command wont/dosen't work. The approach is based on reviewers commen

[R] Online literature database search in R - RCurl?

2011-09-13 Thread Brian Oney
Hi All, theoretically, I should be able (with proper proxy setting and IP address) to connect and perform literature database searches, such as in the ISI Web of Science, with R. I can imagine this working nicely with the RCurl package. I am pestering the list to ask if I missed a certain pa

Re: [R] import csv file into R, strange problem

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 12:58 PM, Greg Snow wrote: It would appear that your file is not a csv file, but we cannot tell for sure without seeing an example of the data that you are reading in, can you cut and paste the first few lines of the file? Or post a link to a copy of the file. If the d

Re: [R] space in directory name

2011-09-13 Thread Jean-Christophe BOUËTTÉ
And a final word of advice: be lazy, and use file.choose() to get the string that R expects as a filename. This will help avoid such typos! JC 2011/9/13 Mikkel Grum : > You are missing \\ between Documents and settings and Administrator. > > > - Original Message - > From: dbonneau > To:

Re: [R] count regression zero count comparison

2011-09-13 Thread Achim Zeileis
On Tue, 13 Sep 2011, jekang wrote: Hello, It might be more of a statistical question than an R question. I was reading http://cran.r-project.org/web/packages/pscl/vignettes/countreg.pdf, and I was wondering why the following functions were used to compare zero counts (observed and predicted),

[R] Reading R Code aloud

2011-09-13 Thread andrewH
Dearfolks-- I have been told by an experienced R programmer and teacher whom I trust that it is easier to understand R code if you read it aloud, as the language that it is. However, she was clear that reading it aloud was not simply reading the marks on the screen: you read A.df[5,] as "the fift

Re: [R] How to calculate the power of Wilcoxon signed rank test

2011-09-13 Thread Tianchan Niu
Thank you David for your reply. But I am still confusing. I am using Wilcoxon test to compare two samples to assess whether their means differ because neither of them is normally distributed. Is it possible to compute the power of the Wilcoxon test similar to that of t test? Or is it just a wron

Re: [R] S4 vs Reference Classes

2011-09-13 Thread Steve Lianoglou
Hi, On Tue, Sep 13, 2011 at 1:54 PM, Joseph Park wrote: > >   Hi, I'm looking for some guidance on whether to use >   S4 or Reference Classes for an analysis application >   I'm developing. >   I'm a C++/Python developer, and like to 'think' in OOD. >   I started my app with S4, thinking that was

Re: [R] x %>% y as an alternative to which( x > y)

2011-09-13 Thread Timothy Bates
Dear Duncan and Hadley, I stumbled across the NA behavior of subset a little while ago and thought it might do the trick. But my common usage case is not getting a subsetting sans NAs, but setting values in the whole dataframe. So I need T/F at each row, not just the list of rows that match th

Re: [R] S4 vs Reference Classes

2011-09-13 Thread Douglas Bates
On Tue, Sep 13, 2011 at 12:54 PM, Joseph Park wrote: >   Hi, I'm looking for some guidance on whether to use >   S4 or Reference Classes for an analysis application >   I'm developing. >   I'm a C++/Python developer, and like to 'think' in OOD. >   I started my app with S4, thinking that was the

Re: [R] S4 vs Reference Classes

2011-09-13 Thread Douglas Bates
On Tue, Sep 13, 2011 at 4:11 PM, Steve Lianoglou wrote: > Hi, > > On Tue, Sep 13, 2011 at 1:54 PM, Joseph Park wrote: >> >>   Hi, I'm looking for some guidance on whether to use >>   S4 or Reference Classes for an analysis application >>   I'm developing. >>   I'm a C++/Python developer, and like

Re: [R] x %>% y as an alternative to which( x > y)

2011-09-13 Thread Duncan Murdoch
On 11-09-13 5:17 PM, Timothy Bates wrote: Dear Duncan and Hadley, I stumbled across the NA behavior of subset a little while ago and thought it might do the trick. But my common usage case is not getting a subsetting sans NAs, but setting values in the whole dataframe. So I need T/F at each r

Re: [R] x %>% y as an alternative to which( x > y)

2011-09-13 Thread William Dunlap
I often use the following function is.true <- function(x) !is.na(x) & x and, less often, is.false <- function(x) !is.na(x) & !x to report if elements of a logical vector are TRUE (not FALSE or NA) or FALSE (not TRUE or NA), respectively. Do your complicated logical expression and apply is.true

[R] stupid lm() question

2011-09-13 Thread Carl Witthoft
I feel bad even asking, but: Rgames> data(OrchardSprays) Rgames> model<-lm(decrease~.,data=OrchardSprays) Rgames> model Call: lm(formula = decrease ~ ., data = OrchardSprays) Coefficients: (Intercept) rowpos colpos treatmentB treatmentC 22.705 -2.784 -1.234

[R] CMYK color space

2011-09-13 Thread ivo welch
dear R experts---I am struggling with the requirements to prepare my files for my printers. I am printing in 2/2 format, which means cyan and black for me, which they take from my color-separated pdf files. R comes into play, because it produces all the figures that are embedded in my book (pdfla

Re: [R] function censReg in panel data setting

2011-09-13 Thread Arne Henningsen
Hi Igors On 13 September 2011 13:27, Igors wrote: > Any success in finding possible solutions for my problem? Somewhat. The calculation of the log-likelihood values is numerically much more robust/stable now. The log-likelihood contributions of some individuals became minus infinity in your mode

Re: [R] How to calculate the power of Wilcoxon signed rank test

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 4:51 PM, Tianchan Niu wrote: Thank you David for your reply. But I am still confusing. I am using Wilcoxon test to compare two samples to assess whether their means differ because neither of them is normally distributed. Is it possible to compute the power of the Wilcoxon

Re: [R] ZOO: Learning to apply it to my data

2011-09-13 Thread Gabor Grothendieck
On Tue, Sep 13, 2011 at 2:07 PM, Rich Shepard wrote: >  I have read ?zoo but am not sure how to relate the parameters (x, > order.by, frequency, and style) to my data.frame. The structure of the > data.frame is > > 'data.frame':   11169 obs. of  4 variables: >  $ stream  : Factor w/ 37 levels "Bur

Re: [R] CMYK color space

2011-09-13 Thread Henrik Bengtsson
Did you ask your question being aware of the following from help("pdf") in R v2.13.1: pdf(file = ifelse(onefile, "Rplots.pdf", "Rplot%03d.pdf"), width, height, onefile, family, title, fonts, version, paper, encoding, bg, fg, pointsize, pagecentre, colormodel, useDingbats, useKerning,

Re: [R] ZOO: Learning to apply it to my data

2011-09-13 Thread Rich Shepard
On Tue, 13 Sep 2011, Gabor Grothendieck wrote: As in ?zoo a zoo object is a numeric matrix, numeric vector or factor together with an ordered time index which is unique. Its not clear that that is what you have; however, if we can assume that for each value of param we have a unique set of dates

Re: [R] CMYK color space

2011-09-13 Thread ivo welch
obviously not. thank you, henrik. going back to square 101. (I just googled, because I usually find stuff faster by googling. cmyk and R brought up nothing.) apologies for the bandwidth, everyone. Ivo Welch (ivo.we...@gmail.com) On Tue, Sep 13, 2011 at 4:24 PM, Henrik Bengtsson wrote:

Re: [R] stupid lm() question

2011-09-13 Thread Rolf Turner
``The only stupid question is the one that is not asked.'' --- Anon. (???) The reason that there is no value for level "A" of treatment is that level "A" is the reference level, under the default ``treatment'' contrasts. See ?contr.treatment. Some insight may be obtained by doing:

Re: [R] ZOO: Learning to apply it to my data

2011-09-13 Thread Gabor Grothendieck
On Tue, Sep 13, 2011 at 7:36 PM, Rich Shepard wrote: > On Tue, 13 Sep 2011, Gabor Grothendieck wrote: > >> As in ?zoo a zoo object is a numeric matrix, numeric vector or factor >> together with an ordered time index which is unique. Its not clear that >> that is what you have; however, if we can a

Re: [R] ZOO: Learning to apply it to my data

2011-09-13 Thread Rich Shepard
On Tue, 13 Sep 2011, Gabor Grothendieck wrote: In the example data there were no duplicate dates for a given param value. Gabor, There are duplicate dates for a given param value, but each is on a different stream. It is the triplet (stream, sampdate, param) that is associated with each val

[R] solving linear equations

2011-09-13 Thread varunshivashankar
I have a dataset X Y1 12001.375 40000.464 1333.33 0.148 444.44 0.047 148.148 0.014 49.383 0.005 16.461 0.004 I have to find a curve fit for the above dataset based on a 4-parameter logistic equation viz. Y1 = d + ((a-d)/(1+(X/cc)^b)), where X and Y1 are the values above. I nee

Re: [R] count regression zero count comparison

2011-09-13 Thread jekang
Thank you, it helped me clear some confusions. Jamie -- View this message in context: http://r.789695.n4.nabble.com/count-regression-zero-count-comparison-tp3810907p3811277.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-pro

Re: [R] drc results for binomial data

2011-09-13 Thread PatGauthier
I've been having trouble making sense of the drc results for my binomial response toxicity data. Firstly, the standard errors are far too large for how well the data fit the log-logistic model, particularly compared to other methods of LC50 estimation (e.g. probit and trimmed spearman-karber). Seco

Re: [R] solving linear equations

2011-09-13 Thread Ben Bolker
varunshivashankar gmail.com> writes: > > I have a dataset > > X Y1 > 1200 1.375 > 4000 0.464 > 1333.33 0.148 > 444.440.047 > 148.148 0.014 > 49.3830.005 > 16.4610.004 > > I have to find a curve fit for the above dataset based on a 4-parameter > logist

Re: [R] Reading R Code aloud

2011-09-13 Thread cwdillon
Well, I've bee reading the R Manual, Nathan Yau's book on Flowing Data, and a few others like this to my kids to get them to go to sleep. I started off doing it the *wrong* way, just to keep it boring but slipped into doing it the *correct* way for my own sanity. It puts the kids to sleep in no tim

Re: [R] S4 vs Reference Classes

2011-09-13 Thread Martin Morgan
On 09/13/2011 10:54 AM, Joseph Park wrote: Hi, I'm looking for some guidance on whether to use S4 or Reference Classes for an analysis application I'm developing. I'm a C++/Python developer, and like to 'think' in OOD. I started my app with S4, thinking that was the best

[R] Difficulty with saving .Rhistory

2011-09-13 Thread T
Current systems     OS: RedHat SE Linux 5.7, fully patched     R: R-2.13.1 I am a new R user. The R installation needs to be very tight because I plan to use it after I am employed on Internet facing systems to analyse firewall and filesystem log data.  Everything to my knowledge regarding

[R] how to get a row with its its p value below 0.01 ??

2011-09-13 Thread anand m t
Hi all, I'm analyzing micro array data.. it has produced a file (to be specific matrix) withdimension of 35556 2. first few lines of the matrix are as below.. probe_name control.fdr.pvals.present 10338001 0.000440001 10338002 0.000583093 10338003 0.000528449 10338004 0.000610362 10338005 0.000151

Re: [R] how to get a row with its its p value below 0.01 ??

2011-09-13 Thread Rolf Turner
On 14/09/11 17:49, anand m t wrote: Hi all, I'm analyzing micro array data.. it has produced a file (to be specific matrix) withdimension of 35556 2. first few lines of the matrix are as below.. probe_name control.fdr.pvals.present 10338001 0.000440001 10338002 0.000583093 10338003 0.000528449

Re: [R] help with hclust and cutree

2011-09-13 Thread Madeleine Seeland
Hello, yes, thanks! Also,how do I get the cutoff similarity value which was used for clustering when I cut the tree at a specific heigt? Best regards, Madeleine On Tue, 2011-09-13 at 09:39 -0700, Mikkel Grum wrote: > plot(clust) > rect.hclust(clust, h = 0.65) > > > > - Original Message

<    1   2