[R] Bayesian PLS-Based Dimension Reduction and Classification Methods for Time to Event Data in R

2014-05-03 Thread kabir opeyemi
Dear all, Please I am currently working on Bayesian Partial Least Squares Dimension Reduction and Classification methods for Modelling survival data in R. Can someone advise on useful R packages and good R books to read. I am also interested in working with experienced colleagues who may be

[R] Construction of Optimal Designs in the Language R

2014-05-03 Thread kabir opeyemi
Dear All, Please can someone help advise on procedures for construction of optimal designs and estimation of parameters in the language R? Guide on useful R libraries, codes and books would be appreciated. Kabir Olorede, Department of Statistics and Mathematical Sciences, Kwara State

Re: [R] matafor package: adding additional information under a forest plot

2014-05-03 Thread Viechtbauer Wolfgang (STAT)
1) There is no function that will do that automatically for you, but it's easy to do with text() and paste(). If you have fitted a random-effects model with rma(), then you can find the Q-statistic in res$QE, the I^2 statistic in res$I2, and probably whatever else you need, and then just

Re: [R] Return TRUE only for first match of values between matrix and vector.

2014-05-03 Thread arun
Hi, This should be little more faster. indx - A==B indx1 - which(indx, arr.ind=TRUE) indx[indx1[duplicated(indx1[,1]),]]- FALSE indx ##Speed comparison ##previous method fun1 - function(mat, vec) {     stopifnot(dim(mat)[1] == length(vec))     indx - mat == vec     t(apply(indx, 1, function(x)

Re: [R] SQL vs R

2014-05-03 Thread Dr Eberhard W Lisse
Thank you very much, Mr Arkell. el On 2014-05-03, 07:11 , Bert Gunter wrote: By making the effort to learn R? See e.g. the Introduction to R tutorial that ships with R. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 Data is not information. Information is

[R] help me

2014-05-03 Thread akram yazdani
hi please help me for below warning . may model is in frailtypack with frailtypenal mod.sha.gap - frailtyPenal(Surv(s,time,event) ~ cluster(id)+sex1+age2+bmi1+whr1+tg1+hdl1+chol1+por+pr2+car1+fib2+fat2+terminal(d),formula.terminalEvent =

[R] auglag store parameters

2014-05-03 Thread Delger Enkhbayar
Hi, I am minimizing a non-linear function subject to linear and nonlinear equalities and inequalities. I managed to understand how to use auglag for this case. Now, I'm using loops to get solutions for different assumptions, however, I'm facing a difficulty in storing just $par part of the answer

[R] auglag store parameters

2014-05-03 Thread Delger Enkhbayar
Hi, I am minimizing a non-linear function subject to linear and nonlinear equalities and inequalities. I managed to understand how to use auglag for this case. Now, I'm using loops to get solutions for different assumptions, however, I'm facing a difficulty in storing just $par part of the answer

[R] Logistic regression help!

2014-05-03 Thread Si Qi L.
Hi guys, I have a trouble to solve the specificity and senstitivity for a logistic regression model. I really need your help, would you please help me out? :) Thank you!! This is the model I constructed: model=glm(Status ~ Gender.compl+ X2.4.times.per.month+

[R] error in r

2014-05-03 Thread thanoon younis
dear all members i have error in the code below Error in Y[i, 9] = 0.9 * XI[i, 2] + eps[9] : subscript out of bounds is there anyone helps me please. many thanks in advance thanoon llibrary(mvtnorm) #Load mvtnorm package library(R2WinBUGS) #Load R2WinBUGS package N=500

Re: [R] error in r

2014-05-03 Thread Duncan Murdoch
On 03/05/2014, 11:49 AM, thanoon younis wrote: dear all members i have error in the code below Error in Y[i, 9] = 0.9 * XI[i, 2] + eps[9] : subscript out of bounds is there anyone helps me please. You created Y with 8 columns, then you refer to column 9. Duncan Murdoch many thanks in

[R] knitr - Highlight code/output

2014-05-03 Thread Axel Urbiz
Hello, I find situations where some segments of the code are displayed in the output but not entirely highlighted. I guess there should be a way to fix this, but I could't find it in the options. Here's an example with knitr/LaTex. \documentclass{article} \begin{document} \section{Example 1}

Re: [R] auglag store parameters

2014-05-03 Thread David Winsemius
On May 3, 2014, at 5:07 AM, Delger Enkhbayar wrote: Hi, I am minimizing a non-linear function subject to linear and nonlinear equalities and inequalities. I managed to understand how to use auglag for this case. Now, I'm using loops to get solutions for different assumptions, however, I'm

Re: [R] Seeking well-commented versions of mgcv source code

2014-05-03 Thread William Dunlap
You can download the source code for the package with download.packages(mgcv, type=source, destdir=/tmp) Bill Dunlap TIBCO Software wdunlap tibco.com On Fri, May 2, 2014 at 7:15 PM, Andrew Crane-Droesch andre...@gmail.com wrote: Dear List, I'm looking for well-commented versions of

[R] mgcv gam bs=re questions

2014-05-03 Thread William Shadish
Dear R-helpers, I am working on a project assessing the prevalence and variance (random effects) of linear and nonlinear trends in a data set that has short time series (each time series identified as PID 1 through 5). I am using mgcv (gam) with the bs=”re” option (more on why not gamm or

Re: [R] SQL vs R

2014-05-03 Thread Dr Eberhard W Lisse
Thanks, will try to figure this out :-)-O el On 2014-05-03, 06:40 , Carlos Ortega wrote: Hi, With the new package dplyr you can create equivalent SQL sintaxt queries like the one you need. You can find examples of how to apply it here:

Re: [R] Seeking well-commented versions of mgcv source code

2014-05-03 Thread Andrew Crane-Droesch
It turns out that commented versions of the source code are available on GitHub: https://github.com/cran/mgcv/blob/master/R/mgcv.r __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] SQL vs R

2014-05-03 Thread Rolf Turner
On 04/05/14 00:05, Dr Eberhard W Lisse wrote: Thank you very much, Mr Arkell. I don't get it. Can anyone explain the (joke? allusion?) ? cheers, Rolf Turner On 2014-05-03, 07:11 , Bert Gunter wrote: By making the effort to learn R? See e.g. the Introduction to R tutorial that ships

Re: [R] SQL vs R

2014-05-03 Thread Sarah Goslee
On Sat, May 3, 2014 at 5:42 PM, Rolf Turner r.tur...@auckland.ac.nz wrote: On 04/05/14 00:05, Dr Eberhard W Lisse wrote: Thank you very much, Mr Arkell. I don't get it. Can anyone explain the (joke? allusion?) ? I believe it's a moderately offensive reply from someone who feels unfairly

Re: [R] SQL vs R

2014-05-03 Thread Dr Eberhard W Lisse
Google Pressdram :-)-O el On 2014-05-03, 23:42 , Rolf Turner wrote: On 04/05/14 00:05, Dr Eberhard W Lisse wrote: Thank you very much, Mr Arkell. I don't get it. Can anyone explain the (joke? allusion?) ? cheers, Rolf Turner __

