Re: [R] extended Kalman filter for survival data

2010-06-12 Thread Christos Argyropoulos
If you mean this paper by Fahrmeir: http://biomet.oxfordjournals.org/cgi/content/abstract/81/2/317 I would recommend BayesX: http://www.stat.uni-muenchen.de/~bayesx/. BayesX interfaces with R and estimates discrete (and continuous) time survival data with penalized regression methods. If

[R] sleep timer resolution on OSX

2010-06-12 Thread ivo welch
I am doing some timing experiments. To test looping performance, I used the Sys.sleep function. I noticed something in the docs that is just a little misleading: The resolution of the time interval is system-dependent, but will normally be down to 0.02 secs or better. (On modern

[R] mob (party package) question

2010-06-12 Thread tudor
Dear useRs: I try to use mob from the party package (thanks Achim and Co.!) to model based recursive partition a data set. The model is a logistic regression specified with model=glinearModel and family=binomial(). Running mob results in a few warnings of the type: In glm.fit ... algorithm did

[R] Logic with regexps

2010-06-12 Thread Ted Harding
Greetings, The following question has come up in an off-list discussion. Is it possible to construct a regular expression 'rex' out of two given regular expressions 'rex1' and 'rex2', such that a character string X matches 'rex' if and only if X matches 'rex1' AND X does not match 'rex2'? The

[R] meta analysis with repeated measure-designs?

2010-06-12 Thread Gerrit Hirschfeld
Dear all, I am trying to run a meta analysis of psycholinguistic reaction-time experiments with the meta package. The problem is that most of the studies have a within-subject designs and use repeated measures ANOVAs to analyze their data. So at present it seems that there are three

Re: [R] Compiling R with multi-threaded BLAS math libraries - why not actually ?

2010-06-12 Thread Tal Galili
Hello Gabor, Matt, Dirk. Thank you all for clarifying the situation. So if I understand correctly then: 1) Changing the BLAST would require specific BLAST per computer configuration (OS/chipset). 2) The advantage would be available only when doing _lots_ of linear algebra So I am left

Re: [R] points marking

2010-06-12 Thread khush ........
Hi, Well Thanks for letting me know that pch is of no use with segments petr. I am using lend as it suits to me more as gregory suggested , but I am not getting imite??? think I try to fix it with some other method also, as I have to deal more with the symbols in this case, But I want to the

[R] help with npcdens function package np

2010-06-12 Thread Steve Friedman
Hello, I'm trying to calculate non-parametric probabilities using the np package and having some difficulties. OS is Windows, R version 2.11.1 Here is what I've done so far. library(np) veg - data.frame(factor(Physiogomy), meanAnnualDepthAve, TP) attach(veg) : for clarification dim(veg)

Re: [R] Logic with regexps

2010-06-12 Thread Prof Brian Ripley
I think you have missed grepl(), e.g. X[grepl(rex1, X) !grepl(rex2, X)] grepl is a fairly recent addition (2.9.0) that is used extensively in R's own text-processing operations (e.g. help files, utilities such as 'R CMD check'). On Sat, 12 Jun 2010, ted.hard...@manchester.ac.uk wrote:

Re: [R] Logic with regexps

2010-06-12 Thread Ted Harding
Thanks, Brian. I had indeed overlooked grepl() (too busy delving into the syntax summary)! That is certainly a useful shortcut of the construction I had used. Your Not in general implies that using grep() twice (in this example; more times in more complex combinations) is inevitable -- which of

Re: [R] extended Kalman filter for survival data

2010-06-12 Thread Christophe Dutang
Thanks Christos. BayesX was strongly recommended also by Fahrmeir and other people of that field. So it is clearly where I need to look at! Christophe iPhone.fan Le 12 juin 2010 à 04:46, Christos Argyropoulos argch...@hotmail.com a écrit : If you mean this paper by Fahrmeir:

[R] Can one get a list of recommended packages?

