search for: macqueen1

...tnr.ma and you will see some NAs. -Don -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062 Lab cell 925-724-7509 ? From: Bill Poling <Bill.Poling at zelis.com> Date: Friday, June 1, 2018 at 8:58 AM To: "MacQueen, Don" <macqueen1 at llnl.gov>, array R-help <r-help at r-project.org> Subject: RE: [R] Time-series moving average question Hello Don, thank you for your response. I appreciate your help. ? I am using the forecast package, originally I found it following a forecasting example on bloggers.com ? https://ww...

...the others? Also, I have 5 series of data I am working with using this script and of course each gave me that warning, but only the one series had the same 8 Point Forecast values, is that coincidental you think? Terrific of you to help, I really appreciate it. WHP From: MacQueen, Don [mailto:macqueen1 at llnl.gov] Sent: Friday, June 01, 2018 12:54 PM To: Bill Poling <Bill.Poling at zelis.com>; r-help (r-help at r-project.org) <r-help at r-project.org> Subject: Re: [R] Time-series moving average question You are right that there are no NAs in the practice data. But there are NAs in t...

...n MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062 Lab cell 925-724-7509 From: Luca Meyer <lucam1968 at gmail.com> Date: Monday, April 30, 2018 at 8:08 AM To: Rui Barradas <ruipbarradas at sapo.pt> Cc: "MacQueen, Don" <macqueen1 at llnl.gov>, array R-help <r-help at r-project.org> Subject: Re: [R] How to visualise what code is processed within a for loop Hi Rui Thank you for your suggestion, I have tested the code suggested by you against that supplied by Don in terms of timing and results are very much aligned:...

Thanks. I thought lm() function is for linear model, such as the correlation below: Y= aX + b On Wed, Jun 14, 2017 at 5:25 PM, MacQueen, Don <macqueen1 at llnl.gov> wrote: > Start with the lm() function; i.e., see > > ?lm > > -Don > > -- > Don MacQueen > > Lawrence Livermore National Laboratory > 7000 East Ave., L-627 > Livermore, CA 94550 > 925-423-1062 > > > On 6/14/17, 3:40 PM, "R-help...

...xpr might indeed be the time consuming piece of code. I might try to optimise this piece of code later on, but for the time being I am working on the following part of building a neural network to try indeed classifying some text. Again, thanks, Luca 2018-04-30 17:25 GMT+02:00 MacQueen, Don <macqueen1 at llnl.gov>: > Luca, > > > > If speed is important, you might improve performance by making d0 into a > true matrix, rather than a data frame (assuming d0 is indeed a data frame > at this point). Although data frames may look like matrices, they aren?t, > and they have...

....org/web/packages/forecast/forecast.pdf Since I created this practice data using the structure() function I am unsure why there would be NA?s as there are none apparently in the structure? No worries though, I am going to reach out to the package author. Cheers. WHP From: MacQueen, Don [mailto:macqueen1 at llnl.gov] Sent: Friday, June 01, 2018 11:24 AM To: Bill Poling <Bill.Poling at zelis.com>; r-help (r-help at r-project.org) <r-help at r-project.org> Subject: Re: [R] Time-series moving average question My guess would be that if you inspect the output from ma(dat3[1:28], order=3) yo...

...> > But X can be some function like - sin, cos, log, exp... > > Cheers > Petr > > > -----Original Message----- > > From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of lily li > > Sent: Thursday, June 15, 2017 1:28 AM > > To: MacQueen, Don <macqueen1 at llnl.gov> > > Cc: R mailing list <r-help at r-project.org> > > Subject: Re: [R] about fitting a regression line > > > > Thanks. I thought lm() function is for linear model, such as the > correlation > > below: > > Y= aX + b > > > > On W...

Is there a package or a command that does window aggregation like select sum(col1) over (partition by col2, col3 order by col4 rows between unbounded preceding and current row) as sum1 from table1 ; the above is Netezza syntax, but Postgre has same capability. Stephen B [[alternative HTML version deleted]]

> On 29 Sep 2017, at 22:43 , MacQueen, Don <macqueen1 at llnl.gov> wrote: > > I used to use SAS a lot, but I don't know what the line > *Yield Champagin; > does. Nothing. It's a comment... -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+4...

Um, maybe just dat1$C <- match(dat1$B, unique(dat1$B)) Indexing 1:k with numbers between 1 and k is a bit of a no-op... AFAICT, this even works without stringsAsFactors=FALSE -pd > On 11 May 2018, at 21:30 , MacQueen, Don <macqueen1 at llnl.gov> wrote: > > dat1$C <- seq(length(unique(dat1$B)))[ match( dat1$B, unique(dat1$B) )] -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Office: A 4.23 Email: pd.mes at cbs.dk P...

