Re: [R] se.fit option to the predict.nls() function

2015-01-05 Thread Rune Haubo
On 5 January 2015 at 21:08, Ben Bolker bbol...@gmail.com wrote: Roger Coppock rcoppock at cox.net writes: When will R implement the se.fit option to the predict.nls() function? Is there some schedule? I think this is unlikely to happen, ever (sorry). The exact method for finding

[R] What is the problem now?

2015-01-05 Thread Levent TERLEMEZ
Hi Dear Users, I abonden ggplot2 adventure for the moment because of urgent results (i will of course will do this with ggplot2 because I would like to use ggplot2 much more than standart plots). But this time I came across with a different problem. Problem is different graphical outputs from

Re: [R] What is the problem now?

2015-01-05 Thread PIKAL Petr
Hi I am not sure what do you want to achieve. I would change your text function {text(min(gecici$Tarih)+0.1,y,paste(z,sigma),cex=.85,pos=3)} and instead of copy paste you can open and close appropriate device (pdf, jpg or emf - from devEMF package) Cheers Petr -Original Message-

Re: [R] sapply function and poisson distribution

2015-01-05 Thread peter dalgaard
On 05 Jan 2015, at 00:21 , Pete Brecknock peter.breckn...@bp.com wrote: n - c(1,2,3,4,5) lambda - c(0.1,0.8,1.2,2.2,4.2) mapply(function(x,y) rpois(x,y), n, lambda) Yes. I'd throw in a SIMPLIFY=FALSE to avoid getting results in a different format if n is constant (then again, sapply()

Re: [R] What is the problem now?