2010-06-12 Thread Dr. David Kirkby
R 2.10.1 is used in the Sage maths project. Several recommended packages (Matrix, class, mgcv, nnet, rpart, spatial, and survival) are failing to build on Solaris 10 (SPARC). We would like to be able to get a list of the recommended packages for R 2.10.1, but ideally via a call to R, so it is

Re: [R] Compiling R with multi-threaded BLAS math libraries - why not actually ?

2010-06-12 Thread Douglas Bates
On Sat, Jun 12, 2010 at 6:18 AM, Tal Galili tal.gal...@gmail.com wrote: Hello Gabor, Matt, Dirk. Thank you all for clarifying the situation. So if I understand correctly then: 1) Changing the BLAST would require specific BLAST per computer configuration (OS/chipset). It's BLAS (Basic

[R] Displaying homogeneous groups in aov post-hoc results ?

2010-06-12 Thread Tal Galili
Hello dear R-help mailing list, A friend of mine teaches a regression and experimental design course and asked me the following question. She is trying to find a way to display the homogeneous groups (after performing tukey test on an aov object). here's an example for what she means by

[R] generating ordered, random decimal fractions

2010-06-12 Thread Kurt_Helf
Greetings How do I do this in R? Checking the Cran site produces a bewildering array of packages that I can't seem to find to load. Surely the main program has this function? Cheers Kurt *** Kurt Lewis Helf, Ph.D. Invertebrate

Re: [R] meta analysis with repeated measure-designs?

2010-06-12 Thread Viechtbauer Wolfgang (STAT)
Dear Gerrit, the most appropriate approach for data of this type would be a proper multivariate meta-analytic model (along the lines of Kalaian Raudenbush, 1996). Since you do not know the correlations of the reaction time measurements across conditions for the within-subject designs, a

Re: [R] Logic with regexps

2010-06-12 Thread Gabor Grothendieck
On Sat, Jun 12, 2010 at 5:38 AM, Ted Harding ted.hard...@manchester.ac.uk wrote: Greetings, The following question has come up in an off-list discussion. Is it possible to construct a regular expression 'rex' out of two given regular expressions 'rex1' and 'rex2', such that a character string

Re: [R] Compiling R with multi-threaded BLAS math libraries - why not actually ?

2010-06-12 Thread Tal Galili
Hello Douglas, Thank you for the BLAST!=BLAS correction (I imagine my slip was due to some working I have done recently with an RNA analysis software called BLAST). Also, thank you for the very interesting posting here and in your reply to David's post. My current conclusion from this thread

Re: [R] Can one get a list of recommended packages?

2010-06-12 Thread Tal Galili
Hello David, I am not sure I understood your question. Are you asking what are the packages that the R release comes with? Or are you asking what recommended packages one should have when installing R? (There is a good list to start with

Re: [R] Displaying homogeneous groups in aov post-hoc results ?

2010-06-12 Thread Hadley Wickham
Try multcompView Hadley On Sat, Jun 12, 2010 at 8:42 AM, Tal Galili tal.gal...@gmail.com wrote: Hello dear R-help mailing list, A friend of mine teaches a regression and experimental design course and asked me the following question. She is trying to find a way to display the homogeneous

[R] nonparametric density and probability methods

2010-06-12 Thread Steve Friedman
Hello, I tried to post this earlier, but it seems that it did not appear on the list. If you've rec'd 2 m I'm trying to calculate non-parametric probabilities using the np package and having some difficulties. OS is Windows, R version 2.11.1 Here is what I've done so far. library(np) veg -

Re: [R] Displaying homogeneous groups in aov post-hoc results ?

2010-06-12 Thread Tal Galili
Thank you very much Hadley, exactly what I was looking for. Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) | www.r-statistics.com

Re: [R] Handling character string

2010-06-12 Thread Megh Dal
Thanks Erik for you reply. You have pointed correctly I want to remove the space at the 1st place (if any). In the mean time I have looked into the function sub() and there seems to be one example that mimics my problem : str - '   Now is the time      ' sub('[[:space:]]+$', '', str)[1]   Now

