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 wrote: > Roger Coppock 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 inte

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

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 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[

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 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]) >>> #

Re: [R] class of strptime() output

2015-01-05 Thread Jinsong Zhao
On 2015/1/5 17:28, Ben Bolker wrote: Jinsong Zhao 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

Re: [R] class of strptime() output

2015-01-05 Thread Ben Bolker
Jinsong Zhao 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 b

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 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)) + > geom_point(

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) + scale_colour_gradient(low="li

[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

[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

[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 0.15515780

[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). I

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

2015-01-05 Thread Ben Bolker
Roger Coppock 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 likelihood p

[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 ar

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 w

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 addi

[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] #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 <- USvote[sample(1:

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 en

[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 http://

[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 1:length(dataset$rstu

[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 CHA

[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 ___

[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 "aB+Bc+D

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, Janu

Re: [R] sapply function and poisson distribution

2015-01-05 Thread peter dalgaard
On 05 Jan 2015, at 00:21 , Pete Brecknock 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() in the original q

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-

[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 sc