2015-01-05 Thread PIKAL Petr
Hi cc to list, others can have further comments substitute is superfluous see results of substitute(paste(z,sigma) paste(z,sigma) substitute(paste(z,sigma) paste(z,sigma) Cheers Petr -Original Message- From: Levent TERLEMEZ [mailto:lterle...@anadolu.edu.tr] Sent: Monday, January

Re: [R-es] Resumen de R-help-es, Vol 71, Envío 5

2015-01-05 Thread Jose Manuel Veiga del Baño
Muchas gracias. Solucionado.  r-help-es-requ...@r-project.org escribió: Envíe los mensajes para la lista R-help-es a         r-help-es@r-project.org Para subscribirse o anular su subscripción a través de la WEB         https://stat.ethz.ch/mailman/listinfo/r-help-es O por correo

Re: [R-es] Resumen de R-help-es, Vol 71, Envío 5

2015-01-05 Thread Jose Manuel Veiga del Baño
Muchas gracias. Solucionado.  r-help-es-requ...@r-project.org escribió: Envíe los mensajes para la lista R-help-es a         r-help-es@r-project.org Para subscribirse o anular su subscripción a través de la WEB         https://stat.ethz.ch/mailman/listinfo/r-help-es O por correo

[R] Subsets of Boolean string model

2015-01-05 Thread Alrik Thiem
Dear R-help list, I have a problem regarding text manipulation in R, where my basic knowledge doesn't suffice anymore. It might be a bigger problem, but any help would be greatly appreciated and acknowledged. As input, I have a character string representing some Boolean function, such as

Re: [R] #library(CHAID) - Cross validation for chaid

2015-01-05 Thread Max Kuhn
You can create your own: http://topepo.github.io/caret/custom_models.html I put a prototype together. Source this file: https://github.com/topepo/caret/blob/master/models/files/chaid.R then try this: library(CHAID) ### fit tree to subsample set.seed(290875) USvoteS -

Re: [R] controlling the gamma and correlation axis tick label of taylor diagram in R

2015-01-05 Thread David Winsemius
On Jan 4, 2015, at 10:15 PM, കുഞ്ഞായി kunjaai wrote: Dear all, I am trying to control the font size of correlation axis and gamma label (0.6, 0.4 , 0.2 labeling of attached plot) of Taylor diagram in R. I am using **plotrix** package for making Taylor diagram. Nothing attached in what

[R] problem in R

2015-01-05 Thread thanoon younis
Dear all R members I have a problem in the R when i wanted to call winbugs as following #Input data set for WinBUGS thd18 -read.table(C:/Users/hp/Desktop/thd18.txt) #Input data set for WinBUGS data-list(N1=2000,N2=2000,P=9,R=Ro,z1=yo1,z2=yo2,thd18) #Call WinBUGS model-bugs

Re: [R] What is the problem now?

2015-01-05 Thread David Winsemius
On Jan 5, 2015, at 2:46 AM, PIKAL Petr wrote: Hi cc to list, others can have further comments substitute is superfluous see results of substitute(paste(z,sigma) paste(z,sigma) substitute(paste(z,sigma) paste(z,sigma) Perhaps this would be more instructive (after adding back the

[R] ggplot2, geom_point with some symbols unfilled

2015-01-05 Thread Roy Mendelssohn - NOAA Federal
Hi All: I have a data frame the first part of which is the following: date lonlat mean chlcol missing 1 2003-04-23 203.899 19.664 0.0730 0.15176637 0 2 2003-04-24 204.151 19.821NaNNaN 1 3 2003-04-30 203.919 20.351 0.0740

[R] glmmADMB model fails when fixed effect interaction term is removed

2015-01-05 Thread cara715
Hello, I am looking at a data set predicting number of plants (FERO) based on fixed effects sowing method (Method) and sowing rate (Rate). I also have a slope random effect for the site and intercept random effects for the nested variables Block and Strip.ID (to account for experimental design).

Re: [R] class of strptime() output

2015-01-05 Thread Ben Bolker
Jinsong Zhao jszhao at yeah.net writes: In the following code snippet, # a - strptime(121114 0510, %m%d%y %H%M) b - data.frame(date = a, res = 1:5) class(a) class(b[1,1]) # I am wondering why the class of a and b[1,1] are not the same. How to make the class of a and b[1,1] to be same?

[R] GSTAT Co-Kriging: How can i fix a non positive definte correlation matrix

2015-01-05 Thread Roberts, Blake
QUESTION: Can anyone give me advice on how to resolve an error message regarding a non positive definite correlation matrix so that the computed negative cross-variogram can be used to perform co-kriging? I am performing co-kriging on negatively correlated parameters. data = dtrd x y

Re: [R] ggplot2, geom_point with some symbols unfilled

2015-01-05 Thread Roy Mendelssohn - NOAA Federal
Perfect. Thanks. -Roy On Jan 5, 2015, at 5:20 PM, Ista Zahn istaz...@gmail.com wrote: You were close, but a) if you want to map shape to something it needs to be inside an aes() call, and b) shape 1 is also an open circle. I think this does what you want: ggplot(temp,aes(x=lon,y=lat)) +

[R] class of strptime() output

2015-01-05 Thread Jinsong Zhao
Hi there, In the following code snippet, # a - strptime(121114 0510, %m%d%y %H%M) b - data.frame(date = a, res = 1:5) class(a) class(b[1,1]) # I am wondering why the class of a and b[1,1] are not the same. How to make the class of a and b[1,1] to be same? I would really appreciate you for any

Re: [R] ggplot2, geom_point with some symbols unfilled

2015-01-05 Thread Ista Zahn
You were close, but a) if you want to map shape to something it needs to be inside an aes() call, and b) shape 1 is also an open circle. I think this does what you want: ggplot(temp,aes(x=lon,y=lat)) + geom_point(aes(colour=chlcol, shape=factor(missing)),size=3) +

Re: [R] class of strptime() output

2015-01-05 Thread Jinsong Zhao
On 2015/1/5 17:28, Ben Bolker wrote: Jinsong Zhao jszhao at yeah.net writes: In the following code snippet, # a - strptime(121114 0510, %m%d%y %H%M) b - data.frame(date = a, res = 1:5) class(a) class(b[1,1]) # I am wondering why the class of a and b[1,1] are not the same. How to make the

Re: [R] class of strptime() output

2015-01-05 Thread Jinsong Zhao
On 2015/1/5 20:51, David Winsemius wrote: On Jan 5, 2015, at 5:47 PM, Jinsong Zhao wrote: On 2015/1/5 17:28, Ben Bolker wrote: Jinsong Zhao jszhao at yeah.net writes: In the following code snippet, # a - strptime(121114 0510, %m%d%y %H%M) b - data.frame(date = a, res = 1:5) class(a)

Re: [R] class of strptime() output

2015-01-05 Thread Jeff Newmiller
Because data.frame converts it. I don't know what rationale was originally used to justify including that transformation automatically, but I think it is much better this way than having POSIXlt in the data frame. POSIXlt is itself like a little data frame, and having data frames inside of data

[R-es] Eventos de una serie de tiempo

2015-01-05 Thread Carlos Andres Perez Angarita
Saludo estimados compañeros y compañeras Tengo una matriz de datos de 51 filas por 160 columnas, Cada columna es una serie de tiempo. Existe alguna libreria q me calcule la posicion de la fila donde se hallan estos maximos y minimos o los eventos de estas curvas Agradezco la atencion CARLOS

Re: [R] class of strptime() output

2015-01-05 Thread David Winsemius
On Jan 5, 2015, at 5:47 PM, Jinsong Zhao wrote: On 2015/1/5 17:28, Ben Bolker wrote: Jinsong Zhao jszhao at yeah.net writes: In the following code snippet, # a - strptime(121114 0510, %m%d%y %H%M) b - data.frame(date = a, res = 1:5) class(a) class(b[1,1]) # I am wondering why the

[R] #library(CHAID) - export associated probabilities and rules for each row

2015-01-05 Thread Rodica Coderie via R-help
Hello, How can I export the sample data frame below -USvoteS with the associated probabilities and rules for each row? library(CHAID) example(chaid, package = CHAID) To be more specific, for each of 1000 rows of the USvoteS data frame I want to see in which node is that row in (what are the

Re: [R] se.fit option to the predict.nls() function

2015-01-05 Thread Ben Bolker
Roger Coppock rcoppock at cox.net writes: When will R implement the se.fit option to the predict.nls() function? Is there some schedule? I think this is unlikely to happen, ever (sorry). The exact method for finding confidence intervals on nonlinear fits would be to compute

[R] #library(CHAID) - Cross validation for chaid

2015-01-05 Thread Rodica Coderie via R-help
Hello, Is there an option of cross validation for CHAID decision tree? An example of CHAID is below: library(CHAID) example(chaid, package = CHAID) How can I use a 10 fold cross-validation for CHAID? I've read that caret package is to cross-validate on many times of models, but model CHAID is

[R] find and remove outlier

2015-01-05 Thread Methekar, Pushpa (GE Transportation, Non-GE)
Hi, this is my function to find rstudent of model which will give me outlier But I wonder if I could find out the exact no. of outlier in data set. Like outlierTest() does from car package . rm.outliers = function(dataset,model){ dataset$rstudent = rstudent(model) for(i in

[R] se.fit option to the predict.nls() function

2015-01-05 Thread Roger Coppock
When will R implement the se.fit option to the predict.nls() function? Is there some schedule? __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] find and remove outlier

2015-01-05 Thread John Fox
Dear Pushpa Methekar, You can use which.max(), as in library(car) mod - lm(prestige~ income + education, data=Duncan) which.max(abs(rstudent(mod))) minister 6 In fact, this is exactly what outlierTest() does, as you can see by looking at car:::outlierTest.lm Beyond that, I'd

[R] controlling the gamma and correlation axis tick label of taylor diagram in R

2015-01-05 Thread കുഞ്ഞായി kunjaai
Dear all, I am trying to control the font size of correlation axis and gamma label (0.6, 0.4 , 0.2 labeling of attached plot) of Taylor diagram in R. I am using **plotrix** package for making Taylor diagram. Thank you all in advance -- DILEEPKUMAR. R J R F, IIT DELHI