Re: [R] Logic with regexps

2010-06-12 Thread Ted Harding
Thanks, Gabor, for the initiation to perly regexps! I've only been used to extended ones till now. A pity, perhaps, that perl=TRUE is not an option for the likes of browseEnv(), help.search(), list.files() and ls() (which take extended regexps), but one can always assign the output and then

Re: [R] Can one get a list of recommended packages?

2010-06-12 Thread Douglas Bates
On Sat, Jun 12, 2010 at 8:37 AM, Dr. David Kirkby david.kir...@onetel.net wrote: R 2.10.1 is used in the Sage maths project. Several recommended packages (Matrix, class, mgcv, nnet, rpart, spatial, and survival) are failing to build on Solaris 10 (SPARC). Have you checked the dependencies for

Re: [R] Problems with BRugs

2010-06-12 Thread Uwe Ligges
With BRugsFit() of the current verison of BRugs, the not given inits are not generated automatically (which may change with the next version). Hence you need to do it manually as in: library(BRugs) modelCheck(test.bug) modelData(bugsData(Data)) modelCompile(numChains=1)

Re: [R] Can one get a list of recommended packages?

2010-06-12 Thread Dr. David Kirkby
On 06/12/10 03:31 PM, Tal Galili wrote: Hello David, I am not sure I understood your question. Sorry, perhaps I should have rephrased it better. Are you asking what are the packages that the R release comes with? Sort of. When R is configured, there is an option

[R] Fast way to compute largest eigenvector

2010-06-12 Thread MInh Tang
Hello all, I was wondering if there is a function in R that only computes the eigenvector corresponding to the largest/smallest eigenvalue of an arbitrary real matrix. Thanks Minh -- Living on Earth may be expensive, but it includes an annual free trip around the Sun.

Re: [R] Transforming simulation data which is spread acrossmanyfiles into a barplot

2010-06-12 Thread Ian Bentley
Thanks Gabor - I was able to use that for my purposes. On 11 June 2010 16:27, Bert Gunter gunter.ber...@gene.com wrote: So two time series? Fair enough. But less is more. Plot them as separates series of points connected by lines, different colors for the two different series. Or as two

Re: [R] generating ordered, random decimal fractions

2010-06-12 Thread Peter Ehlers
I must be missing something; can't you just use, e.g., one of: sort(round(runif(10), 3)) sort(sample(1000, 10, TRUE)/1000) ? -Peter Ehlers On 2010-06-12 7:49, kurt_h...@nps.gov wrote: Greetings How do I do this in R? Checking the Cran site produces a bewildering array of packages that

Re: [R] Handling character string

2010-06-12 Thread jim holtman
This is probably what you want: sub(^[[:space:]]*, , Now is the time) [1] Now is the time You need to anchor it at the beginning with '^' On Sat, Jun 12, 2010 at 10:29 AM, Megh Dal megh700...@yahoo.com wrote: Thanks Erik for you reply. You have pointed correctly I want to remove the

Re: [R] Compiling R with multi-threaded BLAS math libraries - why not actually ?

2010-06-12 Thread Barry Rowlingson
On Sat, Jun 12, 2010 at 3:16 PM, Tal Galili tal.gal...@gmail.com wrote: Hello Douglas, Thank you for the BLAST!=BLAS correction (I imagine my slip was due to some working I have done recently with an RNA analysis software called BLAST). Also, thank you for the very interesting posting here

[R] who know how to program Bartlett lewis model using R

2010-06-12 Thread Huda Ibrahim
I am a graduate student working on disaggregation rainfall using Bartlett Lewis model but i have a problem on how can I make a program to estimate the model parameters and using it for the simulation. who can help of those people who are famous in R.   thanks for all in advance

Re: [R] Handling character string

