[R] specifications windows pc

2008-11-28 Thread [EMAIL PROTECTED]
Hello, I am about to order a new workstation at my university that will be used for R (and other research related tasks). I would appreciate any feedback on the specifications of a very fast machine. The machine should run windows (XP probably better than vista). Which chip, memory size and spec

[R] optimization problem

2008-11-28 Thread tedzzx
Hi, all I am facing an optimization problem. I am using the function optim(par,fun), but I find that every time I give different original guess parameter, I can get different result. For example I have a data frame named data: head(data) price s x t 1 1678.0 12817 11200 0.14959

[R] question involving loops from intro level R programming class

2008-11-28 Thread Heidi Wong
a. Write a R function zerdiag.v1(m) using loop to output a square matrix whose diagonal elements are zero and the other elements are filled in by consecutive integers from 1 to m row-wise. For example, zerdiag.v1(6) = [0, 1, 2] [3, 0, 4] [5, 6,

[R] Question involving loops

2008-11-28 Thread hwong1
a. Write a R function zerdiag.v1(m) using loop to output a square matrix whose diagonal elements are zero and the other elements are filled in by consecutive integers from 1 to m row-wise. For example, zerdiag.v1(6) = [0, 1, 2] [3, 0, 4] [5, 6, 0

[R] Passing bigint (8-byte) values from Postgres to R using RPostgreSQL

2008-11-28 Thread psql2555
Hi - I have a simple R script that inserts values into a Postgres table to generate a record that autogenerates a bigint unique ID (basically, a customized bigserial value). When I query the table to retrieve the bigint value, the value returned is not the same as the value created in the databa

[R] regexp help needed

2008-11-28 Thread Lauri Nikkinen
Hello, I have a vector of dates and I would like to grep the year component from this vector (= all digits after the last punctuation character) dates <- c("28.7.08","28.7.2008","28/7/08", "28/7/2008", "28/07/2008", "28-07-2008", "28-07-08") the resulting vector should look like "08" "2008" "08

Re: [R] regexp help needed

2008-11-28 Thread Kornelius Rohmeyer
Lauri Nikkinen: > Hello, > > I have a vector of dates and I would like to grep the year component > from this vector (= all digits > after the last punctuation character) > > dates <- c("28.7.08","28.7.2008","28/7/08", "28/7/2008", "28/07/2008", > "28-07-2008", "28-07-08") > > the resulting vect

[R] R help

2008-11-28 Thread amit nirmalkar
Dear Sir I am Amit Nirmalkar, M. Sc. Statistics. I generally use SPSS package for analysis. I wanted to run Conditional Logistic Regression in SPSS but this technique is not available. So I thought of R s/w. I am totally blank about R. Can you please send me few of R. Also how to run Conditi

Re: [R] question involving loops from intro level R programming class

2008-11-28 Thread Richard . Cotton
> a. Write a R function zerdiag.v1(m) using loop to output a square matrix > whose diagonal elements are zero and the other elements are filled in by > consecutive integers from 1 to m row-wise. > > For example, > zerdiag.v1(6) = [0, 1, 2] > [3, 0, 4] >

Re: [R] regexp help needed

2008-11-28 Thread Peter Dalgaard
Lauri Nikkinen wrote: > Hello, > > I have a vector of dates and I would like to grep the year component > from this vector (= all digits > after the last punctuation character) > > dates <- c("28.7.08","28.7.2008","28/7/08", "28/7/2008", "28/07/2008", > "28-07-2008", "28-07-08") > > the resultin

[R] editor for MacOS X

2008-11-28 Thread Bunny, lautloscrew.com
Hi all, just wondered again, if there is some R editor for Mac OS X comparable to TINN-R on windows. thx in advance.. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-pr

Re: [R] Question involving loops

2008-11-28 Thread Tom Backer Johnsen
hwong1 wrote: a. Write a R function zerdiag.v1(m) using loop to output a square matrix whose diagonal elements are zero and the other elements are filled in by consecutive integers from 1 to m row-wise. For example, zerdiag.v1(6) = [0, 1, 2] [3, 0, 4]

Re: [R] editor for MacOS X

2008-11-28 Thread Philippe Grosjean
Hello, There is also Komodo Edit + SciViews-K (http://www.sciviews.org/SciViews-K), rather close to Tinn-R functionnalities. Please; note that SciViews-K is not compatible yet with Komodo 5.0 and you must install version 4.4. Best, Philippe Grosjean ..

Re: [R] editor for MacOS X

2008-11-28 Thread Mike Lawrence
I personally have found the built-in editor rather buggy and have switched to TextMate (which I also use for my python programming). On Fri, Nov 28, 2008 at 7:17 AM, Michael Zak <[EMAIL PROTECTED]> wrote: > Why don't you use the bult-in editor from R under OS X? It's much better as > the windows

Re: [R] editor for MacOS X

2008-11-28 Thread Michael Zak
Why don't you use the bult-in editor from R under OS X? It's much better as the windows pendant. Take care, Michael On 28.11.2008, at 12:15, Bunny, lautloscrew.com wrote: Hi all, just wondered again, if there is some R editor for Mac OS X comparable to TINN-R on windows. thx in advance.

Re: [R] specifications windows pc

2008-11-28 Thread seanpor
Good morning Ruud, What sort of tasks are you going to be doing in R? Some tasks will be faster on a single core extreme type processor, and other tasks can benefit from a multi-core processor (which run at slower clock speeds than extreme single-core). If you're working with large matrices, th

Re: [R] optimization problem

2008-11-28 Thread axionator
Hi, I guess your function has several local minima and depending on where you start, i.e. your initial variables, you get into another mimimum. HTH Armin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

Re: [R] regexp help needed

2008-11-28 Thread Gabor Grothendieck
On Fri, Nov 28, 2008 at 5:51 AM, Peter Dalgaard <[EMAIL PROTECTED]> wrote: > Lauri Nikkinen wrote: >> Hello, >> >> I have a vector of dates and I would like to grep the year component >> from this vector (= all digits >> after the last punctuation character) >> >> dates <- c("28.7.08","28.7.2008","

Re: [R] question involving loops from intro level R programming class

2008-11-28 Thread Gustaf Rydevik
On Fri, Nov 28, 2008 at 8:35 AM, Heidi Wong <[EMAIL PROTECTED]> wrote: > a. Write a R function zerdiag.v1(m) using loop to output a square matrix > whose diagonal elements are zero and the other elements are filled in by > consecutive integers from 1 to m row-wise. > > For example, > zerdiag.v1(6)

Re: [R] editor for MacOS X

2008-11-28 Thread John Fox
Dear Bunny, I've been using Eclipse with the StatET plug-in under both Windows and Mac OS X. Eclipse + StatET provides much more than a code editor, such as the ability to check and build packages and to interact with an svn archive. On the downside, it requires

Re: [R] 1-Pearson's R Distance

2008-11-28 Thread Rodrigo Aluizio
Claudia, Thomas and William thank you all so much. It's exactly what we needed. Works perfectly. Rodrigo. -Mensagem original- De: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Em nome de Claudia Beleites Enviada em: quinta-feira, 27 de novembro de 2008 14:01 Para: r-help@r-project.org Assun

[R] confidence interval for glm

2008-11-28 Thread Gerard M. Keogh
Hi all, simple Q: how do I extract the upper and lower CI for predicted probabilities directly for a glm - I'm sure there's a one line to do it but I can't find it. the predicted values I get with the predict (.. "response") Thanks Gerard *

Re: [R] More than doubling performance with snow

2008-11-28 Thread Luke Tierney
Hi Markus, I'm happy to participate in this, as I think I said previously. I won't have time to look carefully at the draft until sometime next week, but I remain puzzled about the high time listed for case 3 with snow/Rmpi. It would be good to understand what is going on there -- the discrepan

[R] [R-pkgs] version 2.0-0 of the effects package

2008-11-28 Thread John Fox
Dear R users, I'd like to announce version 2.0-0 of the effects package, which is now on CRAN. The package constructs graphical and tabular "effect displays," e.g., of interactions, for several kinds of models. These displays are generalizations of "adjusted" means in linear models. The major ad

[R] Transpose and replicate

2008-11-28 Thread Alex99
Hi all, I used the "replicate" function to make 2 samples from file test with 5 records in each sample and saved the output in "result". result=replicate(2,test[,sample(colnames(test),5,replace = FALSE)],simplify=FALSE) now I need to transpose each sample, but when I use: t(result) or res

Re: [R] R help

2008-11-28 Thread C.H.
Actually, SPSS do have condition LR which can be simulate using Cox Proportional Hazards regression. For R, Conditional LR is clogit inside survival package. ?clogit may help. Regards, C On Fri, Nov 28, 2008 at 5:43 PM, amit nirmalkar <[EMAIL PROTECTED]> wrote: > Dear Sir > > I am Amit Nirmal

[R] Examples of advanced data visualization

2008-11-28 Thread Hans W. Borchers
Dear R-help, I am looking for ideas and presentations of new and advanced data visualization methods. As an example of what I am searching for, the 'Many Eyes' pages at http://manyeyes.alphaworks.ibm.com/manyeyes/ may provide a good paradigm. I would be interested even if it will not be easy

Re: [R] editor for MacOS X

2008-11-28 Thread Blanchette, Marco
Carbon emacs ( http://homepage.mac.com/zenitani/emacs-e.html) using the ess-mode ( http://ess.r-project.org/). Amazingly good integration of different buffer types for different tasks. You can have your R session running in a buffer, a .R buffer where you edit your functions/sources and with ver

[R] problem with function plm

2008-11-28 Thread valeria pedrina
Hi all, I want to use the estimation method fixed effect on my panel data. I have found plm function but it give me a problem,my model is: 1/5*log(z_{i,t}/r_{i,t})=a+b*r_{i,t}+u_{i,t} , I have 9 individuals and 5 time for each individuals so I have 45 observation this is what I do: > p <- c("ital

[R] Contrast SS

2008-11-28 Thread Ista Zahn
Hi all, Hopefully this is an easy question and I've overlooked something simple. I want to present specific contrasts in an ANOVA table. I can of course get the significance tests with summary.lm, but my adviser wants the results in ANOVA format. Here is a minimal example: > x <- factor(c(rep("con

Re: [R] optimization problem

2008-11-28 Thread John C Nash
tedzxx asked about apparent multiple optima. See below. Users should be aware that optim() does local optimization. The default Nelder-Mead approach is fairly robust at finding such a local minimum, though it may halt if it is on a flat area of the loss function surface. I would recommend tryi

Re: [R] Request for Assistance in R with NonMem

2008-11-28 Thread Douglas Bates
On Wed, Nov 26, 2008 at 12:46 PM, Michael White <[EMAIL PROTECTED]> wrote: > Hi > I am having some problems running a covariate analysis with my > colleage using R with the NonMem program we are using for a graduate > school project. R and NonMem run fine without adding in the > covariates, but the

Re: [R] Contrast SS

2008-11-28 Thread Gabor Grothendieck
Is your objective to break down the sum of squares for x into one for the group 1 vs group 3 contrast and another for the group 2 vs group 3 contrast? If that's it try this. We create the model matrix, make the two columns generated from x into separate explanatory variables, x1 and x2, and then

Re: [R] Examples of advanced data visualization

2008-11-28 Thread David Winsemius
The gapminder.org project ; http://graphs.gapminder.org/world/ Hans Rosling at TED: http://www.ted.com/index.php/talks/hans_rosling_shows_the_best_stats_you_ve_ever_seen.html http://hdr.undp.org/external/gapminder/2004/hdr2004.html And was given Google support last year: http://googleblog.blo

Re: [R] Contrast SS

2008-11-28 Thread Richard M. Heiberger
I believe you are looking for the summary.lm() function and its split= argument. Look at example(summary.aov) for several examples. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http:/

Re: [R] Contrast SS

2008-11-28 Thread Ista Zahn
Thank you! > m.1 <- aov(y ~ x) > summary(m.1, split = list(x = list("groups 1 and 2 vs group 3" = 1, "group 1 > vs group 2" = 2))) does what I want. On Fri, Nov 28, 2008 at 12:02 PM, Richard M. Heiberger <[EMAIL PROTECTED]> wrote: > I believe you are looking for the summary.lm() function and i

[R] side by side boxplots

2008-11-28 Thread Phillip Porter
Good Morning, I am trying to get side by side boxplots of two groups on the same variable. The last item under ?boxplot led me to some useful code. I use "boxwex" to make the boxes narrower, "at" to shift them over and "add" to draw them both on the same graph. Something along the lines of

[R] ranking the results of a questionnaire

2008-11-28 Thread Wolfgang Lindner
dear experts, I reproduced an experiment (questionnaire) some times. The result of the experiment is a vector of 5 factors, say (A,B,C,D,E). In the original article the result is given in 5 pairs of mean and stDev for A .. E, e.g. mean_A=37.4 and sd_A=8.1. The interval for A,B,C,D,E values is 0..5

[R] Help: "mean" and "replicate"

2008-11-28 Thread Alex99
Hi all, I used the "replicate" function to make 2 samples from file "test" with 5 records in each sample and saved the output in "result". result=replicate(2,test[,sample(colnames(test),5,replace =FALSE)],simplify=FALSE) now I need to calculate the mean of each cell; i.e. I need to add "the nu

Re: [R] optimization problem

2008-11-28 Thread tedzzx
If I want to find out the globle minia, how shoul I change my code? Thanks a lot Armin Meier wrote: > > Hi, > I guess your function has several local minima and depending on where > you start, i.e. your initial variables, you get into another mimimum. > > HTH > Armin > > __

[R] Regarding posting a package to R-forge (with one of the dependent packages not in CRAN)

2008-11-28 Thread Marcelo Perlin
Hi Guys, Recently I wrote a package for dealing with Markov Switching Regressions in R and it is included in the Rmetrics project. https://r-forge.r-project.org/projects/rmetrics/ Everything works fine when I use it in computer. But, the package depends on the use of optimization functions from

[R] R window lines limit

2008-11-28 Thread Rthoughts
Hi, Is the number of lines that can be displayed in R GUI limited? When I do a Fourier Tranform of 4000 data points, I can only see the last 300 or so of it. I need to see all of it. How is this possible? Thank you if you can help! -- View this message in context: http://www.nabble.com/R-wi

Re: [R] R window lines limit

2008-11-28 Thread hadley wickham
On Fri, Nov 28, 2008 at 9:00 AM, Rthoughts <[EMAIL PROTECTED]> wrote: > > Hi, > > Is the number of lines that can be displayed in R GUI limited? > > When I do a Fourier Tranform of 4000 data points, I can only see the last > 300 or so of it. I need to see all of it. How is this possible? Perhaps y

Re: [R] R window lines limit

2008-11-28 Thread Gabor Grothendieck
Try ?View On Fri, Nov 28, 2008 at 10:00 AM, Rthoughts <[EMAIL PROTECTED]> wrote: > > Hi, > > Is the number of lines that can be displayed in R GUI limited? > > When I do a Fourier Tranform of 4000 data points, I can only see the last > 300 or so of it. I need to see all of it. How is this possible

Re: [R] side by side boxplots

2008-11-28 Thread Dylan Beaudette
On Fri, Nov 28, 2008 at 10:55 AM, Phillip Porter <[EMAIL PROTECTED]> wrote: > Good Morning, >I am trying to get side by side boxplots of two groups on the same > variable. The last item under ?boxplot led me to some useful code. >I use "boxwex" to make the boxes narrower, "at" to s

Re: [R] R window lines limit

2008-11-28 Thread Duncan Murdoch
On 28/11/2008 10:00 AM, Rthoughts wrote: Hi, Is the number of lines that can be displayed in R GUI limited? Which GUI? In Windows, the console buffer size is configurable. The defaults are 250,000 characters, 8000 lines; it doesn't sound as though you are hitting those limits. But if you

Re: [R] Examples of advanced data visualization

2008-11-28 Thread Tom Backer Johnsen
Hans W. Borchers wrote: Dear R-help, I am looking for ideas and presentations of new and advanced data visualization methods. As an example of what I am searching for, the 'Many Eyes' pages at http://manyeyes.alphaworks.ibm.com/manyeyes/ may provide a good paradigm. I would be interested

[R] save single image when using plot(lm(y~x))

2008-11-28 Thread Schreiber, Stefan
Dear List! When I'm using "plot(lm(y~x))" R shows me 4 different plots. They pop-up individually by clicking on the graph. How can I save for example the 2nd plot? I've found a camera button that says "copy to the clipboard as a metafile" once you place the mouse pointer on it. I tried, though it

[R] NA and logical indexes

2008-11-28 Thread Sebastian P. Luque
Hi, I vaguely remember this issue being discussed at some length in the past, but am having trouble relocating the proper thread (defining an adequate search string to do so): ---<---cut here---start-->--- R> foo <- data.frame(A=gl(2, 5, labels=letters[1:2]), X

Re: [R] save single image when using plot(lm(y~x))

2008-11-28 Thread Juliet Hannah
Try: mylm <- plot(lm(y~x)) plot(mylm,which=2) On Fri, Nov 28, 2008 at 4:12 PM, Schreiber, Stefan <[EMAIL PROTECTED]> wrote: > Dear List! > > When I'm using "plot(lm(y~x))" R shows me 4 different plots. They pop-up > individually by clicking on the graph. How can I save for example the > 2nd plot

Re: [R] save single image when using plot(lm(y~x))

2008-11-28 Thread Schreiber, Stefan
Woohoo! Thanks a lot Juliet! Stefan -Original Message- From: Juliet Hannah [mailto:[EMAIL PROTECTED] Sent: Friday, November 28, 2008 2:42 PM To: Schreiber, Stefan Cc: r-help@r-project.org Subject: Re: [R] save single image when using plot(lm(y~x)) Try: mylm <- plot(lm(y~x)) plot(mylm,

Re: [R] NA and logical indexes

2008-11-28 Thread Ted Harding
On 28-Nov-08 21:25:36, Sebastian P. Luque wrote: > Hi, > I vaguely remember this issue being discussed at some length in the > past, but am having trouble relocating the proper thread (defining an > adequate search string to do so): > > ---<---cut here---start--

[R] RODBC - problems connecting to oracle through linux

2008-11-28 Thread Simon Collins
Hi I'm presently trying to connect to Oracle through RODBC / UnixODBC on linux (Red Hat 5). I've successfully connected through the easysoft drivers but the trial license has run out and I'm trying to connect via the default Oracle 11G driver. However I'm getting the following error li

Re: [R] NA and logical indexes

2008-11-28 Thread Sebastian P. Luque
On Fri, 28 Nov 2008 22:01:15 - (GMT), (Ted Harding) <[EMAIL PROTECTED]> wrote: [...] > It is not! Oops, of course, I meant to say "return" rather than "set" :-/ [...] > NA can seem to have a bewildering logic, but it all becomes clear if > you interpret NA as "value unkown". > You asked

Re: [R] save single image when using plot(lm(y~x))

2008-11-28 Thread Paul Hiemstra
Hi, The plot function in R is a so called generic function, an appropriate function is chosen based on the input object provided by the user. In this case an lm-object (use the class() command to find this). The plot function that R calls is then called plot.lm. ?plot.lm gives you the help fi

[R] error message when using strptime

2008-11-28 Thread Kara Przeczek
I am new to R and have been struggling with the following problem. I apologize if there is an obvious answer that I have missed. I bring in a .txt file (pyr_h.txt) that has a DateTime column in the format "-mm-dd HH:MM" When it comes into R, R does not recognize it as a date/time format.

Re: [R] error message when using strptime

2008-11-28 Thread jim holtman
strptime has output of a POSIXlt object and you want an POSIXct object, so use: as.POSIXct(strptime( On Fri, Nov 28, 2008 at 7:12 PM, Kara Przeczek <[EMAIL PROTECTED]> wrote: > I am new to R and have been struggling with the following problem. I > apologize if there is an obvious answer that

Re: [R] error message when using strptime

2008-11-28 Thread Gabor Grothendieck
You want to convert the datetime to POSIXct, not POSIXlt: > as.POSIXct("2005-12-20 12:45") [1] "2005-12-20 12:45:00 EST" See R News 4/1 for details. On Fri, Nov 28, 2008 at 7:12 PM, Kara Przeczek <[EMAIL PROTECTED]> wrote: > I am new to R and have been struggling with the following problem. I >

Re: [R] save single image when using plot(lm(y~x))

2008-11-28 Thread Schreiber, Stefan
Thanks Paul!! I wasn't aware of the command "plot.lm". Therefore, I couldn't find the help-file either. Cheers, Stefan -Original Message- From: Paul Hiemstra [mailto:[EMAIL PROTECTED] Sent: Friday, November 28, 2008 5:10 PM To: Juliet Hannah Cc: Schreiber, Stefan; r-help@r-project.org

Re: [R] Help processing large data

2008-11-28 Thread jim holtman
Is this what you want: > x <- read.table(textConnection('"read" "no" "length" + 2 2 144 + 7 7 47490 + 9 9 310944 + 11 11 10089 + 14 14 13152 + 17 17 27363 '), header=TRUE) > closeAllConnections() > result <- lapply(1:nrow(x), function(.indx){ + data.frame(read=paste(x$read[.indx], seq(x$length

Re: [R] editor for MacOS X

2008-11-28 Thread C.H.
Aquamacs preinstalled ess-mode. http://aquamacs.org/ On Fri, Nov 28, 2008 at 11:57 PM, Blanchette, Marco <[EMAIL PROTECTED]> wrote: > Carbon emacs ( http://homepage.mac.com/zenitani/emacs-e.html) using the > ess-mode ( http://ess.r-project.org/). Amazingly good integration of > different buffe

Re: [R] NA and logical indexes

2008-11-28 Thread Berwin A Turlach
G'day Sebastian, On Fri, 28 Nov 2008 16:35:13 -0600 "Sebastian P. Luque" <[EMAIL PROTECTED]> wrote: > On Fri, 28 Nov 2008 22:01:15 - (GMT), > (Ted Harding) <[EMAIL PROTECTED]> wrote: > > > NA can seem to have a bewildering logic, but it all becomes clear if > > you interpret NA as "value unko