Re: [R] Question on Cluster Package, agnes() function

2014-05-03 Thread Martin Maechler
Anna F... on Thu, 1 May 2014 22:09:28 + writes: Hi Martin, I am a statistician at National Jewish Health in Colorado, and I have been working on clustering a dataset using Ward's minimum variance. When plotting the dendrogram, the y-axis is labeled as 'height'. Can you

Re: [R] SQL vs R

2014-05-03 Thread Rolf Turner
On 04/05/14 10:16, Dr Eberhard W Lisse wrote: Google Pressdram :-)-O el On 2014-05-03, 23:42 , Rolf Turner wrote: On 04/05/14 00:05, Dr Eberhard W Lisse wrote: Thank you very much, Mr Arkell. I don't get it. Can anyone explain the (joke? allusion?) ? Thank you. cheers, Rolf Turner

Re: [R] SQL vs R

2014-05-03 Thread Rolf Turner
On 04/05/14 09:58, Sarah Goslee wrote: SNIP Personally, I'm not sure how much work someone who appears to have not read the posting guide should really expect the list to do on his behalf. But snarky replies to reasonable requests to read the documentation are easier than doing one's own work.

[R] panel data: variable selection AND multicollinearity

2014-05-03 Thread phil
Dear all, I am running a regression on panel data (127 observations). The Hausman Test indicats that the random-effects model is superior to the fixed-effects model. However, for a linear regression lm() there is a simple command in R to carry out a variable selection procedure: step(). Is there

Re: [R] SQL vs R

2014-05-03 Thread Satish Anupindi Rao
By making the effort to learn R?? very constructive and not condescending at all. We, lesser beings, are indebted to you, sir. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Bert Gunter Sent: Saturday, May 03, 2014 1:12 AM To:

Re: [R] SQL vs R

2014-05-03 Thread Jeff Newmiller
?table ?aggregate Also, packages plyr, data.table, and dplyr. You might consider reading [1], but if your interests are really as simple as your examples then the table function should be sufficient. That function is discussed in the Introduction to R document that you really should have read

Re: [R] legendre quadrature

2014-05-03 Thread Berend Hasselman
On 01-05-2014, at 09:43, pari hesabi statistic...@hotmail.com wrote: Hello everybody I need to approximate the amount of integral by using legendre quadrature. I have written a program which doesn't give me a logical answer; Can anybody help me and send the correct program? For example

Re: [R] SQL vs R

2014-05-03 Thread David Winsemius
On May 3, 2014, at 9:10 PM, Satish Anupindi Rao wrote: By making the effort to learn R?? very constructive and not condescending at all. We, lesser beings, are indebted to you, sir. For Pete's sake. The OP didn't even express his original request in natural language or offer a working

[R-es] Duda sobre como unificar datos en una oración

2014-05-03 Thread Daliana Ramos Garcia
Hola a todos, tengo una duda sobre como unificar en R, datos de una oración Por ejemplo tengo la siguiente oración Siempre que, el atributo A es alto y el atributo B es media y el atributo C es alto y el atributo X es bajo; entonces el atributo Z es alto y quiero que me quede de la