Hi But X can be some function like - sin, cos, log, exp... Cheers Petr > -----Original Message----- > From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of lily li > Sent: Thursday, June 15, 2017 1:28 AM > To: MacQueen, Don <macqueen1 at llnl.gov> > Cc: R mailing list <r-help at r-project.org> > Subject: Re: [R] about fitting a regression line > > Thanks. I thought lm() function is for linear model, such as the correlation > below: > Y= aX + b > > On Wed, Jun 14, 2017 at 5:25 PM, MacQueen, Don &l...

...nterval 01-01-01 and 9999-12-31, showing numerical value. From: Gabe Becker <becker.gabe at gene.com> Date: Monday, 11 June 2018 at 23:59 To: Emil Bode <emil.bode at dans.knaw.nl> Cc: Joris Meys <jorismeys at gmail.com>, Werner Grundlingh <wgrundlingh at gmail.com>, "macqueen1 at llnl.gov" <macqueen1 at llnl.gov>, r-devel <r-devel at r-project.org> Subject: Re: [Rd] Date class shows Inf as NA; this confuses the use of is.na() format.Date <- function (x, ...) { xx <- format(as.POSIXlt(x), ...) names(xx) <- names(x) xx[is.na<http://is.na...

Emil et al., On Mon, Jun 11, 2018 at 1:08 AM, Emil Bode <emil.bode at dans.knaw.nl> wrote: > I don't think there's much wrong with is.na(as_date(Inf, > origin='1970-01-01'))==FALSE, as there still is some "non-NA-ness" about > the value (as difftime shows), but that the output when printing is > confusing. The way cat is treating it is clearer: it

...> d0[[nm]] <- ifelse( regexpr(d1[i,1], d0$X0) > 0, 1, 0) >>> } >>> >>> is exaclty what I needed. >>> >>> Best regards, >>> >>> Luca >>> >>> >>> 2018-04-25 23:03 GMT+02:00 MacQueen, Don <macqueen1 at llnl.gov>: >>> >>> Your code doesn't make sense to me in a couple of ways. >>>> >>>> Inside the loop, the first line assigns a value to an object named "t". >>>> Then, the second line does the same thing, assigns a value to a...

...? nm <- paste0("V", i) >> ?????? d0[[nm]] <- ifelse( regexpr(d1[i,1], d0$X0) > 0, 1, 0) >> ???? } >> >> is exaclty what I needed. >> >> Best regards, >> >> Luca >> >> >> 2018-04-25 23:03 GMT+02:00 MacQueen, Don <macqueen1 at llnl.gov>: >> >>> Your code doesn't make sense to me in a couple of ways. >>> >>> Inside the loop, the first line assigns a value to an object named "t". >>> Then, the second line does the same thing, assigns a value to an object >&gt...

My guess would be that if you inspect the output from ma(dat3[1:28], order=3) you will find some NAs in it. And then forecast() doesn't like NAs. But I can't check, because I can't find the ma() and forecast() functions. I assume they come from some package you installed; it would be helpful to say which package. -Don -- Don MacQueen Lawrence Livermore National Laboratory 7000

...1])){ newtst_new[c(j, j + 1)] <- newtst[c(i, i + 1)] i <- i + 2 }else{ newtst_new[c(j, j + 1)] <- c(newtst[c(i)], "Fval: ") i <- i + 1 } j <- j + 2 } newtst_new which is also not very pretty. HTH Ulrik On Thu, 13 Jul 2017 at 16:48 MacQueen, Don <macqueen1 at llnl.gov> wrote: > Using Ulrik?s example data (and assuming I understand what is wanted), > here is what I would do: > > ex.dat <- c("FName: fname1", "Fval: Fval1.name1", "Fval: ", "FName: > fname2", "Fval: Fval2.name2", &qu...

I use the method, df$Time = as.POSIXct(df$Time), but it has the warning message: Error in as.POSIXlt.character(x, tz, ...) : character string is not in a standard unambiguous format On Thu, Dec 14, 2017 at 1:31 PM, MacQueen, Don <macqueen1 at llnl.gov> wrote: > In addition to which, I would recommend > > df <- read.table("DATAM", header = TRUE, fill = TRUE, > stringsAsFactors=FALSE) > > and then converting the Time column to POSIXct date-time values using > as.POSIXct() > specifying the form...

I have found that order() fails in a rather arcane circumstance, as in this example: > foo <- I( c('x','\265g') ) > order(foo) Error in if (xi > xj) 1L else -1L : missing value where TRUE/FALSE needed > foo <-c('x','\265g') > order(foo) [1] 1 2 > sessionInfo() R version 3.1.1 (2014-07-10) Platform: x86_64-apple-darwin13.1.0 (64-bit)

> On Jan 26, 2018, at 9:51 AM, MacQueen, Don <macqueen1 at llnl.gov> wrote: > > What Dave said, plus here's a hint. Try this example (which uses base graphics): > > plot(1:5) > plot(1:5, cex.lab=2) > > Then look at the help page for par > help('par') > or > ?par > to search for other graphics parame...