2010-06-12 Thread Megh Dal
Thanks Jim for this reply. This is the way what I was looking for. However would you please explain me the meaning of ^[[:space:]]* or '[[:space:]]+$'? When should I use ^ or * or +$? Thanks for your time. --- On Sat, 6/12/10, jim holtman jholt...@gmail.com wrote: From: jim holtman

Re: [R] Clustering algorithms don't find obvious clusters

2010-06-12 Thread Henrik Aldberg
Dave, I used daisy with the default settings (daisy(M) where M is the matrix). Henrik On 11 June 2010 21:57, Dave Roberts dvr...@ecology.msu.montana.edu wrote: Henrik, The clustering algorithms you refer to (and almost all others) expect the matrix to be symmetric. They do not seek a

[R] calling a function with new inputs every 1 minute

2010-06-12 Thread KstuS
I have inputs to a function which are changing all the time - I pull these values from the internet. I then apply a function to the values. What I'd like to do is automate the process so it runs every one minute and adds the output of the function as a new element of a vector. Pseudo code: at

Re: [R] Clustering algorithms don't find obvious clusters

2010-06-12 Thread Dave Roberts
Henrik, Given your initial matrix, that should tell you which authors are similar/dissimilar to which other authors in terms of which authors they cite. In this case authors 1 and 3 are most similar because they both cite authors 2 and 4. Authors 2 and 3 are most different because they

[R] Configuration of host address for database access

2010-06-12 Thread julia . jacobson
Hello everybody, Whatever I use as the host address of my PostgreSQL database server, I always get a connection if I'm accessing the database from R on the same machine and none if I'm trying to do so from a remote client. PostgreSQL is running on a windows xp machine as the RDBMS, JDBC, R, DBI

Re: [R] Handling character string

2010-06-12 Thread jim holtman
^[[:space:]]* indicates that you want to match as many (zero or more) spaces at the beginning of the string. The '^' anchors the match at the beginning. [[:space:]]*$ says to match as many (zero or more that is what * asks for) at the end of the string; $ anchors the search at the end of the

Re: [R] mob (party package) question

2010-06-12 Thread Achim Zeileis
On Fri, 11 Jun 2010, tudor wrote: Dear useRs: I try to use mob from the party package (thanks Achim and Co.!) to model based recursive partition a data set. The model is a logistic regression specified with model=glinearModel and family=binomial(). Running mob results in a few warnings of

Re: [R] Can one get a list of recommended packages?

2010-06-12 Thread Dr. David Kirkby
On 06/12/10 05:27 PM, Douglas Bates wrote: On Sat, Jun 12, 2010 at 8:37 AM, Dr. David Kirkby david.kir...@onetel.net wrote: R 2.10.1 is used in the Sage maths project. Several recommended packages (Matrix, class, mgcv, nnet, rpart, spatial, and survival) are failing to build on Solaris 10

Re: [R] Fast way to compute largest eigenvector

2010-06-12 Thread Ravi Varadhan
You can use the power method for computing the dominant eigenvector. A more sophisticated approach (for large matrices) is the Lancsoz algorithm for Hermitian matrices, which is based on the power method. The `arpack' function in the igraph package uses the more general Arnoldi iteration,

[R] help with R

2010-06-12 Thread li li
Hi all, I want to solve the following equation for x with rho - 0.5 pnorm(-x)*pnorm((rho*dnorm(x)/pnorm(x)-x)/sqrt(1-rho^2))==0.05 Is there a function in R to do this? Thank you very much! Hannah [[alternative HTML version

[R] Sim function

2010-06-12 Thread Noah Silverman
I'm reading Gellman's book Data Analysis Using Regression and Multilevel-Hierarchical Models In Chapter 7 (and later), he makes frequent referent to a function names sim. I can't find the function anywhere, not in my standard R install, or in any of the packages. Doe anyone have a

Re: [R] help with R

2010-06-12 Thread John Fox
Dear Hannah, If I understand you correctly, you want a solution when rho is 0.5; if so, f - function(x){ + pnorm(-x)*pnorm((0.5*dnorm(x)/pnorm(x)-x)/sqrt(1-0.5^2)) - 0.05 + } uniroot(f, c(-3, 3)) $root [1] 0.8031289 $f.root [1] -1.565857e-06 $iter [1] 11 $estim.prec [1] 6.103516e-05 I

[R] Scope and sapply

2010-06-12 Thread Worik R
I am puzzled by the scope rules that apply with sapply. If I want to modify a vector with sapply I tried... N - 10 vec - vector(mode=numeric, length=N) test - function(i){ vec[i] - i } sapply(1:N, test) vec but it not work. How can this be done? Worik [[alternative HTML version

Re: [R] Scope and sapply

2010-06-12 Thread Sarah Goslee
Hi, On Sat, Jun 12, 2010 at 9:22 PM, Worik R wor...@gmail.com wrote: I am puzzled by the scope rules that apply with sapply. If I want to modify a vector with sapply I tried... N - 10 vec - vector(mode=numeric, length=N) test - function(i){  vec[i] - i } sapply(1:N, test) vec You

Re: [R] mob (party package) question

2010-06-12 Thread tudor bodea
Achim: Thanks a lot for your suggestions. Tudor Date: Sun, 13 Jun 2010 00:26:08 +0200 From: achim.zeil...@uibk.ac.at To: tudor_bo...@hotmail.com CC: r-help@r-project.org Subject: Re: [R] mob (party package) question On Fri, 11 Jun 2010, tudor wrote: Dear useRs: I try

Re: [R] Scope and sapply

2010-06-12 Thread Gabor Grothendieck
On Sat, Jun 12, 2010 at 9:22 PM, Worik R wor...@gmail.com wrote: I am puzzled by the scope rules that apply with sapply. If I want to modify a vector with sapply I tried... N - 10 vec - vector(mode=numeric, length=N) test - function(i){  vec[i] - i } sapply(1:N, test) vec but it not

Re: [R] Rgui crashed on Windows XP Home

2010-06-12 Thread Jinsong Zhao
On 2010-6-12 5:58, Duncan Murdoch wrote: Jinsong Zhao wrote: Hi there, I just installed R 2.11.1 on my PC, which runs a Windows XP Home. The installation is successful, however, when I double click on the R icon, I get the following error message: R for Windows GUI front-end has encountered

Re: [R] mob (party package) question

2010-06-12 Thread tudor
Achim - Thank you so much for your suggestions. Tudor -- View this message in context: http://r.789695.n4.nabble.com/mob-party-package-question-tp2252500p2253141.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org

[R] Break in the y-axis

2010-06-12 Thread beloitstudent
Hello all, I have been having trouble getting a break in my y-axis. All of my data points are up around 100-200, but the graph has to start at zero, so i would like to remove all the white space using a break symbol. I have been able to get the break and labels to be correct, however, I can't

Re: [R] Sim function

2010-06-12 Thread Dennis Murphy
Hi: The book has an accompanying R package called arm. Within that package is a function named sim, which could be what you're looking for... HTH, Dennis On Sat, Jun 12, 2010 at 6:02 PM, Noah Silverman n...@smartmediacorp.comwrote: I'm reading Gellman's book Data Analysis Using Regression and

[R] using latticeExtra plotting confidence intervals

2010-06-12 Thread Joe P King
I am wanting to plot a 95% confidence band using segplot, yet I am wanting to have groups. For example if I have males and females, and then I have them in different races, I want the racial groups in different panels. I have this minor code, completely made up but gets at what I am wanting, 4

Re: [R] Scope and sapply

2010-06-12 Thread Worik R
I was careless. Here is a better example of what I am trying to. With the '-' you offered. ?- That was exactly what I needed, thankyou. cheers Worik N - 10 ## x simulate a return series x - runif(N)-.5 ## Build an array of cumulative returns of a portfolio starting with $1 as it changes